Correlation of respiratory oscillometry with CT image analysis in a prospective cohort of idiopathic pulmonary fibrosis.

BACKGROUND: Markers of idiopathic pulmonary fibrosis (IPF) severity are based on measurements of forced vital capacity (FVC), diffusing capacity (DLCO) and CT. The pulmonary vessel volume (PVV) is a novel quantitative and independent prognostic structural indicator derived from automated CT analysis. The current prospective cross-sectional study investigated whether respiratory oscillometry provides complementary data to pulmonary function tests (PFTs) and is correlated with PVV.
METHODS: From September 2019 to March 2020, we enrolled 89 patients with IPF diagnosed according to international guidelines. We performed standard spectral (5-37 Hz) and novel intrabreath tracking (10 Hz) oscillometry followed by PFTs. Patients were characterised with the gender-age-physiology (GAP) score. CT images within 6 months of oscillometry were analysed in a subgroup (26 patients) using automated lung texture analysis. Correlations between PFTs, oscillometry and imaging variables were investigated using different regression models.
FINDINGS: The cohort (29F/60M; age=71.7±7.8 years) had mild IPF (%FVC=70±17, %DLCO=62±17). Spectral oscillometry revealed normal respiratory resistance, low reactance, especially during inspiration at 5 Hz (X5in), elevated reactance area and resonance frequency. Intrabreath oscillometry identified markedly low reactance at end-inspiration (XeI). XeI and X5in strongly correlated with FVC (r2=0.499 and 0.435) while XeI was highly (p=0.004) and uniquely correlated with the GAP score. XeI and PVV exhibited the strongest structural-functional relationship (r2=0.690), which remained significant after adjusting for %FVC, %DLCO and GAP score.
INTERPRETATION: XeI is an independent marker of IPF severity that offers additional information to standard PFTs. The data provide a cogent rationale for adding oscillometry in IPF assessment. © Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

The pulmonary function test parameters of forced vital capacity and diffusing capacity are currently used as metrics of disease severity in idiopathic pulmonary fibrosis (IPF) and are key factors in the gender-age-physiology IPF score.

Our study revealed that respiratory oscillometry, in particular the intrabreath reactance assessment at end-inspiration (XeI) was uniquely correlated with the gender-age-physiology score and strongly associated with pulmonary vessel volume, a new quantitative metric derived from automated CT image analysis.

Our data suggest that the novel intrabreath oscillometry provides additional information to standard pulmonary function tests and may be a simple surrogate metric of CT derived structural indices in IPF.

Introduction

Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease (ILD) of unknown aetiology. Affecting older adults, it has an exceedingly poor prognosis with a 5-year survival of 45%1–3 and confounded by the fact that it is often recognised late due to the non-specific nature of the clinical presentation. The joint European Respiratory Society (ERS)/American Thoracic Society (ATS)/Japanese Respiratory Society (JRS)/Latin American Thoracic Association (ALAT) guidelines recommend that IPF be diagnosed by consensus following a multidisciplinary review of the clinical presentation, pulmonary function tests (PFTs), high resolution CT and serology with/without surgical biopsy.2 4

Markers of IPF severity are based on a combination of demographics, imaging and physiological measurements that include forced vital capacity (FVC), diffusing capacity (DLCO) and hypoxemia during the 6 min walk (6MW) test. The GAP (gender, age, physiology) score is preferentially weighted by FVC and DLCO.5 CT imaging analyses are increasingly focused on developing quantitative metrics of disease severity and fibrosis. A computer-aided lung informatics for pathology evaluation and rating (CALIPER) tool has identified pulmonary vessel volume (PVV) to be the most predictive of IPF progression and mortality that is independent of PFT parameters.6–9

Oscillometry is a PFT that is highly sensitive to the lung periphery, a region not well assessed by conventional PFTs. It measures respiratory impedance (Zrs), expressed by resistance (Rrs) and reactance (Xrs) grossly reflecting the dynamic resistive and elastic properties, respectively, of the respiratory system.10 11 The difference between the mean inspiratory (X5in) and expiratory (X5ex) reactance at 5 Hz, was reported to distinguish ILD from asthma and chronic obstructive lung disease12 and combined pulmonary fibrosis-emphysema (CPFE) from ILD-only.13 14 Two recent studies found that whole-breath reactance at 5 Hz (X5) and X5in correlated with spirometry and the GAP score,14 and that the reactance values are correlated with CT metrics of structural abnormalities.15 While these data suggest that oscillometry could be useful to assess disease severity, it is not known whether oscillometry, a more easily applied diagnostic tool, can be used in place of or whether it provides additional information beyond spirometry.

The goals of the current study are to characterise IPF with the usual (spectral) and intrabreath oscillometry16 17 and to determine if oscillometry can provide additional information beyond conventional PFT and CT imaging.

Methods

Only patients (n=89) who met the ATS/ERS/JRS/ALAT criteria for IPF2 were eligible and prospectively recruited from the ILD Clinic at UHN, a tertiary referral centre. Written informed consent was obtained from participants prior to oscillometry; conventional PFTs were conducted as part of usual clinical care. We excluded patients with CPFE or who were already enrolled in a clinical trial due to existing protocol guidelines. The diagnosis of CPFE was made in a patient with a compatible history of smoking, evidence of emphysema on CT imaging, spirometry and clinical review by a multidisciplinary team. Healthy subjects were recruited from the UHN PFT Laboratory Biological control procedures and an unrelated study with oscillometry to screen for lung disease.

Oscillometry was conducted according to ERS guidelines11 following published quality control standards,18 using the tremoflo (model C-100, Thorasys, Montreal, Canada) device. Following spectral oscillometry with standard 5–37 Hz multifrequency measurements, monofrequency (10 Hz) recordings of 30 s were conducted for intrabreath analysis of Zrs at a 0.1 s temporal resolution as previously described,17 by using the research modality of the tremoflo device. Spirometry, body plethysmography and DLCO were performed after oscillometry, following international guidelines.19–21 All testing was conducted at the Toronto General Pulmonary Function Laboratory.

Anthropometric and clinical data were extracted from patients’ electronic medical records. Available high resolution CT scans of the chest done within 6 months of paired oscillometry—conventional PFTs were analysed by a fellowship trained thoracic radiologist using Lung Texture Analysis (Imbio, Minneapolis, Minnesota, USA) which is based on the CALIPER technology. In brief, lung texture analysis is a machine learning programme that automatically segments the lung from the chest and removes the central airways and central vasculature. It is trained to classify each voxel of lung as normal, hyperlucent, ground glass, reticular or honeycombing texture using a large previously described dataset.22 The parameters such as the ILD extent, which is the sum of volumes of ground glass, reticular and honeycombing volumes and the PVV, volume of vessel related structures, are expressed a per cent of the total lung capacity as defined by CT (TLCCT).6

Statistical analyses were performed using RStudio (The R Foundation, Boston, Massachusetts, USA). The absolute and % predicted values are presented as mean±SD or median and IQR, as appropriate. Correlations between PFTs, oscillometry and imaging variables were investigated using linear, polynomial with 2 degrees and exponential regression. The coefficient of determination R2 was used to compare goodness-of-fit in different models. Multivariate linear regression models with natural logarithm of PVV or ILD extent as outcomes were used to assess the correlation with XeI.

The study was conducted without patient and public involvement beyond peer-reviewed of the funding proposal.

Results

Demographics and conventional PFT

From 6 September 2019 to 12 March 2020, 120 patients with IPF were identified from the UHN ILD clinic. Seven patients were excluded as they were in a clinical drug trial and nine declined enrolment. Of the 104 consented patients, 10 were subsequently excluded when clinical review identified mixed lung diseases that did not meet inclusion criteria, for example, CPFE. Three patients did not complete oscillometry testing for various reasons (did not feel well; late for appointments). Two additional patients were excluded when post-test review of the PFT and/or oscillometry data failed to meet quality control.18

The remaining 89 participants were predominantly older males (67%) and overweight (mean body mass index=27.2±4.2; table 1). The group was similar to previously reported study cohorts of IPF.5 14 23 All patients performed paired oscillometry and spirometry during the same visit to the PFT laboratory. Fifteen patients of the 89 were unable to complete lung volume and DLCO measurements due to dyspnoea and/or fatigue. Similarly, four patients could not complete the 6MW test. The study group had mild to moderate impairment in conventional pulmonary function measures (mean % predicted FVC=69.6±17.0; % predicted DLCO=62.3±16.6). Six of the 85 patients who completed the 6MW test did so while on supplemental O2; the mean distance was normal with minimum O2 saturations being 87.7%±7.1%. There was no difference in anthropometric or PFT parameters between the overall cohort and the subgroup of 26 patients in whom CT scans were available (table 1). Thirty-eight per cent of the cohort were non-smokers and the others were equally distributed as ex-smokers of greater and less than 20 pack years (table 1). One-way analysis of variance revealed no statistical differences among the three groups with respect to the standard PFTs nor oscillometry with the exception of higher absolute FVC in the >20 pack-year group (online supplemental table). This was likely due to differences in height as %predicted FVC was similar among all groups. Of the oscillometry metrics, only the reactance values were different, but paradoxically, the worse values (more negative X5 and XeE, and higher AX) were observed in the non-smokers (online supplemental table S2).

Table 1

Summary of patient demographics, PFTs and 6MW parameters

	Overall cohort (n=89)		Subgroup with chest CT(n=26)
Sex (M/F)	60/29		18/8
Age (years)	71.7±7.8		72.9±6.3
Height (cm)	167.4±9.6		167.2±8.7
Weight (kg)	76.3±14.4		74.2±11.2
BMI (kg/m2)	27.2±4.2		26.5±3.4
Smoking Status
Never smoked	34		8
Ex-smoker <20 pk years	28		10
Ex-smoker ≥ 20 pk years	27		8
	Absolute values	% predicted	Absolute values	% predicted
FVC (L)	2.6±0.8	70±17	2.6±0.8	71±16
FEV1(L)	2.1±0.6	74±18	2.1±0.7	76±18
FEV1/FVC (%)	80.2±5.5	106±7	80.3±6.2	106±9
TLC (L)	4.3±1.0*	70±13*	4.8±1.0	72±13
RV (L)	1.6±0.4*	71±15*	1.7±0.5	74±19
RV/TLC (%)	38.6±7.2*	92±16*	39.2±9.1	91±21
DLCO (mL/min/mm Hg)	10.9±3.8*	62±17*	11.0±3.6	62±15
6MW min O2 sat (%)	87.7±7.1†	--	89.0±5.7	--
6MW distance (m)	453.8±101.4†	103±22†	460.1±82.5	107±18

Data are shown as mean±SD. The overall cohort and the subgroup were similar.

*n=74.

†6 patients walked with supplemental O2.

BMI, body mass index; DLCO, diffusing capacity for carbon monoxide; FEV1, forced expiratory volume in one second; FVC, forced vital capacity; 6MW, 6-min walk; PFT, pulmonary function test; RV, residual volume; TLC, total lung capacity.

Oscillometry of the IPF cohort

The spectral oscillometry graphs from an patient with IPF (figure 1A) and a healthy subject (figure 1D) showed alterations in oscillation mechanics in IPF with marked changes in the reactance versus frequency plots, with a rightward shift of Xrs resulting in elevated resonance frequency (Fres) and reactance area (AX). The intrabreath impedance-volume (figure 1B, E) and impedance-flow loops (figure 1C, F) from the same two subjects reflect an overall decrease in Xrs with an enhanced difference between end-expiratory (XeE) and end-inspiratory (XeI) values in IPF. This is manifested in a larger area of the Xrs vs flow loop and a distinct clockwise rotation.

Figure 1

Examples in a patient with IPF (A–C) and a healthy subject (D–F) of multifrequency (spectral) oscillometry (A,D) and intrabreath oscillometry at 10 Hz (B, C, E, F). Respiratory resistance (Rrs) and reactance (Xrs) values are indicated by closed and open circles, respectively. Fres denotes resonance frequency and AX denotes reactance area between 5 Hz and Fres. XeE and XeI, respectively, mark the end-expiratory and end-inspiratory reactance values; The solid arrows mark the beginning of inspiration and dotted arrows indicate beginning of expiration. IPF, idiopathic pulmonary fibrosis.

Comparison of the oscillometry for the whole cohort and the subgroup with chest CT measurements found the two groups to be similar (table 2). While Rrs at the lowest frequency (R5) is only slightly higher than its predicted value in the IPF group, this difference disappears by 19 Hz; this results in a considerably higher frequency dependence (R5–19) than that predicted and observed in the group of healthy controls (online supplemental table S1). X5 and its inspiratory and expiratory mean values (X5in and X5ex, respectively) were much lower than predicted (with a negligible difference between the latter). As a consequence, the reactance indices, AX and Fres were significantly higher than both the predicted values24 and the healthy group (online supplemental table S1).

Table 2

Summary of spectral and Intrabreath oscillometry measurements

	Total cohort (n=89)		Subgroup with chest CT (n=26)
	Absolute values	% predicted	Absolute values	% predicted
R5 (cm H2O s/L)	3.41 (1.26)	114 (36)	3.13 (0.72)	110 (15)
R19 (cm H2O s/L)	2.75 (0.96)	99 (30)	2.63 (0.64)	96.99 (20)
R5–19 (cm H2O s/L)	0.62 (0.63)	245 (409)	0.61 (0.45)	253 (280)
X5 (cm H2O s/L)	−2.05 (1.15)	172 (88)	−1.85 (1.17)	159 (47)
X5in (cm H2O s/L)	−2.10 (1.09)		−2.06 (1.18)
X5ex (cm H2O s/L)	−1.96 (1.40)		−1.61 (1.18)
ΔX5 (cm H2O s/L)	−0.25 (0.50)		−0.30 (0.37)
AX (cm H2O/L)	12.80 (12.80)	413 (408)	10.30 (11.70)	365 (229)
Fres (Hz)	21.29 (5.0)	167 (48)	20.10 (6.30)	168 (37)
ReE (cm H2O s/L)	2.98 (1.20)		2.64 (0.96)
ReI (cm H2O s/L)	2.17 (0.76)		2.07 (0.53)
ReE-ReI (cm H2O s/L)	0.69 (0.82)		0.50 (0.82)
XeE (cm H2O s/L)	−0.24 (0.64)		−0.15 (0.60)
XeI (cm H2O s/L)	−0.66 (0.74)		−0.62 (0.63)
XeE-XeI (cm H2O s/L)	0.38 (0.46)		0.47 (0.51)

Data are shown as median (IQR).

AX, reactance area between 5 Hz and Fres; ex, during expiration phase; Fres, resonance frequency; in, during inspiration phase; R5, resistance at 5 Hz; R19, resistance at 19 Hz, R5–19: R5–R19; ReE, resistance at end-expiration; ReI, resistance at end-inspiration; X5, reactance at 5 Hz; XeE, reactance at end-expiration; XeI, reactance at end-inspiration; ΔX5, difference of X5in and X5ex.

The intrabreath metrics reflected slightly elevated Rrs at resting lung volume (ReE) as compared with the healthy controls (online supplemental table S1); however, this difference disappeared at end-inspiration. Accordingly, the tidal change in Rrs, (ReE-ReI), was higher in the IPF group and more pronounced when the changes were normalised by VT. The zero-flow Xrs metrics (XeE and XeI) and the tidal changes (XeE-XeI and (XeE-XeI)/VT) were significantly different from the healthy group (online supplemental table S1). The contrast between the subtle alterations in Rrs and the marked changes in Xrs in IPF is also illustrated by the rank order plots of ReI and XeI; while the IPF and healthy subjects are intermingled in the whole range of ReI values, the IPF group was distinctly clustered at low values of XeI (figure 2).

Figure 2

Rank order plot of end-inspiratory reactance (XeI (A)) and resistance (ReI (B)) measured in patients with IPF (n=89; grey bars) and healthy subjects (n=45; cyan bars). IPF, idiopathic pulmonary fibrosis.

We had complete data to calculate the GAP score in 74 patients (table 3). Of these, 48 patients were in GAP Stage I, 22 in GAP Stage II and 5 were in GAP Stage III. No differences in age or sex were found among the three groups, while those in Stage III had significantly lower body mass index compared with GAP Stages I and II. Smoking status was different among the groups. We also observed an elevation in XeE-XeI and XeE-XeI/VT with increasing GAP stage although this was not statistically significant. Comparison of the oscillometry among the three GAP categories was notable for the significant progressive decrease in XeI with increasing GAP stages, accompanied by no change in XeE. Of all the oscillometry parameters, XeI was the only metric that was significantly different among the GAP stages (p<0.004).

Table 3

Patient demographics and oscillometry parameters according to the GAP score

	GAP I (n=48)	Gap II (n=21)	Gap III (n=5)	P value
Sex (M/F)	31/17	15/6	5/0	0.26
Age (yrs)	71.2±8.8	70.5±4.9	70.7±8.4	0.94
BMI (kg/m2)	27.2±4.0	27.7±3.6	21.2±3.9	0.004
Smoking status
Never smoked	17	6	2	<0.001*
Ex-smoker <20 pk years	19	5	2	<0.001*
Ex-smoker ≥20 pk years	12	10	1	<0.001*
6MW distance (% predicted)	111±17α	98±25β	86±20	0.005
R5 (cm H2O s/L)	3.8±1.2	3.8±1.3	2.7±0.7	0.15
R5-19 (cm H2O s/L)	0.8±0.7	0.7±0.6	0.8±0.5	0.76
X5 (cm H2O s/L)	−2.1±1.1	−2.2±0.7	−2.2±0.3	0.93
X5in (cm H2O s/L)	−2.1±1.0	−2.3±0.8	−2.5±0.5	0.67
X5ex (cm H2O s/L)	−2.1±1.3	−2.1±0.9	−2.0±0.3	0.98
ΔX5 (cm H2O s/L)	−0.0±0.9	−0.2±0.7	−0.5±0.5	0.48
AX (cm H2O/L)	14.6±10.8	15.0±7.8	14.9±4.6	0.99
Fres (Hz)	20.1±4.6	20.9±3.5	21.3±3.1	0.68
ReE-ReI (cm H2O s/L)	0.8±0.6	1.0±0.7	0.4±0.2	0.20
XeE (cm H2O s/L)	−0.4±0.8	−0.4±0.6	−0.5±0.4	0.89
XeI (cm H2O s/L)	−0.6±0.4	−0.8±0.5	−1.1±0.4	0.004
XeE-XeI (cm H2O s/L)	0.2±0.7	0.4±0.6	0.6±0.6	0.24

Data are shown as mean±SD. Statistics were performed using analysis of variance. The bolded values indicate statisical significance.

α (n=47), β (n=20).

The bolded values indicate statisical significance.

*Kruskal-Wallis one-way analysis of variance.

BMI, body mass index.

Quantitative analyses of CT images and correlation with oscillometry

Thirty participants had a non-contrast high-resolution CT of the chest within 6 months of paired oscillometry-conventional PFT. Segmentation and texture analysis of the lung in a patient with IPF using the Lung Texture Analysis tool are shown in online supplemental figure S1. Segmentation failed in four cases after manual review of the output (one case of pneumomediastinum, one case of pacemaker, one case of shoulder arthroplasty and one case of failed central airway segmentation). The demographics, conventional PFT and oscillometry measurements of the remaining 26 patients were not significantly different from the overall cohort (tables 1 and 2). CT texture analysis found that TLCCT was lower than the PFT-measured TLC (3.5±1.1 vs 4.3±1.0 L), with 76.9%±15.9% normal lung volume, 20.8%±16.3% ILD extent and 5.3%±2.4% PVV (online supplemental table S3).

Comparison of spectral and intrabreath oscillometry data with automated CT parameters revealed the highest correlations for XeI with almost all CT-derived parameters (table 4). The highest coefficients of determination for XeI were with PVV and ILD extent (R2=0.69 and 0.62, respectively). For spectral oscillometry, moderate correlations (R2<0.20) were found with X5, X5in, X5ex and AX with the PVV and ILD extent. After adjusting for %FVC, %DLCO, GAP score, XeI remained significantly correlated with PVV and ILD extent (p<0.001 and p=0.002, respectively). The adjusted R2 values indicate that the correlation of XeI is higher for PVV than ILD extent (table 5).

Table 4

Correlation of Oscillometry with CT-derived metrics

	AX	Fres	X5	X5in	X5ex	XeE	XeI	XeE-XeI
TLCCT	0.387	0.216	0.387	0.356	0.389	0.303	0.536	0.095
PVV	0.271	0.116	0.168	0.116	0.178	0.148	0.690	0.314
ILD extent	0.230	0.135	0.141	0.133	0.156	0.119	0.615	0.265
% ground glass	0.298	0.197	0.217	0.208	0.228	0.152	0.538	0.222
% reticular	0.000	0.000	0.000	0.000	0.000	0.000	0.504	0.215
% honeycomb	0.000	0.014	0.000	0.000	0.000	0.000	0.038	0.098

Data are shown as the coefficient of determination, and R2 was derived using polynomial with 2 degrees regression. Similar observations (not shown) were made using linear and exponential regression.

ILD extent=sum of %ground glass, %reticular and %honeycomb.

%, per cent predicted; PVV, pulmonary vessel volume; TLCCT, total lung capacity measured by CT.

Table 5

Results of multivariate analysis

Parameter	Adjusted R2
	ILD extent	PVV
%FVC+%DLCO	0.000	0.106
%FVC+%DLCO+XeI	0.485	0.620
%FVC+%DLCO+XeI+GAP stage	0.448	0.593

%DLCO, % predicted diffusing capacity for carbon monoxide; %FVC, % predicted forced vital capacity; GAP, gender-age- physiology; XeI, reactance at end-inspiration.

Correlations of conventional PFTs with oscillometry

Analysis of intrabreath oscillometry measurements (table 6) revealed that XeI was most highly correlated with TLC and fit better than the spectral oscillometry measurement of X5in (r2=0.58 vs 0.45). Similarly, XeI correlated better than X5in with FVC (R2=0.50 vs 0.44). Of the main physiological parameters in the GAP index, FVC was most significantly correlated with XeI, followed by X5in, X5, AX and Fres, while DLCO was most correlated with X5in and X5.

Table 6

Correlation of oscillometry with conventional PFTs

	AX	Fres	R5	X5	X5in	X5ex	XeE	XeI
%FVC	0.004	0.066	0.000	0.003	0.032	0.000	0.001	0.165
FVC	0.257	0.234	0.076	0.339	0.435	0.215	0.098	0.499
%FEV1	0.085	0.150	0.023	0.094	0.130	0.050	0.032	0.241
FEV1	0.352	0.291	0.132	0.435	0.506	0.305	0.184	0.491
%TLC	0.000	0.012	0.000	0.000	0.011	0.000	0.023	0.186
TLC	0.225	0.126	0.036	0.333	0.450	0.189	0.088	0.579
%DLCO	0.000	0.000	0.065	0.000	0.000	0.000	0.000	0.066
DLCO	0.152	0.074	0.033	0.239	0.360	0.116	0.012	0.285

Data are shown as the coefficient of determination, and R2 was derived using polynomial with 2 degrees regression.

DLCO, diffusing capacity for carbon monoxide; FEV1, forced expiratory volume in one second; FVC, forced vital capacity; PFT, pulmonary function test; TLC, total lung capacity.

Discussion

This is the first study to apply both conventional (spectral) and intrabreath oscillometry for characterisation of patients with IPF who meet the consensus guidelines definition, and to evaluate the correlations of oscillometry with CT-derived metrics of lung volume and lung injury. Intrabreath oscillometry is a novel modality that evaluates dynamic changes in the airways and respiratory tissues during the breathing cycle using a tracking signal with sufficient temporal resolution (0.1 s in the current study).17 25 We observed excellent correlations of XeI with ILD extent and PVV, which remained significant after adjusting for %FVC, %DLCO and GAP score. These novel findings suggest that intrabreath oscillometry identified a physiological-structure signal that was undetectable by conventional PFTs nor standard oscillometry and revealed only by automated CT chest image analysis.6

PVV is a biomarker unique to CALIPER-based analyses that has no visual correlate.6 It is a quantification of vessel-related structures, as determined by the machine learning algorithm when assessing for ground-glass opacities, reticulation and honeycombing. While previously unrecognised, PVV has been shown to be the most important CT metric in predicting disease progression and mortality in patients with IPF.6–9 26 The strong and unique correlation of PVV with XeI in particular indicates that intrabreath oscillometry provides additional information regarding IPF severity that complements conventional PFT and routine visual CT analyses.

Various theories have been proposed to explain the correlation of PVV with IPF severity. Some relate to the intuitive abnormalities of blood-flow diversion, vascular collateralization and pulmonary hypertension6 as well as other potential vascular abnormalities in patients with IPF.27 It has been proposed that increased PVV is secondary to markedly negative intrathoracic pressures generated by the low-compliance, fibrotic lung at end-inspiration causing expansion of the pulmonary vasculature.6 The fact that the strongest correlation of oscillometry with PVV was exhibited by the end-inspiratory reactance (XeI) supports the theory that PVV is a marker of highly negative intrathoracic pressures in patients with IPF at end-inspiration. The marked increase in the retracting forces transmitted to the bronchial wall also explains the typical absence of airway obstruction in IPF; in fact, the ReI values were as low as that in the healthy group and the other Rrs measures were not much higher either (online supplemental table S1).

XeI proved to be the single oscillometry parameter that was the most highly correlated with FVC (r2=0.50) and TLC (r2=0.58); the correlations were higher than the spectral parameter of X5in. XeI was also the only oscillometry parameter that correlated with the GAP Stage. We observed a trend to increasingly large changes in the reactance during a tidal breath as reflected by higher XeE-XeI values with increasing GAP stage (table 3). The intrabreath metrics, XeI and XeE, are exquisitely sensitive to the changes in lung elastance as they are estimated at zero flow, that is, at the endpoints of a tidal breath. The possibility that XeI, XeE-XeI and other intrabreath reactance indices could provide useful information regarding IPF severity will be explored with longitudinal follow-up. As the majority of our patients had mild disease (mostly GAP stage I) at enrolment, we will have the opportunity to evaluate XeI in tracking disease progression and/or drug efficacy if treatment is initiated during follow-up.

Spectral oscillometry showed a distinct pattern of basically normal respiratory resistance with a rightward shift of the reactance-frequency curve (figure 1), leading to increased AX and Fres when compared with the upper limits of normal of 9 cm H2O/L and 12 Hz, respectively.24 The good correlations of X5 and X5in with FVC and DLCO are in keeping with previous observations in two similarly-sized IPF cohorts14 15 although we found higher correlation coefficients between X5 with FVC14 15 and DLCO.14 Our findings also differ from the previous studies which reported normal AX (3.9 and 7.0 cmH2O/L) and Fres values (10.4 and 12.1 Hz).14 15 Our IPF cohort had lower X5 (−2.1 cmH2O s/L), higher AX (12.8 cmH2O/L) and Fres (21.3 Hz). While it is possible that these discrepancies are due to differences in disease severity between the groups (mean %predicted FVC=69.6 vs 95.114 and 83.3,15 it is more likely that they relate to intrinsic difference of the oscillometry devices. Comparisons of measurements made with the currently available commercial oscillometers28 show a large systematic bias with lower AX and Fres and higher (ie, less negative) reactance values with the device used in these studies;14 15 therefore, our present oscillometry results cannot be interpreted numerically in the context of their findings.

We found good correlations of AX and Fres with ILD extent. Our findings corroborate previous observations made by Yamamoto et al,15 the only published study to-date that has addressed the relationship between oscillometry and radiographic metrics. However, direct comparisons cannot be made as the Yamamoto study used a visual scoring system of quantifying structural abnormalities,15 which does not offer the ability to quantify PVV. Furthermore, automated CT image analysis has been shown to be superior to visual scoring of IPF severity.29 As previously reported,6 we also observed high correlation of % ground glass, ILD extent and PVV with conventional PFTs, particularly for %predicted TLC (r2=0.81, 0.72 and 0.60, respectively; online supplemental figure S2). Notably, correlation of PVV was higher for XeI than for %predicted TLC (r2=0.69 vs 0.60).

While all 89 patients completed spirometry and oscillometry, 15 were unable to perform full PFTs and 4 were too symptomatic to complete the 6MW test. This is not unusual in IPF as respiratory symptoms often limit patients’ ability to breath hold or to be off oxygen for the timed needed to perform lung volume measurements. In contrast, oscillometry is conducted during normal tidal breathing, eliminating the need for forced expiratory manoeuvres. A full oscillometry inclusive of spectral and monofrequency recordings can be completed in less than 10 min, offering an advantage over spirometry and conventional PFTs.

Interpretation of the novel intrabreath oscillometry measures is burdened by the lack of normative values, which are available for spectral oscillometry and conventional PFTs. Therefore, our comparison with a modest-size and not age-matched healthy group is only of illustrative value. Our study is also limited to a single centre, with the majority of patients being older men with mild IPF. With increased enrolment and follow-up, the proportion of patients with moderate and severe disease is expected to grow. As such, the strength of the associations between the oscillometry parameters with CT-derived metrics and disease severity is expected to improve. Furthermore, only a subset of patients had imaging in the appropriate time-frame for analysis.

In conclusion, the intrabreath parameter of XeI corresponds to the end-tidal elastance of the respiratory tissues that is not available from conventional PFTs. The high correlation of XeI with CT-derived metrics of PVV and ILD extent and the finding that these correlations remain significant after accounting for %FVC, %DLCO and the GAP score suggest that XeI provides complementary information regarding IPF severity, particularly in settings where access to subspecialty thoracic radiology expertise and/or automated quantitative CT analysis is limited. While quality control is important for accurate measurements, training of technicians to conduct oscillometry is relatively simple.18 Given the ease of conduct of oscillometry during normal tidal breathing, high patient satisfaction18 30 and the small environmental footprint of oscillometry devices, our findings support the addition of oscillometry and in particular, its intrabreath technique as an adjunct to routine PFTs in the management of patients with IPF.