- Original research article
- Open Access
Assessment of survival in patients with idiopathic pulmonary fibrosis using quantitative HRCT indexes
Multidisciplinary Respiratory Medicinevolume 13, Article number: 43 (2018)
The assessment of Idiopathic Pulmonary Fibrosis (IPF) using HRCT requires great experience and is limited by a significant inter-observer variability, even between trained radiologists. The evaluation of HRCT through automated quantitative analysis may hopefully solve this problem. The accuracy of CT-histogram derived indexes in the assessment of survival in IPF patients has been poorly studied.
Forty-two patients with a diagnosis of IPF and a follow up time of 3 years were retrospectively collected; HRCT and Pulmonary Function Tests (PFTs) performed at diagnosis time were analysed; the extent of fibrotic disease was quantified on HRCT using kurtosis, skewness, Mean Lung Density (MLD), High attenuation areas (HAA%) and Fibrotic Areas (FA%). Univariate Cox regression was performed to assess hazard ratios for the explored variables and a multivariate model considering skewness, FVC, DLCO and age was created to test their prognostic value in assessing survival. Through ROC analysis, threshold values demonstrating the best sensitivity and specificity in predicting mortality were identified. They were used as cut-off points to graph Kaplan-Meier curves specific for the CT-indexes.
Kurtosis, skewness, MLD, HAA% and FA% were good predictors of mortality (HR 0.44, 0.74, 1.01, 1.12, 1.06; p = 0.03, p = 0.01, p = 0.02, p = 0.02 and p = 0.017 respectively). Skewness demonstrated the lowest Akaike’s information criterion value (55.52), proving to be the best CT variable for prediction of mortality. Significant survival differences considering proposed cut-off points were also demonstrated according to kurtosis (p = 0.02), skewness (p = 0.005), MLD (p = 0.003), HAA% (p = 0.009) and FA% (p = 0.02) – obtained from quantitative HRCT analysis at diagnosis time.
CT-histogram derived indexes may provide an accurate estimation of survival in IPF patients. They demonstrate a correlation with PFTs, highlighting their possible use in clinical practice.
Idiopathic pulmonary fibrosis (IPF) is a rare pulmonary disease characterized by a progressive fibrosis of the lungs that leads to the inexorable worsening of lung function . IPF is also characterized by a very poor survival rate, ranging from 3 to 5 years from diagnosis, which is worse than in many cancers [2,3,4,5].
ILD assessment using HRCT requires great experience and is limited by a significant inter-observer variability, even between trained radiologists . In the longitudinal evaluation of disease progression, some studies are based on a “visual” semi-quantitative analysis of HRCT; however, this method has been considered “tedious and not reproducible” .
A more objective evaluation of HRCT through automated quantitative analysis may hopefully solve this problem. Recent studies started exploring this second possibility using, in most cases, complex algorithms, able to perform texture analysis. However, although very promising, texture analysis requires high computational steps, currently reserving its use for a research setting. To date, only a few studies have focused on CT-histogram evaluation, which is a different quantitative analysis method requiring fewer computational steps than texture analysis. These characteristics of objectivity and relative simplicity make this analysis potentially usable and reproducible in a realistic clinical setting.
Therefore, the aim of this study is to evaluate in a real-life work setting, the prognostic value of histogram-based quantitative HRCT variables in predicting mortality in a cohort of patients affected by IPF.
A retrospective analysis of our interstitial lung disease (ILD) database identified a total of 210 patients who had received a multidisciplinary team (MDT) diagnosis of IPF according to 2011 ATS/ERS/JRS/ALAT IPF guidelines over a period of three years (January 2012 to August 2015) . Only patients with an unenhanced, supine, volumetric thin-section CT exam performed in our departmental CT-scanner a well conducted clinical and functional follow up were included in the analysis. A median follow up time of 3 years was also required to assess survival. A total of 42 consecutive IPF patients fulfilling these criteria were finally included in the analysis (Fig. 1).
This study received approval the ethical committee of Policlinico-Vittorio Emanuele Hospital of Catania.
All clinical information was obtained from medical records. PFTs (Forced Vital Capacity (FVC) and Diffusing Capacity for carbon monoxide (DLCO)) were performed according to the ATS/ERS guidelines using Vmax Sensor Medics . Results were expressed as percentages of predicted values. PFTs were matched with HRCT examinations (mean distance between PFTs and HRCTs was 9.59 days ±44). All patients were treated, as recommended, with antifibrotic therapies (pirfenidone or nintedanib). No difference in terms of treatment dose was observed between deceased and non (p = 0.19).
HRCT examinations were acquired according to the following technical parameters: thickness ranged between 0.625–1.25 mm; sharp kernel imaging reconstruction, contiguous or overlap images; no contrast media administration. A user-friendly DICOM (digital imaging and communications in medicine) - based image processing software for automatic lung parenchyma segmentation was used to calculate HRCT quantitative parameters. This tool has been developed in the MATLAB R2016a simulation environment (The MathWorks, Natick, MA, USA), running on iMac (3,5 GHz Intel Core i7 processor, 16 GB memory random-access memory; Apple Computer, Cupertino, CA, USA) with Mac Operating System OS X El Capitan.
HRCT attenuation histogram was extracted after automatic lung parenchyma segmentation using region growing algorithm to include voxels between − 200 Hounsfield units (HU) and − 1.024 HU. kurtosis, Mean Lung Density (MLD), skewness, High Attenuation Areas (HAA) % and Fibrotic Areas (FA) % were subsequently calculated (Fig. 2).
It is known that kurtosis, skewness and MLD may provide prognostic estimation of disease severity in subjects with IPF and other ILDs. In particular, low kurtosis, low skewness and high MLD values may indicate severe IPF .
HAA may reflect the presence of parenchymal lesions, such as ground-glass opacity and reticulation. HAA% was calculated as percentages of the extracted whole lung volume with attenuation values ranging from − 250 HU to − 600 HU. In the same way, FA% was calculated considering values raging from − 1.024 HU to − 700 HU .
Characteristics of the study population were expressed as mean (SD), median (interquartile range) or as percentage of the relative frequency as appropriated. Wilcoxon test for unpaired data was used to assess differences between baseline and follow up time. Reported data were adjusted using Holm’s correction. Spearman’s rank correlation was used to assess relation between CT-histogram based variables. Univariate and multivariate Cox proportional hazard regression analysis were performed to assess Hazard Ratios for predictors of survival. Akaike’s information criterion (AIC) was used to estimate the relative quality of the statistical models. Kaplan-Meier survival analysis was used to assess overall survival while ROC analysis was used to assess sensitivity and specificity. All the statistical analyses were performed using STATA/IC 14.2 version. A p less than 0.05 was considered significant.
The study evaluated 42 consecutive subjects with an MDT diagnosis of IPF. Thirty-five were male (83.33%) and 7 were female (16.67%). There were 33 (78.57%) former/current smokers and 9 (21.43%) never-smokers. The mean age at diagnosis was 68.26 ± 6.34 years. The mean follow up time was 1037 ± 430 days. Thirty-height (90.48%) patients were treated only with pirfenidone while 4 (9.52%) switched from pirfenidone to nintedanib because of drug-related adverse events. A total of 10 patients (23.33%) died during the follow up period. All demographic, PFTs and CT-derived variables at baseline are summarized in Table 1.
A subset of 23 patients underwent a follow up HRCT within 12 ± 3 months from baseline. No significant differences were found in the variation of kurtosis, skewness, MLD, HAA%, FA%, FVC and GAP index (p = 0.25, p = 0.25, p = 0.40, p = 0.40, p = 0.68, p = 0.88, p = 0.21 respectively) (see Table 2).
Univariate Cox proportional hazard regression analysis demonstrated, as expected, that both FVC% and DLCO% at baseline were significant predictors of survival (HR 0.93 and 0.93; p = 0.019 and p = 0.005). Similarly, also kurtosis, skewness, MLD, HAA% and FA% were good predictors of mortality (HR 0.44, 0.74, 1.01, 1.12, 1.06; p = 0.03, p = 0.01, p = 0.02, p = 0.02 and p = 0.017 respectively). Results of univariate Cox hazard regression are summarized in Table 3.
Akaike’s information criterion (AIC) was calculated to estimate relative quality of the statistical model. Skewness demonstrated the lowest AIC value (AIC = 55.52) (see Table 4).
A multivariate model was created considering FVC, DLCO, age and skewness, which was considered the best CT-derived variable because of the lowest AIC value. This model demonstrated an HR of 0.98, 0.94, 0.98 and 0.28 respectively for the four considered variables (AIC = 56.86) (see Table 5). ROC analysis shown in Fig. 3, demonstrated a sensitivity of 70% and a specificity of 71.9% for a threshold value of ≤0.76 of kurtosis (AUC = 0.70; p = 0.05), a sensitivity of 80% and a specificity of 75% for a threshold value of ≤1.15 of skewness (AUC = 0,731; p = 0.021), a sensitivity of 70% and a specificity of 78% for a threshold value of > − 762.97 of MLD (AUC = 0.722; p = 0.032), a sensitivity of 80% and a specificity of 62.5% for a threshold value of > 14.76 of HAA% (AUC = 0.681; p = 0.104) and a sensitivity of 80% and a specificity of 62.5% for a threshold value of > 24.26 of FA% (AUC = 0.706; p = 0.045) in predicting mortality.
Each threshold value demonstrating the best sensitivity and specificity in predicting mortality was used as cut-off point to graph Kaplan-Meier curves specific for the CT indexes as shown in Fig. 4.
In recent years, HRCT has become a cornerstone in the diagnosis of IPF. Furthermore, the less-invasiveness of the exam compared to lung biopsy and the possibility of evaluating lung damage using quantitative analysis, has promoted HRCT as a reasonable marker of disease extension [10, 12,13,14,15,16,17,18]. To our knowledge, only three studies, evaluating the predictive value of CT histogram indexes in assessing survival, have been published [11, 19, 20].
Analysing these studies, we found some differences in terms of CT quantitative analysis protocols, treatments and duration of follow up. In particular, Tanizawa et al. and Ash et al. analysed scans (1 or 2 mm thick) every 10 mm of intervals while Best et al. analysed (1 mm thick) scans every 20 mm . Best et al. considered patients enrolled in a double-blind placebo-controlled clinical trial of interferon β-1a while Tanizawa et al. evaluated a non-homogeneous cohort of IPF and other fibrotic ILDs that received multiple drugs (corticosteroids, immunosuppressive agents or pirfenidone) . No information about treatment is reported in the study of Ash et al. . These studies also differ regarding duration of follow up. Tanizawa et al. reported a median follow up time of 1,140 days while Best et al. and Ash et al. reported 547.5 and 465 days respectively.
For this reason, in our study, we tried to re-evaluate the predictive value of CT histogram indexes in assessing survival, analysing 1 mm thick scans without interval in a homogeneous cohort of consecutive IPF patients treated with antifibrotics and with a follow up time of 3 years.
Univariate Cox proportional analysis, as in previous studies, demonstrated that kurtosis, skewness, MLD, HAA% and FA% are significant predictors of mortality. Of all of these, skewness demonstrated the lowest AIC value and, according to our analysis, should be considered as the best CT histogram derived variable in predicting mortality. We also evaluated sensitivity and specificity of these variables with respect to mortality. Based on this, we identified possible cut-off points that we used to graph Kaplan-Meier curves. Although these cut-off points were extracted from a relatively small number of patients and need a validation in larger-group cohorts, we demonstrated how CT-variables and cut-off points may help in identifying patients with different life expectations. This finding could greatly help the clinician in the management of patients, especially if they are not able to correctly perform PFTs.
We also tested the role of skewness, FVC, DLCO and age in predicting mortality in a Multivariate model. As shown in Table 5, although skewness is able to correctly predict mortality, it demonstrates in our model a less significant weight than DLCO, which was confirmed as a strong predictor of mortality.
Our study also evaluated changes after 12 ± 3 from baseline in a small subgroup of 23 patients. Differently from the study of Best et al., no difference was demonstrated, neither for CT quantitative variables nor for FVC changes. A possible explanation of this finding could reside in the treatment the patients received. In our study, subjects homogeneously received antifibrotic therapy resulting in a relative functional and radiological stability after one year. To our knowledge this represents the first study analysing patients treated with pirfenidone and nintedanib with this approach.
This study has some strengths. We analysed the entire volume of the lungs taking into account every thin section in the final calculation of the indexes; every exam was performed with the same protocol and in the same CT-scanner both at baseline and follow up; all patients were treated with antifibrotic drugs reproducing a real clinical setting; patients were followed for a long enough period to assess survival.
This study, however, also has several limitations. The retrospective design and the single-centre nature of the study represent the first great limitation. We also analysed a small number of patients, even if only the study of Best et al. considered a truly greater number of patients with this method. Our analysis considered HRCT scans performed with the same protocol and in the same CT-scanner. This is the reason for the exclusion of a large number of newly-diagnosed IPF patients over the considered period. Moreover, as is known, this kind of analysis may be greatly influenced by the level of inspiration of the patients. An incorrect manoeuvre or the presence of artifacts may erroneously affect quantitative CT indexes. Given the retrospective nature of the study, we did not used a spirometric control during acquisition to ensure an optimal inspiration. However, all CT-scans were reviewed and approved for the inclusion in the study by our radiologists. Another notable and previously mentioned limitation of the study includes the need for some computational steps required for the image processing to get the CT-indexes.
Based on our results, baseline HRCT indexes are able to provide a prediction of survival and can be used as a surrogate marker for prognosis. Given the characteristics of simplicity, objectivity and reproducibility, this method could actually help the clinician and the radiologist in the evaluation of sub-clinical changes and the prognosis of patients in daily clinical practice. Further studies with larger cohorts and an external validation are needed to confirm these data.
Akaike’s information criterion
Diffusing Capacity for carbon monoxide
Forced vital capacity
High attenuation areas
High-resolution computed tomography
Mean Lung Density
Raghu G, Collard HR, Egan JJ, Martinez FJ, Behr J, Brown KK, et al. An Official ATS/ERS/JRS/ALAT Statement: Idiopathic pulmonary fibrosis: Evidence-based guidelines for diagnosis and management. Am J Respir Crit Care Med. 2011;183:788–824.
Palmucci S, Roccasalva F, Puglisi S, Torrisi SE, Vindigni V, Mauro LA, et al. Clinical and radiological features of idiopathic interstitial pneumonias (IIPs): a pictorial review. Insights Imaging. 2014;5:347–64.
Puglisi S, Torrisi SE, Vindigni V, Giuliano R, Palmucci S, Mulè M, Vancheri C. New perspectives on management of idiopathic pulmonary fibrosis. Ther Adv Chronic Dis. 2016;7:108–20.
Puglisi S, Torrisi SE, Giuliano R, Vindigni V, Vancheri C. What we know about the pathogenesis of idiopathic pulmonary fibrosis. Semin Respir Crit Care Med. 2016;37:358–67.
Torrisi SE, Pavone M, Vancheri A, Vancheri C. When to start and when to stop antifibrotic therapies. Eur Respir Rev. 2017;26(145):170053. https://doi.org/10.1183/16000617.0053-2017.
Walsh SL, Calandriello L, Sverzellati N, Wells AU, Hansell DM, On behalf of the UIP observer consort. Interobserver agreement for the ATS/ERS/JRS/ALAT criteria for a UIP pattern on CT. Thorax. 2016;71:45–51.
Bartholmai BJ, Raghunath S, Karwoski RA, Moua T, Rajagopalan S, Maldonado F, et al. Quantitative computed tomography imaging of interstitial lung disease. J Thorac Imaging. 2013;28:298–307.
Kim HJ, Brown MS, Chong D, Gjertson DW, Lu P, Kim HJ, et al. Comparison of the quantitative CT imaging biomarkers of idiopathic pulmonary fibrosis at baseline and early change with an interval of 7 months. Acad Radiol. 2015;22:70–80.
Salaffi F, Carotti M, Di Donato E, Di Carlo M, Ceccarelli L, Giuseppetti G. Computer-aided tomographic analysis of interstitial lung disease (ILD) in patients with systemic sclerosis (SSc). Correlation with pulmonary physiologic tests and patient-centred measures of perceived dyspnea and functional disability. PLoS One. 2016;11:1–13.
Lynch DA, Godwin JD, Safrin S, Starko KM, Hormel P, Brown KK, et al. High-resolution computed tomography in idiopathic pulmonary fibrosis: Diagnosis and prognosis. Am J Respir Crit Care Med. 2005;172:488–93.
Best AC, Meng J, Lynch AM, Bozic CM, Miller D, Grunwald GK, Lynch DA. Idiopathic pulmonary fibrosis: physiologic tests, quantitative CT indexes, and CT visual scores as predictors of mortality. Radiology. 2008;246:935–40.
Jacob J, Bartholmai BJ, Rajagopalan S, Kokosi M, Nair A, Karwoski R, et al. Automated quantitative computed tomography versus visual computed tomography scoring in idiopathic pulmonary fibrosis. J Thorac Imaging. 2016;31:304–11.
Jacob J, Bartholmai BJ, Rajagopalan S, Kokosi M, Nair A, Karwoski R, et al. Mortality prediction in idiopathic pulmonary fibrosis: evaluation of computer-based CT analysis with conventional severity measures. Eur Respir J. 2017;49(1):1601011. https://doi.org/10.1183/13993003.01011-2016.
Humphries SM, Yagihashi K, Huckleberry J, Rho B-H, Schroeder JD, Strand M, et al. Idiopathic pulmonary fibrosis: data-driven textural analysis of extent of fibrosis at baseline and 15-month follow-up. Radiology. 2017;285:270–8.
Iwasawa T, Asakura A, Sakai F, Kanauchi T, Gotoh T, Ogura T, et al. Assessment of prognosis of patients with idiopathic pulmonary fibrosis by computer-aided analysis of CT images. J Thorac Imaging. 2009;24:216–22.
Maldonado F, Moua T, Rajagopalan S, Karwoski RA, Raghunath S, Decker PA, et al. Automated quantification of radiological patterns predicts survival in idiopathic pulmonary fibrosis. Eur Respir J. 2014;43:204–12.
Lee SM, Seo JB, Oh SY, Kim TH, Song JW, Lee SM, Kim N. Prediction of survival by texture-based automated quantitative assessment of regional disease patterns on CT in idiopathic pulmonary fibrosis. Eur Radiol. 2018;28:1293–300.
Jacob J, Bartholmai BJ, Rajagopalan S, Kokosi M, Egashira R, Brun AL, et al. Serial automated quantitative CT analysis in idiopathic pulmonary fibrosis: functional correlations and comparison with changes in visual CT scores. Eur Radiol. 2018;28:1318–27.
Tanizawa K, Handa T, Nagai S, Hirai T, Kubo T, Oguma T, et al. Clinical impact of high-attenuation and cystic areas on computed tomography in fibrotic idiopathic interstitial pneumonias. BMC Pulm Med. 2015;15:74. https://doi.org/10.1186/s12890-015-0069-0.
Ash SY, Harmouche R, Vallejo DLL, Villalba JA, Ostridge K, Gunville R, et al. Densitometric and local histogram based analysis of computed tomography images in patients with idiopathic pulmonary fibrosis. Respir Res. 2017;18(1):45. https://doi.org/10.1186/s12931-017-0527-8.
Departmental funding from the “Progetto Dipartimentale #A” has been used for this study. The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The data are not publicly available due to the containing information that could compromise research participant consent.
SET and SP contributed equally.
Ethics approval and consent to participate
This study (protocol QH- IPF) has been approved by the ethical committee of Policlinico-Vittorio Emanuele Hospital of Catania with letter number 0029883 released on June 2018 by Ethical Committee n. 114/2018/PO.
Consent for publication
SET has received honoraria for lectures from Boehringer Ingelheim and F. Hoffmann-La Roche Ltd. SP has received honoraria for lectures and a financial support for an educational program from Boehringer Ingelheim and F. Hoffmann-La Roche Ltd. CV is part of the F. Hoffmann-La Roche Ltd. Scientific board. He has received consulting fees and/or speaker fees from AstraZeneca, Boehringer Ingelheim, Chiesi, F. Hoffmann-La Roche Ltd. and Menarini. The authors AS, GR, AGT, DF, MG, MP, AV, GS, DS, LAM, EG and AB declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.