Journal of Applied Physiology AJP: Advances in Physiology Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 91: 1687-1693, 2001;
8750-7587/01 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
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 Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Frey, U.
Right arrow Articles by Suki, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Frey, U.
Right arrow Articles by Suki, B.
Vol. 91, Issue 4, 1687-1693, October 2001

Analysis of the harmonic content of the tidal flow waveforms in infants

Urs Frey1, Michael Silverman2, and Bela Suki3

1 Department of Paediatrics, University Hospital of Berne, 3010 Berne, Switzerland; 2 Department of Child Health, Leicester University, Leicester LE2 7LX, United Kingdom; and 3 Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The aim of this study was to examine whether the spectral characteristics of tidal flow waveform reflect the interaction between the control of breathing and lung mechanics in 10 healthy infants (H), 10 infants with a history of wheezing disorders (W), and 10 infants with chronic lung disease (CLD). From the flow waveform, we calculated a shape index, the harmonic distortion (kd), which quantifies the extent to which a periodic signal deviates from a sine wave. The kd of the entire tidal flow waveform did not significantly discriminate between diagnostic groups. However, kd was sensitive to maturation: it increased from 0.26 at 1 mo to 0.37 at 6 mo of age (P < 0.002). Furthermore, the frequency (f) spectra of the flow (V) amplitudes between 0.13 and 10 Hz followed a power law: V(f) ~ f-s, where s (slope) is the exponent in the power law. The exponent of the healthy infants s(H) was 4.24 [95% confidence interval (CI) = 0.2] at 1 mo, 4.39 (CI = 0.16) at 6 mo, and 4.35 (CI = 0.19) at 12 mo and not significantly changing with age. The mean value of s(W) was marginally lower (4.09 ± 0.28; P < 0.05) than that of s(H). The mean s(CLD) was significantly lower (3.04 ± 0.31; P < 0.001). Lower values of s and higher values of kd indicate an increased complexity of the feedback mechanisms determining tidal flow waveform and may be associated with disease.

respiratory function tests; control of breathing; respiratory mechanics


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

IT HAS BEEN PROPOSED THAT simple tidal flow indexes, such as the ratio of time to peak tidal expiratory flow (TPTEF) to expiratory time (TE) (TPTEF/TE) (13), which can be assessed noninvasively (e.g., Refs. 6, 10, 13), may contain information regarding airway mechanics in infants. However, the early phase of expiration is not only determined by respiratory mechanics but, to a large extent, by the postinspiratory activity of the respiratory muscle drive (19). Ueda et al. (18) demonstrated a close correlation between TPTEF/TE and the pressure 100 ms after airway occlusion at onset of inspiration, a measure of inspiratory drive in infants. Furthermore, Thomas et al. (16, 17) have recently shown in infants that inspiratory and expiratory phases are related to respiratory drive, peripheral chemoreceptor sensitivity, and sleeping patterns in infants. Thus inspiratory and expiratory waveforms are closely related, and the expiratory waveform is not independent of breathing frequency, inspiratory waveform, and inspiratory off-switch mechanisms (3, 8, 11, 12). Thus it seems more appropriate to consider the whole flow waveform as an integrated output of the neural respiratory oscillator in the brain stem, reflecting its interaction with the various chemo- and stretch-receptor feedback mechanisms and the passive mechanical properties of the lung and the chest wall.

Therefore, instead of analyzing the expiratory flow waveform using simple indexes (e.g., TPTEF/TE) based on limited information, one should define and use parameters that are more suitable to characterize the tidal flow wave as the output of an integrated control system. The neurorespiratory control system is a highly complex feedback system that maintains a certain dynamic equilibrium in health. In disease, this dynamic equilibrium may become disturbed if one of its components, for example the mechanical properties of the airways or the chemo- or stretch-receptor feedback, is altered. We hypothesize that such a disturbance to the system leads to a concomitant change in the shape of the periodic flow waveform. Changes in the shape of a periodic waveform must in turn be reflected in the spectral characteristics of the waveform.

Accordingly, the aims of this study were 1) to develop simple indexes that are suitable to characterize the shape and hence the spectral features of the tidal flow waveform in the frequency domain; 2) to assess whether these new parameters are sensitive to alterations in respiratory control and/or lung mechanics because of developmental changes and disease; and 3) to test whether these indexes are more reliable and robust descriptors of the changes in respiratory control and/or airway mechanics than TPTEF/TE, which is one of the most commonly used tidal flow indexes.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The study includes analysis of data obtained from both new experiments and newly analyzed tidal flow data sets from previously published studies on TPTEF/TE (1, 4, 5). First, we examined tidal flow waveforms of 10 healthy infants during development at the ages of 1, 6, and 12 mo. Based on these data, we introduced two new indexes derived from the tidal flow waveform in the frequency domain. We then compared the breath-to-breath variabilities of the new parameters with that of TPTEF/TE calculated from the same series of tidal flow data.

Second, we compared the new indexes obtained from healthy infants with measurements in infants with wheezing disorders (10 infants) and chronic lung disease (CLD) of prematurity (10 infants). To quantify the degree of airway obstruction, we performed forced flow-volume loops using the rapid thoracic compression (RTC) technique in each group.

Third, to assess whether any of the tidal flow waveform parameters were associated with altered airway mechanics, we reanalyzed data from previously published bronchial challenge studies (1, 5). Histamine challenge tests were performed as previously described in healthy infants (1, 5), and the data were reanalyzed using our new indexes.

Experimental Study

Subjects. We analyzed tidal flow waveforms in 10 healthy infants, each measured at the ages of 1, 6, and 12 mo; 10 infants with a known history of wheezing disorders [mean age, 13.9 ± 1.6 (SD) mo] asymptomatic at the time of measurement; and 10 infants with CLD of prematurity [mean postnatal age, 8.1 ± 7.2 (SD) mo and gestational age, 27.7 ± 2.4 wk]. The mean inspired O2 fraction of the CLD group was 29.1 ± 6.1%, whereas in all other groups it was 21%. The characteristics of these infants are given in Table 1. In the 10 healthy infants, 14 histamine challenge tests were performed (in 4 infants at 2 different ages). The studies were authorized by the Ethics Committee of the Hammersmith Hospital, London, UK, where the measurements were performed. Parental consent was obtained for all studies.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Physical characteristics of subjects

Measurements. All infants were sedated by using a maximal dose of 150 mg/kg triclofos sodium. Partial expiratory flow at functional residual capacity (Vmax FRC) was measured by using the RTC technique (5). The supine infants wore an inflatable polythene thoracoabdominal jacket (Medical Engineering Department, Royal Postgraduate Medical School, London, UK) with arms out. Flow was measured by using a face mask (size 1, Rendell Baker Soucek, Ambu International, Bath, Avon, UK) and a Fleisch no. 1 (Gould) pneumotachograph equipped with a differential pressure transducer (MP45, Validyne, Northridge, CA). The linearity of the pneumotachograph was found to be accurate within 2% up to flows of 1 l/s. The frequency response of the flow transducer was flat up to 12 Hz. The flow signals were low-pass filtered (hardware filter, cutoff at 20 Hz), sampled at 100 Hz using a 12-bit analog-to-digital converter, and processed by using the RASP software (Physiologic, Newbury, UK).

In all infants, tidal flow sequences of 8-10 s and forced expiratory flows were recorded by using the RTC technique described in detail previously (5). After a series of initial RTC measurements to determine the optimal jacket pressure, 10 RTC recordings were carried out at end-tidal inspiration using this optimal jacket pressure, and the corresponding Vmax FRC was calculated. The mean value of all technically satisfactory Vmax FRC was determined. Transcutaneous PO2 (TMC3, Radiometer, Copenhagen, Denmark) and transcutaneous arterial O2 saturation (Biox 3740, Omeda, Omaha, NE) were measured. The head position was kept fixed between forced flow measurements and tidal breathing measurements. For both techniques, we used the same face mask and putty ring.

Bronchial challenge tests. As previously described (5), after the baseline measurements in the healthy infants, we administered histamine in doubling concentrations until Vmax FRC decreased by at least 30%. Sequences of tidal breathing were recorded and analyzed before and after the bronchial challenge test.

Analysis of tidal flow signals. The primary aim of the analysis was to develop a method to characterize the waveform of the flow signal. Because the flow signal in the time domain does not appear to follow a simple shape (e.g., a sinusoid), we analyzed the shape in the frequency domain. In particular, we examined the relationship between breathing frequency and its higher harmonics, which are the frequency components of the signal at integer multiples of the breathing rate. The harmonics can be calculated by using the fast Fourier transform (FFT) algorithm applied to one or more periods of the signal. The FFT is invariably applied to several individual breathing cycles, and the individual spectra are then ensemble averaged. Unfortunately, the time period of the breathing cycle can vary by 20-50% from cycle to cycle. One solution is to find the longest cycle and zero pad all other cycles to obtain blocks of data having the same number of time points. When the spectra of these blocks are ensemble averaged, the resulting spectrum will be a smeared version of the individual spectra, because frequency of breathing with the highest energy content would vary from block to block. If, on the other hand, the FFT is applied to the entire signal, spectral leakage will occur, as the breathing frequency will not be the same in the individual cycles. In both cases, the frequency dependence of the higher harmonics will be masked by these effects. To avoid these problems, we first isolated all individual respiratory flow cycles using a zero-crossing detection algorithm that identifies the beginning and the end of a cycle. The number of data points (Ni, where i is the cycle number) was different for the different cycles. Because we were interested in the shape of the spectra (rather than the absolute frequencies), instead of zero padding the time domain data and using FFT of the same length for all cycles, we applied the discrete Fourier transform of window length Ni (MATLAB, 5.0, MATHWORKS) to the individual flow cycles. We then analyzed the power spectra of the individual respiratory cycles based on the two methods described below. No filtering or windowing was used, as these operations affect the shape of the spectra. However, to smooth the spectra, we applied a two-point box car averaging technique. The above analysis was then repeated for each respiratory cycle of a single recording. This method allowed us to define parameters (see below) that can characterize the shape of the power spectra on a cycle-by-cycle basis and then calculate the average and variability of these parameters.

Interpretation of Flow Signal Harmonics

We examined two methods of characterizing the shape of the flow signal as expressed by the frequency dependence of the higher harmonics in the flow signal. These methods are the 1) power law analysis (e.g., Refs. 7, 9) and the 2) harmonic distortion method (15, 20).

Power law analysis. As an example, Fig. 1 shows the tidal waveform of a healthy 6-mo-old infant. The corresponding power spectra of the tidal flow signals (Fig. 1) of the infant are shown in Fig. 2. The spectra apparently decrease linearly on the log-log graph. Thus the logarithm of the harmonics of the flow signals appears to be linearly related to the logarithm of frequency. Indeed, a linear regression provides good fit to the data, as shown in Fig. 2. This can be formally written as
log [<A><AC>V</AC><AC>˙</AC></A>(<IT>f</IT>)]<IT>=c−s ∗ </IT>log (<IT>f</IT>) (1)
where V is flow, f is frequency, c is the y-intercept, and s is the negative slope of the best-fit linear equation. Taking the inverse logarithm of Eq. 1, we obtain
<A><AC>V</AC><AC>˙</AC></A>(<IT>f</IT>)<IT>=A ∗ f</IT><SUP>−<IT>s</IT></SUP> (2)
where A is log (c), another constant. Eq. 2 demonstrates that the flow harmonics follow a power law in the frequency domain. The parameter s is called the exponent in the power law. A decrease in s represents a slower rate of decrease of flow harmonics and, hence, may signify a stronger contribution of higher harmonics to the tidal flow shape. We also note that, because a sine wave has no higher order harmonics, the stronger the harmonics are, the more nonsinusoidal the tidal waveform. Thus, in the power law analysis method, an increase in s (steeper slope) is associated with a lower content of harmonics at higher frequencies, resulting in a more sinusoidal shape, whereas a smaller s (flatter slope) represents higher amplitudes of the harmonics and, hence, a less sinusoidal shape.


View larger version (8K):
[in this window]
[in a new window]
 
Fig. 1.   Example of a tidal flow waveform including the ratio of time to peak expiratory flow to total expiratory time in a healthy infant 6 mo of age.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 2.   Example of the flow power spectrum of the infant with chronic lung disease (CLD) in Fig. 1 in a log-log presentation demonstrating a linear decrease of the amplitude of the harmonics of the fundamental breathing frequency. Slope s represents the slope of the power law.

Harmonic distortion. The extent to which a periodic signal deviates from a single sine wave can also be characterized by the so-called harmonic distortion index (kd), defined as the square root of the relative power in the signal above the fundamental frequency (15, 20). In other words, kd is the square root of the ratio of the sum of the squares of the higher harmonic amplitudes above the fundamental frequency F and the sum of the squares of all amplitudes. The kd is usually specified in percent. For example, the harmonic distortion of a sine wave is obviously zero, whereas that of a triangular wave is 37%. Higher values of kd indicate the presence of stronger harmonic components in the signal and, hence, a shape in the time domain that deviates more from an ideal sinusoidal wave.

Statistical Analysis

In individual infants, s, kd, and TPTEF/TE were calculated for each respiratory cycle and averaged (mean) for the whole flow signal. In healthy infants, the group means (and 95% confidence interval) of the parameters s, kd, and TPTEF/TE were calculated for the 1-, 6-, and 12-mo-old groups. Because the healthy infant groups at 1, 6, and 12 mo consisted of longitudinal data sets, we tested age dependence of the paramaters s, kd, and TPTEF/TE by using paired t-tests. Because there were only three age groups and the main change was in the first 6 mo of age, we used paired t-tests instead of ANOVA statistics.

We compared the infants with CLD and the infants with wheezing disorders with the healthy infants of similar ages. The s, kd, and TPTEF/TE in the infants with CLD were compared with those in the healthy infants at the age of 6 mo, whereas the wheezy infants were compared with healthy infants at the age of 12 mo (see Tables 1-3). The mean postconceptional as well as postnatal ages of the compared groups were not statistically different.

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Parameters


                              
View this table:
[in this window]
[in a new window]
 
Table 3.   Comparison of respiratory cycle variabilities expressed by its coefficient of variations

To compare the breath-to-breath variability of the parameters s, kd, and TPTEF/TE, the coefficient of variation (CV) for all respiratory cycles around the mean value was calculated and expressed in percentage of the mean CV for each individual subject. For each group, the mean CV was calculated for each parameter. Within these particular groups, the CVs of s, kd, and TPTEF/TE were compared by using a paired t-test because the mean CVs of the various physiological and diagnostic groups were not expected to be similar. However, to account for the effect of multiple comparisons, we set the limits of the P value to 0.01 instead of the usual 0.05.

Posthistamine parameters of s, kd, and TPTEF/TE from the bronchial challenge tests were compared with the corresponding baseline values, using a paired t-test for the whole group of healthy infants.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

A power law behavior of flow harmonics was found in all healthy infants between 0.13 and 10 Hz. An example is shown in Fig. 2, and the mean (±SD) power spectra of all infants at the age of 6 mo is shown in Fig. 3. The respiratory cycle variability (RCV%) was found to be significantly lower in s and kd than in TPTEF/TE (except for kd in the healthy 1-mo-old infants) (Tables 2 and 3). Longitudinal measurements of s, kd, and TPTEF/TE in healthy infants at the ages of 1, 6, and 12 mo revealed no change in s but a significant increase in kd (P < 0.0025, power > 0.8) and decrease in TPTEF/TE between 1 and 6 mo (P < 0.006, power > 0.8) (Table 2). There was an inverse relationship between kd and TPTEF/TE in healthy infants [log(TPTEF/TE) = -0.017 kd - 0.12; r2 = 0.51]. Although s and kd appeared to be weakly correlated, this relationship did not reach a statistically significant level.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 3.   Mean (±SD) of the power spectra of 10 healthy infants at the age of 6 mo (black-triangle), all exhibiting a power law with similar slope s (see also Table 2). In comparison, the mean (±SD) power spectra of the group of 10 infants with CLD (open circle ) differed from that of the healthy group of infants of a similar age.

Comparing the different groups of infants, we found only a trend difference in the mean (95% confidence interval) values of s between infants with a history of wheezing disorders (4.35 ± 0.19) and the 12-mo-old healthy infants (4.09 ± 0.28) (P < 0.05, power < 0.8; Table 2). There was a tendency for Vmax FRC to be lower in the group with wheezing disorders (Table 1), but it was not statistically significant (P = 0.07). However, s was highly significantly lower in infants with CLD (3.84 ± 0.31; P < 0.001, power > 0.8) compared with the group of healthy 6-mo-old infants (4.39 ± 0.16) (Figs. 3 and 4). The corresponding Vmax FRC was also statistically lower in the group with CLD [58.4 ± 47.3 (SD) ml/s] than in the normal infants at age 6 mo [299 ± 158 (SD) ml/s] (P < 0.0002, power > 0.8). The kd and TPTEF/TE were not statistically different among these diagnostic groups.


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 4.   Slope s of the CLD group (box plot: median; box: 25th and 75th; errors bars: 10th and 90th; circles: 5th and 95th percentile intervals) was significantly different from that of the 6-mo-old healthy infants, whereas s in the wheezy group was not different from that of their age-related healthy group.

Despite the fact that Vmax FRC changed by >30% in all 14 healthy infants, there was no significant change in s, kd, or TPTEF/TE after histamine challenge, indicating that tidal breathing parameters are not sensitive to acute changes in airway mechanics.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The primary finding of this study is that the flow harmonics of all cases studied follow a power law functional form. The exponent s in the power law was sensitive to disease. We also applied another shape index, kd, which was not sensitive to disease but was sensitive to maturation. Additionally, these two parameters were more robust than a widely used conventional time domain descriptor of the tidal expiratory flow waveform TPTEF/TE.

The kd and Power Law Harmonics in Tidal Flow Waveforms

The kd was originally introduced by Suki et al. (15) and extended by Zhang et al. (20). They defined kd as the square root of the relative power in the signal above the fundamental frequency. The kd characterizes the power of the harmonics in the spectrum, regardless of the distribution of the harmonics in the flow spectrum and, hence, the shape of the spectrum in the frequency domain. If, for example, the harmonic energy in the flow moves to higher frequencies because of an altered complex feedback mechanism, kd would not reflect this change as long as the total power remains the same. Thus a disadvantage of using kd is that it is insensitive to the redistribution of harmonic energy in the spectrum. On the other hand, it has the advantage that it does not assume any particular shape for the spectrum, such as a power law form.

The second method, the power law analysis, is based on the primary finding of this study. The decrease in amplitude of the harmonic components of the tidal flow waveform follows a power law, as evidenced by the linear decrease of log harmonic amplitudes vs. log frequency (e.g., Fig. 2). The negative slope of the best linear fit provides an estimate of the exponent s in the power law. The lower the amplitude of these harmonics compared with the amplitude of the fundamental frequency, the steeper the linear decrease of the harmonics and the higher the value of s. Thus s is a simple index related to the degree of the nonsinusoidal behavior present in the tidal flow signal. In our study, we found only a trend but not a significant inverse correlation between s and kd in healthy infants. However, in theory, s and kd are not completely independent. If the flow spectrum follows a perfect power law, kd can be calculated from s. Thus s better characterizes the shape of the flow waveform. Nevertheless, in situations in which the shape deviates significantly from a power law, s loses its meaning, whereas one can still use kd to characterize the complexity of the waveform. The power law behavior was characteristic of the flow signals in all infants studied, raising two fundamental questions: what determines the flow harmonics and why do the harmonics follow a power law function?

The origin of the power law behavior of the tidal flow waveform is not entirely clear. TE profiles have been reported to follow one or at most the sum of two exponents (2). The power spectrum of an exponent has an infinite number of harmonics. This spectrum on a log-log graph is flat for low frequencies and turns into a power law with a fixed exponent of s = 2 above the corner frequency corresponding to the time constant of the exponent. Our power spectra, however, do not show a broad flat portion for low frequencies. Additionally, s varies among infants, ranging from 3 to 5. Therefore, the expiratory flow and hence the passive mechanical properties of the respiratory system alone cannot produce a power law behavior of the flow spectra. It is likely that the output of the respiratory oscillator (as described in Refs. 3, 11, 12) is strongly influencing the power law behavior of the flow waveform.

Changes in s and kd may be interpreted from a system point of view as follows. If a linear system is driven with a sine wave, the output will also be sine wave but, in general, with a different amplitude and phase. Thus, whenever the output of the system that is driven by a sine wave is nonsinusoidal, internal nonlinearities in the system must also contribute to the output. With regard to the tidal flow signal, unfortunately, the input to the system that generates the flow signal is not known. It is very likely that the input to the respiratory oscillator is not an ideal sine wave. Therefore, the presence of a power law spectrum or a kd > 0 in itself may not allow us to conclude too much about the saturation type of nonlinearities in the respiratory oscillator, because the power law may be a combination of the input and the internal nonlinearities of the system. Thus a decrease in s (stronger harmonics) or an increase in kd signifies that either the input to the oscillator has become more complex or the internal nonlinearities of the system have become stronger. In either case, the internal complexity of the entire system generating the flow tidal signal would increase.

Variability of Tidal Flow Indexes

To compare breath-to-breath variability of the parameters s, kd, and TPTEF/TE, we calculated the CV of the individual values for s, kd, and TPTEF/TE from cycle to cycle in the entire tidal flow signal. We found that RCV% of s and kd was significantly lower than that of TPTEF/TE for all groups of infants except the healthy infants at 1 mo of age (kd). This indicates that the parameters s and kd are more stable indexes to characterize tidal flow waveforms than TPTEF/TE. We can think of two possible explanations. The variability of TPTEF/TE is large, partly because its calculation is based on only three points of the expiratory part of the flow signal. From a mathematical point of view, the more data points of a complex wave are analyzed, the better the waveform can be characterized. From the point of view of signal-to-noise ratio, the expected noise on TPTEF/TE based on three points from the cycle can be much higher than the noise on a parameter that is estimated via some smoothing operation, such as a linear regression. Additionally, expiratory flow waveform is influenced by both the respiratory oscillator and passive mechanics of the lung and chest wall, whereas the inspiratory portion is dominated by the respiratory oscillator. Therefore, TPTEF/TE may carry mixed information about passive mechanics and respiratory control. In contrast, the calculation of s and kd involves the whole cycle, and, hence, they will be less noisy and, therefore, more reliable. Additionally, the kd and power law indexes may be more closely related to the true physiological origin of the tidal flow waveform.

Sensitivity of Tidal Flow Indexes to Maturation

We found changes in kd and TPTEF/TE between 1 and 6 mo of age in longitudinal data sets of healthy infants but no changes in s. Theoretically, kd and s should be correlated; however, in our data set, we only found a nonsignificant trend. We believe that this might be a potential explanation for the inconsistent behavior of s and kd during maturation. An increase of kd during development would be consistent with a higher degree of nonsinusoidal behavior of the mechanism determining tidal flow waveform, maybe with increased complexity of the neurorespiratory regulatory system. Developmental changes of kd and TPTEF/TE could be explained by changes in respiratory control during the first 6 mo of life and/or by changes in lung mechanics. However, our results from the histamine challenge in these healthy infants (see below) suggest that airway mechanics do not play a major role in any of these indexes because we did not find a consistent change in any of these tidal flow parameters. Therefore, it is more likely that changes in control of breathing during development might be the dominant factor responsible for the increase in kd. Additionally, based on this interpretation, the parameters s, kd, and TPTEF/TE are likely to be dependent on sleep stage, which we did not test in this paper. Because all measurements were made during quiet sleep, sleep stage is unlikely to explain any differences.

Tidal Flow Indexes in Disease

In disease, the situation becomes more complex. We only found marginally decreased s in infants with wheezing disorders, whereas we found s to be significantly lower in infants with CLD of prematurity. It is possible that the sensitivity of s is not enough to detect changes in mild lung disease, and larger data sets are necessary to find differences. However, s was very sensitive to the presence of severe lung disease such as CLD. The Vmax FRC in the CLD group was statistically significantly lower, indicating increased lung resistance, but the infants with wheezing disorders, who were asymptomatic at the time of measurement, had only a nonsignificant trend to decreased forced expiratory flows. Another possibility is as follows. Because most of these infants with CLD still needed oxygen at the time of measurements, it could be that the patients with CLD had additional changes in control of breathing. The indexes kd and TPTEF/TE were not significantly different in these groups of infants, indicating that these parameters are less sensitive to changes in the disease state.

Tidal Flow Indexes During Histamine Challenge

None of the tidal flow indexes changed consistently after histamine challenge, despite an acute change in airway obstruction as demonstrated by the fall of Vmax FRC by >30%. In contrast, in a chronic state of lung disease, such as in CLD, the power law slope s decreased, consistent with an increased level of nonsinusoidal behavior of the mechanisms determining flow waveform. Thus either lung mechanics were much more severely affected by CLD than after bronchial challenge or CLD has a direct effect on respiratory control. Nevertheless, these data support the argument that the tidal flow waveform is strongly dependent on neurorespiratory control and adaptive mechanisms rather than lung mechanics.

Summary

We proposed two new indexes (kd, s) derived from tidal flow waveform that can characterize the complexity of the neuromechanical respiratory system by quantifying the extent to which the tidal flow waveform deviates from the shape of a single sinusoid. We found that the tidal flow harmonics follow a power law that can be characterized by a single number, the exponent s. This new index and kd also have a significantly smaller breath-to-breath variability than TPTEF/TE, making them more robust descriptors of tidal flow waveforms. We also found them to be significantly different in disease but not, however, during a histamine challenge. The results from the histamine challenge indicate that the tidal flow waveform is not significantly affected by acute changes in airway mechanics and that it is, therefore, more likely that the neurorespiratory control plays a key role in determining flow waveform. These findings are consistent with previous work (16-19). To study the separate influence of control of breathing and lung mechanics on tidal flow and these parameters, future studies should employ hypoxic and/or hypercapnic challenges. Furthermore, it is important to consider that infant breathing may also inherit long-range correlations, as it has been shown by Small et al. (14). Thus the parameters s and kd will have to be studied in a very long time series to get more information on the control of breathing in infants. Because s and kd can be determined noninvasively, they might be useful in monitoring CLD in infants, assessing therapeutic interventions, and/or giving evidence of changes in neurorespiratory control during development and growth.


    ACKNOWLEDGEMENTS

U. Frey was supported by the Swiss National Science Foundation (no. 32-51974.97), the British National Asthma Campaign, and the British Society for the Protection of Infants Life.


    FOOTNOTES

Address for reprint requests and other correspondence: U. Frey, Dept. of Paediatrics, Univ. Hospital of Berne, Inselspital 3010, Berne, Switzerland (E-mail: urs.frey{at}insel.ch).

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.

Received 25 August 2000; accepted in final form 15 June 2001.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Aston, H, Clarke J, and Silverman M. Are tidal breathing indices useful in infant bronchial challenge tests? Pediatr Pulmonol 17: 225-230, 1994[ISI][Medline].

2.   Bates, JH, Shardonofsky F, and Stewart DE. The low-frequency dependence of respiratory system resistance and elastance in normal dogs. Respir Physiol 78: 369-382, 1989[ISI][Medline].

3.   Botros, SM, and Bruce EN. Neural network implementation of a three phase model of respiratory rhythm generation. Biol Cybern 63: 143-153, 1990[ISI][Medline].

4.   Clarke, JR, Aston H, and Silverman M. Evaluation of tidal expiratory flow index in healthy and diseased infants. Pediatr Pulmonol 17: 285-290, 1994[ISI][Medline].

5.   Clarke, JR, Salmon B, and Silverman M. Bronchial responsiveness in the neonatal period as a risk factor for wheezing in infancy. Am J Respir Crit Care Med 151: 1434-1440, 1995[Abstract].

6.   Duffy, P, Spriet L, Bryan MH, and Bryan AC. Respiratory inductive plethysmography (Respitrace): an evaluation of its use in the infant. Am Rev Respir Dis 123: 542-546, 1981[ISI][Medline].

7.   Frey, U, Silverman M, Barabasi AL, and Suki B. Irregularities and power law distribution in the breathing pattern of preterm and term infants. J Appl Physiol 86: 789-797, 1998[Abstract/Free Full Text].

8.   Milic-Emili, J, and Zin WA. Relationship between neuromuscular respiratory drive and ventilatory output. In: Handbook of Physiology. The Respiratory System. Mechanics of Breathing. Bethesda, MD: Am. Physiol. Soc, 1986, sect. 3, vol. III, pt. 2, chapt. 35, p. 631-646.

9.   Montroll, EW, and Schlesinger MF. On 1/f noise and distributions with long tails. Proc Natl Acad Sci USA 79: 3380-3383, 1982[Abstract/Free Full Text].

10.   Morris, MJ, and Lane DJ. Tidal expiratory flow pattern in airflow obstruction. Thorax 36: 135-142, 1981[Abstract].

11.   Ogilvie, MD, Gottschalk A, Andersm K, Richter DW, and Pack AI. A network model of respiratory rythmogenesis. Am J Physiol Regulatory Integrative Comp Physiol 263: R962-R975, 1992[Abstract/Free Full Text].

12.   Richter, DW, and Ballantyne D. A three phase theory about the basic respiratory pattern generator. In: Central Neuron Environment, edited by Schlafke ME, Koepchen HP, and See WR.. Berlin: Springer Verlag, 1983, p. 165-174.

13.   Sackner, MA, Watson H, Belsito AS, Feinerman D, Suarez M, Gonzalez G, Bisouski F, and Krieger B. Calibration of respiratory induced plethysmograph during natural breathing. J Appl Physiol 66: 410-420, 1989[Abstract/Free Full Text].

14.   Small, M, Judd K, Lowe M, and Stick SJ. Is breathing in infants chaotic? Dimension estimates for respiratory patterns during quiet sleep. J Appl Physiol 86: 359-376, 1999[Abstract/Free Full Text].

15.   Suki, B, Hantos Z, Daroczy B, Alkaysi G, and Nagy S. Nonlinearity and harmonic distortion of dog lungs measured by low frequency forced oscillations. J Appl Physiol 71: 69-75, 1991[Abstract/Free Full Text].

16.   Thomas, DA, Poole K, McArdle EK, Beardsmore CS, Goodenough PC, Thomson J, and Simpson H. The effect of sleep deprivation on sleep states, breathing events, peripheral chemoresponsiveness and arousal propensity in healthy 3 month old infants. Eur Respir J 9: 932-938, 1996[Abstract].

17.   Thomas, DA, Poole K, McArdle EK, Beardsmore CS, Goodenough PC, Thomson J, and Simpson H. Peripheral chemoresponses of infants measured by a minimally invasive method utilizing two-breath alternations in FIO2. Eur Respir J 9: 1261-1268, 1996[Abstract].

18.   Ueda, Y, Stick SM, Hall G, and Sly PD. Control of breathing in infants born to smoking mothers. J Pediatr 135: 226-232, 1999[ISI][Medline].

19.   Van der Ent, CK, van der Grinten CP, Meessen NE, Luijendijk SC, Mulder PG, and Bogaard JM. Time to peak tidal expiratory flow and the neuromuscular control of expiration. Eur Respir J 12: 646-652, 1998[Abstract].

20.   Zhang, Q, Suki B, and Lutchen KR. Harmonic distortion from nonlinear systems with broadband applications to lung mechanics. Ann Biomed Eng 23: 672-681, 1995[ISI][Medline].


J APPL PHYSIOL 91(4):1687-1693
8750-7587/01 $5.00 Copyright © 2001 the American Physiological Society



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
R. H. Habib, K. H. Pyon, S. E. Courtney, and Z. H. Aghai
Spectral characteristics of airway opening and chest wall tidal flows in spontaneously breathing preterm infants
J Appl Physiol, May 1, 2003; 94(5): 1933 - 1940.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
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 Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Frey, U.
Right arrow Articles by Suki, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Frey, U.
Right arrow Articles by Suki, B.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online