3.1. Demographic and clinical characteristics of the study population
The study involved a group of 48 IPF individuals, with an average age of 72±2 years. Of the participants, 33 out of 48 (69%) were male. The average age at which IPF was diagnosed was 70±7 years. Among the group, 63% were ex-smokers, 33% were non-smokers, and only 6% were current smokers, with an average of 47±38 pack years. Interestingly, a higher proportion of males were ex-smokers compared to females. About 54% of the participants were under antifibrotic treatment for IPF, with 48% receiving pirfenidone, and the remainder receiving nintedanib. The main characteristics of the sample and according to gender are presented in Table 1.
Males were significantly more likely to have a history of smoking compared to females, indicating a higher prevalence of smoking in the male population. Additionally, disease severity, as measured by the GAP index, was greater in males, suggesting that male patients may have a worse prognosis. CAD score, was also significantly higher in males than in females, highlighting an increased risk of cardiovascular comorbidities in male IPF patients. These findings emphasize the influence of gender on smoking history, disease progression, and comorbidities in IPF.
Table 1. The main characteristics of the study population (N=48), and comparisons between genders
Parameter
|
Sum (N=48)
|
Males (n=33)
|
Females (n=15)
|
p-value
|
Age, years
|
72±2
|
73±8
|
71±5
|
0.541
|
Age of IPF diagnosis, years
|
70±7
|
70±8
|
70±6
|
0.953
|
Current smokers, n (%)
|
3 (6)
|
2 (6)
|
1 (7)
|
0.780
|
Ever smokers, n (%)
|
30 (63)
|
25 (76)
|
5 (33)
|
0.019
|
Nonsmokers, n (%)
|
15 (31)
|
6 (18)
|
9 (6)
|
0.890
|
Pys
|
47±38
|
43±22
|
66±77
|
0.180
|
Under antifibrotic treatment, n (%)
|
26 (54)
|
20 (61)
|
6 (40)
|
0.155
|
Pirfenidone, n (%)
|
13 (27)
|
10 (30)
|
3 (20)
|
0.862
|
Nintedanib, n (%)
|
13 (27)
|
10 (30)
|
3 (20)
|
0.862
|
MRC score
|
2±1
|
2±1
|
2±1
|
0.636
|
GAP index
|
4±1
|
4±1
|
3±1
|
0.016
|
GAP stage I
|
26 (54)
|
15 (45)
|
11 (73)
|
0.106
|
GAP stage II
|
16 (33)
|
12(37)
|
4 (27)
|
0.106
|
GAP stage III
|
6 (13)
|
6 (18)
|
0
|
0.106
|
CAD score
|
18±9
|
21±8
|
10±4
|
0.001
|
Note: The data are presented as mean value ± SD or frequencies (percentages)
Abbreviations: CAD score, Coronary Artery Disease score; IPF, idiopathic pulmonary fibrosis; GAP index, Gender, Age, Physiology index; MRC score, Medical Research Council score; Pys, pack per years
3.2. Lung function parameters of the study population
The lung function parameters of the study population are presented in Table 2. 52% (25/48) of the patients had restrictive pattern (TLC<70%). 6% (3/48) of the participants had obstructive pattern (FEV1/FVC<0.70). The mean DLCO/VA% for the study population was 51±16% predicted. Pulmonary hyperinflation was detected in IPF patients as the mean RV/TLC was above 40%.
Abnormal FOT and IOS values of AX, R5-R20 and Fres were found, indicating the presence of SAD in the study population.
Table 2. Lung function parameters of the study population and comparisons between genders
Parameters
|
Sum (N=48)
|
Males (n=33)
|
Females (n=15)
|
p-value
|
FVC, ml
|
2.7±0.8
|
2.9±0.8
|
2.1±0.6
|
<0.001
|
FVC%
|
83±22
|
80±22
|
89±20
|
0.171
|
FEV1, ml
|
2.2±0.7
|
2.5±0.7
|
1.8±0.4
|
<0.001
|
FEV1%
|
88±22
|
86±23
|
92±18
|
0.382
|
FEV1/FVC
|
82±7
|
81±7
|
84±6
|
0.147
|
FEV1/FVC<70%
|
3 (6)
|
3 (0.9)
|
0
|
0.315
|
PEF, ml
|
6.4±2.3
|
7±2.5
|
6.1±1.6
|
0.006
|
PEF%
|
90±28
|
89±29
|
92±25
|
0.756
|
FEF25-75, ml
|
2.3±1.2
|
2.7±1.3
|
2.2±0.8
|
0.161
|
FEF25-75%
|
94±44
|
94±44
|
94±47
|
0.948
|
FEF25-75<70%
|
15 (31)
|
9 (27)
|
6 (40)
|
0.289
|
TLC, ml
|
4.1±1.1
|
4.4±1.2
|
3.7±0.6
|
0.024
|
TLC%
|
70±17
|
67±17
|
77±14
|
0.067
|
TLC<70%
|
25 (52)
|
19 (58)
|
6 (40)
|
0.177
|
RV, ml
|
1.7±0.6
|
1.7±0.6
|
1.7±0.8
|
0.754
|
RV%
|
70±28
|
64±23
|
83±34
|
0.031
|
RV/TLC%
|
66±28
|
64±30
|
73±26
|
0.359
|
RV/TLC>40%
|
35 (72)
|
23 (69)
|
12 (80)
|
0.484
|
DLCO, ml/min/mmHg/L
|
4.1±1.4
|
4.4±1.4
|
3.5±1.2
|
0.072
|
DLCO/VA
|
51±16
|
50±17
|
50±16
|
0.966
|
DLCO<80%
|
48 (100)
|
33 (100)
|
15 (100)
|
0.978
|
FOT parameters
|
|
|
|
|
R5, kPa s L−1
|
0.3±0.1
|
0.28±0.09
|
0.34±0.09
|
0.038
|
R5 z-score
|
-0.2±1.2
|
-0.1±1.2
|
-0.6±1.2
|
0.132
|
R5%
|
100.1±34.6
|
105±35
|
89±30
|
0.146
|
R20 kPa, s L−1
|
0.26±0.06
|
0.24±0.05
|
0.29±0.06
|
0.100
|
R5-20, kPa s L−1
|
0.05±0.04
|
0.05±0.03
|
0.05±0.04
|
0.611
|
R5-20 >0.03 kPa s L−1
|
31 (65)
|
20 (60)
|
11 (73)
|
0.302
|
R5-20 z-score
|
0.4±1.6
|
0.4±1.6
|
0.4±1.7
|
0.950
|
X5, cmH₂O/L/s
|
-0.19±0.08
|
-0.16±0.07
|
-0.23±0.08
|
0.004
|
AX, KpA/L
|
1.2±0.8
|
1.2±0.8
|
1.5±0.7
|
0.133
|
AX>0.5 KpA/L
|
48 (100)
|
33 (100)
|
15 (100)
|
0.978
|
AX z-score
|
1.5±0.9
|
1.6±1
|
1.3±0.8
|
0.279
|
AX%
|
380±285
|
426±326
|
284±138
|
0.116
|
Fres, Hz
|
20±4.3
|
20±5
|
20±4
|
0.641
|
Fres >12 Hz
|
45 (94)
|
31 (94)
|
14 (93)
|
0.685
|
IOS parameters
|
|
|
|
|
R5, kPa s L−1
|
0.3±0.05
|
0.3±0.05
|
0.3±0.06
|
0.523
|
R5%
|
95±19
|
98±18
|
79±17
|
0.112
|
R20, kPa s L−1
|
0.25±0.05
|
0.2±0.04
|
0.3±0.08
|
0.141
|
R20%
|
90±18
|
91±18
|
83±25
|
0.516
|
R5-20, kPa s L−1
|
0.1±0.02
|
0.1±0.02
|
0.03±0.02
|
0.183
|
R5-R20 ≥0. 1 kPa s L−1
|
29 (60)
|
23 (70)
|
6 (40)
|
0.048
|
X5, kPa s L−1
|
-0.07±0.05
|
-0.06±0.03
|
-0.09±0.01
|
0.122
|
X5%
|
271±204
|
303±205
|
89±34
|
0.096
|
AX, kPa s L−1
|
0.38±0.22
|
0.39±0.23
|
0.32±0.12
|
0.661
|
AX>0.1 kPa s L−1
|
48 (100)
|
33 (100)
|
15 (100)
|
0.999
|
Fres, Hz
|
16±3
|
17±3
|
14±4
|
0.219
|
Fres> 12 Hz
|
43 (90)
|
29 (87)
|
14 (93
|
0.284
|
Note: The data are presented as mean value ± SD or frequencies (percentages)
Abbreviations: AX, area under the curve; between X5 Hz and cross-point of the curve to horizontal line; DLCO/VA, diffusing capacity for carbon monoxide per liter of lung volume corrected for alveolar transfer; DLCO/VA%, predicted diffusing capacity for carbon monoxide per liter of lung volume corrected for alveolar transfer; FEF25-75, forced mid-expiratory flow; FEF25-75%, predicted forced mid-expiratory flow; FEV1%, forced expiratory volume in 1 s of FVC test; FVC%, forced vital capacity; PEF, peak expiratory flow, PEF, predicted peak expiratory flow; R5, actual measure of resistance at frequency 5 Hz; R5%, predicted resistance at frequency 5 Hz; R20, actual measure of resistance at frequency 20 Hz; R20%, predicted resistance at frequency 20 Hz; R5–20, actual lung resistance when subtraction resistance at frequency 5 Hz minus resistance at frequency 20 Hz; RV, residual volume; RV%, predicted residual volume; TLC, total lung capacity; TLC%, predicted total lung capacity; X5%, predicted lung reactance at frequency 5 Hz; X5, actual measure of lung reactance at frequency 5 Hz.
Abnormal IOS values: R5 ≥ 0.5 kPa s L−1, R5% >150%, R20 ≥ 0.3 kPa s L−1, R20%>150%, Di5-20 ≥ 0.1 kPa s L−1, X5≤-0.1 kPa s L−1, Fres>12 Hz, AX ≥ 0.1 kPa s L−1. [14,15]
Normal Tremoflo values: R5: 0.2 – 0.4 kPA/L/s; R20: 0.2 – 0.4 KpA /L/s; R5-R20 ≤ 0.3 cmHO2.s/L; X5: -0.1 to -0.4 KpA/L/s; AX:< 0.5 KpA/L; AX z-score: <1.64; R5-R20 z-score: <1.64, R5 z-score <1.64 [16,17]
52% of patients with IPF exhibited reduced lung volumes, indicative of restrictive patterns, while a smaller subset (6%) presented with airflow limitation (FEV1/FVC<70%), suggesting obstructive conditions. Impaired gas exchange was noted across all individuals. Many participants also showed signs of small airway obstruction, evidenced by reduced mid-expiratory flow rates. Additionally, more than half of the participants demonstrated air trapping, reflected in elevated RV/TLC ratios.
The findings from both FOT and IOS further support these observations, with many participants exhibiting increased small airway resistance, reduced lung compliance, and heightened lung stiffness.
The R5-R20 parameter (in both FOT and IOS) measures resistance in the peripheral or small airways, and an elevated R5-R20 value is a direct marker of small airway dysfunction. In this study, a substantial 65% of participants had an abnormal R5-R20 value (greater than 0.03 kPa·s/L) in FOT, indicating increased small airway resistance. Similarly, in IOS, 60% of participants had R5-R20 values exceeding the abnormal threshold (0.1 kPa·s/L), underscoring widespread small airway dysfunction in this population.
Additionally, Fres proved to be highly sensitive in detecting lung stiffness and small airway obstruction, with 94% of participants in FOT and 90% in IOS showing values above 12 Hz. This further highlights the prevalence of small airway dysfunction and increased lung stiffness. While R5-R20 specifically assesses small airway resistance, Fres is more sensitive to broader changes in lung mechanics, encompassing both lung stiffness and small airway obstruction.
The AX parameter, which reflects reduced lung compliance and increased small airway resistance, was elevated in all participants, with 100% of participants in both FOT and IOS showing AX values above the abnormal threshold. In FOT, AX values exceeded 0.5 kPa·s/L in all participants, while in IOS, AX exceeded 0.1 kPa·s/L in all cases. This universal abnormality in AX highlights the severity and prevalence of small airway dysfunction and reduced lung compliance across the study group, reinforcing its role as a critical marker in detecting restrictive and obstructive lung pathologies.
The FEF25-75% parameter indicated SAD in 31% of participants, with values below 70%, making it a less sensitive marker for small airway dysfunction compared to AX, R5-R20 and Fres in FOT and IOS.
Lastly, the RV/TLC ratio proved particularly effective in identifying air trapping, with 72% of participants displaying ratios above 40%, a hallmark of obstructive lung diseases. However, it was slightly less sensitive than AX, Fres in detecting overall lung dysfunction.
Collectively, these findings point to a combination of restrictive and obstructive lung pathologies, with SAD being particularly prevalent in this cohort.
3.3. Correlations between CAD Score and lung function parameters: The role of reactance and lung stiffness in cardiovascular risk
There was a positive correlation (r =0.283, p = 0.045) between CAD score and AX (Figure 1). This suggests that higher AX (indicating increased lung stiffness or small airway obstruction) was associated with a higher CAD score, suggesting that patients with more severe lung stiffness or small airway obstruction were more likely to have coronary artery disease. There was a negative correlation (r =-0.314, p = 0.037) of CAD with X5, which relates to lung compliance (Figure 2). This suggests that as X5 becomes less negative (indicating better lung compliance or less stiffness), the CAD score decreases, meaning patients with less lung stiffness were less likely to have coronary artery disease.
3.4. Correlations Between GAP Score and MRC and Lung Function Parameters: The Role of Reactance and Lung Stiffness in IPF Severity
No significant correlations were observed between the MRC or GAP Index and the IOS or FOT parameters in this dataset, indicating that while SAD is prevalent in IPF patients, these oscillometry parameters may not correlate strongly with disease severity as measured by the MRC score and GAP index in this cohort.
3.5. Other correlations between lung function, dyspnea, severity and cardiovascular risk in IPF patients
Correlation analyses showed that there was a significant negative correlation (r = -0.522, p = 0.001) between DLCO% and MRC score, indicating that as DLCO% decreases, MRC score increase, which reflects worsening dyspnea as gas exchange becomes impaired. There was an even stronger negative correlation (r = -0.711, p < 0.001) between DLCO% and GAP score, reinforcing that as DLCO% decreases, the GAP index increases, reflecting worsening lung function and prognosis. There was a positive correlation (r = 0.474, p = 0.001) between GAP score and MRC, suggesting that as GAP index points increase, which indicates a worse prognosis in IPF, the MRC score also increases, indicating higher levels of dyspnea. There was a positive correlation (r = 0.255, p =0.046), suggesting that a higher GAP index (worse prognosis in IPF) is associated with a higher CAD score, indicating that coronary artery disease risk may increase as IPF severity worsens. There was a positive correlation (r = 0.428, p = 0.002) between RV/TLC and CAD score, indicating that as RV/TLC increases (reflecting more air trapping), the CAD score increases, which could indicate more coronary artery disease risk in patients with more severe lung dysfunction.
There was a positive correlation between age and GAP index (r = 0.381, p = 0.008), suggesting that older age at diagnosis is associated with a higher GAP index, indicating a worse prognosis. There was a positive correlation (r = 0.506, p < 0.001), indicating that older age at diagnosis is associated with a higher CAD score, reflecting that older individuals are more likely to have coronary artery disease.
There was no significant association between IOS/FOT parameters and IPF patients’ age, medication or smoking status.
3.6. Correlation Between DLCO% and FOT or IOS Parameters
There was a significant positive correlation between DLCO% and X5 (r = 0.377, p = 0.014), indicating that as DLCO% improves (reflecting better gas exchange), X5 becomes less negative. This suggests that better lung function was associated with improved lung compliance and less lung stiffness. Conversely, DLCO% showed a strong negative correlation with Fres (r = -0.705, p < 0.001), meaning that as gas exchange worsens, Fres increases, which is indicative of increased lung stiffness and small airway obstruction (Figure 3). Additionally, there was a significant negative correlation between DLCO% and AX (r = -0.502, p < 0.001), demonstrating that decreased gas exchange was associated with greater airway obstruction and worse lung compliance (Figure 4). Together, these findings showed that as lung function declines, stiffness and resistance in the airways increase.
3.7. Correlation Between Disease Duration and Lung Function Decline in FOT Parameters
The only statistically significant correlation found was between disease duration and X5, suggesting that as the disease progresses, lung compliance decreases, leading to increased stiffness. Other trends, such as increasing small airway obstruction (AX) and lung stiffness (Fres) with disease duration, were observed but did not reach statistical significance.