Several specific high-resolution computed tomography patterns correlate with survival in patients with idiopathic pulmonary fibrosis.

BACKGROUND: Evidence of honeycombing in high-resolution computed tomography (HRCT) is a recognized risk factor for shortened survival in patients with idiopathic pulmonary fibrosis (IPF), but few studies have evaluated the feasibility of exploiting other specific patterns for predicting survival. The aim of this study was to examine the extent of specific HRCT patterns in IPF and determine whether they correlate with clinical features, pulmonary function tests (PFT), and survival.
METHODS: Both the presence and extent of specific HRCT patterns, such as traction bronchiectasis, honeycombing, architectural distortion, reticulation, emphysema, and ground glass opacity, in 129 HRCT examinations were scored semi-quantitatively in three zones of each lung. HRCT examinations were also re-classified according to the 2011 and 2018 international statements. Correlations were calculated between the scores of specific HRCT patterns, clinical features, PFT, and patient survival.
RESULTS: The extent of traction bronchiectasis was found to be an independent risk factor of shortened survival (HR 1.227, P=0.001). Patients with a possible usual interstitial pneumonia (UIP) pattern had a better median survival than the patients with a definite UIP pattern (61 vs. 37 months, P=0.026). The extents of traction bronchiectasis, honeycombing, and architectural distortion displayed an inverse correlation with all PFT values at the time of diagnosis. There were few differences between the radiological classifications of the 2011 and 2018 international statements.
CONCLUSIONS: We conclude that several specific HRCT patterns displayed a correlation with shortened survival in IPF; these may help in evaluating the risk of death in IPF patients. 2021 Journal of Thoracic Disease. All rights reserved.

Introduction

Idiopathic pulmonary fibrosis (IPF) is a heterogeneous disease that has a variable course, which makes predicting survival challenging (1). High-resolution computed tomography (HRCT) is the most important noninvasive diagnostic method in the diagnosis of IPF, but the feasibility of exploiting HRCT for evaluating the course of the disease remains unclear. In addition, several previous studies investigating specific HRCT patterns in IPF, diagnosed according to the 2011 international statement, have mainly concentrated on diagnosis (2-4).

Honeycombing is considered as one of the core findings when the UIP pattern is diagnosed in HRCT (5). Many of the studies comparing disease progression with HRCT findings have concentrated on the honeycombing pattern, revealing its association with shortened survival in IPF (6-8). Patients with a honeycombing score higher than 15% in HRCT have been shown to have a shorter survival than those with a honeycombing score less than 15% (6). In a recent study of 350 IPF patients, honeycombing was associated with poorer survival (7) and it has also been reported that honeycombing was associated with an increased risk of death (8).

The previous studies investigating either the presence or extent of traction bronchiectasis for predicting survival in IPF patients who have been diagnosed according to 2011 international IPF statement, have shown discrepant results (7-9). In the study of Arcadu et al., the presence of traction bronchiectasis was not associated with shortened survival, unlike the global assessment of interstitial abnormalities (7). Salisbury et al. found that the fibrosis score (number of lobes where traction bronchiectasis was present) was associated with an increased risk of death (8). Recently, Sverzellati et al. found that the extent of traction bronchiectasis was associated with shortened survival (9).

Since relatively little is known about the correlation between specific HRCT patterns (except honeycombing) and survival in IPF patients diagnosed according to the 2011 international statement, we aimed to evaluate both the presence and extent of specific HRCT patterns and to examine their correlation with clinical features, such as age, gender, smoking history, dyspnea, cough, long term oxygen therapy, and the number of hospitalizations due to respiratory reasons, as well as data from pulmonary function tests (PFT), and finally with patient survival. Furthermore, we have classified the HRCT examinations according to both the 2011 and 2018 international statements in order to determine whether the application of different guidelines would change the radiological categorization (1,5). We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/jtd-20-1957).

Methods

Patient selection and data collection

The subjects of this study were obtained from Kuopio University Hospital pulmonology clinic medical records (digitized from 2004 onwards) from January 2002 to December 2012 using the International Classification Codes (ICD-10) J84.1, J84.8 and J84.9. From the 610 patients discovered, 132 IPF-patients were identified from their clinical information, radiology, and histology when available, and accepted after reaching a multidisciplinary consensus. Patients with pulmonary fibrosis caused by a known reason, for example, asbestosis or connective tissue disease, were excluded. Smoking history, clinical information, and PFT were retrospectively gathered. Death or lung transplantation was considered as the outcome. The collection of this data has been described in more detail in three previous publications (10-12).

Because of the retrospective nature of this study and the fact that most of the study subjects were already deceased, no informed consents for study inclusion were collected. The study protocol was approved by the Research Ethics Committee of the Northern Savo Hospital District (statement 17/2013) and by the National Institute for Health and Welfare (Dnro THL/1052/5.05.01/2013). This study was conducted in compliance with the Declaration of Helsinki (as revised in 2013).

The CT protocol

Digital HRCT scans and analogue film images, at the time of diagnosis, were available from 99 of 132 patients (75%) and 28 patients (21%), respectively. Other comparable CT examinations, such as CT pulmonary angiography, were available from five patients. All HRCT scans were spaced axial scans. The radiologists evaluated the material other than digital HRCT scans eligible for further analysis. A total of 86 out of 132 images (65%) were evaluated as good quality, 43 as suboptimal quality, and three as poor quality/not diagnostic. The poor-quality images were excluded, and thus the final number of study subjects was 129, of which 97 (75%) had undergone digital HRCT scanning.

Re-classifying of CT

Two radiologists independently re-classified the HRCT scans according to the 2011 international statement into three classes: usual interstitial pneumonia (UIP), possible UIP, and inconsistent with UIP (1). The first radiologist (10 years of expertise in thoracic radiology) re-classified the 129 HRCT examinations in 2014. In 2018, the second radiologist (2 years of expertise in thoracic radiology) re-classified 121 of the original HRCT examinations available according to both the 2011 and 2018 international statements i.e., as UIP, probable UIP, indeterminate with UIP, and alternative diagnosis (5). There were no differences in the categorical classification of the HRCT scans between good and suboptimal quality HRCTs.

Scoring of specific HRCT patterns

The first radiologist assessed the 129 HRCT examinations in detail using a form designed for a previous study (13). The presence and extent of several specific HRCT patterns such as honeycombing, traction bronchiectasis, emphysema, architectural distortion, reticulation, ground-glass opacity (GGO), pleural plaques, consolidation, bronchiectasis, and crazy paving were assessed separately in three zones of each lung. The extents of GGO, reticulation, and honeycombing were semi-quantitatively graded on a scale of 0 to 4 (0= absent, 1= minor peripheral scattered changes, 2= uniform peripheral or minor central scattered changes, 3= substantial peripheral changes that penetrated deeply into the lung parenchyma, 4= very abundant peripheral and central changes). Instead, emphysema, traction bronchiectasis, and architectural distortion were scored using a scale of 0 to 3 (0= absent, 1= single scattered changes, 2= larger single or several minor changes, 3= uniform or substantial changes) (). Scores from all six zones were summed to obtain the total score for each of the patterns, the maximum scores being 24 or 18, respectively. The scoring system has been described in more detail in a previous publication (13). To assess the repeatability of the scoring, the second radiologist scored 26 of the original HCRT examinations, of which 96% were good and only 4% were of suboptimal quality. There were mainly no differences in specific HRCT patterns observed between good and suboptimal quality HRCT examinations, except for traction bronchiectasis. There were more HRCT examinations where traction bronchiectasis was observed among the good quality HRCTs than among the suboptimal quality HRCTs (91% vs. 67%, respectively, P=0.001).

Figure 1

Representative HRCT images from four UIP patients with emphysema and fibrosis, categorized according to the 2018 international statement. (A) HRCT shows uniform peripheral scattered changes (score 2) of honeycombing in the right upper zone and substantial signs of peripheral honeycombing that penetrated deeply into the lung parenchyma (score 3) in the left upper zone; (B) HRCT shows minor peripheral scattered changes (score 1) of reticulation in the middle zones and uniform peripheral scattered changes (score 2) of honeycombing on the right side, as well as minor peripheral scattered changes (score 1) of honeycombing on the left side; (C) HRCT shows uniform peripheral scattered changes (score 2) of reticulation in the lower zones and larger single scattered changes (score 2) of traction bronchiectasis (arrows), as well as minor peripheral scattered changes (score 1) of honeycombing only on the right side; (D) HRCT reveals larger single scattered changes (score 2) of paraseptal and bullous emphysema in the right lower zone, as well as larger single scattered changes (score 2) of traction bronchiectasis (arrows) and honeycombing. HRCT, high-resolution computed tomography; UIP, usual interstitial pneumonia.

Pulmonary function tests

The PFT values at the time of diagnosis are presented as percent predicted. The change in forced vital capacity (FVC) and diffusion capacity to carbon monoxide (DLCO) during six and 12 months’ follow-up have been calculated as follows: (FVC or DLCO at 6 or 12 months—FVC or DLCO at the time of diagnosis)/FVC or DLCO at the time of diagnosis ×100, with the result expressed in percentages. Furthermore, the FVC change during 6 and 12 months was classified into two categories: <10% decline or improvement of FVC and ≥10% decline in FVC. Similarly, DLCO was classified into two categories using a 15% decline as a cut-off point (14).

Statistical analysis

Data were expressed as means and standard deviations with ranges or frequencies with percentages. The presence and extent of specific HRCT patterns, clinical features, PFT and survival were analyzed in the different HRCT categories according to the 2011 international statement (1). Independent-samples Kruskal-Wallis or Mann-Whitney U tests were used for continuous variables, with Chi-square test also being used when appropriate. Correlations between variables were investigated with the Spearman rank correlation coefficient. The Kaplan-Meier estimator and Cox regression analysis were used to evaluate the association between variables and survival. Clinical features such as age, gender, PFT, FVC and DLCO absolute change for 12 months, oxygen therapy, dyspnea, cough, and smoking status, as well as specific HRCT patterns were analyzed in a univariate Cox regression analysis. Some of the data of the clinical features has been published previously (12). Thereafter, a stepwise Cox regression model was built using variables with univariate P values <0.20, entering variables at the 0.05 significance level and removing them at the 0.10 level. Inter-observer agreements of the HRCT examinations according to the 2011 international statement and the specific HRCT patterns are presented as a kappa (κ) value: good agreement κ=0.61–0.80, moderate agreement κ=0.41–0.60 and fair agreement κ=0.21–0.40. P values <0.05 were considered statistically significant. Missing data was excluded. IBM statistics SPSS software, version 22.0, was used in the statistical analysis of the whole data.

Results

HRCT classification according to the 2011 and 2018 international IPF statements

Sixty-five HRCT examinations were classified as UIP (50%), 42 as possible UIP (32%), and 22 as inconsistent with UIP (17%) according to the 2011 guidelines. There were only subtle differences between the radiological classifications of the 2011 and 2018 international statements, since the UIP category did not change, whereas the HRCT examinations in the possible UIP category were divided mostly into the probable and the indeterminate with UIP categories ().

Figure 2

Radiological categorization of HRCT examinations in patients with idiopathic pulmonary fibrosis (IPF) (A) according to the 2011 IPF statement (N=129 and N=121, respectively) and (B) according to the 2018 IPF statement (N=121). Numbers are presented as N. HRCT, high-resolution computed tomography; UIP, usual interstitial pneumonia.

Inter-observer agreement of the radiologists

The inter-observer agreement of the radiologists was moderate (κ=0.429) in the re-classification of the patterns into three categories, i.e., UIP, possible UIP, and inconsistent with UIP. When the UIP pattern was compared with the non-UIP pattern (the possible and inconsistent with UIP patterns combined) inter-observer agreement was higher (κ=0.551). Inter-observer agreement was good when evaluating honeycombing (κ=0.783) and emphysema (κ=0.623).

Patient characteristics and specific HRCT patterns

The majority of the patients were men (73%) with a mean age of 71 years. The characteristics of the patients and the different radiological categories are presented in . Some of the data have been published previously (10-12). Reticulation was the most predominant HRCT pattern found in 99% of cases, followed by traction bronchiectasis (83%). The extent of traction bronchiectasis was higher in the UIP group than in the other two groups. present the specific HRCT patterns in the whole cohort and in the different categories.

Table 1

Characteristics of the patients with idiopathic pulmonary fibrosis categorized according to the 2011 IPF statement

Feature	All (N=129)	UIP (N=65)	Possible UIP (N=42)	Inconsistent with UIP (N=22)
Age, mean (SD)	70.8 (9.5)	70.7 (10.1)	72.2 (8.9)	68.1 (8.4)
FVC%, mean (SD)	76.7 (18.6)	75.7 (19.5)	79.7 (18.4)	74.2 (16.6)
FEV1%, mean (SD)	77.1 (17.1)	75.9 (18.1)	78.9 (16.3)	76.8 (16.1)
DLCO%, mean (SD)	55.9 (17.2)	53.0 (14.9)†	63.0 (19.3)§	50.7 (15.9)
Male	94 [73]	53 [82]‡	31 [74]§	10 [46]
Female	35 [27]	12 [18]‡	11 [26]§	12 [54]
Ever-smoker	82/126 [65]	46/63 [73]‡	25/41 [61]	11 [50]
Dyspnea	96 [74]	48 [74]	29 [69]	19 [86]
Cough	77 [60]	34 [52]†	30 [71]	13 [59]
Oxygen therapy	45 [35]	19 [29]‡	13 [31]§	13 [59]
ReHosp	87 [67]	38 [59]	32 [76]	17 [77]

Numbers are presented as mean (SD) or N [%]. P value <0.05 †UIP vs. possible UIP, ‡UIP vs. inconsistent with UIP, §possible UIP vs. inconsistent with UIP. UIP, usual interstitial pneumonia; FVC%, forced vital capacity percent predicted; FEV1%, forced expiratory volume in one second percent predicted; DLCO%, diffusion capacity to carbon monoxide percent predicted; ReHosP, hospitalization due to respiratory reasons; N, number; SD, standard deviation.

Table 2

Specific HRCT patterns of patients with idiopathic pulmonary fibrosis categorized according to the 2011 IPF statement

Radiological pattern	All (N=129)	UIP (N=65)	Possible UIP (N=42)	Inconsistent with UIP (N=22)
Reticulation	128 [99]	64 [99]	42 [100]	22 [100]
Traction bronchiectasis	106/128 [83]	59 [91]†	30/41 [73]	17 [77]
Honeycombing	91/122 [75]	64 [99]†‡	18/37 [47]	9/21 [43]
Architectural distortion	63 [49]	44 [68]†‡	12 [29]	7 [32]
GGO	61/117 [52]	16/59 [27]†‡	29/38 [76]	16/20 [80]
Emphysema	34 [26]	24 [37]†	6 [14]	4 [18]
Plaques	36/128 [28]	–	–	–
Consolidation	3 [2]	–	–	–
Bronchiectasis	4 [3]	–	–	–
Crazy paving	1 [1]	–	–	–

Numbers are presented as N [%]. P value <0.05 †UIP vs. possible UIP, ‡UIP vs. inconsistent with UIP. HRCT, high-resolution computed tomography; UIP, usual interstitial pneumonia; GGO, ground glass opacity; N, number.

Table 3

The extent of specific HRCT patterns in patients with idiopathic pulmonary fibrosis categorized according to the 2011 IPF statement

Radiological pattern§	All (N=129)	UIP (N=65)	Possible UIP (N=42)	Inconsistent with UIP (N=22)
Reticulation	8.9±3.1 (0–19)	8.39±3.30 (0–19)‡	9.14±2.95 (3–16)	10.19±2.44 (5–14)
Traction bronchiectasis	3.5±2.9 (0–13)	4.49±2.94 (0–13)†‡	2.22±2.29 (0–11)	2.82±2.58 (0–10)
Honeycombing	3.6±3.5 (0–14)	5.77±3.22 (1–14)†‡	1.11±1.61 (0–8)	1.52±2.58 (0–10)
Architectural distortion	2.0±3.0 (0–14)	3.06±3.44 (0–14)†‡	0.83±2.07 (0–12)	0.86±1.36 (0–4)
GGO	3.3±4.4 (0–18)	1.17±2.39 (0–10)†‡	4.39±3.67 (0–13)	7.58±6.25 (0–18)
Emphysema	1.4±3.2 (0–14)	1.97±3.47 (0–14)†	0.57±1.95 (0–11)	1.55±3.90 (0–14)

§GGO, reticulation and honeycombing score max =24, traction bronchiectasis, architectural distortion and emphysema score max =18, score numbers are presented as mean ± SD (min–max). P value <0.05 †UIP vs. possible UIP, ‡UIP vs. inconsistent with UIP. HRCT, high-resolution computed tomography; UIP, usual interstitial pneumonia; GGO, ground glass opacity; N, number.

Correlation of the HRCT patterns with survival

The mean follow-up time was 50 months (range, 0–188 months). One hundred and twenty-three patients (95%) died during the follow-up and two patients underwent lung transplantation. The extents of traction bronchiectasis, honeycombing, architectural distortion, and emphysema were associated with shortened survival in a univariate Cox regression analysis (). A univariate Cox regression analysis of the clinical features is presented in the Supplementary file (Table S1). The extent of traction bronchiectasis and the presence of emphysema independently associated with shortened survival in the multivariate Cox regression analysis (). The extent of emphysema did not associate with shortened survival in a multivariate analysis (data not shown). The possible UIP group had a longer median survival than the UIP patients in the unadjusted Kaplan-Meyer estimator analysis. The survival was also longer in patients without honeycombing, traction bronchiectasis, architectural distortion or emphysema ( and Table S2).

Table 4

Univariate Cox regression analysis of the presence and extent of specific HRCT patterns in patients with idiopathic pulmonary fibrosis

Radiological pattern	Presence				Extent
	HR	95% CI	P		HR	95% CI	P
Honeycombing	2.017	1.297–3.137	0.002		1.122	1.065–1.182	<0.001
Traction bronchiectasis	1.757	1.078–2.863	0.024		1.212	1.130–1.299	<0.001
Architectural distortion	2.220	1.521–3.240	<0.001		1.160	1.094–1.231	<0.001
Emphysema	1.765	1.169–2.666	0.007		1.055	1.000–1.114	0.051
GGO	0.800	0.549–1.164	0.243		0.981	0.937–1.027	0.412
Reticulation	0.753	0.104–5.424	0.778		1.023	0.954–1.098	0.521
UIP—possible UIP, UIP	Reference	–	–		–	–	–
Possible UIP	0.644	0.428–0.970	0.035		–	–	–

P value <0.05 regarded as significant. HRCT, high-resolution computed tomography; CI, confidence interval; HR, hazard ratio; GGO, ground glass opacity; UIP, usual interstitial pneumonia.

Table 5

Multivariate Cox regression analysis of the extent of the specific HRCT patterns and clinical features in patients with idiopathic pulmonary fibrosis

Variable	HR	95% CI	P
Age	1.058	1.030–1.087	<0.001
Gender, female	Reference	–	–
Male	1.946	1.087–3.483	0.025
DLCO%	0.956	0.938–0.973	<0.001
DLCO absolute change in 12 months	0.970	0.958–0.983	<0.001
The presence of emphysema	1.971	1.134–3.428	0.016
The extent of traction bronchiectasis	1.227	1.088–1.385	0.001
The extent of architectural distortion	0.984	0.871–1.112	0.799

P value <0.05 regarded as significant. HRCT, high-resolution computed tomography; HR, hazard ratio; CI, confidence interval; DLCO%, diffusion capacity to carbon monoxide percent predicted.

Figure 3

Unadjusted Kaplan-Meier survival curves in patients with idiopathic pulmonary fibrosis and the presence of specific HRCT patterns stratified by (A) traction bronchiectasis, (B) honeycombing, (C) emphysema, (D) architectural distortion, and (E) UIP vs. possible UIP pattern. Survival curves and P values represent the time to death or lung transplantation. HRCT, high-resolution computed tomography; UIP, usual interstitial pneumonia.

Correlation of the HRCT patterns with clinical features and PFT

The extents of traction bronchiectasis, honeycombing, and architectural distortion had an inverse correlation with all PFT values at the time of diagnosis. The UIP pattern correlated with the extent of traction bronchiectasis, honeycombing, architectural distortion, and emphysema (r=0.401, P<0.001; r=0.742, P<0.001; r=0.491, P<0.001; r=0.233, P<0.01, respectively) showing an inverse correlation with GGO and reticulation (r=−0.522, P<0.001; r=−0.202, P<0.05, respectively). In addition, the extent of reticulation correlated with dyspnea (r=0.195, P<0.05) and the extent of emphysema inversely correlated with cough (r=−0.306, P<0.01) (Table S3). Cough was more common in the possible UIP patients than in the UIP group. The number of hospitalizations due to respiratory reasons between the different radiological categories did not differ in the various groups (). However, more hospital visits due to respiratory reasons were observed among the non-UIP patients than in the UIP patients [49/64 (78%) vs. 38/65 (59%), respectively, P=0.03]. In addition, the hospital visits due to respiratory reasons correlated inversely with the extent of emphysema (Table S3). There were more ever-smokers among men than women (78% vs. 26%, P<0.001) as well as among the patients with emphysema than those without emphysema (88% vs. 55%, P=0.001). The men had more extensive honeycombing, architectural distortion, and emphysema than women, as did the ever-smokers in comparison to the non-smokers (Table S4).

Discussion

In this real-life study, we demonstrated that the extent of traction bronchiectasis correlated with shortened survival in IPF patients, even when adjusted for clinical features and another indicator of fibrosis. The presence of emphysema seemed to have a similar correlation as a categorical variable, whereas the emphysema score did not associate with survival. Furthermore, we found that the possible UIP patients had a longer survival in comparison to the UIP patients and that honeycombing and architectural distortion were associated with shortened survival. In addition, we made three novel findings. Firstly, the extent of several specific HRCT patterns, such as traction bronchiectasis, honeycombing, architectural distortion, and emphysema positively correlated with the UIP pattern. Secondly, the extent of reticulation correlated with dyspnea and the extent of emphysema inversely correlated with cough. Thirdly, there were a similar number of hospitalizations due to respiratory reasons in the different radiological categories. The radiological categorizations between the 2011 and 2018 international guidelines varied only marginally.

Although several previous studies have investigated specific HRCT-patterns, many of them have included fibrotic idiopathic interstitial pneumonia (IIP) patients rather than focusing on IPF patients only (15-17) or IPF patients were diagnosed according to an outdated international statement (18). Several studies focusing on IPF diagnosed according to the 2011 international statement, have found similar results regarding the correlation of survival with traction bronchiectasis, whereas the results of emphysema have been partly different (7-9,19-21). The two previous studies that have evaluated the correlation of specific HRCT patterns with survival in IPF analyzed fewer variables and moreover, the quantification of HRCT patterns was more modest than in our study (7,8). In a study of 58 IPF patients, Sverzellati et al. discovered that the visually evaluated extent of traction bronchiectasis was associated with shortened survival in a univariate Cox regression analysis at baseline and in two years’ follow up (9). Salisbury et al. found that the number of lung lobes with traction bronchiectasis was associated with an increased risk of death (8). Arcadu et al. did not detect any correlation between the presence of traction bronchiectasis and patient survival in a study of 350 IPF patients. However, in Arcadu’s study, the extent of traction bronchiectasis was not analyzed, which could explain the contradictory result (7). Although Jacob et al. evaluated several HRCT patterns using automated CT analysis, traction bronchiectasis was evaluated visually. It was claimed that visually evaluated traction bronchiectasis was associated with shortened survival and a FVC decline during 12 months (19,22). In more recent studies using quantitative CT-analysis, CT histogram values and vessel-related structures were associated with survival in IPF (23-25). Since traction bronchiectasis has not yet been studied with automated CT analysis methods, we believe that this present work is so far the most detailed study of traction bronchiectasis in IPF.

Several previous studies evaluating the correlation between pulmonary fibrosis, emphysema, and survival have included also patients with other interstitial lung diseases (ILDs) in addition to IPF, or they have not diagnosed IPF patients according to the current international IPF statement (26-28). In a study focusing on combined pulmonary fibrosis and emphysema (CPFE), there was no difference in survival between the patients with ≥10% emphysema and less fibrosis and those with more fibrosis and <10% emphysema (20). Jacob et al. demonstrated that emphysema was associated with increased mortality in IPF patients, but independently the extent of emphysema did not predict mortality (21). Similarly, we did not find that the emphysema score would associate independently with survival even though the presence of emphysema seemed to correlate with shortened survival.

The pattern of architectural distortion described previously as “abnormal displacement of bronchi, vessels, fissures, or septa caused by diffuse or localized lung disease, particularly interstitial fibrosis” has not been studied previously in cases diagnosed according to the 2011 international IPF statement (29). We demonstrated that architectural distortion was associated with shortened survival and furthermore, it correlated with PFT values and several clinical features. Two previous studies with different study protocols have reported results that are in line with our findings. Arcadu et al. observed that the extent of “global assessment of interstitial abnormalities” was associated with increased mortality and correlated inversely with FVC and DLCO percent predicted (7). In a study of fibrotic idiopathic interstitial pneumonia (IIP) patients, the definite UIP patients with the overall extent of parenchymal abnormalities >45% had shortened survival in comparison to those patients with an overall CT score <45% (17).

Our result of the survival difference between the possible and definite UIP patterns is in line with previous studies (7,8). Salisbury et al. showed that the possible UIP pattern was associated with longer survival in comparison to the definite UIP pattern and that the difference remained after adjustments (8). Arcadu et al. reported that the definite UIP pattern was associated with shortened survival in comparison to the possible and the inconsistent with UIP patterns (7). Controversially, a study of 42 IPF patients did not detect a survival difference between the definite UIP and the possible UIP patients, a finding which may be partly due to the small size of the cohort (6).

As in the previous studies, honeycombing was associated with shortened survival also in our study (6-8). Novelli et al. found that IPF patients with a honeycombing score ≥15% had worse survival than those who had a honeycombing score <15% (6). In two more recent studies, honeycombing was associated with an increased risk of death in univariate Cox regression analyses (7,8). In a quantitative CT-analysis study, the area of honeycombing was also associated with worse survival in IPF (30). Controversially, in a Japanese study, 43 IPF patients did not have honeycombing in HRCT at baseline and the emergence of honeycombing in the three years’ follow-up did not affect patient survival. The small cohort size and relatively young mean age (64.5 years) of subjects examined in that study could partly explain the contradictory result (31).

The limitations of our study are those typically and often encountered in real-life studies. The study was retrospective, and the cohort size was relatively small. In addition, the CT protocol was variable, and the study was performed in one institution. The inter-observer agreement between the two radiologists for the UIP classifications in our study was moderate, which is similar to previous studies (32-34). However, the amount of suboptimal quality HRCT examinations (31%) may have affected the radiologists’ ability to interpret fine differences when evaluating the radiological classifications. In addition, there was a rather wide difference between the radiologists’ experience in thoracic radiology (10 and 2 years of expertise, respectively). Nevertheless, in our study, the inter-observer agreement for specific HRCT patterns was good, which is also in line with previous studies (6,9,15,16). Furthermore, not all of the material was available for the second HRCT analysis in 2018.

Conclusions

In conclusion, several specific HRCT patterns were associated with shortened survival in IPF. The extent of traction bronchiectasis was one of the strongest patterns to be individually associated with shortened survival, which suggests that this pattern could be an important prognostic factor in IPF. These results may help in evaluating the risk of death in IPF patients.

The article’s supplementary files as