Skip to main content

Combined forced oscillation and fractional-order modeling in patients with work-related asthma: a case–control study analyzing respiratory biomechanics and diagnostic accuracy

Abstract

Background

Fractional-order (FrOr) models have a high potential to improve pulmonary science. These models could be useful for biomechanical studies and diagnostic purposes, offering accurate models with an improved ability to describe nature. This paper evaluates the performance of the Forced Oscillation (FO) associated with integer (InOr) and FrOr models in the analysis of respiratory alterations in work-related asthma (WRA).

Methods

Sixty-two individuals were evaluated: 31 healthy and 31 with WRA with mild obstruction. Patients were analyzed pre- and post-bronchodilation. The diagnostic accuracy was evaluated using the area under the receiver operating characteristic curve (AUC). To evaluate how well do the studied models correspond to observed data, we analyzed the mean square root of the sum (MSEt) and the relative distance (Rd) of the estimated model values to the measured resistance and reactance measured values.

Results and discussion

Initially, the use of InOr and FrOr models increased our understanding of the WRA physiopathology, showing increased peripheral resistance, damping, and hysteresivity. The FrOr model (AUC = 0.970) outperformed standard FO (AUC = 0.929), as well as InOr modeling (AUC = 0.838) in the diagnosis of respiratory changes, achieving high accuracy. FrOr improved the curve fitting (MSEt = 0.156 ± 0.340; Rd = 3.026 ± 1.072) in comparison with the InOr model (MSEt = 0.367 ± 0.991; Rd = 3.363 ± 1.098). Finally, we demonstrated that bronchodilator use increased dynamic compliance, as well as reduced damping and peripheral resistance.

Conclusions

Taken together, these results show clear evidence of the utility of FO associated with fractional-order modeling in patients with WRA, improving our knowledge of the biomechanical abnormalities and the diagnostic accuracy in this disease.

Background

Asthma is an umbrella label for various conditions characterized by chronic airway and/or lung disease. This condition includes several different phenotypes and is likely to have several different underlying mechanisms [1]. It is a treatable chronic airway disease that affects all age groups and has a high prevalence, morbidity, and mortality worldwide [2, 3].

Previous studies showed that 10 to 15% of adult-onset asthma cases are directly caused by occupational factors, while another 10% result from worsening pre-existing asthma due to workplace conditions. Work-related asthma (WRA) is characterized by obstructive airways and hyperreactivity due to conditions in the workplace rather than stimuli from outside the workplace. It is divided into two categories: occupational asthma (OA), attributed to the particular causes and conditions of the work environment, and exacerbated occupational asthma (EOA), referred to as pre-existing or concurrent asthma that is aggravated by occupational exposures. OA is one of the most prevalent occupational respiratory diseases in industrialized countries [3]. Occupational and environmental pollution in the form of dust, fumes, vapors, and toxic gases are essential risk factors for this disease.

Asthma functional assessment through spirometry establishes the diagnosis, documents the severity of airflow obstruction, and monitors the course of the disease and changes resulting from treatment [1]. However, this technique requires great patient cooperation in performing respiratory maneuvers, which may limit its use in the elderly or individuals with cognitive impairment [4]. Also, forced maneuver subjects the bronchi to stress, which can alter the bronchial tone, as well as exhaustion by repetition [5]. Whole-body plethysmography allows the measurement of lung volumes, capacities, and resistances [6, 7]. However, this method demands a cooperation maneuver similar to spirometry.

The Forced Oscillation Technique (FOT), also known as Oscillometry, is a system identification method used to evaluate the respiratory system resistance and reactance [8]. In this method, a sinusoidal pressure variation is applied in the opening of the airway through a mouthpiece, overlapping spontaneous ventilation. Thus, it requires only passive cooperation and no forced expiratory maneuvers are required [9]. Due to several recent technical enhancements, FOT currently represents state of the art in lung function evaluation [10]. Several authors have argued that it has the potential to improve diagnosis and monitor the treatment of respiratory diseases and that further studies are needed in this area [11, 12]. In this context, the FOT and the associated traditional parameters have been used to simplify the routine evaluation and to improve our understanding of the pathophysiology of several respiratory diseases [8]. A recent consensus statement pointed out that the use of FOT in occupational diseases requires further research [8].

A further improvement in the respiratory system evaluation is obtained by the use of inverse modeling. These models are based on electrical components analogous to the respiratory system properties of resistance, compliance, and inertance [13, 14]. It allows us to gain additional insight into the anatomical or pathophysiological changes that occur in respiratory diseases by obtaining detailed mechanical information about the respiratory system [10, 15, 16].

Recently, models based upon fractional derivatives and integrals have gained considerable popularity in the analysis of biological systems because they are more appropriate to describe the dynamic response of living systems than models based upon classical or integer-order derivatives and integrals [17,18,19]. In the particular case of the respiratory system, the ability of fractional-order models to effectively describe fractional power laws, hysteresis, and system memory is of pivotal importance. Parameters from these models can describe alterations in the distribution of gas flow within the system (spatial inhomogeneities), which are associated with peripheral abnormalities in the diseased lung. These parameters are also able to describe parenchymal deformity and elastance as well as changes in lung structure. Thanks to these characteristics, the fractional-order models have been the focus of much investigation in the search for improving respiratory clinical science and practice [10, 15, 16, 20].

In recent years, various fractional-order models have been developed and introduced [15, 21,22,23,24] to increase our understanding of cystic fibrosis [24], asthma in children [23], asbestos-exposed workers [25], chronic obstructive pulmonary disease [21, 22, 26], and sickle cell anemia [27]. Thus, there is a growing body of literature that recognizes the importance of fractional-order models in the interpretation of respiratory system functioning. The parameters obtained from these models hold the promise of improving the diagnosis and treatment of respiratory diseases.

Currently, however, these models show several unclear points, which resulted in much-debated topics, including the association of the model parameters and lung pathology [21], as well as the use of these parameters in the diagnostic of lung abnormalities [22]. The exact clarification of these points demands studies in well-defined groups of patients with different respiratory abnormalities [22]. The use of these models in asthma is limited to a preliminary study, including all asthma phenotypes [28]. This study, however, was not specific since it included all phenotypes and did not investigate the bronchodilator effect.

In this paper, we introduce integer and fractional-order models in the evaluation of the respiratory changes in work-related asthma and, for the first time, apply it to evaluate the bronchodilator response in these patients. This paper is organized as follows: first, the use of the new perspective of the FrOr model to interpret the changes in respiratory mechanics due to WRA will be evaluated. Then, we turn to the association between this model and the bronchodilator response in these patients. Finally, we evaluate the diagnostic accuracy of FrOr parameters and the best-performing parameters in the identification of the abnormal respiratory effects in WRA, as well as in the evaluation of the bronchodilator response in these patients.

Results

Thirty-one asthmatics and 31 healthy volunteers were selected for this study. Their anthropometric data are presented in Table 1. Among the 31 individuals with WRA, 25 never smoked, 3 were former smokers (average smoking load = 3.2 pack-years), and 3 were smokers (average smoking load = 8.3 pack-years).

Table 1 Biometric characteristics (mean ± SD) of control individuals and patients with work-related asthma

The spirometric and plethysmographic characteristics of the studied subjects are described in Table 2. Considering the severity of airway obstruction, according to spirometric parameters [29], 19 patients were classified as mild (61%), nine were moderate (29%), while 3 (10%) presented severe obstruction. The analysis of obstruction severity by the plethysmographic resistance [6] showed that 45% of the patients were classified as normal, 26% had mild and moderate obstruction, while 29% presented severe obstruction.

Table 2 Spirometric and plethysmographic characteristics of control individuals and patients with work-related asthma

Figure 3 describes the results of the mean respiratory resistance (a) and reactance (b) curves in controls and patients with WRA pre- and post-bronchodilator (BD) use. WRA increased resistance values (Fig. 3a) and introduced more negative reactance values (Fig. 3b). The BD use resulted in a reduction in resistance (Fig. 3a) and less negative reactance values (Fig. 3b).

Table 3 shows the effects of WRA and BD use on the resistive and reactive properties of the respiratory system. A detailed graphical description of these results is provided in Additional file 1: Figure S1 and Additional file 2: Figure S2. Considering the resistive parameters, R0 and R4 increased significantly in WRA (p < 0.05), while S and R20 − R4 also presented significant changes (p < 0.0001). All of the reactive parameters showed significant changes in the comparison between the control group and the WRA pre-BD. Interestingly, the use of BD in WRA introduced significant changes (p < 0.01) in all of the studied reactive parameters.

Table 3 Oscillometric characteristics of control individuals and patients with work-related asthma and the effect of bronchodilator use

Figure 4 displays the results obtained from the extended resistance–inertance–compliance (eRIC) model. The presence of abnormalities associated with WRA introduced significant increases in peripheral resistance (Rp, p < 0.0001) and total resistance (Rt, p < 0.01), as well as significant decreases in inertance (I, p < 0.001) and compliance (C, p < 0.01). In this figure, it is also apparent that the use of bronchodilator resulted in a significant decrease in all of the resistive parameters [central resistance (R), p < 0.0001; Rp, p < 0.05; Rt, p < 0.001] and a significant increase in C (p < 0.0001).

The parameters obtained from the FrOr modeling are described in Fig. 5. As can be seen from this figure, G and η increased significantly in patients with WRA (Fig. 5a, p < 0.01 and Fig. 1c, p < 0.0001, respectively), while H showed a significant decrease (p < 0.0001, Fig. 5c). Besides, the use of bronchodilator resulted in a significant decrease in G (p < 0.05).

Fig. 1
figure 1

Extended RIC model used to interpret respiratory impedance. R reflects the central airway resistance, and Rp is related to peripheral resistance, while I and C are associated with respiratory inertance and compliance, respectively

The errors analysis in the studied models is displayed in Table 4, which stands out the significant difference in Rd observed in patients before BD use.

Table 4 Mean square errors (MSE) and relative distances (Rd) in the integer and fractional-order models studied in control individuals and patients with work-related asthma pre- and post-bronchodilator use

Table 5 shows that R4 − R20 presented the best correlations among spirometric and traditional FOT resistive parameters (R = − 0.66), while fr presented the highest correlation among the reactive parameters (R = − 0.66). Among the integer-order model parameters, Rp presented the highest association with spirometry (R = − 0.62), while among the FrOr model parameters, η showed the closest association (R = − 0.40). The degrees of these associations were reasonable to good.

Table 5 Correlation analysis among traditional, eRIC, and fractional-order parameters, and the spirometric results

The results of the correlation analysis among FOT, eRIC, and FrOr parameters and the plethysmographic measurements are described in Table 6. Reasonable to good associations were observed. S showed the closest correlations among plethysmographic and traditional FOT resistive parameters (R = − 0.57), while Xm presented the best correlation among the reactive parameters (R = − 0.58). Among the integer-order model parameters, Rp presented the best association with plethysmography (R = − 0.50), while among the FrOr model parameters, η and G showed the highest correlation (R = − 0.53).

Table 6 Association among traditional, eRIC and fractional-order parameters and plethysmographic analysis

Three of the studied traditional FOT parameters presented adequate diagnostic accuracy (0.8 < AUC < 0.9) to identify respiratory abnormalities in patients with WRA (Table 7: S, Xm, Axt, and Axi). This initial ROC analysis also showed that fr was able to achieve high diagnostic accuracy (AUC > 0.9).

Table 7 Diagnostic accuracy (mean and 95% confidence interval) of the traditional FOT parameters in the diagnostic of respiratory abnormalities in patients with work-related asthma

Table 8 provides the results of the ROC analysis performed using the eRIC and FrOr models. Rp was the only parameter obtained using the eRIC model able to achieve adequate diagnostic accuracy. Considering the parameters obtained by the FrOr model, H achieved adequate diagnostic performance, while η showed high diagnostic accuracy.

Table 8 Diagnostic accuracy (mean and 95% confidence interval) of the eRIC and fractional-order parameters in the detection of respiratory alterations in patients with work-related asthma

Figure 6a shows the performance of the best parameters obtained from traditional FOT analysis, eRIC, and FrOr modeling to identify respiratory changes in WRA. The AUCs of fr and η were significantly higher than Rp (p < 0.04 and p < 0.02, respectively), while the AUCs of fr and η were not significantly different. The results of the LOOCV analysis in the most discriminative parameters observed in the detection of respiratory abnormalities in WRA using traditional parameters (Table 7), the eRIC, and the FrOr model (Table 8), are described in Fig. 6b. Fr achieved adequate diagnostic value, while η obtained a high accuracy.

The ability of the traditional FOT parameters to detect the changes due to bronchodilator use in patients with WRA is described in Table 9. Similar analyses for the eRIC and FrOr models are presented in Table 10.

Table 9 Area under the receiving operator curve (mean and 95% confidence interval) of the traditional FOT parameters in the identification of the bronchodilator response in patients with work-related asthma
Table 10 Diagnostic accuracy (mean and 95% confidence interval) of the eRIC and fractional-order parameters in the detection of respiratory effects of bronchodilator use in patients with work-related asthma

The results of the ROC analysis performed in the best parameters obtained from traditional FOT analysis (Table 9), and eRIC and FrOr modeling (Table 10) to identify bronchodilator responses in WRA are described in Fig. 2. LOOCV analysis considering these parameters was not able to reach adequate diagnostic values (AUC ≥ 0.7).

Fig. 2
figure 2

Receiver operator characteristic analysis and associated area under the curve in the three best parameters observed in the evaluation of the bronchodilator responses. R4: resistance in 4 Hz; Cdyn: dynamic compliance and Z4: impedance modulus in 4 Hz

Discussion

This is the first study that has quantitatively evaluated the performance of the eRIC model and a FrOr model in the analysis of airway obstruction and the bronchodilator response in work-related asthma. The most exciting findings of this study were as follows: (1) patients with WRA showed increased peripheral resistance, damping, and hysteresivity when compared with controls; (2) fractional-order analysis outperformed standard FOT, as well as integer-order modeling in the diagnosis of respiratory changes in these patients; (3) the bronchodilator use in WRA resulted in increased dynamic compliance and reduced damping and peripheral resistance; and (4) standard FOT analysis outperformed integer and fractional-order modeling in the identification of the bronchodilator effects in these patients.

Table 1 shows that the two studied groups were of comparable age, height, and gender distribution. Although there were slight differences between the groups related to body mass, this parameter is not determinants in terms of alterations in respiratory impedance. We should highlight the fact that the main parameter that has a significant impact on the impedance—subject height [30]—was quite similar between the groups under study. The changes in volumes and flows observed in WRA patients before and after BD use (Table 2) were consistent with the involved physiopathology [31].

There is a consensus in the literature that FOT is in state of the art in the analysis of pulmonary function, contributing to increasing our knowledge about respiratory diseases, as well as in its diagnosis. However, only one study used this method to investigate WRA [3]. This study focused only on the evaluation of methacholine challenge and was also limited by the evaluation to just R0. Variables relating to respiratory reactance and respiratory modeling were not investigated. Supporting and adding new information to these previous results, Fig. 3a shows an increase in the respiratory obstruction in WRA. This increase was more discriminating in the 4–16 Hz range, which resulted in increased values of R0 and R4 (Table 3, Additional file 1: Figure S1), as well as highly significant reductions in ventilation homogeneity (S and R4 − R20). Considering the reactive changes (Fig. 3b), WRA introduced more negative values in reactance in comparison with the control group, which resulted in significant changes in all of the reactive parameters (Table 3). These results are consistent with preliminary results in a smaller group [32] and can be explained by the presence of bronchoconstriction and inflammatory processes in asthmatics. These abnormalities reduce the diameter of the internal airways introducing increased airways resistance and changes in the time constants in the ventilatory process of these patients.

Fig. 3
figure 3

a Mean respiratory resistance and b reactance curves in controls and patients with work-related asthma (WRA) pre- and post-bronchodilator (BD) use

One interesting finding was that S, Xm, fr, Axt, and Axi presented fair values of AUC for clinical application (Table 7). These observations were in line with the usual interpretation of these parameters as being related to small airways disease [33,34,35] and the pattern of predominantly peripheral airway abnormalities in patients with mild obstruction, as may be characterized by the studied group.

Concerning the evaluation of the best methodology for calculating Ax, the correlations of Axi with spirometry (Table 5) and plethysmography (Table 6) were slightly higher than that obtained by Axt. Interestingly, the diagnostic accuracy of Axt in the identification of WRA respiratory abnormalities was higher than that observed in Axi (Table 7). It may be explained by the fact that, as can be seen in Fig. 3b, the reactance curve is not a perfect triangle. The approximation of the reactance area by a triangle amplifies the differences observed among the curves, improving the performance of Axt. Accordingly, Table 9 shows that the performance of Axt in the identification of the changes due to BD use was also higher than that obtained by Axi. Thus, although using a less accurate method to estimate area, Axt is more accurate than Axi in terms of clinical use. This counterintuitive finding may help elucidate the debate about the proper methodology for calculating Ax [8].

The bronchodilator use introduced a reduction in the resistance values and associated parameters (Fig. 3a, Table 3, and Additional file 1: Figure S1), as well as less negative values of respiratory reactance and parameters (Fig. 3b, Table 3). These results are in close agreement with the reduction in airway obstruction and the improvement in ventilation homogeneity usually observed after BD use in these patients [31].

The effects of WRA and BD use in the parameters associated with the eRIC model are described in Fig. 4. WRA does not introduce alterations in R, which indicates small changes in the central airways of the studies patients. This result is consistent with the data obtained using spirometry and plethysmography (Table 2), which described the presence of a predominantly small or moderate obstruction in the studied sample. BD use resulted in a reduction in R (Fig. 4a). A possible explanation for this result may be the typical smooth muscle relaxation that occurs in these individuals. The resulting mean R values were smaller than that measured in controls. These findings are also in line with the predominantly small or moderate obstruction observed in the studied WRA population (Table 2).

Fig. 4
figure 4

Influence of work-related asthma (WRA) and bronchodilator (BD) use on parameter values estimated from the eRIC model. Central airway resistance (R; a), peripheral resistance (Rp; b), total resistance (Rt; c), lung inertance (I; d), and alveolar compliance (C; 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

Peripheral resistance increased in WRA (Fig. 4a), which could be attributed to inflammation and airway wall remodeling. These results are also likely to be related to airway smooth muscle shortening, which introduces peripheral constriction. It was interesting to note that, even after the BD use and a reduction in its value, the Rp of WRA patients remained higher than that observed in controls. This result is in contrast to the reduction of R to values smaller than those measured in controls (Fig. 4a), which is probably related to the inflammatory effect of the disease.

Figure 4c shows that WRA introduced increased values of Rt. Considering that Rt = R + Rp [14] and that R was not increased in patients (Fig. 4a), we can speculate that this increase was associated with the increase observed in the peripheral resistance (Fig. 4b). This result reflects the fact that airway changes in asthma usually begin at the peripheral airways, and that the studied patients with WRA presented predominantly small or moderate obstruction (Table 2) so that they can be considered as in the early stages of the disease.

Respiratory inertance primarily describes the mass of gas that is moved during tidal breathing. It may be interpreted as an index related to pressure losses, as well, mostly due to the acceleration of the gas column in the central airways [8]. Respiratory inertance was reduced in WRA (Fig. 4d), which can be explained by the concepts of choke points [36] and apparent inertance [37]. Usually, inertance integrates the inertial characteristics of the whole respiratory system. As the respiratory obstruction advances, the oscillatory signal used by FOT to assess the impedance is prevented from passing through the choke points. It precludes FOT from considering the lung beyond the choke point so that the measured inertance reflects the airways proximal to the choke points. As a result, we observed a reduction in the apparent mass of the gas measured by the FOT, in the associated pressure necessary for the acceleration of the gas, and consequently, in the measured inertance. This process is similar to that observed in the apparent compliance and results in an apparent inertance. In line with this interpretation, direct associations were observed between inertance and spirometric indexes of peripheral airway obstruction (Table 5). Further additional supports to this hypothesis was provided by the inverse relationship obtained between inertance and the Raw and the direct association observed with SGaw (Table 6). Bronchodilation does not cause a significant change in I (Fig. 4d). Similar to S (Additional file 1: Figure S1C), R4 − R20 (Additional file 1: Figure S1G), and Rp (Fig. 4b), I remained distinct from the results obtained in the control group after BD use. These results probably reflect the irreversible inflammatory effect of the disease.

WRA introduced a decrease in C (Fig. 4e). This finding is consistent with the work of Bhatawadekar et al. [38], which used a single compartment model fit to estimate Ers (1/Crs). These authors pointed out that Ers is associated with small airways, and potentially a very clinically useful measure in asthma. This parameter includes the lungs and bronchial wall compliances, the compliance of the chest wall/abdomen compartment, and the thoracic gas compression. Thus, this result may be related to airway remodeling and frequency dependence of dynamic compliance due to non-uniform ventilation. The deformation of the thoracic wall associated with lung hyperinflation also needs to be considered since it introduces an essential restrictive factor in the interaction between the lung and thoracic wall. In Fig. 4e, it is also apparent that the use of bronchodilator resulted in a significant increase in C, which became similar to that presented in regular patients. These results further support the idea of the reduction in airway obstruction and the improvement in ventilation homogeneity after BD use in these patients.

Considering the diagnostic use of eRIC parameters, only Rp reached an adequate value for clinical use (Table 8). This finding is in close agreement with the interpretation of this parameter as reflecting peripheral airway resistance, and the presence of peripheral changes in our studied patients, which shows mainly mild obstruction (Table 2).

Recently, the concept of FrOr modeling of the respiratory system has received significant interest in the research community [23, 28, 39,40,41]. Theoretically, these emerging models present an improved sensitivity to pathologic changes, due to an improved ability to capture the characteristics of respiratory mechanics. In reviewing the literature, however, no data were found on the question of FrOr modeling in patients with WRA. The current study found increased values of G in WRA, presenting a significant reduction after BD use (Fig. 5a). These findings broadly support the interpretation linking WRA with increased energy dissipation in the respiratory system [15], which may be explained by the increased airway obstruction and reduced respiratory compliance. This finding also supports evidence from clinical observations reporting increased respiratory work and dyspnea on small efforts in these patients. The reduction after BD use is also consistent with the reduction of dyspnea usually observed after BD use in these patients [1].

Fig. 5
figure 5

Impact of work-related asthma (WRA) and bronchodilator (BD) use on parameter values estimated from the fractional-order model. Respiratory damping (G; a), elastance (H; b), and hysteresivity (η; c). 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 current study observed values of H in controls (Fig. 5b) similar to that obtained previously [25, 26]. In contrast with the results previously reported in non-specific asthma [28], mild obstruction in patients with WRA introduced a highly significant reduction in H. This provides evidence that asthma resulting from occupational exposure results in more aggressive changes in terms of the elastic properties of the respiratory system than in the average of asthmatics. These findings may have essential implications in the development of objective methods for the differential diagnosis between WRA and non-work-related asthma.

It is also interesting to point out that H was also reduced in mild patients with other obstructive diseases, including chronic obstructive pulmonary disease [26], and silicosis [42], but not in asbestos-exposed workers with mild abnormalities [25]. This difference may be attributed to the restrictive nature of the asbestosis. It is a compelling finding since it provides another evidence that H may help in the differential diagnosis of work-related respiratory diseases.

Perhaps the most interesting finding in FrOr analysis was the increase in η values observed in patients with WRA (Fig. 5c). It is in close agreement with the involved physiology, reflecting chronic airway inflammation and remodeling, which predisposes the lung to a more heterogeneous pattern of peripheral airway constriction. A comparison between the present results and those of a preliminary study, including all asthma phenotypes [28], confirms the ability of this parameter to describe the presence of heterogeneous peripheral ventilation in the specific phenotype of WRA. Additional supports of this interpretation are provided by other studies performed recently in patients with sickle cell anemia [27], chronic obstructive pulmonary disease [22, 26], and asbestos-exposed workers [25]. The hysteresivity increases with the hysteresis area of the pressure–volume loop [43], which associates this parameter with the work of breathing [15, 21]. Correlation analysis was consistent with this interpretation, describing inverse associations with spirometric indexes of airway obstruction, and direct associations with Raw (Tables 5 and 6, respectively). These findings indicate that η clearly describes the respiratory abnormalities in WRA, which are characterized by increased respiratory work [1].

Another interesting finding was the absence of changes in η values as a consequence of BD use (Fig. 5c). This result is in contrast with the reduction observed in G after BD use (Fig. 5a) and provided additional evidence of the association between η and peripheral abnormalities. Among FrOr parameters, η presented the highest correlations with spirometric and plethysmographic parameters of airway obstruction (Tables 5 and 6, respectively).

The range of measured values in asthmatics was reduced after BD use (Figs. 4 and 5). Patients with mild airway obstruction (61%) mainly compose the studied group of asthmatics. However, there are also 29% of patients with moderate and 10% with severe obstruction. This may explain the observed large range of measured values in asthmatics before BD use. After BD use, the airway obstruction tends to be reduced, and the respiratory system properties of the asthmatics tends to be closer to normal, reducing the range of measured values in these patients (Figs. 4 and 5).

It is now well established that fractional-order dynamic behavior may be linked to fractal structure, implying that properties of both function and structure are fundamentally interconnected [44]. It has been shown that recurrent fractal geometry may result in fractional-order terms [15]. In the particular case of the bronchial tree of normal subjects, a highly complex fractal structure is observed, in which the presence of self-similarity in its spatial structure is closely linked to a healthy lung function. In contrast, diseased lung presents asymmetry as a result of inhomogeneities due to the physiopathological process. The bronchial tree of a patient with WRA shows progressive loss of complexity in its spatial structure related to inflammation, airway remodeling, bronchoconstriction, edema, and fluid accumulation in the airways [31]. In line with these principles, previous studies from our laboratory demonstrate a consistent reduction in respiratory impedance complexity with increased airflow obstruction in a preliminary group, including all asthma phenotypes [45]. Further studies in similar groups of asthmatic subjects showed a significant increase in η and G with airway obstruction [28], indicating that these parameters are inversely related to respiratory complexity in these patients. It was hypothesized that the increase observed in η and G may be explained, at least partially, by the reduction in the complexity of the spatial structure of the airway tree of patients with asthma. Figure 5a, c provides further support to this hypothesis, extending this evidence to the specific case of WRA.

On the question of diagnostic use, this study found that η, obtained from FrOr modeling, reached a high diagnostic accuracy in identifying WRA abnormalities (Table 8). The comparison of the more accurate parameters obtained in traditional analysis, eRIC, and FrOr modeling showed that η was more accurate than Rp (Fig. 6). These results corroborate the propositions of previous authors, who suggested that FrOr models have the potential to improve respiratory clinical science and practice [10, 15, 16]. Also in line with this proposition and the observed results, it is apparent from the data in Table 4 that the FrOr model provided an improved description of the measured impedance. Following the present results, previous studies have demonstrated that FrOr models provide a more suitable fitting than integer-order models [26, 28]. It could be associated with the nature of the FrOr models, whose flexibility allows these models to adjust to fractional values of 20 dB/decade. Integer-order models, however, are only able to adjust to integer multiples of 20 dB/decade.

Fig. 6
figure 6

Receiver operator characteristic curves of the best parameters obtained from traditional FOT analysis, eRIC, and FrOr modeling to identify respiratory changes in WRA (a). Fr: resonance frequency; Rp: peripheral resistance; η: hysteresivity coefficient, AUC: area under the receiver operating characteristic curve; Se: sensitivity, Sp: specificity. More restrictive leave-one-out cross-validation (LOOCV) analyses performed on the same parameters (b)

These results are in keeping with previous studies, in which the detailed analysis offered by FrOr modeling improved our knowledge about several biomedical areas, including the properties of the arterial wall in brain aneurysms [48], the description of the red blood cell membrane mechanics [49], and the blood flow in the cranial network [50]. Similar improvements were also observed in modeling the viscoelastic non-linear compressible properties of the lung [46], the blood ethanol concentration [47], and improving the chemotherapy used in cancer treatment [18].

The evidence presented in Table 8 and Fig. 6 supports the notion that FrOr models may be useful in clinical use. The increase in diagnostic accuracy obtained in the present study (Fig. 6) is in close agreement with improvements observed in the differentiation between malignant and benign breast lesions detected on X-ray screening mammography [48], cancer detection [49], screening for hemodialysis patients [50], differentiation of low- and high-grade pediatric brain tumors [51], and Parkinson’s disease severity assessment [52].

Regarding these aspects, an initially unexpected result was observed from data in Tables 9 and 10 and Fig. 2. In contrast to the results observed in the identification of changes in respiratory mechanics in WRA, there was no evidence that the parameters obtained from the FrOr model were more accurate than the traditional FOT parameters to identify the respiratory effects of BD use. When these results are analyzed more carefully, it can be observed that the highest accuracy observed among the traditional FOT parameters (Table 9) was obtained by R4, which is related to central airway obstruction. Among model parameters (Table 10), the most accurate was R. This eRIC model parameter also reflects mainly central airway resistance. Both results are in close agreement with the spirometric and plethysmographic changes observed after BD use in the present study (Table 2), which described changes associated mainly to central airways. These results are in line with the recent work of Bhatawadekar et al. [38] investigating the bronchodilator response in asthma. Thus, FrOr parameters were not the most adequate to diagnose these changes because these parameters are more related to peripheral airways, while the observed BD responses are involved with more central airways.

The findings in this study are subject to at least three limitations. First, the present work is limited to patients with WRA. This focus allowed us to investigate this specific phenotype, clarifying the use of FOT and respiratory modeling in this critical disease. However, many other types of asthma exhibit different features. Therefore, further studies are needed to assess these specific disorders.

Secondly, one could argue that the bronchodilator analysis was limited to evaluate the adequacy of the studied parameters to reflect the respiratory changes due to BD use. A comparative analysis of a group of BD responders and non-responders using spirometry as a reference and including the administration of placebo and BD medication in a large sample of patients could establish cut-off points for changes in parameters derived from FOT models. It has important practical application and should, therefore, be addressed in future studies.

Finally, the present study investigated a relatively small sample size. Although this limitation was minimized using the LOOCV method, it is still a limitation, and additional studies, including a more significant number of subjects, are necessary. The present analysis, however, significantly contributes to the essential debates in the literature concerning the proper methodology for calculating Ax [8], the use of FOT in occupational health [8], particularly in WRA [3], as well as introduced respiratory modeling in this disease.

Conclusion

The present study provided clear evidence that patients with WRA show increased peripheral resistance, damping, and hysteresivity when compared with controls. These results provide new physiological insight into the effects of WRA on respiratory biomechanics. It was demonstrated that a combination of FOT and fractional-order modeling outperformed standard FOT, as well as integer-order modeling in the diagnosis of respiratory abnormalities in these patients, leading to high diagnostic accuracy. It was also shown that the use of bronchodilator in WRA resulted in increased dynamic compliance and reduced damping and peripheral resistance. FOT parameters may adequately identify these changes. Taken together, these results show the utility of the FOT associated with fractional-order modeling in the analysis of the respiratory abnormalities in patients with WRA.

Materials and methods

The ethical investigation clearance was obtained from the Research Ethics Committee of the Pedro Ernesto University Hospital (protocol 456-CEP/HUPE). The study obeys the Declaration of Helsinki. Before data collection, the participants received an explanation of the project. Upon obtaining written informed consent from patients, the respiratory analysis was carried out.

Study design, volunteers, and inclusion and exclusion criteria

This study utilized a cross-sectional design involving two groups of subjects: a group of patients with WRA and a group of controls. The asthmatics were diagnosed according to the GINA criteria [1]. These patients were not allowed to use the bronchodilator for at least 12 h before the test. They have not smoked for at least 2 h, nor have drunk coffee or alcohol for 6 h before the tests. The control group had no present or previous cardiorespiratory disease or medication, had no respiratory symptoms, and the findings in flow-volume spirometry were normal. For both groups, additional exclusion criteria were tuberculosis, inability to perform the tests, chemotherapeutic or radiotherapeutic treatment, and respiratory infections in the last 30 days.

The WRA group measurements were analyzed for two moments: pre-bronchodilator (WRA Pre-BD) and post-bronchodilator (WRA Post-BD). The use of bronchodilator medication (400 μg salbutamol sulfate spray) was performed immediately after the first set of pulmonary function tests. The second group of tests was performed 15 min later.

Forced oscillation

FOT was evaluated using small amplitude pressure oscillations (≤ 2 cmH2O) generated by a speaker applied during spontaneous breathing at the entrance of the airway through the oral cavity. Three tests were conducted, each 16 s long, with the mean score being adopted as the final result. The test was considered acceptable if the volunteers presented stable tidal volumes and rate and free of pauses. Common artifacts, such as swallows, coughs, and leaks, were identified by the evaluation of flow and pressure signals. The acquisition was repeated until three stable and free of artifact measurements were obtained. The forced pseudo-random noise used in this study was composed of a frequency range between 4 and 32 Hz. To reduce the influence of the spontaneous breathing signal in the lowest frequency range, the minimum coherence function (CF) used for the acceptance of results was 0.9. The exams are repeated until all analyzed frequencies have this minimum CF value. To exclude outlying values, the coefficient of variability at the lowest oscillation frequency (4 Hz) of the three used tests was ≤ 10%. The analyses were performed using an OSCILAB 2.0 impedance analyzer developed at our laboratory [53].

The real part of impedance was submitted to linear regression analysis in the 4–16 Hz range, which yielded respiratory resistance extrapolated at 0 Hz (R0), and frequency dependence of Rrs was expressed as the slope (S) of the linear relationship between the resistive impedance and frequency. The mean resistance (Rm) in this frequency range was also evaluated. R0 is related to the low-frequency range. Newtonian resistance of the respiratory system, including the airways and resistance of tissue originating from the lung and chest wall, along with the effect of gas redistribution (pendelluft) [54]. S describes the resistance change with frequency and is related to respiratory system non-homogeneities [55]. These non-homogeneities are associated with increased peripheral resistance. The wide range of frequencies used in FOT allows for independent assessment of both proximal and peripheral airways. Lower frequencies (e.g., 4 Hz) penetrate deeply in the lung structure and are related to the sum of the proximal and peripheral airways. As the frequency rises, the measurement signal penetrates less and less into the lung structure, so that at higher frequencies (e.g., 20 Hz) is related only to proximal airways. Consequently, the resistance change with frequency has been used as a surrogate marker to describe peripheral resistance. Rm describes the mid-frequency range resistance, which is sensitive to the airway caliber, reflecting resistance in the central airways [56]. Other usual indexes of respiratory resistance were also studied; the resistances in 4 Hz (R4), 12 Hz (R12), and 20 Hz (R20), representing the low-, mid-, and high-frequency spectra, respectively, and the frequency dependence of resistance, which was represented as the difference between R4 and R20 (R4 − R20).

The results associated with the imaginary part of the impedance were interpreted using five parameters: mean reactance (Xm), resonance frequency (fr), respiratory system dynamic compliance (Cdyn), the impedance module, and the area under the reactance curve. Mean reactance is generally related to respiratory system inhomogeneity, and in this study, it was calculated using the 4 to 32 Hz frequency range. The frequency at which Xrs becomes zero is known as the resonance frequency [57]. Cdyn is related to the total compliance of the respiratory system, comprising pulmonary compliance, chest wall compliance, and airway compliance. This parameter is also related to the homogeneity of the respiratory system [56]. Cdyn was calculated based on the reactance at 4 Hz (Cdyn = 1/2πfX4). The 4 Hz impedance module (Z4) was also studied, which reflects the total mechanical load of the respiratory system [58].

The area under the reactance curve from the lowest frequency to fr (Ax) is also reported. There is a debate about the proper methodology for calculating Ax that has gained recent prominence [8], which includes the frequency resolution and numerical integration method as essential points to be clarified. To help address these research gaps, we evaluated Ax using two different methods. First, we used the approximation of the area by a triangle (AXt) defined by the lowest frequency studied (4 Hz), the corresponding reactance value (X4), and the resonance frequency. The second method uses the integral based on the trapezoidal rule to obtain a more accurate area estimate (Axi).

Spirometry and plethysmography

After FOT, spirometry was performed and interpreted according to the recommendations of the American Thoracic Society/European Respiratory Society [59]. Plethysmographic exams were conducted using the reference values described by Ref. [60] and a constant volume and variable pressure plethysmograph (HD CPL nSpire Health Ltd., Hertford, UK). The studied parameters were the airway resistance (Raw), the total lung capacity (TLC), functional residual capacity (FRC) and residual volume (RV), as well as their relationships (RV/TLC and FRC/TLC).

Integer-order modeling

The extended resistance–inertance–compliance (eRIC) compartmental model was used to interpret the changes in the respiratory system impedance due to the presence of WRA (Fig. 1). The cited model was introduced as an improved version of the basic RIC model [10].

This model offers a detailed description of the respiratory system properties, using R to describe the central airway resistance, while the peripheral resistance (Rp) is associated with the small airways. Rp allows for the frequency dependence of resistance values typically observed in patients, which is beyond the basic RIC capability. Reactive properties are described in this model by the respiratory inertance (I) and compliance (C) [14].

Fractional-order modeling

A recently described [23] fractional model of the respiratory impedance (ZFrOr) was used, according to Eq. (1):

$$ Z_{\text{FrOr}} \left( {j\omega } \right) = {\text{FrL}}\left( {j\omega } \right)^{\alpha } + \frac{1}{{{\text{FrC}}\left( {j\omega } \right)^{\beta } }}, $$
(1)

where FrL represents a frequency-dependent fractional-order inertance, associated with a fractional inertance coefficient (0 ≤ α ≤ 1), while FrC represents the constant-phase fractional-order compliance, which in turn is associated with a fractional compliance coefficient (0 ≤ β ≤ 1). Figure 7 depicts this model.

Fig. 7
figure 7

Fractional-order model based on a series association of a constant phase inertance (CPL) and constant phase compliance (CPC). These elements are composed of a frequency-dependent fractional inertia (FrL) and a frequency-dependent fractional compliance (FrC) elements, both related to their associated fractional exponents α and β

Figure 7 shows that the coupling of airway structures with dissipative energy (resistive), as well as storage energy properties (inertial), is described by FrL. The associated α coefficient modulates the influence of FrL in the frequency domain. As α approaches to zero, the influence of FrL in the airway resistance increases, with a concomitant reduction of its influence in airway inertance. Considering the impedance curves, it describes increased resistance values with frequency and more negative values of reactance in higher frequencies.

As described in Fig. 7, more peripheral structures of the respiratory system were described as the constant-phase fractional-order compliance (FrC). This parameter reflects structures that elastically store energy, and that is also simultaneously coupled to those that dissipate energy. The influence of FrC in resistance and reactance is modulated by β, where lower values are related to increased resistance and reduced compliance. In impedance traces, it results in increased resistance values, as well as more negative values of reactance in the low-frequency range.

The cited parameters will change according to the physiological relations of the respiratory system, involving morphology and geometry [15]. The damping factor (G) is usual in the interpretation of these physiological relations, representing the energy dissipation in the respiratory system [15]. This parameter is defined as follows:

$$ G = \frac{1}{C}\cos \left( {\frac{\pi }{2}\beta } \right). $$
(2)

The elastance (H) is a measure of potential elastic energy accumulation, according to Eq. (3):

$$ H = \frac{1}{C}\sin \left( {\frac{\pi }{2}\beta } \right). $$
(3)

Another widely used parameter is the hysteresivity coefficient (η), which is proportional to the heterogeneity of ventilation in the lung [15] and defined as

$$ \eta = \frac{G}{H}. $$
(4)

The fitting of the selected parameters was implemented using the Levenberg–Marquardt algorithm, determining the set of parameters of the model that best represents the input dataset in terms of least squares. The measured FOT data in the frequency range between 4 and 32 Hz was used as input dataset. The software for curve fitting was developed in the LABVIEW™ 12.0 (National Instruments, Austin, TX) environment, providing as a measure of the goodness of fit of the model the total error value (MSEt). This parameter was calculated as the square root of the sum of the real (MSEr) and imaginary (MSEx) impedance estimation errors. Following the procedure used by Oostveen et al. [61], a further error analysis was performed using the mean relative distance from the model and measured resistance and reactance values (Rd).

Statistical analysis

A commercial software (Origin® 8.0, Microcal Software Inc., Northampton, Massachusetts, United States) was used to assess normality (Shapiro–Wilk test) and to perform statistical tests (t test or Mann–Whitney). The results are presented as the mean ± SD, and statistical significance was considered when p ≤ 0.05. The association of model parameters and pulmonary function was investigated using the correlation analysis. Pearson or Spearman correlation was used depending on whether data distribution is normal or not. The following categorization of the strength of these associations was used [62]:

  • Small or no correlation (− 0.25 to 0.25);

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

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

  • Very good to excellent correlation: (0.75 to 1, or − 0.75 to − 1).

It was used a correction in the significance level to minimize the chances of making a Type I error (modified Bonferroni) due to the computation of several correlations [63]. This correction was performed dividing the p-value by an estimate of the effective number of independent correlations used. Usually, four independent variables are observed in traditional pulmonary function exams, while two independent variables are associated with the resistive and reactive properties of the FOT. Thus, a corrected significance level of 0.0063 (0.05/8), associated with eight independent correlations, was used.

Receiver operating characteristic (ROC) curves were used to evaluate the clinical potential of the FOT indexes. The values of sensitivity, specificity, and area under the curve (AUC) were obtained based on the optimal cut-off point, as determined by the ROC curve analysis. MedCalc 12 (MedCalc Software, Mariakerke, Belgium) was used in these analyses, which followed the STARD [64] requirements for studies of diagnostic accuracy. Leave-one-out cross-validation (LOOCV) was employed to minimize the statistical problem of finite patients. These evaluations were conducted as described in Ref. [65].

The minimum value of the AUC considered adequate for the identification of the changes due to WRA was 0.8 [66, 67]. Following previous studies concerning the use of bronchodilators in asthmatics [68, 69], 0.7 was considered to be a reasonable cut-off value for the diagnosis of respiratory changes due to BD use.

The sample size was calculated using MedCalc 12.3 (MedCalc Software, Mariakerke, Belgium). The used criterium was the comparison of the area under a ROC curve with a null hypothesis value [70]. Here, the objective was to show that adequate diagnostic accuracy (AUC = 0.8) [66, 67] was significantly different from the null hypothesis (AUC = 0.5), which indicates no clinical diagnostic value. The initial results obtained in a pilot study, including 17 controls and 12 patients [32], were used. Considering adequate type I and type II errors of 0.10, this analysis resulted in a minimum of 29 volunteers per group.

Availability of data and materials

The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

α :

Fractional inertance coefficient

β :

Fractional compliance coefficient

η :

Hysteresivity coefficient

AUC:

Area under the receiver operating characteristic curve

Axi:

Area under the reactance curve using the integral based on the trapezoidal rule

Axt:

Area under the reactance curve using the approximation by a triangle

BD:

Bronchodilator

C :

Compliance

C dyn :

Dynamic compliance

CF:

Coherence function

EOA:

Exacerbated occupational asthma

eRIC:

Extended resistance–inertance–compliance model

FEV1 :

Forced expiratory volume in the first second

FOT:

Forced oscillation technique

fr:

Resonant frequency

FrC:

Fractional-order compliance

FRC:

Functional residual capacity

FrOr:

Fractional-order

FrL:

Fractional-order inertance

FVC:

Forced vital capacity

G :

Damping factor

H :

Elastance

I :

Respiratory inertance

LOOCV:

Leave-one-out cross-validation

MSEr:

Real mean square error value

MSEt:

Total mean square error value

MSEx:

Imaginary mean square error value

OA:

Occupational asthma

R :

Central airway resistance

R p :

Peripheral resistance

R 0 :

Intercept resistance

R 4 :

Respiratory resistance in 4 Hz

R 12 :

Respiratory resistance in 12 Hz

R 20 :

Respiratory resistance in 20 Hz

R 4 − R 20 :

Difference between the respiratory resistance in 4 Hz and 20 Hz

R aw :

Airway resistance

R d :

Mean relative distance from the model and measured resistance and reactance

R m :

Mean resistance in the 4–16 Hz range

ROC:

Receiver operating characteristic curve

R rs :

Respiratory resistance

RV:

Residual volume

S :

Angular coefficient of resistance in the 4–16 Hz range

Se:

Sensitivity

Sp:

Specificity

TLC:

Total lung capacity

WRA:

Work-related asthma

X 4 :

Respiratory reactance at 4 Hz

X m :

Mean reactance evaluated considering the 4 to 32 Hz frequency range

X rs :

Respiratory reactance

Z 4 :

Impedance module in 4 Hz

Z FrOr :

Fractional-order respiratory impedance

References

  1. World Health Organization. GINA—Global initiative for asthma. Geneva: WHO; 2015.

    Google Scholar 

  2. Vila-Rigat R, Valls RP, Huet EH, Maristany JS, Prat XB, Muñoz-Ortiz L, Monserrat PT, Santacana VR. Prevalence of work-related asthma and its impact in primary health care. Archivos de Bronconeumología. 2015;51:449–55 (English Edition).

    Article  Google Scholar 

  3. Descatha A, Fromageot C, Ameille J, Lejaille M, Falaize L, Louis A, Lofaso F. Is forced oscillation technique useful in the diagnosis of occupational asthma? J Occup Environ Med. 2005;47:847–53.

    Article  Google Scholar 

  4. Carvalhaes-Neto N, Lorino H, Gallinari C, Escolano S, Mallet A, Zerah F, Harf A, Macquin-Mavier I. Cognitive function and assessment of lung function in the elderly. Am J Respir Crit Care Med. 1995;152:1611–5.

    Article  Google Scholar 

  5. Amaral JL, Lopes AJ, Veiga J, Faria AC, Melo PL. High-accuracy detection of airway obstruction in asthma using machine learning algorithms and forced oscillation measurements. Comput Methods Programs Biomed. 2017;144:113–25.

    Article  Google Scholar 

  6. Madama VC. Pulmonary function testing and cardiopulmonary stress testing. 2nd ed. Albany: Delmar Publishers; 1998.

    Google Scholar 

  7. de Mir Messa I, Sardon Prado O, Larramona H, Salcedo Posadas A, Villa Asensi JR. Body plethysmography (I): standardisation and quality criteria. An Pediatr. 2015;83(136):e131–7.

    Google Scholar 

  8. King GG, Bates J, Berger KI, Calverley P, de Melo PL, Dellaca RL, Farre R, Hall GL, Ioan I, Irvin CG, et al. Technical standards for respiratory oscillometry. Eur Respir J. 2020;55:1900753.

    Article  Google Scholar 

  9. LA Silva KKD, Jansen JM, Melo PL. Total inspiratory and expiratory impedance in patients with severe chronic obstructive pulmonary disease. Clinics. 2011;66(12):2085–91.

    Article  Google Scholar 

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

    Google Scholar 

  11. Brusasco V, Barisione G, Crimi E. Pulmonary physiology: future directions for lung function testing in COPD. Respirology. 2015;20:209–18.

    Article  Google Scholar 

  12. MacIntyre NR. The future of pulmonary function testing. Respir Care. 2012;57:154–61 (discussion 161–154).

    Article  Google Scholar 

  13. Meraz EG, Nazeran H, Ramos CD, Nava P, Diong B, Goldman MD, Goldman CA. 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 

  14. 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 

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

    Book  MATH  Google Scholar 

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

    Book  Google Scholar 

  17. Eke A, Herman P, Kocsis L, Kozak LR. Fractal characterization of complexity in temporal physiological signals. Physiol Meas. 2002;23:R1–38.

    Article  Google Scholar 

  18. Pachauri N, Yadav J, Rani A, Singh V. Modified fractional order IMC design based drug scheduling for cancer treatment. Comput Biol Med. 2019;109:121–37.

    Article  Google Scholar 

  19. Magin RL. Fractional calculus models of complex dynamics in biological tissues. Comput Math Appl. 2010;59:1586–93.

    Article  MathSciNet  MATH  Google Scholar 

  20. Kaczka DW, Ingenito EP, Suki B, Lutchen KR. Partitioning airway and lung tissue resistances in humans: effects of bronchoconstriction. J Appl Physiol. 1997;82:1531–41.

    Article  Google Scholar 

  21. Ionescu CM, De Keyser R. Relations between fractional-order model parameters and lung pathology in chronic obstructive pulmonary disease. IEEE Trans Biomed Eng. 2009;56:978–87.

    Article  Google Scholar 

  22. Ionescu C, Derom E, De Keyser R. Assessment of respiratory mechanical properties with constant-phase models in healthy and COPD lungs. Comput Methods Programs Biomed. 2010;97:78–85.

    Article  Google Scholar 

  23. Ionescu C, Desager K, De Keyser R. Fractional order model parameters for the respiratory input impedance in healthy and in asthmatic children. Comput Methods Programs Biomed. 2011;101:315–23.

    Article  Google Scholar 

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

    Article  Google Scholar 

  25. de Sa PM, Castro HA, Lopes AJ, Melo PL. Early diagnosis of respiratory abnormalities in asbestos-exposed workers by the forced oscillation technique. PLoS ONE. 2016;11:e0161981.

    Article  Google Scholar 

  26. Ribeiro CO, Faria AC, Lopes AJ, Melo PL. Forced oscillation technique for early detection of the effects of smoking and chronic obstructive pulmonary disease: contribution of fractional-order modeling. Int J COPD. 2018;13:3281–95.

    Article  Google Scholar 

  27. Marinho CL, Maioli MCP, do Amaral JLM, Lopes AJ, Melo PL. Respiratory resistance and reactance in adults with sickle cell anemia: Part 2—fractional-order modeling and a clinical decision support system for the diagnostic of respiratory disorders. PLoS ONE. 2018;14:e0213257.

    Article  Google Scholar 

  28. Faria AC, Veiga J, Lopes AJ, Melo PL. Forced oscillation, integer and fractional-order modeling in asthma. Comput Methods Programs Biomed. 2016;128:12–26.

    Article  Google Scholar 

  29. Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, Casaburi R, Coates A, van der Grinten CP, Gustafsson P, Hankinson J, et al. Interpretative strategies for lung function tests. Eur Respir J. 2005;26:948–68.

    Article  Google Scholar 

  30. Landser FJ, Clement J, Van de Woestijne KP. Normal values of total respiratory resistance and reactance determined by forced oscillations: influence of smoking. Chest. 1982;81:586–91.

    Article  Google Scholar 

  31. Chan-Yeung M, Malo JL, Tarlo SM, Bernstein L, Gautrin D, Mapp C, Newman-Taylor A, Swanson MC, Perrault G, Jaques L, et al. Proceedings of the first Jack Pepys occupational asthma symposium. Am J Respir Crit Care Med. 2003;167:450–71.

    Article  Google Scholar 

  32. Tuza FAD, Sá PM, Lopes AJ, Castro HA, Melo PL. Evaluation of the forced oscillation technique and the extended RIC model in the analysis of individuals with work-related asthma. Book oscillation mechanics and response to bronchodilation of individuals with work-related asthma, vol. 70., The International Federation for Medical and Biological Engineering (IFMBE) Proceedings book seriesNew York: Springer; 2018. p. 37–43.

    Google Scholar 

  33. Foy BH, Soares M, Bordas R, Richardson M, Bell A, Singapuri A, Hargadon B, Brightling C, Burrowes K, Kay D, et al. Lung computational models and the role of the small airways in asthma. Am J Respir Crit Care Med. 2019;200:982–91.

    Article  Google Scholar 

  34. Postma DS, Brightling C, Baldi S, Van den Berge M, Fabbri LM, Gagnatelli A, Papi A, Van der Molen T, Rabe KF, Siddiqui S, et al. Exploring the relevance and extent of small airways dysfunction in asthma (ATLANTIS): baseline data from a prospective cohort study. Lancet Respir Med. 2019;7:402–16.

    Article  Google Scholar 

  35. Marinho CL, Maioli MCP, do Amaral JLM, Lopes AJ, Melo PL. Respiratory resistance and reactance in adults with sickle cell anemia: correlation with functional exercise capacity and diagnostic use. PLoS ONE. 2017;12:e0187833.

    Article  Google Scholar 

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

    Article  Google Scholar 

  37. Lima AN, Faria AC, Lopes AJ, Jansen JM, Melo PL. Forced oscillations and respiratory system modeling in adults with cystic fibrosis. Biomed Eng Online. 2015;14:11.

    Article  Google Scholar 

  38. Bhatawadekar SA, Leary D, de Lange V, Peters U, Fulton S, Hernandez P, McParland C, Maksym GN. Reactance and elastance as measures of small airways response to bronchodilator in asthma. J Appl Physiol. 2019;127:1772–81.

    Article  Google Scholar 

  39. Beaulieu A, Bosse D, Micheau P, Avoine O, Praud JP, Walti H. Measurement of fractional order model parameters of respiratory mechanical impedance in total liquid ventilation. IEEE Trans Biomed Eng. 2012;59:323–31.

    Article  Google Scholar 

  40. Copot D, De Keyser R, Derom E, Ionescu C. Structural changes in the COPD lung and related heterogeneity. PLoS ONE. 2017;12:e0177969.

    Article  Google Scholar 

  41. Ionescu C, Lopes A, Copot D, Machado JAT, Bates JHT. The role of fractional calculus in modeling biological phenomena: a review. Commun Nonlinear Sci Numer Simul. 2017;51:141–59.

    Article  MathSciNet  MATH  Google Scholar 

  42. Faria ACD, Carvalho ARS, Guimaraes ARM, Lopes AJ, Melo PL. Association of respiratory integer and fractional-order models with structural abnormalities in silicosis. Comput Methods Programs Biomed. 2019;172:53–63.

    Article  Google Scholar 

  43. Fredberg JJ, Stamenovic D. On the imperfect elasticity of lung-tissue. J Appl Physiol. 1989;67:2408–19.

    Article  Google Scholar 

  44. Thamrin C, Frey U, Kaminsky DA, Reddel HK, Seely AJ, Suki B, Sterk PJ. Systems biology and clinical practice in respiratory medicine. The Twain Shall Meet. Am J Respir Crit Care Med. 2016;194:1053–61.

    Article  Google Scholar 

  45. Veiga J, Lopes AJ, Jansen JM, Melo PL. Fluctuation analysis of respiratory impedance waveform in asthmatic patients: effect of airway obstruction. Med Biol Eng Comput. 2012;50:1249–59.

    Article  Google Scholar 

  46. Birzle AM, Wall WA. A viscoelastic nonlinear compressible material model of lung parenchyma—experiments and numerical identification. J Mech Behav Biomed Mater. 2019;94:164–75.

    Article  Google Scholar 

  47. Qureshi S, Yusuf A, Shaikh AA, Inc M, Baleanu D. Fractional modeling of blood ethanol concentration system with real data application. Chaos. 2019;29:013143.

    Article  MathSciNet  MATH  Google Scholar 

  48. Bickelhaupt S, Steudle F, Paech D, Mlynarska A, Kuder TA, Lederer W, Daniel H, Freitag M, Delorme S, Schlemmer HP, Laun FB. On a fractional-order calculus model in diffusion-weighted breast imaging to differentiate between malignant and benign breast lesions detected on X-ray screening mammography. PLoS ONE. 2017;12:e0176077.

    Article  Google Scholar 

  49. Sinkus R, Siegmann K, Xydeas T, Tanter M, Claussen C, Fink M. MR elastography of breast lesions: understanding the solid/liquid duality can improve the specificity of contrast-enhanced MR mammography. Magn Reson Med. 2007;58:1135–44.

    Article  Google Scholar 

  50. Wu JX, Li CM, Chen GC, Ho YR, Lin CH. Peripheral arterial disease screening for hemodialysis patients using a fractional-order integrator and transition probability decision-making model. IET Syst Biol. 2017;11:69–76.

    Article  Google Scholar 

  51. Sui Y, Wang H, Liu G, Damen FW, Wanamaker C, Li Y, Zhou XJ. Differentiation of low- and high-grade pediatric brain tumors with high b-value diffusion-weighted MR imaging and a fractional order calculus model. Radiology. 2015;277:489–96.

    Article  Google Scholar 

  52. Jan Mucha J, Mekyska J, Galaz Z, Faundez-Zanuy M, Lopez-de-Ipina K, Vojtech Zvoncak V, Kiska T, Smekal Z, Brabenec L, Rektorova I. Identification and monitoring of Parkinson’s disease dysgraphia based on fractional-order derivatives of online handwriting. Appl Sci. 2018;8:2566–84.

    Article  Google Scholar 

  53. 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 

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

    Article  Google Scholar 

  55. Brochard L, Pelle G, de Palmas J, Brochard P, Carre A, Lorino H, Harf A. Density and frequency dependence of resistance in early airway obstruction. Am Rev Respir Dis. 1987;135:579–84.

    Google Scholar 

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

    Article  Google Scholar 

  57. 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 

  58. 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.

    Article  Google Scholar 

  59. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, Van Der Grinten CPM, Gustafsson P, et al. Standardisation of spirometry. Eur Respir J. 2005;26(2):319–38.

    Article  Google Scholar 

  60. Neder JA, Andreoni S, Castelo A, Nery LE. Reference values for lung function tests. I. Static volumes. Braz J Med Biol Res. 1999;32:703–17.

    Article  Google Scholar 

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

    Google Scholar 

  62. Dawson B, Trapp RG. Basic & clinical biostatistics. 4th ed. New York: Lange Medical Books-McGraw-Hill, Medical Pub. Division; 2004.

    Google Scholar 

  63. Conneely KN, Boehnke M. So many correlated tests, so little time! Rapid adjustment of P values for multiple correlated tests. Am J Hum Genet. 2007;81:1158–68.

    Article  Google Scholar 

  64. Bossuyt PM, Reitsma JB, Bruns DE, Gatsonis CA, Glasziou PP, Irwig L, Lijmer JG, Moher D, Rennie D, de Vet HC, et al. STARD 2015: an updated list of essential items for reporting diagnostic accuracy studies. BMJ. 2015;351:h5527.

    Article  Google Scholar 

  65. Witten IH, Frank E, Hall MA, Pal CJ. Data mining: practical machine learning tools and techniques. Burlington: Morgan Kaufmann; 2016.

    Google Scholar 

  66. Swets JA. Measuring the accuracy of diagnostic systems. Science. 1988;240:1285–93.

    Article  MathSciNet  MATH  Google Scholar 

  67. Golpe R, Jimenez A, Carpizo R, Cifrian JM. Utility of home oximetry as a screening test for patients with moderate to severe symptoms of obstructive sleep apnea. Sleep. 1999;22:932–7.

    Google Scholar 

  68. Goedhart DM, Zanen P, Kerstjens HA, Lammers JW. Discriminating asthma and COPD based on bronchodilator data: an improvement of the methods. Physiol Meas. 2005;26:1115–23.

    Article  Google Scholar 

  69. Dundas I, Chan EY, Bridge PD, McKenzie SA. Diagnostic accuracy of bronchodilator responsiveness in wheezy children. Thorax. 2005;60:13–6.

    Article  Google Scholar 

  70. Hanley JA, Mcneil BJ. The meaning and use of the area under a receiver operating characteristic (Roc) curve. Radiology. 1982;143:29–36.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the Brazilian Council for Scientific and Technological Development (CNPq), the Rio de Janeiro State Research Supporting Foundation (FAPERJ), FIOTEC (Fundação para o Desenvolvimento Científico e Tecnológico em Saúde), and in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brasil (CAPES)-Finance Code 001.

Author information

Authors and Affiliations

Authors

Contributions

FAAT collected data regarding pulmonary function. FAAT and PMS analyzed the data, and drafted the manuscript. HAC provided funding, data, and subject identification. AJL and PLM mentored FAAT, participated in the data analysis process, and helped draft the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pedro Lopes de Melo.

Ethics declarations

Ethics approval and consent to participate

The Research Ethics Committee of the Pedro Ernesto University Hospital granted approval for this study, and the collection and use of the data analyzed in this study.

Consent for publication

All the authors consent for the publication of this manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Comparative analysis of the traditional FOT resistive parameters obtained in controls and patients with work-related asthma (WRA) pre and post bronchodilator use: resistance in 0 Hz (R0; A), mean resistance (Rm; B), slope of the resistance values (S; C), resistance in 4 Hz (R4; D), resistance in 12 Hz (R12; E), resistance in 20 Hz (R20; F) and difference in the resistance 4 Hz and 20 Hz (R4 − R20; G).

Additional file 2: Figure S2.

Comparative analysis of the traditional FOT reactive parameters obtained in controls and patients with work-related asthma (WRA) pre and post bronchodilator (BD) use: mean reactance (Xm; A), resonant frequency (Fr; B), dynamic compliance (Cdyn; C), reactance area approximated by a triangle (Axt: D), by an integral (Axi; E), impedance module in 4 Hz (Z4; F).

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://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

Tuza, F.A.d., de Sá, P.M., Castro, H.A. et al. Combined forced oscillation and fractional-order modeling in patients with work-related asthma: a case–control study analyzing respiratory biomechanics and diagnostic accuracy. BioMed Eng OnLine 19, 93 (2020). https://doi.org/10.1186/s12938-020-00836-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12938-020-00836-6

Keywords