Characteristics of study subjects
A total of 138 eligible patients with IPF were enrolled in this study (Figure1). The median age of the IPF patient population was 63.0 years (IQR:59.0, 69.0). 73.2%(n=101) enrolled patients were male. Most of the patients (n=84, 60.9%) were smokers. Regarding comorbidities, there were 33(23.9%) patients complicated with emphysema, 42(30.4%) patients with hypertension, 24(17.4%) with diabetes, and 4(2.9%) with lung cancer. Among the total 138 IPF patients, 53 patients (38.4%) were diagnosed with spirometry-SAD.
Comparison of small airway function between IPF patients and controls
Before PSM, compared with the control group composed of 186 participants without IPF, IPF patients were more likely to be smokers(p<0.001). In addition, the differences in age (63.0y vs. 59.0y, p<0.001), sex (male:73.2% vs. 36.0%, p<0.001) and BMI (24.4kg/m2 vs. 23.6 kg/m2, p=0.032) were significant. After a 1:1 PSM analysis, 79 matched pairs were obtained, and all clinical characteristics were balanced between the two groups. The pulmonary function parameters differed significantly between IPF patients and controls after PSM practice. In the PSM groups, there were significant differences in FVC (75.9% vs. 85.3%, p=0.001) and DLco (69.2% vs. 94.1%, p<0.001) between the two groups (Table 1). IPF patients also showed lower FEV1/FVC, RV, and TLC, RV/TLC. In addition, compared with controls, IPF patients tended to manifest poorer functional status on indicators which reflect small airway function: MEF25(68.3% vs. 90.4%, p= 0.001), MMEF (69.5% vs. 86.8%, p<0.001). According to the criteria of SAD mentioned earlier, 30(38.0%) IPF patients were diagnosed with spirometry-SAD among 79 match pairs, while only 12(15.2%) controls were diagnosed with spirometry-SAD. The differences were statistically significant (p=0.001).
Table 1 Clinical characteristics and pulmonary function of patients with and without IPF
|
Before PSM
|
After PSM
|
IPF patients
(n=138)
|
Control
(n=186)
|
p
|
IPF patients
(n=79)
|
Control
(n=79)
|
p
|
Age, y
|
63.0(59.0,69.0)
|
59.0(50.0,66.3)
|
<0.001
|
62.4±9.0
|
62.1±9.0
|
0.806
|
Gender: male, n (%)
|
101(73.2)
|
67(36.0)
|
<0.001
|
44(55.7)
|
36(45.6)
|
0.203
|
Smoker, n (%)
|
84(60.9)
|
28(15.1)
|
<0.001
|
26(32.9)
|
26(32.9)
|
1.000
|
BMI, kg/m2
|
24.4±3.1
|
23.6±3.2
|
0.032
|
24.0±3.4
|
24.8±3.2
|
0.136
|
Pulmonary function test
|
|
|
|
|
|
|
FVC (% predicted)
|
76.7(59.2,87.5)
|
83.8(77.7,91.2)
|
<0.001
|
75.9±19.9
|
85.3±12.2
|
0.001
|
FEV1/FVC (% predicted)
|
104.2(97.2,110.4)
|
92.2(85.8,102.0)
|
<0.001
|
103.3(97.0,110.0)
|
91.5(86.0,110.8)
|
0.005
|
RV (% predicted)
|
88.8(78.3,110.5)
|
116.4(103.2,131.8)
|
<0.001
|
94.5(78.5,116.1)
|
114.7(99.7,132.2)
|
<0.001
|
TLC (% predicted)
|
79.6(68.3,95.6)
|
93.9(85.1,103.1)
|
<0.001
|
83.2±18.6
|
92.1±10.9
|
<0.001
|
RV/TLC (% predicted)
|
113.5(103.0,126.2)
|
122.0(113.6,132.7)
|
<0.001
|
115.7(103.0,128.0)
|
121.7(109.3,133.8)
|
0.046
|
DLco (% predicted)
|
64.4(47.4,78.9)
|
100.1(85.1,111.9)
|
<0.001
|
69.2±24.6
|
94.1±19.3
|
<0.001
|
MEF50 (% predicted)
|
87.9(64.6,116.0)
|
85.2(68.7,101.2)
|
0.268
|
92.8±34.7
|
90.7±26.4
|
0.675
|
MEF25 (% predicted)
|
70.4(45.8,101.6)
|
80.2(61.4,108.6)
|
0.011
|
68.3(43.1,95.5)
|
90.4(70.3,117.0)
|
<0.001
|
MMEF (% predicted)
|
70.8(51.8,99.1)
|
83.3(66.8,99.9)
|
0.003
|
69.5(52.8,99.0)
|
86.8(70.7,109.6)
|
0.001
|
MEF50(%)<65%(n,%)
|
34(24.6)
|
37(19.9)
|
0.307
|
18(22.8)
|
12(15.2)
|
0.224
|
MEF25(%)<65%(n,%)
|
60(43.5)
|
53(28.5)
|
0.005
|
34(43.0)
|
17(21.5)
|
0.004
|
MMEF(%)<65%(n,%)
|
57(41.3)
|
39(21.0)
|
<0.001
|
33(41.8)
|
13(16.5)
|
<0.001
|
Spirometry-SAD, n (%)
|
53(38.4)
|
39(21.0)
|
0.001
|
30(38.0)
|
12(15.2)
|
0.001
|
The data are presented as mean ± SD, median (interquartile range)
Abbreviations: IPF=Idiopathic pulmonary fibrosis, PSM=propensity score match, BMI=Body-mass index, FVC=forced vital capacity, FEV1=forced expiratory volume in 1 second, RV=residual volume, TLC total lung capacity, DLco=the diffusing capacity for carbon monoxide, MEF=maximum expiratory flow, MMEF=maximal mid-expiratory flow, SAD=small airway dysfunction.
Small airway function of patients of varying disease severity
To explore the small airway status of IPF patients with varying levels of disease severity, we divided them into three groups according to FVC decline degree: mild (FVC%>75%), moderate (55%≤FVC%≤75%), and severe(FVC%<55%). (Table2) With increasing severity of FVC, values of FEV1/FVC (p=0.004), RV (p<0.001), TLC (p<0.001), and RV/TLC (p<0.001) significantly decreased, showing progressive restrictive pulmonary ventilation disorder. The value of Rtot (p=0.005) and Reff (p=0.018) were higher in the more severe group. Lung diffusing capacity was also altered in three groups(p<0.001). More than half of patients with severe IPF had spirometry-SAD (n=16, 61.5%). Even in the mild disease group, 30.6% (n=22) of patients were diagnosed with spirometry-SAD. Furthermore, the more severe the FVC decline was, the worse the small airway function was. The incidences of SAD between mild, moderate, and severe groups were significantly different(30.6% vs. 37.5% vs. 61.5, p=0.021). Pairwise contrasts of the three group is showed in Figure 2.
Table 2 Pulmonary function data of different disease severity level
|
Mild
FVC%>75%
(n=72)
|
Moderate
55%≤FVC%≤75%
(n=40)
|
Severe
FVC%<55%
(n=26)
|
p
|
Age, y
|
63.0(59.0,69.0)
|
63.5(58.0,71.5)
|
62.0(57.0,66.0)
|
0.732
|
Gender: male, n (%)
|
59(81.9)
|
27(67.5)
|
15(57.7)
|
0.036
|
Smoker, n (%)
|
48(66.7)
|
21(52.5)
|
15(57.7)
|
0.316
|
BMI, kg/m2
|
24.0±3.0
|
24.8±3.0
|
24.8±3.6
|
0.254
|
Pulmonary function test
|
|
|
|
|
FVC (% predicted)
|
87.3(82.4,100.3)
|
63.6(59.7,68.7)
|
51.2(45.3,54.0)
|
<0.001
|
FEV1/FVC (% predicted)
|
100.6±9.7
|
106.4±9.3
|
106.0±9.2
|
0.004
|
RV (% predicted)
|
103.6(84.2,117.1)
|
81.5(67.8,99.2)
|
78.5(66.5,104.7)
|
<0.001
|
TLC (% predicted)
|
94.8(83.8,103.3)
|
70.6(63.8,77.4)
|
60.9(55.1,71.4)
|
<0.001
|
RV/TLC (% predicted)
|
109.2(100.0,116.0)
|
114.4(103.6,127.2)
|
126.0(115.9,146.9)
|
<0.001
|
Rtot
|
112.7(95.9,143.0)
|
139.0(109.1,175.8)
|
148.6(106.5,172.8)
|
0.005
|
Reff
|
101.2(84.7,127.9)
|
120.4(86.6,155.4)
|
134.1(97.8,167.1)
|
0.018
|
DLco (% predicted)
|
78.2±22.1
|
56.7±14.8
|
43.5±11.1
|
<0.001
|
MEF50 (% predicted)
|
98.2(69.8,123.1)
|
87.4(71.3,115.9)
|
52.4(48.6,83.0)
|
0.001
|
MEF25 (% predicted)
|
74.2(48.9,107.1)
|
71.9(42.5,102.4)
|
61.3(34.6,71.2)
|
0.062
|
MMEF (% predicted)
|
77.2(54.5,104.7)
|
74.6(54.6,99.4)
|
44.6(38.9,73.7)
|
0.002
|
MEF50(%)<65%(n,%)
|
12(16.7)
|
8(20.0)
|
14(53.8)
|
0.001
|
MEF25(%)<65%(n,%)
|
28(38.9)
|
16(40.0)
|
16(61.5)
|
0.119
|
MMEF(%)<65%(n,%)
|
23(31.9)
|
16(40.0)
|
18(69.2)
|
0.004
|
Spirometry-SAD, n (%)
|
22(30.6)
|
15(37.5)
|
16(61.5)
|
0.021
|
The data are presented as mean ± SD, median (interquartile range)
Abbreviations: BMI = Body-mass index; FVC = forced vital capacity, FEV1 = forced expiratory volume in 1 second, RV = residual volume, TLC = total lung capacity, Rtot = total airway resistance, Reff = effective airway resistance, DLco = the diffusing capacity for carbon monoxide, MEF = maximum expiratory flow, MMEF = maximal mid-expiratory flow, SAD = small airway dysfunction.
The diagnostic value of IOS in identifing SAD
Although spirometry was widely used to define SAD, IOS may also be an effort-independent method that is convenient and more sensitive to detecting SAD in IPF patients. In IPF patients, the values of MEF50 and MMEF were significantly correlated with many IOS parameters (Figure 3), although the correlation was very limited.
In our study, there were 60 patients available for the advanced analysis of IOS. IOS-SAD ratios were 68.3%, which was higher than those of SAD diagnosed by spirometry (n=25, 41.7%). The Cohens kappa agreement test was used to check the agreement and concordance between IOS and spirometry in the diagnosis of SAD. As described in Table 3, the concurrence of spirometry and IOS was generally low in the assessment of SAD (Kappa=-0.068, p=0.542). The results of McNemar's test remained consistent (p=0.009). Taking spirometry as standard, the sensitivity and specificity of IOS were 64.0% and 28.9%, respectively. The area under the curve (AUC) was 0.463 (0.313-0.612). For 35 patients diagnosed with non-SAD by spirometry, 25 patients were diagnosed with SAD by IOS. For patients with normal FVC (FVC%>80%), the number of IOS-only SAD was twice that of spirometry. Moreover, for patients with decreased FVC(FVC%<80%), IOS-only SAD were three to four time more than spirometry-only SAD. IOS criteria exhibited enhanced sensitivity to spirometry, especially in those patients with FVC less than 80% predicted(Table 3). The diagnosis values of other IOS parameters except for R5-R20 were listed in supplementary statement table 1(Additional file 1).
Table 3 The consistency test between spirometry and IOS parameters. (N=60)
|
Impulse Oscillometry
|
Spirometry
|
Kappa test
|
McNemar’s test
|
Negative (-)
|
Positive (+)
|
Kappa value
|
p
|
p
|
Patients with normal FVC
(n=11)
|
Negative (-)
|
6
|
4
|
-0.125
|
0.537
|
0.388
|
Positive (+)
|
8
|
3
|
Patients with decreased FVC
(n=39)
|
Negative (-)
|
4
|
5
|
-0.083
|
0.519
|
0.017
|
Positive (+)
|
17
|
13
|
Total Patients
|
Negative (-)
|
10
|
9
|
-0.068
|
0.542
|
0.009
|
(n=60)
|
Positive (+)
|
25
|
16
|
Abbreviations: FVC = forced vital capacity
IOS parameters related to peripheral airways resistance, including R5-R20[0.08kPa/(L/S) vs. 0.11kPa/(L/S), p=0.019], X5[-0.11kPa/(L/S) vs. -0.17 kPa/(L/S), p<0.001], Fres (15.72Hz vs. 17.59 Hz, p=0.001), were all different between IPF patients and controls. According to diagnostic criteria of SAD based on IOS(IOS-SAD), 41(68.3%) of 60 IPF patients were with R5-R20 greater than 0.07 kPa/(L/S). Compared with controls, IPF patients were more susceptible to IOS-SAD (52.2% vs. 68.3%, p=0.090) (Table 4).
Table 4 Differences in IOS parameters between IPF patients and controls.
IOS parameter
|
Control
(n=46)
|
IPF patients
(n=60)
|
p
|
R5, [kPa/(L/S)]
|
0.36(0.27,0.44)
|
0.35(0.26,0.44)
|
0.779
|
R20, [kPa/(L/S)]
|
0.28(0.22,0.34)
|
0.23(0.19,0.31)
|
0.059
|
R5-R20, [kPa/(L/S)]
|
0.08(0.04,0.12)
|
0.11(0.05,0.15)
|
0.019
|
X5, [kPa/(L/S)]
|
-0.11(-0.14,-0.08)
|
-0.17(-0.25,-0.12)
|
<0.001
|
Fres, (Hz)
|
15.72(12.00,18.5)
|
17.59(15.4,21.5)
|
0.001
|
R5–R20>0.07kPa/(L/S)
|
24(52.2)
|
41(68.3)
|
0.090
|
X5<−0.12 kPa/(L/S)
|
19(41.3)
|
43(71.7)
|
0.002
|
Fres> 14.14Hz
|
29(63.0)
|
56(93.3)
|
<0.001
|
The data are presented as median (interquartile range)
Abbreviations: IPF=Idiopathic pulmonary fibrosis, IOS=impulse oscillometry, R5=resistances at 5 Hz, R20=resistances at 20 Hz, R5-R20=resistances at 5 and 20 Hz, X5=reactance at 5 Hz, Fres=resonant frequency in Hz.
The IOS parameters also altered in three FVC groups(Figure 4). The value of R5-R20, X5, Fres were different between mild and moderate groups(R5-R20:p=0.002, X5:p=0.065, Fres:p=0.001).
Clinical features and PFT parameters of IPF patients with or without SAD
Compared with patients without spirometry-SAD, patients with spirometry-SAD were more likely to be smokers (58.8% vs. 64.2%, p=0.533). More patients with SAD were in GAP stage Ⅱ-Ⅲ than patients without SAD (23.5% vs. 39.6%, p=0.044). There were significant differences in FVC (80.0% vs. 68.3%, p=0.034) and RV/TLC (110.5% vs.116.0%, p=0.011) as well as Rtot (118.2% vs. 145.1%, p=0.003) and Reff (102.1% vs. 131.8%, p=0.001) between IPF patients with SAD and without SAD.(Table 5)
Table 5 Clinical characteristic and pulmonary function differences between patients with SAD and without SAD.
|
Non-SAD
(n=85)
|
SAD
(n=53)
|
p
|
Age, y
|
64.6±8.6
|
62.6±8.7
|
0.140
|
Gender:male, n (%)
|
65(76.5)
|
36(67.9)
|
0.270
|
Smoker, n (%)
|
50(58.8)
|
34(64.2)
|
0.533
|
BMI, kg/m2
|
24.2±3.1
|
24.8±3.1
|
0.471
|
Emphysema, n (%)
|
19(22.4)
|
14(26.4)
|
0.586
|
GAP score
|
|
|
0.044
|
GAP Ⅰ, n (%)
|
65(76.5)
|
32(60.4)
|
|
GAP Ⅱ-Ⅲ, n (%)
|
20(23.5)
|
21(39.6)
|
|
Pulmonary function test
|
|
|
|
FVC (% predicted)
|
80.0(64.0,90.3)
|
68.3(55.1,87.3)
|
0.034
|
FEV1/FVC (% predicted)
|
107.7±7.7
|
96.3±8.8
|
<0.001
|
RV (% predicted)
|
88.5(78.5,110.0)
|
100.4(81.2,113.9)
|
0.281
|
TLC (% predicted)
|
82.6±18.8
|
79.8±18.3
|
0.397
|
RV/TLC (% predicted)
|
110.5(100.2,120.4)
|
116.0(106.8,127.6)
|
0.011
|
Rtot (% predicted)
|
118.2(91.6,147.2)
|
145.1(105.3,172.8)
|
0.003
|
Reff (% predicted)
|
102.1(78.3,131.2)
|
131.8(95.8,165.5)
|
0.001
|
DLco (% predicted)
|
65.9(52.2,80.0)
|
62.5(44.6,75.8)
|
0.219
|
The data are presented as mean ± SD, median (interquartile range)
Abbreviations: SAD=small airway dysfunction, BMI=Body-mass index; FVC = forced vital capacity, FEV1 = forced expiratory volume in 1 second, RV = residual volume, TLC = total lung capacity, Rtot = total airway resistance, Reff = effective airway resistance, DLco = the diffusing capacity for carbon monoxide.
Risk factor analysis for IPF patients with SAD
Table 6 summarizes the ORs of SAD. Multiple logistics regression revealed that smoker was an independent risk factor for SAD (p=0.049). The risk of SAD was approximately threefold higher in those smokers than in non-smokers.
Table 6 Risk factors for SAD in IPF patients. (N=138)
|
Univariate
|
|
Multivariate
|
OR
|
95% CI
|
p
|
OR
|
95% CI
|
p
|
Age, y
|
0.973
|
0.934-1.013
|
0.187
|
|
0.981
|
0.939-1.024
|
0.374
|
Male, n (%)
|
0.652
|
0.303-1.399
|
0.272
|
0.296
|
0.083-1.061
|
0.062
|
Smoker, n (%)
|
1.253
|
0.617-2.544
|
0.533
|
3.354
|
1.005-11.193
|
0.049
|
BMI, kg/m2
|
1.061
|
0.949-1.186
|
0.302
|
1.049
|
0.931-1.182
|
0.433
|
FVC (% predicted)
|
0.980
|
0.963-0.998
|
0.026
|
0.978
|
0.956-1.001
|
0.056
|
DLco (% predicted)
|
0.993
|
0.979-1.007
|
0.343
|
1.008
|
0.989-1.028
|
0.413
|
Abbreviations: OR = odds ratio, CI = confidence interval, FVC = forced vital capacity, DLco = the diffusing capacity for carbon monoxide
Small airway function in different GAP groups
To explore the relationship between SAD and prognosis, we divided study subjects into two groups according to the GAP stage, namely GAP Ⅰ and GAP Ⅱ-Ⅲ (Table 7). Supplementary statement table 2(Additional file 1) showed GAP score distribution. In our study, 97 patients were in GAP stage Ⅰ, and 41 were in GAP stage Ⅱ-Ⅲ. With increasing stage of GAP, values of FVC (84.5% vs. 55.1%, p<0.001), FEV1/FVC (102.1% vs. 106.2%, p=0.028), RV (96.4% vs. 78.5%, p<0.001), TLC (88.0% vs. 66.0%, p<0.001), DLco (73.1% vs. 41.3%, p<0.001) significantly decreased, while RV/TLC (111.3% vs. 120.4%, p=0.019) increased. Apart from that, patients in GAP stage Ⅱ-Ⅲ showed lower MMEF50 (91.0% vs. 76.8%, p=0.040), MEF25 (72.2% vs. 62.1%,p=0.148) and MMEF (71.1% vs. 62.1%, p=0.070), indicating descending small airway function. Among patients in GAP stage Ⅱ-Ⅲ, more than half patients had spirometry-SAD (n=21, 51.2%). Compared with GAP Ⅰ patients, more GAP Ⅱ-Ⅲ patients were diagnosed with spiratory-SAD (33.0% vs. 51.2%), and the difference in SAD proportion between the two groups reached statistical significance(p=0.044). As showed in Supplementary statement table 3(Additional file 3), patients in GAP stage Ⅱ-Ⅲ showed worse R5-R20, X5 and Fres compared with controls and patients in GAP stage Ⅰ (p=0.008; p<0.001; p=0.001).
Table 7 Differences in spirometry parameters between GAP subgroups. (N=138)
Variable
|
GAP points
|
GAP Stage I
(n=97)
|
GAP Stage II-III
(n=41)
|
p
|
Gender, n (%)
|
|
|
|
0.676
|
Female
|
0
|
27(27.8)
|
10(24.4)
|
|
Male
|
1
|
70(72.2)
|
31(75.6)
|
|
Age(y), n (%)
|
|
|
|
<0.001
|
≤60
|
0
|
44(45.4)
|
6(14.6)
|
|
61-65
|
1
|
23(23.7)
|
8(19.5)
|
|
>65
|
2
|
30(30.9)
|
27(65.9)
|
|
Physiology
|
|
|
|
|
FVC (% predicted)
|
|
|
|
<0.001
|
≥75
|
0
|
70(72.2)
|
2(4.9)
|
|
50-75
|
1
|
26(26.8)
|
25(61.0)
|
|
<50
|
2
|
1(1.0)
|
14(34.1)
|
|
DLco (% predicted)
|
|
|
|
<0.001
|
>55
|
0
|
81(83.5)
|
8(19.5)
|
|
36-55
|
1
|
16(16.5)
|
24(58.5)
|
|
≤50
|
2
|
0(0.0)
|
9(22.0)
|
|
Pulmonary function test
|
|
|
|
FVC (% predicted)
|
|
84.5(68.1,96.8)
|
55.1(47.4,67.1)
|
<0.001
|
FEV1/FVC (% predicted)
|
|
102.1±9.1
|
106.2±11.0
|
0.028
|
RV (% predicted)
|
|
96.4(81.5,115.4)
|
78.5(64.5,91.3)
|
<0.001
|
TLC (% predicted)
|
|
88.0±16.2
|
66.0±14.4
|
<0.001
|
RV/TLC (% predicted)
|
|
111.3(102.4,122.8)
|
120.4(105.6,133.0)
|
0.019
|
DLco (% predicted)
|
|
73.1(58.3,86.1)
|
41.3(35.4,51.1)
|
<0.001
|
MEF50(% predicted)
|
|
91.0(69.9,116.3)
|
76.8(47.2,119.1)
|
0.040
|
MEF25(% predicted)
|
|
72.2(49.5,102.9)
|
62.1(37.2,91.5)
|
0.148
|
MMEF (%predicted)
|
|
71.1(54.9,98.2)
|
62.1(39.4,99.4)
|
0.070
|
MEF50(%)<65%(n,%)
|
|
16(16.5)
|
18(43.9)
|
0.001
|
MEF25(%)<65%(n,%)
|
|
39(40.2)
|
21(51.2)
|
0.233
|
MMEF(%)<65%(n,%)
|
|
36(37.1)
|
21(51.2)
|
0.124
|
Spirometry-SAD, n (%)
|
|
32(33.0)
|
21(51.2)
|
0.044
|
The data are presented as mean ± SD, median (interquartile range)
Abbreviations: FVC=forced vital capacity, FEV1=forced expiratory volume in 1 second, RV=residual volume, TLC=total lung capacity, DLco=the diffusing capacity for carbon monoxide, MEF =maximum expiratory flow, MMEF =maximal mid-expiratory flow, SAD=small airway dysfunction.