Quantitative CT analysis of honeycombing area predicts mortality in idiopathic pulmonary fibrosis with definite usual interstitial pneumonia pattern: A retrospective cohort study

Abstract

Background

Honeycombing on high-resolution computed tomography (HRCT) images is a key finding in idiopathic pulmonary fibrosis (IPF). In IPF, honeycombing area determined by quantitative CT analysis is correlated with pulmonary function test findings. We hypothesized that quantitative CT-derived honeycombing area (HA) might predict mortality in patients with IPF.

Materials and methods

Chest HRCT images of 52 IPF patients with definite usual interstitial pneumonia (UIP) pattern were retrospectively evaluated. Mortality data up to July 31, 2016, were recorded. Using a computer-aided system, HA and percentage of HA (%HA) were measured quantitatively. Predictors of 3-year mortality were evaluated using logistic regression models.

Results

The median %HA, %predicted forced vital capacity (FVC) and composite physiologic index (CPI) were 3.8%, 83.6%, and 33.6, respectively. According to GAP (gender, age, and physiology) stage, 20, 14, and 5 patients were classified under stages I-II-III, respectively. Percentage of HA was significantly correlated with %FVC, CPI, and GAP stage (all, p < 0.001). In univariate analysis, %HA, %FVC, and CPI were statistically significant predictors of mortality. In multivariate analysis using the stepwise regression method, only %HA (odds ratio [OR], 1.27; p = 0.011) was a significant independent predictors of mortality. Patients with %HA ≥ 4.8% had significantly lower survival rates than those with lesser %HA (median survival time, 1.3 vs 5.0 years; log-rank test; p < 0.001).

Conclusion

Quantitative CT-derived HA might be an important and independent predictor of mortality in IPF patients with definite UIP pattern.

About the Authors:

Hiroaki Nakagawa

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

ORCID logo http://orcid.org/0000-0003-3830-765X

Emiko Ogawa

Roles Resources, Writing – review & editing

Affiliation: Health Administration Center, Shiga University of Medical Science, Shiga, Japan

Kentaro Fukunaga

Roles Resources, Writing – review & editing

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

Daisuke Kinose

Roles Resources, Writing – review & editing

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

Masafumi Yamaguchi

Roles Resources, Writing – review & editing

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

Taishi Nagao

Roles Resources, Writing – review & editing

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

Sachiko Tanaka-Mizuno

Roles Formal analysis, Resources, Writing – review & editing

Affiliation: Department of Medical Statistics, Shiga University of Medical Science, Shiga, Japan

Yasutaka Nakano

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliation: Division of Respiratory Medicine, Department of Internal Medicine, Shiga University of Medical Science, Shiga, Japan

ORCID logo http://orcid.org/0000-0001-7516-8690

Introduction

Idiopathic pulmonary fibrosis (IPF) is a specific form of chronic and progressive fibrosing interstitial pneumonia of unknown cause [1–3]. Honeycombing is one of the key findings on high-resolution computed tomography (HRCT) images in IPF with definite usual interstitial pneumonia (UIP) pattern. The extent of honeycombing area on HRCT images is an independent predictor of mortality in IPF [1, 4, 5]. Therefore, quantifying the extent of honeycombing on HRCT images allows better management of patients with IPF. Experts have evaluated the extent of honeycombing using CT visual scores [6–9]. However, this evaluation process is not uniform, yielding only fair to moderate agreement regarding honeycombing area among expert radiologists and pulmonologists [10]. The 2011 international statement also mentions the lack of uniformity in approach towards quantification of honeycombing area [4]. Recently, there are some reports about computer-based quantitative CT analyses [11–18]. These methods may be the useful and reliable for objective assessment of HRCT findings including honeycombing. However, the algorithms of their methods remain unclear, and not everyone can use these methods.

We recently developed a computer-aided system for quantitative CT analysis of honeycombing area in IPF [19]. This system has no interobserver errors and could, therefore, address the issues associated with current quantification methods for honeycombing area. Furthermore, everyone can use this system through the public-domain computer program, ImageJ.

Pulmonary function tests (PFTs) are a standard approach for objective monitoring of patients with IPF. Both forced vital capacity (FVC) and diffusing capacity of the lungs for carbon monoxide (DLCO) are strongly correlated with disease extent determined histologically or estimated visually using CT images [1]. However, in patients with IPF, interpretation of PFT findings is confounded by coexisting emphysema [20, 21]. Therefore, the composite physiologic index (CPI) was developed to improve on previous prognostic measures in IPF by adjusting for emphysema [20]. There exists another composite index, the GAP stage, which uses gender (G), age (A), and physiology (P) (FVC and DLCO) variables for predicting mortality in IPF [22]. Although PFTs are very important in both these indexes, one should be aware that patients with severe IPF may not be able to perform PFTs.

The aim of this study was to assess the relationship between honeycombing area determined by quantitative CT analysis and mortality in patients with IPF. We hypothesized that quantitative CT-derived honeycombing area (HA) might predict survival in IPF with definite UIP pattern.

Materials and methods

Study design and patients

All IPF patients with definite UIP pattern—diagnosed in accordance with the 2011 IPF guidelines—who visited the outpatient clinic of our hospital between April 2011 and March 2014 were enrolled in this study. Patients with secondary interstitial pneumonitis (such as collagen vascular disease, pneumonia caused by occupational or environmental exposure, chronic hypersensitivity pneumonitis, and drug-induced pneumonia) were excluded. Mortality and survival data up to July 31, 2016, were recorded. In order to evaluate 3-year mortality, patients with follow-up period < 3 years were excluded. For comparison of clinical features, patients with IPF were classified as “alive” or “deceased” on the basis of survival status at the end of the 3-year follow-up period. The study protocol conformed to the Declaration of Helsinki and was approved by the ethics committee of Shiga University of Medical Science (approval number 24–182; May 1, 2013), with a waiver of the need for informed consent because of the retrospective study design.

Chest HRCT images at the first visit were acquired using Toshiba Aquillion ONE (Toshiba Medical Systems Corp., Otawara, Tochigi, Japan) or SOMATOM Sensation Cardiac (Siemens Japan K.K., Shinagawa, Tokyo, Japan). The technical parameter for Toshiba Aquillion ONE was as follows: slice thickness 0.5–1.0 mm, x-ray voltage 120kVp, x-ray median current 200mA (range 80-350mA), rotation speed 0.5 sec, and reconstruction algorithm FC52. The technical parameter for SOMATOM Sensation Cardiac: slice thickness 0.5–1.0 mm, x-ray voltage 120kVp, x-ray median current 288mA (range 226-481mA), rotation speed 0.375 sec, and reconstruction algorithm B70f. The images used for the analysis were selected with 10-mm intervals. Among the enrolled patients, only PFT data acquired within 3 months of the selected CT images were analyzed. For these patients, CPI and GAP stage were determined as previously described [20, 22]. The PFTs had been performed in accordance with the Japanese Respiratory Society guidelines [23].

A computer-aided method for quantitative CT analysis of honeycombing area was used to automatically measure honeycombing areas on whole HRCT images, as previously described [19]. Briefly, quantitative CT analysis of honeycombing area involved the following process: (1) low-attenuation areas enclosed by thick walls were extracted as honeycombing areas; (2) the sum of honeycombing areas in all slices was calculated as the HA. The values of HA and percentage of HA to total lung area (%HA) were measured for each patient (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Honeycombing area detected by quantitative CT analysis.

A 76-year-old man with idiopathic pulmonary fibrosis. (A) Original CT image was acquired at the level of the lower lungs using the lung window setting. (B) Honeycombing area was detected by quantitative CT analysis at the same level.

https://doi.org/10.1371/journal.pone.0214278.g001

Statistical analysis

Data are expressed as median and interquartile range (IQR). The correlation of HA with PFT findings and patient characteristics was evaluated by Spearman’s correlation coefficient analysis.

For assessment of 3-year mortality, the clinical features of alive and deceased patients were compared using the Mann–Whitney U test. The effect of each variable on 3-year mortality was evaluated by univariate logistic regression analysis adjusted by age, sex, body mass index (BMI), and pack-years. The model with the greatest predictive utility was determined by multivariable logistic regression analysis adjusted by age, sex, BMI, and pack-years using the stepwise method. The cutoff value for %HA was determined on the basis of the area under the receiver operating characteristic curve and Youden’s index [24]. Survival analysis was performed in accordance with the methods of Kaplan–Meier, with the endpoint being death. Differences in mortality were assessed by the log-rank test. All statistical analyses were performed using JMP version 9.0.2 (SAS Institute, Cary, NC). For all tests, p values < 0.05 were considered statistically significant.

Results

Patient characteristics

Of the 72 consecutive IPF patients with definite UIP pattern who visited the outpatient clinic of our hospital between April 2011 and March 2014, one was excluded for the lack of HRCT data at the first visit; additionally, 19 patients with follow-up period < 3 years were excluded. Thus, 52 patients with IPF were finally enrolled in this study. Among the 52 enrolled patients, 45 had corresponding PFT data acquired within 3 months of the selected CT images. The demographic, clinical, and physiologic characteristics of these 52 patients are summarized in Table 1.

[Figure omitted. See PDF.]

Table 1. Demographic, clinical, and physiologic characteristics of patients with idiopathic pulmonary fibrosis.

https://doi.org/10.1371/journal.pone.0214278.t001

CT quantitation of honeycombing area

The median (IQR) %HA was 3.8% (2.2–5.5%). The relationships between %HA and patient characteristics are shown in Table 2; %HA was significantly correlated with FVC, DLCO, and CPI (ρ, -0.53, p < 0.001; ρ, -0.46, p = 0.004; and ρ, 0.60, p < 0.001, respectively). Furthermore, %HA was significantly correlated with GAP stage (ρ, 0.51; p < 0.001).

[Figure omitted. See PDF.]

Table 2. Relationship between patient characteristics and %HA determined by quantitative CT analysis.

https://doi.org/10.1371/journal.pone.0214278.t002

3-year mortality and baseline parameters

Of the 52 patients, 23 were deceased and 29 alive at the end of the 3-year follow-up period. Deceased group had the lower BMI and %FVC (p = 0.028, and 0.035, respectively; Table 3) and the higher GAP stage and %HA (p = 0.005, and 0.004, respectively; Table 3). Table 4 presents the results of univariate logistic regression analysis of correlation between baseline variables and 3-year mortality. Among the evaluated clinical, physiologic, and radiographic characteristics, %HA, %FVC, CPI, and GAP stage were statistically significant predictors of mortality. Although %HA correlated with CPI and GAP stage (Table 2), multivariate logistic regression analysis revealed that only %HA (odds ratio [OR], 1.27; 95% confidence interval [CI], 1.05–1.62; p = 0.011) was a significant predictor of 3-year mortality.

[Figure omitted. See PDF.]

Table 3. Comparison of clinical features between alive and deceased patients.

https://doi.org/10.1371/journal.pone.0214278.t003

[Figure omitted. See PDF.]

Table 4. Results of univariate logistic regression analysis for predictors of mortality.

https://doi.org/10.1371/journal.pone.0214278.t004

Prediction of mortality by %HA

From the results of receiver operating characteristic analysis, the Youden’s index and accuracy of %HA were estimated to be 4.8% and 73.1%, respectively (sensitivity, 86.2%; specificity, 56.5%; S1 Fig). The Kaplan–Meier survival curve demonstrated a significant difference in mortality, which were defined by a cutoff point of 4.8% for %HA (Fig 2). The median (IQR) survival times were 1.3 (0.8–2.8) and 5.0 (2.1–5.9) years, respectively, among patients with %HA greater and lesser than 4.8% (log-rank test; p < 0.001).

[Figure omitted. See PDF.]

Fig 2. Kaplan–Meier plot of survival probability.

The median survival times (interquartile ranges) were 1.3 (0.8–2.8) and 5.0 (2.1–5.9) years among patients with greater (≥ 4.8%; blue line) and smaller (< 4.8%; red line) %HA (log-rank test; p < 0.001), relative to the cutoff. %HA = computed-tomography-derived %honeycombing area.

https://doi.org/10.1371/journal.pone.0214278.g002

Discussion

The 2011 international statement mentions that extent of honeycombing on HRCT images is an independent predictor of mortality in IPF, although there is no uniformity in approach to quantitation [4]. We hypothesized that honeycombing area measured by quantitative CT analysis might be associated with an increased risk of mortality in IPF. The present results demonstrated that %HA was significantly correlated with 3-year mortality (Table 4). Furthermore, %HA was an independent predictor of 3-year mortality as the result of multivariate analysis using the stepwise regression method. These findings show that %HA might be a quantitative prognostic factor in IPF.

Both FVC and DLCO are considered the gold standards for assessment of IPF. However, it is important to consider that, in patients with IPF, FVC and DLCO findings are confounded by coexisting emphysema [20, 21]. Emphysema is present in approximately one-third of patients with IPF [25]. FVC measurements may not be appropriate for monitoring disease progression in patients with IPF and severe emphysema [26]. In patients with IPF with coexistent emphysema, CPI is thought to be a more accurate prognostic determinant than any individual parameter of PFTs [20]. Additionally, in IPF, GAP stage is thought to be a better prognostic determinant than any individual parameter of PFTs [22, 27]. In the present study, %HA was significantly correlated with not only FVC and DLCO but also CPI and GAP stage (Table 2). This is the retrospective cohort study of the consecutive IPF patients with HRCT scan data. Unfortunately, 7 patients did not have PFT within 3 months of the selected HRCT. It is inaccurate to discuss the relationship between selected HRCT and PFT performed at different time. Therefore, PFT of 7 patients was treated as missing data.

In order to identify better prognostic factors in IPF, we performed logistic regression analysis of correlation of 3-year mortality with %HA and other variables, including PFT findings. We selected the 3-year time point for mortality because several retrospective longitudinal studies have reported that the median survival time from the time of IPF diagnosis is between 2 and 3 years [4, 27–31]. Among the evaluated variables, only %HA was found to be significantly correlated with 3-year mortality upon multivariate logistic regression analysis using the stepwise method. This result showed that %HA was an independent and useful variable for prediction of 3-year mortality in the present study population. Among FVC, forced expiratory volume in 1 s (FEV1), DLCO, CPI, GAP stage, and %HA, all variables exhibited strong correlations with each other. We, therefore, performed bivariate logistic regression analysis for the relationship of FVC, FEV1, DLCO, CPI, and GAP stage with %HA (S1 Table) and found that %HA was a significant predictor of mortality in combination with each of the variables except CPI and GAP stage. On the other hand, in combination with each other, neither %HA nor CPI or GAP stage was a significant predictor of mortality. These results seem to be derived from the strong correlation between CPI and %HA, GAP stage and %HA. These results show that in addition to CPI and GAP stage, %HA was a useful variable for prediction of 3-year mortality.

In evaluating the mortality in IPF, the analysis for time to death is also helpful. In univariate Cox proportional hazards models, %HA, %FVC, %FEV1, and CPI were statistically significant predictors of mortality (S2 Table). In multivariate Cox proportional hazards models including all variables using the stepwise regression method, only %HA (hazard ratio [HR], 1.13; 95% CI, 1.04–1.21; p = 0.006) was significant predictor of mortality. The results of bivariate Cox proportional hazards models (S3 Table) were almost similar to those of bivariate logistic regression analysis.

Because we wanted to evaluate honeycombing area by quantitative CT analysis, the present study population only included patients with IPF with definite usual interstitial pneumonia pattern on HRCT images. Although our study population was not representative of the entire spectrum of IPF, the relationship between mortality and GAP stage in our population was almost similar to that reported in previous studies [22, 32] on IPF (S4 Table). In evaluation of prognosis, we believe that the present results might be comparable to the trend among the general IPF population.

According to the review of computer-based quantitative CT image analyses in IPF, there were 3 types of quantitative CT analyses; the density histogram analysis, the density mask analysis like our method, and texture classification analysis [33]. Including semi-quantitative analysis by visual scoring, each analysis has its merits that can be used in clinical settings. The texture classification analysis may be the most useful and reliable method for objective management of IPF. However, the specific CT index and the algorithm remain unclear [33]. It is desirable to establish a method which is transparent and which is freely available for widespread use, such as the low attenuation area (LAA) for chronic obstructive pulmonary disease (COPD) [34, 35].

The present method of quantitative CT analysis of percentage of honeycombing area has important advantages relative to previously reported methods. First, since this method does not require PFT data (especially DLCO), honeycombing area can be quantified in almost all cases of IPF. Second, this method does not involve interobserver or intraobserver errors and can be used by anybody as a public-domain computer program. Third, it allows quantitative assessment of prognosis in IPF. Our method using 10-mm interval images is not perfect. However, we think it is not inferior to the past methods using 3–4 CT images of fixed position such as aortic arch, carina, or right inferior pulmonary vein.

There are several limitations to the present study. First, this was a retrospective study with a small sample size, conducted at a single institution. Second, our patients did not receive uniform treatment for IPF. The rates of anti-fibrotic drugs usage were not different between the alive group and deceased group (Table 3). Third, we analyzed only IPF patients with definite UIP pattern to investigate whether the extent of honeycombing area was a prognostic predictor or not. In other words, we analyzed only the IPF patients with honeycombing area. We did not analyzed the patients without honeycombing area. A prospective multicenter study on quantitative analysis of honeycombing area is needed to confirm the usefulness of %HA in management of patients with IPF.

Conclusions

Relative extension of honeycombing in comparison to lung volume measured by quantitative CT analysis might be an important and independent predictor of mortality in IPF patients with definite UIP pattern. Furthermore, this analysis might provide new insights into disease progression in the future.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Results of logistic regression analysis for predictors of mortality and relationship between %HA and other variables.

https://doi.org/10.1371/journal.pone.0214278.s001

(DOCX)

S2 Table. Results of univariate Cox proportional hazards models for predictors of mortality.

https://doi.org/10.1371/journal.pone.0214278.s002

(DOCX)

S3 Table. Results of Cox proportional hazards models for predictors of mortality and relationship between %HA and other variables.

https://doi.org/10.1371/journal.pone.0214278.s003

(DOCX)

S4 Table. Comparison of predicted or observed mortality according to GAP stage.

https://doi.org/10.1371/journal.pone.0214278.s004

(DOCX)

S1 Fig. The receiver operating characteristic curve with the relative Youden’s index.

The area under the curve with 4.8% cutoff point of %HA was 0.735. The sensitivity, the specificity, and the accuracy were 86.2%, 56.5%, and 73.1%, respectively). %HA = computed-tomography-derived %honeycombing area.

https://doi.org/10.1371/journal.pone.0214278.s005

(DOCX)

Citation: Nakagawa H, Ogawa E, Fukunaga K, Kinose D, Yamaguchi M, Nagao T, et al. (2019) Quantitative CT analysis of honeycombing area predicts mortality in idiopathic pulmonary fibrosis with definite usual interstitial pneumonia pattern: A retrospective cohort study. PLoS ONE 14(3): e0214278. https://doi.org/10.1371/journal.pone.0214278

References

1. 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–493. pmid:15894598

2. Nishimura K, Kitaichi M, Izumi T, Nagai S, Kanaoka M and Itoh H. Usual interstitial pneumonia: histologic correlation with high-resolution CT. Radiology. 1992; 182: 337–342. pmid:1732946

3. American Thoracic Society. Idiopathic pulmonary fibrosis: diagnosis and treatment. International consensus statement. American Thoracic Society (ATS), and the European Respiratory Society (ERS). Am J Respir Crit Care Med. 2000; 161: 646–664. pmid:10673212

4. 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. pmid:21471066

5. Sumikawa H, Johkoh T, Colby TV, Ichikado K, Suga M, Taniguchi H, et al. Computed tomography findings in pathological usual interstitial pneumonia: relationship to survival. Am J Respir Crit Care Med. 2008; 177: 433–439. pmid:17975197

6. Oda K, Ishimoto H, Yatera K, Naito K, Ogoshi T, Yamasaki K, et al. High-resolution CT scoring system-based grading scale predicts the clinical outcomes in patients with idiopathic pulmonary fibrosis. Respir Res. 2014; 15: 10. pmid:24479411

7. Best AC, Meng J, Lynch AM, Bozic CM, Miller D, Grunwald GK, et al. Idiopathic pulmonary fibrosis: physiologic tests, quantitative CT indexes, and CT visual scores as predictors of mortality. Radiology. 2008; 246: 935–940. pmid:18235106

8. Rosas IO, Yao J, Avila NA, Chow CK, Gahl WA and Gochuico BR. Automated quantification of high-resolution CT scan findings in individuals at risk for pulmonary fibrosis. Chest. 2011; 140: 1590–1597. pmid:21622544

9. Kazerooni EA, Martinez FJ, Flint A, Jamadar DA, Gross BH, Spizarny DL, et al. Thin-section CT obtained at 10-mm increments versus limited three-level thin-section CT for idiopathic pulmonary fibrosis: correlation with pathologic scoring. AJR Am J Roentgenol. 1997; 169: 977–983. pmid:9308447

10. Watadani T, Sakai F, Johkoh T, Noma S, Akira M, Fujimoto K, et al. Interobserver variability in the CT assessment of honeycombing in the lungs. Radiology. 2013; 266: 936–944. pmid:23220902

11. 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–222. pmid:19704326

12. Iwasawa T, Ogura T, Sakai F, Kanauchi T, Komagata T, Baba T, et al. CT analysis of the effect of pirfenidone in patients with idiopathic pulmonary fibrosis. Eur J Radiol. 2014; 83: 32–38. pmid:22465123

13. 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–212. pmid:23563264

14. Park SO, Seo JB, Kim N, Lee YK, Lee J and Kim DS. Comparison of usual interstitial pneumonia and nonspecific interstitial pneumonia: quantification of disease severity and discrimination between two diseases on HRCT using a texture-based automated system. Korean J Radiol. 2011; 12: 297–307. pmid:21603289

15. Yoon RG, Seo JB, Kim N, Lee HJ, Lee SM, Lee YK, et al. Quantitative assessment of change in regional disease patterns on serial HRCT of fibrotic interstitial pneumonia with texture-based automated quantification system. Eur Radiol. 2013; 23: 692–701. pmid:22918563

16. 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.

17. 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: Validation Against Pulmonary Function. J Thorac Imaging. 2016; 31: 304–311. pmid:27262146

18. Colombi D, Dinkel J, Weinheimer O, Obermayer B, Buzan T, Nabers D, et al. Visual vs Fully Automatic Histogram-Based Assessment of Idiopathic Pulmonary Fibrosis (IPF) Progression Using Sequential Multidetector Computed Tomography (MDCT). PLoS One. 2015; 10: e0130653. pmid:26110421

19. Nakagawa H, Nagatani Y, Takahashi M, Ogawa E, Tho NV, Ryujin Y, et al. Quantitative CT analysis of honeycombing area in idiopathic pulmonary fibrosis: Correlations with pulmonary function tests. Eur J Radiol. 2016; 85: 125–130. pmid:26724656

20. Wells AU, Desai SR, Rubens MB, Goh NS, Cramer D, Nicholson AG, et al. Idiopathic pulmonary fibrosis: a composite physiologic index derived from disease extent observed by computed tomography. Am J Respir Crit Care Med. 2003; 167: 962–969. pmid:12663338

21. Bodlet A, Maury G, Jamart J and Dahlqvist C. Influence of radiological emphysema on lung function test in idiopathic pulmonary fibrosis. Respir Med. 2013; 107: 1781–1788. pmid:24051272

22. Ley B, Ryerson CJ, Vittinghoff E, Ryu JH, Tomassetti S, Lee JS, et al. A multidimensional index and staging system for idiopathic pulmonary fibrosis. Ann Intern Med. 2012; 156: 684–691. pmid:22586007

23. [Guideline of respiratory function tests—spirometry, flow-volume curve, diffusion capacity of the lung]. Nihon Kokyuki Gakkai Zasshi. 2004; Suppl: 1–56.

24. Schisterman EF, Perkins NJ, Liu A and Bondell H. Optimal cut-point and its corresponding Youden Index to discriminate individuals using pooled blood samples. Epidemiology. 2005; 16: 73–81. pmid:15613948

25. Mejia M, Carrillo G, Rojas-Serrano J, Estrada A, Suarez T, Alonso D, et al. Idiopathic pulmonary fibrosis and emphysema: decreased survival associated with severe pulmonary arterial hypertension. Chest. 2009; 136: 10–15. pmid:19225068

26. Cottin V, Hansell DM, Sverzellati N, Weycker D, Antoniou KM, Atwood M, et al. Effect of Emphysema Extent on Serial Lung Function in Patients with Idiopathic Pulmonary Fibrosis. Am J Respir Crit Care Med. 2017; 196: 1162–1171. pmid:28657784

27. Bjoraker JA, Ryu JH, Edwin MK, Myers JL, Tazelaar HD, Schroeder DR, et al. Prognostic significance of histopathologic subsets in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 1998; 157: 199–203. pmid:9445300

28. King TE Jr., Schwarz MI, Brown K, Tooze JA, Colby TV, Waldron JA Jr., et al. Idiopathic pulmonary fibrosis: relationship between histopathologic features and mortality. Am J Respir Crit Care Med. 2001; 164: 1025–1032. pmid:11587991

29. Nicholson AG, Colby TV, du Bois RM, Hansell DM and Wells AU. The prognostic significance of the histologic pattern of interstitial pneumonia in patients presenting with the clinical entity of cryptogenic fibrosing alveolitis. Am J Respir Crit Care Med. 2000; 162: 2213–2217. pmid:11112140

30. Rudd RM, Prescott RJ, Chalmers JC and Johnston ID. British Thoracic Society Study on cryptogenic fibrosing alveolitis: Response to treatment and survival. Thorax. 2007; 62: 62–66. pmid:16769717

31. Flaherty KR, Toews GB, Travis WD, Colby TV, Kazerooni EA, Gross BH, et al. Clinical significance of histological classification of idiopathic interstitial pneumonia. Eur Respir J. 2002; 19: 275–283. pmid:11866008

32. Kim ES, Choi SM, Lee J, Park YS, Lee CH, Yim JJ, et al. Validation of the GAP score in Korean patients with idiopathic pulmonary fibrosis. Chest. 2015; 147: 430–437. pmid:25211374

33. Ohkubo H, Nakagawa H and Niimi A. Computer-based quantitative computed tomography image analysis in idiopathic pulmonary fibrosis: A mini review. Respir Investig. 2018; 56: 5–13. pmid:29325682

34. Sakai N, Mishima M, Nishimura K, Itoh H and Kuno K. An automated method to assess the distribution of low attenuation areas on chest CT scans in chronic pulmonary emphysema patients. Chest. 1994; 106: 1319–1325. pmid:7956377

35. Nakano Y, Muro S, Sakai H, Hirai T, Chin K, Tsukino M, et al. Computed tomographic measurements of airway dimensions and emphysema in smokers. Correlation with lung function. Am J Respir Crit Care Med. 2000; 162: 1102–1108. pmid:10988137