Article Text

Download PDFPDF

Acute exacerbation in rheumatoid arthritis-associated interstitial lung disease: a retrospective case control study
  1. Hironao Hozumi1,
  2. Yutaro Nakamura1,
  3. Takeshi Johkoh2,
  4. Hiromitsu Sumikawa3,
  5. Thomas V Colby4,
  6. Masato Kono1,
  7. Dai Hashimoto1,
  8. Noriyuki Enomoto1,
  9. Tomoyuki Fujisawa1,
  10. Naoki Inui1,
  11. Takafumi Suda1,
  12. Kingo Chida1
  1. 1Second Division, Department of Internal Medicine, Hamamatsu University School of Medicine, Hamamatsu, Japan
  2. 2Department of Radiology, Kinki Central Hospital of Mutual Aid Association of Public School Teachers, Itami-City, Hyogo, Japan
  3. 3Diagnostic and Interventional Radiology, Osaka University Graduate School of Medicine, Suita, Osaka, Japan
  4. 4Department of Laboratory Medicine and Pathology, Mayo Clinic Arizona, Scottsdale, Arizona, USA
  1. Correspondence to Dr Yutaro Nakamura; nakayuta{at}hama-med.ac.jp

Abstract

Objectives To investigate the risk factors and prognosis associated with acute exacerbation (AE) in patients with rheumatoid arthritis-associated interstitial lung disease (RA-ILD).

Design A retrospective case–control study.

Setting A single academic hospital.

Participants 51 consecutive patients diagnosed with RA-ILD between 1995 and 2012. All patients fulfilled the diagnostic criteria of the American College of Rheumatology for RA. ILD was diagnosed on the basis of clinical presentation, pulmonary function tests, high-resolution CT (HRCT) findings and lung biopsy findings.

Main outcome measures Overall survival and cumulative AE incidence were analysed using Kaplan-Meier method. Cox hazards analysis was used to determine significant variables associated with AE occurrence and survival status.

Results A total of 11 patients (22%) developed AE, with an overall 1-year incidence of 2.8%. Univariate analysis revealed that older age at ILD diagnosis (HR 1.11; 95% CI 1.02 to 1.21; p=0.01), usual interstitial pneumonia (UIP) pattern on HRCT (HR 1.95; 95% CI 1.07 to 3.63; p=0.03) and methotrexate usage (HR 3.04; 95% CI 1.62 to 6.02; p=0.001) were associated with AE. Of 11 patients who developed AE during observation period, 7 (64%) died of initial AE. In survival, AE was a prognostic factor for poor outcome (HR 2.47; 95% CI 1.39 to 4.56; p=0.003).

Conclusions In patients with RA-ILD, older age at ILD diagnosis, UIP pattern on HRCT and methotrexate usage are associated with the development of AE. Furthermore, AE has a serious impact on their survival.

  • Rheumatology

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Article summary

Strengths and limitations of this study

  • Acute exacerbation has a serious impact on the survival of patients with rheumatoid arthritis-associated interstitial lung disease.

  • Given its retrospective study design, it is subject to several possible biases. Therefore, prospective studies are necessary to confirm our results.

Introduction

Acute exacerbation (AE) is a recently established and an increasingly recognised occurrence in idiopathic pulmonary fibrosis (IPF).1 AE is characterised by acute deterioration in respiratory status, with newly developed bilateral ground-glass opacities and/or consolidations on chest radiographs or CT scans. It should be in the absence of other alternative causes such as infection, left heart failure, pulmonary embolism or an identifiable cause of lung injury. AE reportedly occurs not only in patients with IPF but also in patients with other interstitial lung diseases (ILDs), including idiopathic non-specific interstitial pneumonia (NSIP), collagen vascular disease-associated ILDs (CVD-ILDs) and other forms of ILD.2–4 The in-hospital mortality associated with AE in patients with CVD-ILD was demonstrated to be as high as that in patients with IPF,1–4 suggesting that AE may be associated with poor prognosis in patients with CVD-ILD.

Rheumatoid arthritis (RA) is a chronic, systemic, autoimmune disorder of unknown aetiology that primarily involves joints.5 ILD is one of the common extra-articular manifestations.6 ,7 The reported prevalence is variable (1–58%) and depends on the detection and diagnostic method, or the selected population.8–13 We and Park et al4 previously reported that RA-associated ILD (RA-ILD) was the most common CVD-ILD associated with AE.3 However, the risk factors and prognosis associated with AE in patients with RA-ILD are not clarified.

In the current study, we attempted to elucidate the cumulative incidence of AE, its risk factors and prognostic factors in patients with RA-ILD.

Materials and methods

Subjects

This study was approved by the Institutional Review Board of the Hamamatsu University School of Medicine (approval number 24–191). Patient approval or informed consent was waived because the study involved a retrospective review of patient records and images. We retrospectively reviewed 51 consecutive patients diagnosed with RA-ILD between 1995 and 2012 at Hamamatsu University Hospital in Japan. All patients fulfilled the diagnostic criteria of the American College of Rheumatology for RA.14 Patients with other coexisting CVD were excluded. In addition, because our aim was to investigate the features of AE among patients with a chronic course of RA-ILD, patients with no evidence of chronic ILD were also excluded. Chronic ILD was defined as the ILD which had been stable for over 3 months.

ILD was diagnosed on the basis of clinical presentation, pulmonary function tests, high-resolution CT (HRCT) findings and lung biopsy findings. HRCT findings, such as bilateral areas with ground-glass attenuation, reticular opacities and honeycomb patterns, were interpreted and defined as ILD by a consensus between radiologists and pulmonologists. All cases underwent transbronchial lung biopsy and bronchoalveolar lavage, and 21 (41%) cases underwent surgical lung biopsy (SLB) to definitively diagnose ILD or rule out other diseases. The patients with environmental exposures, suspected of drug-induced pneumonia (ILD developed within 1 year after initiation of new drug) or with other known causes of ILD were excluded after considering exposure history and the findings of appropriate tests and histopathological examinations.

Acute exacerbation

AE was defined using recently proposed criteria1 that were slightly modified for adaptation to RA-ILD: previous diagnosis of RA-ILD; unexplained worsening or development of dyspnoea within 30 days of onset; new bilateral ground-glass abnormalities and/or consolidation superimposed on a reticular or honeycomb pattern on HRCT; no evidence of pulmonary infection on negative respiratory culture, including endotracheal aspirate or bronchoalveolar lavage and serological test results for respiratory pathogens; exclusion of alternative causes such as left heart failure, pulmonary embolism or an identifiable cause of lung injury. Patients with AE were required to meet all five criteria. In our institution, cultures of sputum, blood, urine and bronchoalveolar lavage fluid, and serological tests examined for mycobacteria, fungi, bacteria and some viruses were routinely performed. Echocardiography and, if necessary, CT scanning with intravenous contrast were performed to rule out left heart failure or pulmonary thromboembolism. The patients who developed acute pneumonitis within 1 year after initiation of drug for RA were excluded because drug-induced pneumonitis was not completely ruled out. Patients who did or did not develop AE during the observation period were classified into the AE and non-AE groups, respectively.

Data collection

All clinical and laboratory data were collected from medical records. The observation period was calculated from the date of diagnosis of RA-ILD to the last visit. The AE-free period was defined as the time elapsed between the date of RA-ILD diagnosis and the first AE occurrence (AE group) or the last visit (non-AE group).

Review of radiographic findings

HRCT images taken at the time of RA-ILD diagnosis were reviewed. These images comprised 1–2.5 mm collimation sections at 10 mm intervals. They were reconstructed by a high spatial frequency algorithm, and were displayed at window settings appropriate for viewing the lung parenchyma (window level, −600 to −800 Hounsfield units; window width, 1200–2000 Hounsfield units). HRCT images were randomised and reviewed independently by two expert chest radiologists (with 23 and 12 years of experience) who were unaware of the related clinical information. Interobserver agreement was classified as follows: poor (κ=0–0.20), fair (κ=0.21–0.40), moderate (κ=0.41–0.60), good (κ=0.61–0.80) and excellent (κ=0.81–1.00).

RA-ILD on HRCT was classified as a usual interstitial pneumonia (UIP) pattern or a non-UIP pattern according to recent guideline with slight modification.15 Briefly, a UIP pattern on HRCT was characterised by the following: subpleural and basal predominance; reticular abnormalities; honeycombing with or without traction bronchiectasis and the absence of features listed as inconsistent with a UIP pattern; including upper or mid-lung predominance, peribronchovascular predominance, extensive ground-glass abnormalities, profuse micronodules, discrete cysts, diffuse mosaic attenuation/air-trapping and consolidation in bronchopulmonary segment(s)/lobe(s). If the HRCT pattern did not meet the criteria, it was interpreted as a non-UIP pattern. Disagreements regarding HRCT interpretation were resolved by a consensus between both radiologists.

Review of histopathological findings

SLB specimens were obtained from at least two sites. A diagnosis of RA-ILD was originally made in all cases on the basis of histological features evaluated by a lung pathologist at our hospital, correlated with the clinical and radiological findings. All SLB specimens were also reviewed by a second lung pathologist with 36 years of experience. The histological classification of interstitial pneumonia was based on the consensus statement criteria for idiopathic interstitial pneumonias,16 and histological patterns that could not be classified according to the criteria were categorised as unclassified interstitial pneumonia.

Statistical analysis

All values are expressed as median (range) or number (%). The Mann-Whitney U test was used for non-parametric comparisons involving continuous data, whereas Fisher's exact test was used for comparing categorical data. Interobserver agreement on HRCT pattern was analysed using the κ statistic test. Survival was evaluated using the Kaplan-Meier method and the log-rank test according to observation period. The cumulative AE incidence was obtained from a Kaplan-Meier survival curve by treating AE as the death variable according to AE-free period. The 1-year AE incidence was based on person-year method and calculated by dividing the number of patients who developed AE by the total of AE-free period (year) of all the patients. Cox hazards analysis was used to determine significant variables associated with survival and AE occurrence. In all analyses, p<0.05 was considered statistically significant. All data were analysed using commercially available software (JMP V.9.0.3a, SAS Institute Inc, Cary, North Carolina, USA).

Results

Patient characteristics

A total of 82 patients with RA-ILD were identified from medical records at Hamamatsu University Hospital. Of these, 11 patients with other coexisting CVDs, 8 with no evidence of chronic ILD and an initial acute/subacute course of ILD and 3 suspected of drug-induced pneumonitis were excluded. In addition, nine patients for whom no initial HRCT images were available for review were also excluded. The remaining 51 patients with RA-ILD were included in this study (table 1 and figure 1).

Table 1

Clinical characteristics

Figure 1

Number of patients included in this study. Eighty-two patients with rheumatoid arthritis-associated interstitial lung disease (RA-ILD) were assessed for eligibility. From these, 31 patients were excluded, and 51 patients with chronic course of RA-ILD were included. During observation period, 11 patients developed acute exacerbation.

The median age at the onset of RA and diagnosis of RA-ILD was 61 (range 28–82) years and 62 (range 31–83) years, respectively. The median observation period and AE-free period was 8.5 (range 1–17) years and 7 (range 1–17) years, respectively.

During the observation period, 11 patients (22%) developed AE. The median age at the onset of AE was 72 (range 60–86) years. There were no statistically significant differences in age at the onset of RA, age at diagnosis of RA-ILD, sex, observation period, AE-free period, smoking habits, predicted forced vital capacity at diagnosis of RA-ILD and arterial oxygen pressure at diagnosis of RA-ILD between the AE and non-AE groups. Also there were no statistically significant differences in articular disease activities (RF and disease activity score 28 C reactive protein (DAS28-CRP)) at the last visit (non-AE group) or the first AE occurrence (AE group).

In the AE group, 7 of 11 patients (64%) died of respiratory failure caused by initial AE during the observation period. In the non-AE group, 2 of 40 patients (5%) died of respiratory failure. Regarding death caused by respiratory failure, there was a statistically significant difference in the mortality rate between the AE and non-AE groups (p=0.02).

Extra-articular manifestations except for ILD had been observed in 8 (16%) of 51 patients (pericarditis in 2 patients, pleuritis in 1, neuropathy in 2, cutaneous vasculitis in 2, glomerulonephritis in 1, rheumatoid nodule in 2). There were three (27%) patients with extra-articular manifestations beyond ILD in AE group while five (13%) in non-AE group. No significant difference in the incidence was observed between AE group and non-AE group (p=0.26).

HRCT pattern

Of the 51 patients with RA-ILD, 14 (27%) exhibited the UIP pattern on HRCT. HRCT images of 8 of 40 patients (20%) in the non-AE group and 6 of 11 (55%) patients in the AE group demonstrated a UIP pattern. Interobserver agreement between both radiologists was moderate (κ statistic test, κ=0.44). The number of patients who exhibited a UIP pattern on HRCT was significantly higher in the AE group than in the non-AE group (table 1).

When HRCT pattern was identified as a UIP pattern, the positive predictive value and negative predictive value for AE occurrence was 42.9% and 86.5%, respectively.

Diagnostic accuracy of HRCT pattern for histopathological pattern

Among the 51 patients, 21 (41%) underwent SLB. Histopathological analysis revealed UIP in 12 patients, NSIP in 7, desquamative interstitial pneumonia in 1 and an unclassified interstitial pneumonia in 1. The diagnostic accuracy of HRCT pattern for histopathological pattern was evaluated in 21 histopathologically confirmed patients. The sensitivity, specificity, positive predictive value and negative predictive value of HRCT for histopathological evidence of UIP is 41.7%, 100%, 100% and 56.3%, respectively (table 2).

Table 2

Diagnostic accuracy of HRCT pattern for histopathological pattern in 21 histopathologically confirmed patients

Treatment for RA

Of the 51 patients, 37 (73%) had received treatment for RA at their final visit (non-AE group) or at AE onset (AE group). Specifically, 29 patients (57%) had received corticosteroids, 17 (33%) received immunosuppressants except for methotrexate (MTX), 10 (20%) received MTX and 17 (33%) received other drugs. There were 6 (55%) patients who had received MTX in AE group while 4 (10%) patients in non-AE group and statistically significant difference was observed among the groups (p=0.001; table 3).

Table 3

Treatment for RA at final visit† or AE onset‡

None of the patients had received MTX earlier in the AE group, whereas two in non-AE group had received MTX during observation period but discontinued prior to final visit. Thus, there were 6 (55%) patients who had experienced MTX treatment during observation period in AE group and 6 (15%) patients in non-AE group. Median cumulative MTX dose were 1952 mg (384–3872) in six patients who had experienced MTX among AE group and 802 mg (32–2496) in six patients among non-AE group. No statistically significant difference was observed between two groups in cumulative MTX dose (p=0.20).

Articular disease activities (RF and DAS28-CRP) were compared between patients who had received MTX at the last visit or at the first AE occurrence, and those who had not. Median RF (range) were 112 (16–334) and 284 IU/mL (6–2666; p=0.25), median DAS28-CRP (range) were 1.94 (1.02–3.87) and 1.97 (1.02–5.13; p=0.78), respectively. There were no statistically significant differences in articular disease activities.

AE incidence

The cumulative AE incidence in patients with RA-ILD is shown in figure 2. The 5-year AE incidence was 11% (95% CI 2% to 21%) in all patients, 33% (95% CI 7% to 60%) in patients with a UIP pattern on HRCT (UIP pattern group) and 3% (95% CI 0% to 9%) in patients with a non-UIP pattern on HRCT (non-UIP pattern group). The AE incidence was significantly higher in UIP pattern group than in non-UIP pattern group (log-rank test, p=0.018). The overall 1-year AE incidence was 2.8%. The 1-year AE incidence was 6.5% in the UIP pattern group and 1.7% in the non-UIP pattern group.

Figure 2

Cumulative acute exacerbation (AE) incidence in patients with rheumatoid arthritis-associated interstitial lung disease. Patients with usual interstitial pneumonia (UIP) pattern on high resolution CT (HRCT) had a significantly higher incidence of AE compared with non-UIP pattern on HRCT (log-rank p=0.018).

Univariate analysis (table 4) revealed that older age at ILD diagnosis (HR 1.11; 95% CI 1.02 to 1.21; p=0.01), UIP pattern on HRCT (HR 1.95; 95% CI 1.07 to 3.63; p=0.03) and MTX usage (HR 3.04; 95% CI 1.62 to 6.02; p=0.001) were associated with the development of AE.

Table 4

Risk factors for AE occurrence according to univariate cox hazard analysis

Survival

The overall survival of patients with RA-ILD and the survival according to HRCT pattern subgroup are as shown in figure 3. The 5-year survival was 90% (95% CI 81% to 98%) in all patients, 70% (95% CI 44% to 94%) in the UIP pattern group and 97% (95% CI 92% to 100%) in the non-UIP pattern group. Survival was significantly poorer in the UIP pattern group than in the non-UIP pattern group (log-rank test, p=0.04).

Figure 3

Overall survival and the survival according to high-resolution CT (HRCT) pattern subgroup. Patients with usual interstitial pneumonia (UIP) pattern on HRCT had a significantly worse survival compared with those with non-UIP pattern on HRCT (log-rank p=0.04).

The survival statistics of AE group and non-AE group are shown in figure 4. Significantly worse survival was demonstrated in the AE group compared with the non-AE group (log-rank test, p=0.001).

Figure 4

Survival of acute exacerbation (AE) group and non-AE group. Patients who developed AE during observation periods (AE group) had a significantly worse survival compared with those who did not (non-AE group; log-rank p=0.001).

Univariate analysis (table 5) revealed that older age at ILD diagnosis (HR 1.08; 95% CI 0.99 to 1.17; p=0.057) and UIP pattern on HRCT (HR 1.74; 95% CI 0.97 to 3.12; p=0.06) showed a trend towards poor outcome, while AE (HR 2.47; 95% CI 1.39 to 4.56; p=0.003) was significantly associated with poor survival.

Table 5

Prognostic factors for survival, univariate cox hazard analysis

Discussion

To our knowledge, the current study is the first to focus on the risk factors and prognosis associated with AE in patients with RA-ILD. The overall 1-year AE incidence was 2.8% among the RA-ILD patients. Univariate analysis identified that a UIP pattern on HRCT, MTX usage and older age at ILD diagnosis were associated with the development of AE. Furthermore, AE was a prognostic factor for poor survival.

The 1-year AE incidence among patients with IPF has been reported to be 5–19%,17 ,18 whereas that among patients with other ILDs was 4.2% for idiopathic NSIP, 3.3% for CVD-ILD and 5.6% for CVD-UIP.4 In the current study, the 1-year AE incidence among patients with a UIP pattern on HRCT was 6.5% and cumulative incidence was significantly higher in UIP pattern group than in non-UIP pattern group. Therefore, a UIP pattern on HRCT was associated with the development of AE. In a previous report of small number of biopsy-proven cases, AE incidence among patients with RA-UIP was reported to be 11.1%, similar to that among patients with IPF.4 The correlation between pathological and radiological UIP patterns has been established in IPF, but not in RA-ILD.19 In the current study, determination of a UIP pattern on HRCT exhibited high specificity and positive predictive values for the detection of pathological UIP pattern. Therefore, a typical UIP pattern on HRCT was assumed to be highly suggestive of pathological UIP pattern in patients with RA-ILD as well as IPF.15 ,19 Collectively, these results suggest that AE incidence in a subgroup with pathological UIP pattern or a UIP pattern on HRCT may be higher than that in other patterns.

The current study also identified that MTX usage was associated with the development of AE. MTX is widely used to treat RA, although its toxicity affects the pulmonary system. Diagnostic criteria for MTX-associated pneumonitis (MTX-pneumonitis) have been proposed previously.20 However, the specificity of the criteria has not been fully examined as yet. MTX-pneumonitis typically occurs with the acute/subacute onset early in the course of MTX treatment,20–23 in particular, mostly within the first year24 and often presents a hypersensitive pneumonitis pattern. Patients with MTX-pneumonitis generally respond to discontinuation of MTX or corticosteroid treatment and have a favourable prognosis.20–23 In the current study, 6 of 11 patients with AE had received MTX, and five of the six patients had been treated for three or more years. The remaining patient had been treated for 1 year. In all six patients, their respiratory conditions progressively deteriorated despite discontinuation of MTX, and they poorly responded to high-dose corticosteroid. Consequently, on the basis of inconsistency in terms of MTX treatment duration and clinical features, we diagnosed these patients with AE of RA-ILD.

Similar to AE in patients with IPF, the aetiology of AE in RA-ILD is unknown; however, several reports have suggested that AE may be a distinct manifestation of the primary disease process, a distinct condition associated with undiagnosed infection, or a subsequent acceleration of a fibroproliferative process caused by acute direct stress to the lung.1 In the current study, DAS28-CRP in all patients with AE showed low disease activity or remission at the time of AE occurrence, and no significant difference in articular disease activities was observed between AE group and non-AE group or regardless of MTX treatment; therefore, RA activity may not be related to AE occurrence, a finding similar to that of a previous study.4 MTX possibly accelerates the fibroproliferative process of RA-ILD. It was reported that MTX treatment was a risk factor for RA-ILD progression.12 For RA-ILD, MTX may be associated not only with newly developing drug-induced pneumonitis but also with deterioration of pre-existing ILD, including AE.

This study found that patients with a UIP pattern on HRCT had significantly poorer survival than those with a non-UIP pattern. In addition, AE was a prognostic factor for poor outcome. AE has a serious impact on the survival of patients with IPF,18 and it may also affect the prognosis in other ILDs.3 ,4 Consistent with these reports, the current study revealed that 7 (64%) of 11 patients with AE died of respiratory failure, giving rise to high in-hospital mortality, and that survival in the AE group was significantly poorer than that in the non-AE group. Therefore, AE is suggested to be a predictor of poor prognosis in patients with RA-ILD as well as patients with IPF.

Other prognostic factors in RA-ILD have been reported, including pathological UIP,25–29 UIP pattern on HRCT30 and age at RA-ILD onset.31 In the current study, pathological UIP could not be identified as prognostic factors because the number of biopsy-proven cases were small. However, we observed similar results regarding UIP pattern on HRCT and age at RA-ILD onset.

In our study, the average time between age at the onset of RA and that of RA-ILD was relatively shorter than some previous reports.7 ,8 However, Gabbay et al11 reported that nearly 60% of the patients have ILD in the recent onset RA. In another point of view, since the current authors’ institution is regional referral centre for ILD, referral bias might have increased the proportion of patients with early-stage ILD, leading above results.

This study had several limitations. First, given its retrospective study design, it is subject to several possible biases. For instance, selection and recall bias may exist. Second, because of the relatively small sample size, it may be difficult to determine the precise incidence of AE and the number of event (AE occurrence or death) was too small to perform multivariate analysis. In this analysis, the incidence of extra-articular manifestations beyond ILD and the proportion of current/former smokers in AE group were higher than those in non-AE group (27% vs 13% and 82% vs 55%, respectively), with no statistically significant differences. Such differences might be statistically significant in a larger study. Therefore, larger studies are necessary to confirm our results. Third, because the presence of a UIP pattern on HRCT had low sensitivity and negative predictive value for the detection of pathologic UIP, it was likely that some patients with a non-UIP pattern on HRCT may have had pathological UIP. HRCT criteria used in our study do not aim to detect pathological UIP but to clinically diagnose IPF and to exclude alternative diagnosis.15 The frequency of honeycombing on HRCT differs between RA-ILD and IPF, and some ILDs by other causes including RA-ILD have some of the findings inconsistent with UIP. Therefore, HRCT pattern in RA-ILD may be more frequently interpreted as non-UIP pattern than that in IPF. In our study, two of five patients with histopathological UIP in the absence of a UIP pattern on HRCT developed AE. Thus, SLB may have to be considered for these patients to make pathological diagnosis and to predict AE. Finally, RA includes a variety of manifestations. Therefore, other comorbidities may have affected the results.

In conclusion, older age at ILD diagnosis, UIP pattern on HRCT and MTX usage are associated with the development of AE in patients with RA-ILD. Furthermore, the mortality associated with AE is high and AE is prognostic factor for poor outcome. Larger prospective studies investigating AEs in RA-ILD are indicated.

Acknowledgments

We appreciate the help of Dr T Tsuchiya in the statistical advice. This study was partly supported by the Diffuse Lung Diseases Research Group from the ministry of Health, Labour and Welfare, Japan.

References

Footnotes

  • Contributors HH, YN, TS and KC designed the research; HH, DH, TF and NI contributed in acquisition of the data; HH, YN, TJ, HS, TC, MK, NE and TS interpreted and analysed the data. HH, YN, TJ and TC wrote the paper; HS, MK, DH, NE, TF, NI, TS and KC revised the paper. All authors approved the final version of the paper.

  • Funding This study was partly supported by a grant to the Diffuse Lung Diseases Research Group from the minisry of Health, Labour and Welfare, Japan.

  • Competing interests None.

  • Ethics approval Institutional Review Board of the Hamamatsu University School of Medicine approved this study (approval number 24–191) and waived patient approval or informed consent because the study involved a retrospective review of patient records and images.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.