Skip to main content

Forced oscillations and respiratory system modeling in adults with cystic fibrosis

Abstract

Background

The Forced Oscillation Technique (FOT) has the potential to increase our knowledge about the biomechanical changes that occur in Cystic Fibrosis (CF). Thus, the aims of this study were to investigate changes in the resistive and reactive properties of the respiratory systems of adults with CF.

Methods

The study was conducted in a group of 27 adults with CF over 18 years old and a control group of 23 healthy individuals, both of which were assessed by the FOT, plethysmography and spirometry. An equivalent electrical circuit model was also used to quantify biomechanical changes and to gain physiological insight.

Results and discussion

The CF adults presented an increased total respiratory resistance (p < 0.0001), increased resistance curve slope (p < 0.0006) and reduced dynamic compliance (p < 0.0001). In close agreement with the physiology of CF, the model analysis showed increased peripheral resistance (p < 0.0005) and reduced compliance (p < 0.0004) and inertance (p < 0.005). Significant reasonable to good correlations were observed between the resistive parameters and spirometric and plethysmographic indexes. Similar associations were observed for the reactive parameters. Peripheral resistance, obtained by the model analysis, presented reasonable (R = 0.35) to good (R = 0.64) relationships with plethysmographic parameters.

Conclusions

The FOT adequately assessed the biomechanical changes associated with CF. The model used provides sensitive indicators of lung function and has the capacity to differentiate between obstructed and non-obstructed airway conditions. The FOT shows great potential for the clinical assessment of respiratory mechanics in adults with CF.

Introduction

Pulmonary disease, with a mortality rate of 85%, is the major cause of morbidity and mortality in individuals with Cystic Fibrosis (CF) [1]. CF care has dramatically improved over the last decades. In 2011, more than 48% of CF patients in the Patient Registry of the CF Foundation were adults, and this number continues to grow [1]. The emergence of a larger population of adults living with CF necessitates further understanding of the overall changes in their lung function. Such studies will provide essential information for clinical care and research. Understanding the mechanisms of lung injury may guide choices in clinical care and the development of new therapies.

The forced expiratory volume in 1 second (FEV1) has been used as a predictor of survival and has been widely adopted in clinical studies of lung function in children, adolescents and adults with CF [1]. However, changes in respiratory mechanics are not always detected by this test [2]. People of all ages with CF have some lung damage, even when their lung function is normal as evaluated by the FEV1, and this measurement has the limitation of being insensitive to early changes in lung disease [1].

The Forced Oscillation Technique (FOT) offers a simple and detailed approach to investigate the mechanical properties of the respiratory system and represents the current state-of-the-art assessment of lung function [3,4]. This method consists of the application of sinusoidal signals during normal respiration by means of an external pressure generator, allowing for the measurement of respiratory system impedance. The main advantage of the FOT is the simplicity of the test administration, which requires little cooperation from the patient and involves the production of parameters describing the resistive and reactive properties of the respiratory system that are complementary to traditional pulmonary assessment methods [3-5]. These characteristic were important to explore lung function in children with CF [6-11]. Recently, this technique was successfully applied in our laboratory to study the early respiratory changes in smokers [12] and patients with sarcoidosis [13] and systemic sclerosis [14]. Compartmental models of the multi-frequency impedance allowed a detailed description of the human respiratory system by equivalent electrical circuits [15,16]. These models allow us to gain additional insight into the anatomical or pathophysiological changes that occur in respiratory diseases. In addition to being useful in furthering our understanding of respiratory biomechanical function, these model parameters could improve the detection, diagnosis, and treatment of different respiratory diseases.

Therefore, the FOT has the potential to increase our knowledge about the biomechanical abnormalities of adults with CF. However, as noted in a recent review, there are no reports in the literature evaluating the FOT in adults with CF, and further research is needed in this area before the FOT can be widely used in clinical practice [17].

In this context, the purpose of the present study was to analyze the changes in the resistive and reactive biomechanical properties of the respiratory system in adults with CF.

Materials and methods

Patients and study design

This study utilized a cross-sectional design involving two groups of subjects: a group of CF adults aged 19–43 years and a group of controls. The diagnosis of CF was based on the presence of at least two of the following criteria: sweat chloride concentration >60 mEq/mL; two clinical features consistent with CF [18]; or genetic testing demonstrating two mutations associated with CF [19].

The CF group was composed of patients who had no reports of smoking, hemoptysis and/or recent hemopysis, thoracic surgery, cold or flu in the last month, or recent worsening respiratory infection (within 2 weeks of the exams). The control group was composed of individuals who had no reports of flu or cold during the last month, respiratory symptoms, lung disease, cognitive impairment, heart disease, history of smoking, hospitalization in the past six months or spirometric abnormalities. Demographic data, including age, height and weight, were obtained from each subject at the time of the procedures.

This study was approved by the Medical Research Ethics Committee of the State University of Rio de Janeiro. The work was performed in accordance with The Declaration of Helsinki. Informed consent was obtained from all volunteers before inclusion in the study.

Spirometry

A closed circuit spirometer (Warren E. Collins, Inc., Braintree, MA, USA) was used to measure the forced vital capacity (FVC), FEV1, FEV1/FVC ratio, expiratory flow between 25% and 75% of the FVC (FEF25%-75%) and FEF25%-75%/FVC ratio. These parameters were expressed as absolute values and as a percentage of the predicted values (% of predicted), and the reference values were obtained from the equations of Knudson et al. [20]. These measurements were conducted according to the recommendations of the American Thoracic Society/European Respiratory Society [21]. An obstructive ventilatory defect was defined as an FEV1/FVC ratio below the normal lower limit, as recommended in previous guidelines [21,22].

Plethysmography

Using a constant volume and variable pressure plethysmograph (Warren E. Collins, Inc., Braintree, MA, USA), we evaluated the total lung capacity (TLC), functional residual capacity (FRC) and residual volume (RV), as well as their relationships (RV/TLC and FRC/TLC). The airway resistance (Raw) and specific airway conductance (SGaw) were also measured [23]. The reference values were obtained from Goldman and Becklake [24].

Forced Oscillation

The respiratory system impedance was measured between 2 and 32 Hz with 2-Hz increments using a pseudorandom noise forced oscillation system that was built in our laboratory [12,14,25]. Pressure oscillations, which were produced by a loudspeaker coupled to the respiratory system by means of a mouthpiece, were applied with amplitude of approximately 1 cmH2O. The resulting flow (V’) and pressure (P) signals were measured near the mouth by means of a pneumotachograph and a pressure transducer, respectively. The common mode rejection ratio of the airflow measurement system (65.2 dB) was adequate for these studies [26]. The frequency response of the pressure measurement system is flat in the used frequency range. The small errors associated to the attenuation and phase shift introduced by the flow measurement system was corrected digitally using calibration curves obtained by means of the analysis of a reference mechanical load [27].

The respiratory system impedance (Zrs(f)) was estimated from the Fourier analysis of these signals, as described by equation (1).

$$ {\mathrm{Z}}_{\mathrm{rs}}\left(\mathrm{f}\right) = \frac{\overline{{\mathrm{G}}_{\mathrm{PV}}}}{\overline{{\mathrm{G}}_{\mathrm{VV}}}} $$
(1)

Where \( \overline{{\mathrm{G}}_{\mathrm{PV}}} \) and \( \overline{{\mathrm{G}}_{\mathrm{VV}}} \) are the mean crosspectrum between pressure and flow and mean autospectrum of flow, respectively. Three FOT exams were performed; each exam lasted approximately 16 seconds, and the exams were separated by one-minute intervals. During the exams, each individual wore a nose clip to prevent air escape and breathed calmly through a silicone mouthpiece. The subjects remained seated with their chests straight and their heads in a neutral position relative to the device. Additionally, the individuals firmly supported their extrathoracic airways to minimize the potential confounding effects of airway wall flow shunting [3-5]. To express the amount of noise on the pressure and flow signal, a coherence function was calculated at each investigated frequency. Only exams with a coherence function between 0.9 and 1.0 were retained [14,28]. The final result of the tests was estimated by calculating the average of three satisfactory exams.

Linear regression analysis of the resistance values in the frequency range between 4 and 16 Hz was used to determine the intercept resistance (R0) and the slope of the resistive component of the impedance (S). R0 is associated with the airway and tissue Newtonian resistance plus the delayed airway resistance resulting from gas redistribution [14]. S reflects both spatial and temporal ventilatory non-homogeneity [28,29]. The value of Rrs measured in 4 Hz (R4) was also analyzed.

The reactance in the frequency range between 4 and 32 Hz was used to measure the resonant frequency (fr), the respiratory system mean reactance (Xm) [30], and the area under the negative part of the reactance curve (Ax) [15]. Fr is defined as the lowest frequency at which the reactance crossed 0 from negative to positive. Fr reflects the airway heterogeneity changes and tissue changes. Xm is usually related to respiratory system non-homogeneity [31]. Ax is related to respiratory compliance, and reflects changes in the degree of peripheral airway obstruction [15]. The dynamic compliance of the respiratory system (Cdyn) was estimated using the respiratory reactance at 4 Hz (X4Hz) and the equation Cdyn = −1/(2πfX4Hz) [32]. The same frequency was used to evaluate the absolute value of the respiratory impedance (Zrs). This parameter is associated with the work performed by the respiratory muscles to overcome resistive and elastic loads, promoting the movement of air in the respiratory system [14,33].

Respiratory impedance model

Compartmental model analysis was performed using the extended RIC (eRIC) model (Figure 1), in which R is analogous to central airway resistance, Rp describes peripheral resistance, I is associated with lung inertance and C is associated with alveolar compliance [16]. This model is proposed as an improvement to the basic RIC model [4,15]. Specifically, the added peripheral resistance Rp allows for observation of the frequency dependence of typical real impedance data, which is beyond the capability of the RIC model. This additional component describes the resistance presented by the small airways of the respiratory system. We also evaluated the total resistance (Rt = R + Rp), which includes the effects of the central and peripheral airways. Model parameters were estimated using the Levenberg-Marquardt algorithm to determine the set of coefficients of the nonlinear model that best represents the input data set in the least squares sense. Along with the corresponding model estimates, this analysis also provided the evaluation of the total error value, which is an overall measure of “goodness of fit” for the model. Herein, this parameter is defined as the square root of the sum of the real and imaginary squared impedance estimation errors. The mean relative distance from the model and measured resistance and reactance values was also measured [34,35].

Figure 1
figure 1

Electrical representation of a two-compartment model used to analyze respiratory impedance. Resistance, inductance and capacitance are the analogs of mechanical resistance, inertance and compliance, respectively. R is analogous to central airway resistance and Rp describes peripheral resistance, I is associated with airway, gas and tissue inertance, while C is related with alveolar compliance. This analysis also evaluated the total resistance (Rt = R + Rp), which included the effects of central and peripheral airways.

Statistical analysis

Statistical analysis was performed using STATISTICA 5.0 (StatSoft, Inc., Tulsa, OK, USA) and Microcal Origin version 6.0 (Microcal Software, Inc. Ostend, Belgium). First, the Shapiro-Wilks W test was used to analyze whether the data were normally distributed. Student’s t test was used to analyze normally distributed data, and the Mann–Whitney U test and Wilcoxon Matched Pairs Test were used when the data were not normally distributed. The results with p values < 0.05 were considered statistically significant.

To measure the overall agreement between the variables related to the FOT, respiratory modeling, spirometry, plethysmography and diffusion, we calculated Pearson’s correlation coefficient for the whole group of studied volunteers. These correlations were classified as follows [36]:

  • Small or no correlation: between 0 and 0.25 (or −0.25);

  • Reasonable correlation: between 0.25 and 0.50 (or −0.25 to −0.50);

  • Moderate to good correlation: between 0.50 and 0.75 (or −0.50 to −0.75);

  • Very good to excellent correlation: greater than 0.75 (or −0.75).

Results

A total of 23 control subjects and 27 CF patients were included in the study. The patient characteristics are shown in Table 1. The biometric characteristics of the two studied groups were well matched, and there were no significant differences between the groups.

Table 1 Anthropometric characteristics of the studied subjects (mean ± standard deviation)

The average time of CF diagnosis was 12.4 years. Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene were distributed as follows: 56 (52.8%) patients had two severe mutations (classes I–III and VI); 28 (26.4%) patients had one severe mutation and one unknown mutation; 11 (10.4%) patients had one severe and one mild mutation (class IV–V); and 11 (10.4%) patients had no known mutations on both alleles. In 5 (4.7%) patients, no mutations were identified, and in 6 (5.7%) patients, mutation-analysis was not available at the time of study inclusion; however, all 11 of these patients had previous repeat pilocarpine-iontophoresis sweat tests with pathological results to establish their CF diagnoses. In this sample, 13 (48.1%) patients were chronically infected by Pseudomonas aeruginosa; 2 (7.4%) patients were chronically infected by the Burkholderia cepacia complex; 4 (14.8%) patients were chronically infected by both P. aeruginosa and B. cepacia complexes; and 6 (22.2%) patients had no chronic pulmonary infection. The prevalence of pancreatic insufficiency and diabetes mellitus was 63% and 11.1%, respectively.

A significant reduction was observed in all of the studied spirometric parameters in the CF group (Table 2). An obstructive ventilatory defect was diagnosed in 23 cases (85.2%). Regarding the plethysmographic parameters, no patient had a TLC below the normal lower limit. There were no significant differences in the TLC and FRC, whereas the RV significantly increased. The Raw was significantly increased in CF patients (p < 0.0001), whereas the SGaw was significantly reduced (p < 0.0001).

Table 2 Spirometric, plethysmografic and diffusion characteristics of the studied subjects

Figure 2 shows the mean curves of the respiratory resistance (Rrs) and reactance (Xrs) as functions of the frequencies in normal and CF subjects. The Rrs curve in CF adults was significantly different from the control curve (Figure 2A; p < 0.0001), whereas the Xrs curves displayed a smaller difference between the groups (Figure 2B; p < 0.03).

Figure 2
figure 2

Mean course of respiratory resistance (A) and reactance (B) as a function of frequency in the normal subjects and in adult patients with cystic fibrosis.

Adults with CF presented an increased R0 (Figure 3A; p < 0.0001) and R4 (p < 0.0001), whereas S was significantly more negative (p < 0.0006). The presence of CF was associated with higher fr values (p < 0.0002), more negative Xm values (p < 0.0001), more positive Ax (p < 0.0003) and lower Cdyn values (p < 0.0001). In contrast, the Zrs values were higher in CF patients (p < 0.0001).

Figure 3
figure 3

Comparative analysis of the resistive parameters obtained from groups of adults with Cystic Fibrosis and control: Respiratory system resistance (R0; Figure A), resistance in 4 Hz (R4; Figure B), slope of respiratory resistance (S; Figure C). Similar analysis for reactive parameters are also described: resonant frequency (fr; Figure D), mean reactance (Xm; Figure E), reactance area (Ax; Figure F), dynamic compliance (Cdyn; Figure G) and the impedance modulus (Zrs; Figure H). The top and the bottom of the box plot represent the 25th- to 75th-percentile values, while the circle represents the mean value, and the bar across the box represents the 50th-percentile value. The whiskers outside the box represent the 10th-to 90th-percentile values.

The changes in the parameters of the eRIC model associated with CF are described in Figure 4. The mean total error in the model estimates was 0.13 cmH2O/L/s and 0.10 cmH2O/L/s in the control group and FC group, respectively. The relative distance from the model and measured resistance and reactance values was 3.5 ± 1.1% in controls and 2.1 ± 0.7% in patients. There was no significant change in the R values, whereas the Rp significantly increased in CF (Figures 4A and B, respectively). A significant increase in the Rt occurred in CF (Figure 4C). Significant reductions in I and C values were observed (Figures 4D and 4E, respectively).

Figure 4
figure 4

Influence of cystic fibrosis in adult patients on parameter values estimated from the model described in Figure 1 : central airway resistance (R; Figure A), peripheral resistance (Rp; Figure B), total resistance (Rt; Figure C), lung inertance (I; Figure D) and alveolar compliance (C; Figure E). The top and the bottom of the box plot represent the 25th- to 75th-percentile values, while the circle represents the mean value, and the bar across the box represents the 50th-percentile value. The whiskers outside the box represent the 10th-to 90th-percentile values.

In general, moderate to good relationships between the resistive FOT parameters (R0, R4 and S) and spirometric indexes (Table 3) were observed. The reactive parameters (fr, Ax and Xm), presented better relationships; specifically, the Cdyn presented excellent correlations with the FVC (L) and FEV1 (L) (Table 3). Considering the associations obtained from compartmental model analysis, it was observed that R was not related with the spirometric parameters, whereas Rp and Rt presented moderate to good relationships. On the other hand, C showed a very good direct correlation with the FVC (L) and FEV1 (L).

Table 3 Correlation coefficient and significance level of the analysis including FOT, compartmental model and spirometric exams

Significant reasonable to good correlations were observed between the resistive parameters (R0, R4 and S) and plethysmographic indexes (Table 4). Similar associations were observed for the reactive parameters (fr, Xm, Ax, Cdyn and Zrs). The strongest relationships were observed for Xm and Zrs with RV/TLC (R = −0.69; p < 0.0001). Via the model analysis, it was observed that R was not correlated with the plethysmographic parameters, whereas Rp presented reasonable (R = 0.35) to good (R = 0.64) relationships with the parameters. C showed a good inverse correlation with the RV/TLC (R = −0.65).

Table 4 Correlation coefficient and significance level of the analysis including FOT, compartmental model and plethysmographic exams

Discussion

This study is the first to systematically evaluate the use of the forced oscillation technique in adults with CF. Two major findings were observed: 1) adult patients with CF presented increased Rrs and more negative Xrs values, and 2) the FOT parameter alterations were correlated with the standard pulmonary function analysis methods.

In agreement with other studies, we observed that the majority of adults with CF had obstructive airway disease [17,37]. In CF, several mechanisms contribute to airflow limitation, including chronic inflammation, increased smooth muscle in the airway walls, mucus hypersecretion with the formation of viscous secretions, and bronchiectasis [38,39]. The average FEV1 in our patients was 65.8%, similar to that observed by Yankaskas et al. [40] in an American, adult CF population. Compared with the controls, the CF patients in the present study demonstrated a significant increase in both the RV and the RV/TLC ratio. Interestingly, the clinical status of adults with CF is negatively correlated with the amount of air trapping in poorly ventilated areas of the lungs [37]. Thus, patients with worse clinical conditions tend to have greater air trapping.

In contrast to the relatively uniform resistance values in the different frequencies observed in healthy subjects, the resistance was increased at lower frequencies and decreased with increasing frequency applied via the FOT in CF patients (Figure 2). Thus, the negative dependence of the resistance on the frequency in patients with airway obstruction may have resulted from a redistribution of alveolar gas due to ventilation heterogeneity or changes in airway compliance [28]. In the present study, the abnormalities in the resistance curve may be related to airway obstruction, which was detected in 85.2% of adults with CF.

Resistances measured between 4–32 Hz are related to the airway and tissue Newtonian resistance plus the delayed airway resistance resulting from gas redistribution. R0 is an extrapolation to the intercept, estimating how the cited properties work at low frequencies. Thus, this parameter, as well as R4, does not include tissue viscoelastic properties, but it is important to evaluate how low-frequency Newtonian and delayed airway resistance changes in disease. R0 is related to the total resistance of the respiratory system, and can be proposed as an index to evaluate the degree of obstruction and to assess the reversibility of the airways [28]. Because there are no other studies of the FOT in adults with CF, it is reasonable to compare our results with those of other studies conducted in patients with airway obstruction due to other causes. The average R0 value (Figure 3A; 4.71 ± 1.52 cmH2O/L/s) was similar to those obtained in adults with asthma and patients with COPD with moderate airway obstruction [28,41,42]. R0 and R4 also include the effect of the central airway resistance. Thus, the significant increase of these parameters may also include the involvement of large airways in adults with CF. Despite the initial damage that occurs in the small airways, the central and hilar bronchi are always involved when bronchiectasis is present [43]. The alterations of R0 and R4 were inversely correlated with forced expiratory flows (Table 3) and directly related with pulmonary volumes and airway obstruction (Table 4). This allowed us to infer that these parameters represent alterations of the central and peripheral airways and are also sensitive to air trapping (associated with the RV and the RV/TLC ratio), which is a typical abnormality found in CF. It is interesting to note that, while the associations of R4 were higher than that observed in R0 with spirometric parameters (Table 3), the associations of R0 were higher than that observed for R4 with plethysmographic values (Table 4).

In adults with CF, the average S value was significantly more negative than in healthy subjects (Figure 3C). This finding is in agreement with a previous study of asthma that also evaluated patients with airflow limitation [41]. S is associated with the non-homogeneity of the respiratory system [28]. Among the causes that can lead to non-homogeneity in the respiratory system in CF, the most important cause is the uneven distribution of pulmonary lesions. In CF patients, there is a preference for regional involvement, which is most pronounced in the upper lung zones [43]. In agreement with this interpretation, Table 3 shows that S was associated with the spirometric parameters, mainly those sensitive to both airway closure and airway narrowing (FEV1), and that representing alterations of the central (FEV1/FVC) and the small (FEF25–75%) airways. Additional support for this interpretation is shown in Table 4, in which S is shown to be correlated with residual volume abnormalities and airway obstruction.

Figure 2B shows that CF in adult patients introduced more negative Xrs values. According to Clement et al. [44], the Xrs may contribute to the discrimination between healthy individuals and patients with respiratory symptoms. As a result of the Xrs changes, adults with CF had higher fr values (Figure 3D) and more negative Xm values (Figure 3E). Similarly, Hellinckx et al. [45] also observed more negative reactance values in children with CF. Whereas S is associated with non-homogeneity in terms of resistance distribution, Xm describes non-homogeneity in terms of the reactive properties of the respiratory system. Ax (Figure 3F) presented smaller discrimination than that presented by Xm and Cdyn (p < 0.0001). This parameter also presented smaller association with spirometric (Table 3) and plethysmographic (Table 4) parameters than Xm. The fr value may be a sensitive marker of incipient airway obstruction [25]. In line with this hypothesis, an inverse association was observed with the spirometric parameters (Table 3), as well as direct relationships with the residual volume and airway obstruction (Table 4). Xm was also associated with the spirometric parameters (Table 3), the residual volume and airway obstruction (Table 4). The association between Xm and RV/TLC is interesting, since this ratio is supposed to be of prognostic value. Higher fr and more negative Xm values in patients with airway obstruction may also be related to reduced lung compliance, which results in a predominance of the negative phase of the reactance curve, as can be observed in Figure 2B.

The viscoelastic behavior of lung tissue dominate over the low-frequency range (below about 2 Hz) [4]. For frequencies below 5 Hz, the effect of the inductance is negligible [46,47]. In the medium-frequency range (2–40 Hz), reactance transitions at the resonance frequency from dominance by the tissue elastic properties at low frequencies to dominance by the inertial properties of the gas in the airways at higher frequencies [4]. In this study, the Cdyn was measured at a frequency of 4 Hz, at which the impedance mainly reflects the elastic properties. Hyatt et al. [23] suggested that the reduced Cdyn values may be associated with decreased lung compliance or increased airway resistance. Considering that dynamic lung compliance is a major component of the Cdyn, the most likely causes for their reduction in CF patients (Figure 3G) include a reduced number of functioning small airways, increased airflow resistance, airway inflammation, increased secretions and parenchymal inflammation/fibrosis [17]. This interpretation is supported by the associations observed with the spirometric (Table 3) and plethysmographic (Table 4) indexes of airway obstruction.

The respiratory impedance module is related to the total mechanical load of the respiratory system and describes the effects of the respiratory system resistance and compliance. This parameter is related to the work performed by the respiratory muscles to overcome resistive and elastic loads to promote the movement of air in the respiratory system [33]. It is also related to fatigue and breathlessness, which are among the most important symptoms in predicting the quality of life in respiratory patients. The increase in Zrs values observed in CF (Figure 3H) may be due to several factors, including the increase in respiratory resistance, associated with the reduction of airway caliber, and the reduced dynamic compliance. The associations observed with airway obstruction (Tables 3 and 4) and residual volume (Table 4) are consistent with this interpretation. Similar to Xm, the correlation observed between Zrs and RV/TLC is clinically interesting, since this ratio is supposed to be of prognostic value.

The compartmental model analysis may contribute to a better understanding of the biomechanical abnormalities of CF in adults. Figure 4 shows that CF does not introduce significant increases in the central resistance values (Figure 4A), whereas both the peripheral (Figure 4B) and total resistance (Figure 4C) values significantly increased.

R was not associated with spirometric or plethysmographic parameters (Tables 3 and 4), whereas Rp presented a significant inverse association with the forced expiratory flow (Table 3), and direct relationships with the Raw and the RV (Table 4). Rt presented similar associations, with slightly reduced values. These results are in line with the pathophysiology of CF [17], providing additional evidence that the harmful effects of CF-related dysfunctions are particularly manifested in the distal airways [17,19,45]. These small airways, by nature of their size and architecture, are particularly susceptible to the effects of airway wall inflammation and obstruction [17,19].

Respiratory inertance mainly reflects the mass of gas that is moved during spontaneous ventilation. Since this parameter is proportional to the ratio of airway length to diameter [48], it might be expected that the CF patients would have a higher inertance. However, the estimates of respiratory inertance using the eRIC model showed reduced values in CF patients (Figure 4D). These results are consistent with Diong et al. [16], which observed a decrease in inertance in patients diagnosed with mild obstructive lung disease. Similar results were also obtained by Michaelson et al. [49] and Hayes et al. [50] in obstructive patients, by Ionescu and de Keyser [51] using several different models to evaluate COPD patients and by Iwatsubo et al. [52] in subjects displaying moderate airway obstruction. It is important to point out that Iwatsubo et al. [52] used the standard and head generator methods, which indicates that the reduction of inertance is not associated with upper airway shunt. The reasons for this behavior may be interpreted using the concepts of choke points [53] and apparent compliance [54]. Under normal conditions, inertance reflects the inertial properties of the entire respiratory system. However, when obstruction is present, the oscillatory signal cannot pass through these choke points and reach the alveoli, which prevent the FOT from being able to probe the lung beyond the choke point where obstruction occurs. This way, FOT measurements reflects the mechanical properties of airways proximal to the choke points. This result in a reduction of the apparent mass of the gas accelerated by the FOT signal, and consequently a reduction in the pressure required for gas acceleration, and in the measured inertance. Thus, we may speculate that the airway obstruction may result in an apparent inertance, similar to the process observed in the apparent compliance. Theoretically, these changes would be associated with retained air within the obstructed lung units, which would increase RV and FRC. This hypothesis was confirmed by the inverse association observed between inertance and the FRC and RV/TLC (Table 4). It is also consistent with the direct association observed between inertance and spirometric indices of obstruction (Table 3).

As can be observed in Figure 4E, the analysis based on the eRIC model also revealed reduced values of compliance. These results clearly reflect the CF physiopathology, which includes airway inflammation, increased secretions and parenchymal inflammation/fibrosis [17]. This accurate description of the CF physiopathology afforded by the Cdyn and C allows us to infer that the FOT might be an alternative method to quantitatively and noninvasively assess the elastic properties of the respiratory system in these individuals, thus avoiding the use of invasive techniques such as those using an esophageal balloon.

One may argue about the redundancy of the evaluation of Cdyn and C. Note that these parameters are estimated by two very different methods. The comparison of the parameters obtained by these two different approaches help to understand what is the better way to evaluate the elastic properties of the respiratory system. The performance of these parameters, described in Tables 3 and 4, indicates that they present very similar results. Considering that the estimative of Cdyn is simpler, do not demanding the computational effort to fit a model, it seems more adequate than C.

Similarly to other functional assessment techniques, the limitations and consequences of the FOT must be recognized. The FOT is subject to the influence of upper airway shunting; this must be carefully controlled for when using this method. The resulting effect is a reduced impedance measurement relative to its actual value. This effect becomes progressively stronger as the respiratory resistance increases, as is the case for highly obstructive patients. Because we are studying patients with relatively low resistance values (Figures 2, 3 and 4) compared with typical obstructive patients [41,42], this limitation did not present a notable problem in this study. In addition, the shunt effect was minimized by asking the patients to firmly support their cheeks and mouth floors [3,4]. A second limitation arises from the process of spontaneous breathing, which introduces both random and systematic errors [3,4]. These errors were minimized in the present study using excitation frequencies (4–32 Hz) at least 20 times higher than those present in the spontaneous ventilation process (approximately 0.2 Hz). These errors can be easily evaluated using the coherence function (γ2) between the pressure and airflow signals [3,4]. In this study, any time the coherence value computed for any of the studied frequencies was smaller than a minimal threshold (0.90), the maneuver was not considered valid, and the examination was repeated [33,41,42].

The multiple breath inert gas washout (MBW) test assesses efficiency of gas mixing processes by following the expiration of inert tracer gas from the lungs during tidal breathing [17]. Since there is evidence that this method may help to assess the uniformity of ventilation distribution and to detect early lung disease, there has been considerable investigation as to the role of MBW in CF management [17,55-57]. The use of MBW as a gold standard may contribute to provide further information concerning the use of FOT in adults with CF. This hypothesis warrants further study.

Fitting mathematical models may help to obtain parameters describing lung structure and function. However, there are numerous models, and different models may be fitted by the same impedance data [4]. It is important to point out that more complicated models would not allow statistically reliable parameter estimates [58,59]. Thus it is not possible to include every known component of a complicated system in a mathematical model. Therefore, the choice of model structure requires a decision about which components of the real system are important for the purpose at hand, and which components can be safely ignored [60]. To minimize the problem of model complexity, a simple four-element compartmental model was used for interpreting forced oscillation measurements [16]. However, we believe that studies focusing on fractional order impedance models [46] and more complex models associated with a wider frequency range [4,46,48] could contribute with new information concerning the effects of CF in adults and should therefore be addressed in future research.

Partial correlation coefficients and multiple step wise regression analysis may help to further clarify the physiological meaning and the association between each of the studied parameters and the physiological changes. We believe that these analyses in a higher number of subjects deserve further studies.

A critical analysis of the limitations of the present study is also necessary. First, the present report was a cross-sectional study of clinically stable patients. Future studies should incorporate the analysis of the FOT relative to the long-term results and their role in exacerbations. Second, the volunteers were recruited at one center, which might hinder the generalizability of the results.

In conclusion, in this study, we present evidence of increased total respiratory resistance and reduced ventilation homogeneity and dynamic compliance in adults with CF. The eRIC model analysis showed abnormally increased peripheral resistance values. This analysis contributed to the elucidation of the biomechanical changes in these patients. The FOT associated with the eRIC model is a non-invasive, simple and radiation-free method for the detection of biomechanical abnormalities. Thus, this method has promising applications as a complementary pulmonary function test in adult CF patients.

References

  1. Foundation CF. Cystic Fibrosis Foundation Patient Registry. Annual data report to the center directors. In: Cystic fibrosis foundation patient registry. 2013.

    Google Scholar 

  2. Coe CI, Watson A, Joyce H, Pride NB. Effects of smoking on changes in respiratory resistance with increasing age. Clin Sci. 1989;76:487–94.

    Google Scholar 

  3. Oostveen E, MacLeod D, Lorino H, Farre R, Hantos Z, Desager K, et al. The forced oscillation technique in clinical practice: methodology, recommendations and future developments. Eur Respir J. 2003;22:1026–41.

    Article  Google Scholar 

  4. Bates JH, Irvin CG, Farre R, Hantos Z. Oscillation mechanics of the respiratory system. Compr Physiol. 2011;1:1233–72.

    Google Scholar 

  5. Kaczka DW, Dellaca RL. Oscillation mechanics of the respiratory system: applications to lung disease. Crit Rev Biomed Eng. 2011;39:337–59.

    Article  Google Scholar 

  6. Ionescu CM, Desager K, Vandersteen G, De Keyser R. Respiratory mechanics in children with cystic fibrosis. Biomed Signal Proces. 2014;11:74–9.

    Article  Google Scholar 

  7. Gangell CL, Horak F, Patterson HJ, Sly PD, Stick SM, Hall GL. Respiratory impedance in children with cystic fibrosis using forced oscillations in clinic. Eur Respir J. 2007;30:892–7.

    Article  Google Scholar 

  8. Lebecque P, Stanescu D. Respiratory resistance by the forced oscillation technique in asthmatic children and cystic fibrosis patients. Eur Respir J. 1997;10:891–5.

    Google Scholar 

  9. Hellinckx J, Demedts M. No paradoxical bronchodilator response with forced oscillation technique in children with cystic fibrosis. Chest. 1998;113:55–9.

    Article  Google Scholar 

  10. Pandit C, Valentin R, De Lima J, Robinson P, Fitzgerald D, van Asperen P, et al. Effect of general anesthesia on pulmonary function and clinical status on children with cystic fibrosis. Pediatr Anesth. 2014;24:164–9.

    Article  Google Scholar 

  11. Brennan S, Hall GL, Horak F, Moeller A, Pitrez PMC, Franzmann A, et al. Correlation of forced oscillation technique in preschool children with cystic fibrosis with pulmonary inflammation. Thorax. 2005;60:159–63.

    Article  Google Scholar 

  12. Faria AC, Lopes AJ, Jansen JM, Melo PL. Evaluating the forced oscillation technique in the detection of early smoking-induced respiratory changes. Biomed Eng Online. 2009;8:22.

    Article  Google Scholar 

  13. Faria AC, Lopes AJ, Jansen JM, Melo PL. Assessment of respiratory mechanics in patients with sarcoidosis using forced oscillation: correlations with spirometric and volumetric measurements and diagnostic accuracy. Respiration. 2009;78:93–104.

    Article  Google Scholar 

  14. Miranda IA, Dias Faria AC, Lopes AJ, Jansen JM, Lopes de Melo P. On the respiratory mechanics measured by forced oscillation technique in patients with systemic sclerosis. PLoS One. 2013;8:e61657.

    Article  Google Scholar 

  15. Meraz EG, Nazeran H, Ramos CD, Nava P, Diong B, Goldman MD, et al. Analysis of impulse oscillometric measures of lung function and respiratory system model parameters in small airway-impaired and healthy children over a 2-year period. Biomed Eng Online. 2011;10:21.

    Article  Google Scholar 

  16. Diong B, Nazeran H, Nava P, Goldman M. Modeling human respiratory impedance. Comparing the best method with the least estimation errors. IEEE Eng Med Biol Mag. 2007;26:48–55.

    Article  Google Scholar 

  17. Horsley A, Siddiqui S. Putting lung function and physiology into perspective: cystic fibrosis in adults. Respirology. 2014;20(1):33–45.

    Article  Google Scholar 

  18. Rosenstein BJ, Cutting GR. The diagnosis of cystic fibrosis: a consensus statement. Cystic fibrosis foundation consensus panel. J Pediatr. 1998;132:589–95.

    Article  Google Scholar 

  19. Farrell PM, Rosenstein BJ, White TB, Accurso FJ, Castellani C, Cutting GR, et al. Guidelines for diagnosis of cystic fibrosis in newborns through older adults: Cystic fibrosis foundation consensus report. J Pediatr. 2008;153:S4–14.

    Article  Google Scholar 

  20. Knudson RJ, Lebowitz MD, Holberg CJ, Burrows B. Changes in the normal maximal expiratory flow-volume curve with growth and aging. Am Rev Respir Dis. 1983;127:725–34.

    Google Scholar 

  21. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38.

    Article  Google Scholar 

  22. Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, Casaburi R, et al. Interpretative strategies for lung function tests. Eur Respir J. 2005;26:948–68.

    Article  Google Scholar 

  23. Hyatt RE, Scandon PD, Nakamura M. Interpretation of pulmonary function tests. Phyladelphia: Lippincott-Raven; 1997.

    Google Scholar 

  24. Goldman HI, Becklake MR. Respiratory function tests - normal values at median altitudes and the prediction of normal results. Am Rev Tuberc Pulm. 1959;79:457–67.

    Google Scholar 

  25. de Melo PL, Werneck MM, Giannella-Neto A. New impedance spectrometer for scientific and clinical studies of the respiratory system. Rev Sci Instrum. 2000;71:2867–72.

    Article  Google Scholar 

  26. Farre R, Navajas D, Peslin R, Rotger M, Duvivier C. A correction procedure for the asymmetry of differential pressure transducers in respiratory impedance measurements. Ieee T Bio-Med Eng. 1989;36:1137–40.

    Article  Google Scholar 

  27. Cauberghs M, Vandewoestijne KP. Calibration procedure of the forced oscillation technique. Eur Respir Rev. 1991;1(3):158–62.

    Google Scholar 

  28. Lorino AM, Zerah F, Mariette C, Harf A, Lorino H. Respiratory resistive impedance in obstructive patients: linear regression analysis vs viscoelastic modelling. Eur Respir J. 1997;10:150–5.

    Article  Google Scholar 

  29. Peslin R, Marchal F, Duvivier C, Ying Y, Gallina C. Evaluation of a modified head generator for respiratory impedance measurements. Eur Respir Rev. 1991;1(3):140–5.

    Google Scholar 

  30. Pride NB. Forced oscillation techniques for measuring mechanical properties of the respiratory system. Thorax. 1992;47:317–20.

    Article  Google Scholar 

  31. Ying Y, Peslin R, Duvivier C, Gallina C, Felicio da Silva J. Respiratory input and transfer mechanical impedances in patients with chronic obstructive pulmonary disease. Eur Respir J. 1990;3:1186–92.

    Google Scholar 

  32. Nagels J, Landser FJ, van der Linden L, Clement J, Van de Woestijne KP. Mechanical properties of lungs and chest wall during spontaneous breathing. J Appl Physiol. 1980;49:408–16.

    Google Scholar 

  33. MacLeod D, Birch M. Respiratory input impedance measurement: forced oscillation methods. Med Biol Eng Comput. 2001;39:505–16.

    Article  Google Scholar 

  34. Lorino AM, Atlan G, Lorino H, Zanditenas D, Harf A. Influence of posture on mechanical parameters derived from respiratory impedance. Eur Respir J. 1992;5:1118–22.

    Google Scholar 

  35. Oostveen E, Peslin R, Gallina C, Zwart A. Flow and volume dependence of respiratory mechanical-properties studied by forced oscillation. J Appl Physiol. 1989;67:2212–8.

    Google Scholar 

  36. Dawson B, Trapp RG. Basic & Clinical Biostatistics. 3rd ed. London, England: McGraw-Hill; 2001.

    Google Scholar 

  37. Andrade EF, Fonseca DLO, Abreue Silva FA, Menna-Barreto SS. Spirometry evolution assessment of cystic fibrosis. J Pneumol. 2001;27:130–1.

    Google Scholar 

  38. Chmiel JF, Konstan MW. Inflammation and anti-inflammatory therapies for cystic fibrosis. Clin Chest Med. 2007;28:331.

    Article  Google Scholar 

  39. Regamey N, Ochs M, Hilliard TN, Muhlfeld C, Cornish N, Fleming L, et al. Increased airway smooth muscle mass in children with asthma, cystic fibrosis, and non-cystic fibrosis bronchiectasis. Am J Respir Crit Care Med. 2008;177:837–43.

    Article  Google Scholar 

  40. Yankaskas JR, Marshall BC, Sufian B, Simon RH, Rodman D. Cystic fibrosis adult care - consensus conference report. Chest. 2004;125:1S–39.

    Article  Google Scholar 

  41. Cavalcanti JV, Lopes AJ, Jansen JM, Melo PL. Detection of changes in respiratory mechanics due to increasing degrees of airway obstruction in asthma by the forced oscillation technique. Respir Med. 2006;100:2207–19.

    Article  Google Scholar 

  42. Di Mango AM, Lopes AJ, Jansen JM, Melo PL. Changes in respiratory mechanics with increasing degrees of airway obstruction in COPD: detection by forced oscillation technique. Respir Med. 2006;100:399–410.

    Article  Google Scholar 

  43. Bhalla M, Turcios N, Aponte V, Jenkins M, Leitman BS, Mccauley DI, et al. Cystic-fibrosis - scoring system with thin-section Ct. Radiology. 1991;179:783–8.

    Article  Google Scholar 

  44. Clement J, Landser FJ, Vandewoestijne KP. Total resistance and reactance in patients with respiratory complaints with and without airways obstruction. Chest. 1983;83:215–20.

    Article  Google Scholar 

  45. Hellinckx J, De Boeck K, Demedts M. No paradoxical bronchodilator response with forced oscillation technique in children with cystic fibrosis. Chest. 1998;113:55–9.

    Article  Google Scholar 

  46. Ionescu CM. The human respiratory system: an analysis of the interplay between anatomy, structure, breathing and fractal dynamics. London: Springer; 2013.

    Book  Google Scholar 

  47. Hantos Z, Daroczy B, Suki B, Nagy S, Fredberg JJ. Input impedance and peripheral inhomogeneity of dog lungs. J Appl Physiol. 1992;72:168–78.

    Article  Google Scholar 

  48. Bates JHT. Lung mechanics: an inverse modeling approach. Cambridge: Cambridge University Press; 2009.

    Book  Google Scholar 

  49. Michaelson ED, Grassman ED, Peters WR. Pulmonary mechanics by spectral analysis of forced random noise. J Clin Invest. 1975;56:1210–30.

    Article  Google Scholar 

  50. Hayes DA, Pimmel RL, Fullton JM, Bromberg PA. Detection of respiratory mechanical dysfunction by forced random noise impedance parameters. Am Rev Respir Dis. 1979;120:1095–100.

    Google Scholar 

  51. Ionescu C, de Keyser R. Parametric models for characterizing respiratory input impedance. J Med Eng Technol. 2008;32:315–24.

    Article  Google Scholar 

  52. Iwatsubo Y, Lorino H, Hubert C, Duvivier C, Peslin R, Pham QT, et al. Measurement of respiratory impedance by forced oscillation - comparison of the standard and head generator methods. Eur Respir J. 1994;7:901–6.

    Google Scholar 

  53. Dellaca RL, Duffy N, Pompilio PP, Aliverti A, Koulouris NG, Pedotti A, et al. Expiratory flow limitation detected by forced oscillation and negative expiratory pressure. Eur Respir J. 2007;29:363–74.

    Article  Google Scholar 

  54. Kelly VJ, Sands SA, Harris RS, Venegas JG, Brown NJ, Stuart-Andrews CR, et al. Respiratory system reactance is an independent determinant of asthma control. J Appl Physiol. 2013;115:1360–9.

    Article  Google Scholar 

  55. Fuchs SI, Gappa M, Eder J, Unsinn KM, Steinkamp G, Ellemunter H. Tracking lung clearance index and chest CT in mild cystic fibrosis lung disease over a period of three years. Respir Med. 2014;108:865–74.

    Article  Google Scholar 

  56. Gustafsson PM, Robinson PD, Gilljam M, Lindblad A, Houltz BK. Slow and fast lung compartments in cystic fibrosis measured by nitrogen multiple-breath washout. J Appl Physiol. 2014;117:720–9.

    Article  Google Scholar 

  57. Vanderhelst E, De Meirleir L, Schuermans D, Malfroot A, Vincken W, Verbanck S. Evidence of an acinar response following treatment for exacerbation in adult patients with cystic fibrosis. Respiration. 2014;87:492–8.

    Article  Google Scholar 

  58. Lutchen KR, Jackson AC. Statistical measures of parameter estimates from models fit to respiratory impedance data: emphasis on joint variabilities. IEEE Trans Biomed Eng. 1986;33:1000–9.

    Article  Google Scholar 

  59. Lutchen KR, Jackson AC. Reliability of parameter estimates from models applied to respiratory impedance data. J Appl Physiol. 1987;62:403–13.

    Google Scholar 

  60. Bates JHT. Lung mechanics: an inverse modeling approach. United States of America: Cambridge University Press, New York; 2009.

    Book  Google Scholar 

Download references

Acknowledgements

The Brazilian Council for Scientific and Technological Development (CNPq) and the Rio de Janeiro State Research Supporting Foundation (FAPERJ) supported this study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Pedro L Melo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ANL conducted the measurements for this study, analyzed the data, and drafted the manuscript. ACDF analyzed the data and contributed to the drafting of the manuscript. AJL provided data and subject identification. JMJ mentored ACDF and participated in the data analysis process. PLM mentored ACDF, organized the study and helped to draft the manuscript. All authors have read and approved this manuscript.

Adma N Lima, Alvaro C D Faria, Agnaldo J Lopes, José M Jansen and Pedro L Melo contributed equally to this work.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lima, A.N., Faria, A.C.D., Lopes, A.J. et al. Forced oscillations and respiratory system modeling in adults with cystic fibrosis. BioMed Eng OnLine 14, 11 (2015). https://0-doi-org.brum.beds.ac.uk/10.1186/s12938-015-0007-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12938-015-0007-7

Keywords