Abstract
Introduction
Interstitial lung disease (ILD) is a heterogeneous group of chronic respiratory diseases characterized by inflammation and progressive pulmonary fibrosis.1,2 ILDs are classified according to their etiology and encompass idiopathic interstitial pneumonias, connective tissue disease (CTD)-related ILD, hypersensitivity pneumonitis, and other forms, such as drug-induced and postinfectious ILD. Among these, idiopathic pulmonary fibrosis (IPF) is the most prevalent, accounting for over 30% of ILD cases. Although each type of ILD has distinct pathophysiological features, clinical presentations, and prognoses, all forms can ultimately lead to irreversible pulmonary fibrosis. 3
In terms of symptoms and presentations, ILD is typically characterized by progressive dyspnea and limited physical activity,1,2 which have also been suggested to be associated with high levels of premature mortality. In ILD, lung function is the most commonly used predictor of mortality and the prognosis after diagnosis. Moreover, measures of exercise capacity and physical activity have also been proposed as predictors of mortality. 4 For example, previous studies demonstrated that reduced physical activity, 5 lower exercise capacity indicated by the six-minute walk test (6MWT),6,7 and decreased peak oxygen uptake (VO2 peak)8,9 are associated with a higher mortality risk in individuals with IPF or other ILDs. 10 However, the sample size was relatively small, and the follow-up duration was insufficient. In addition, these studies did not consider changes in exercise capacity and primarily evaluated patients with IPF. Furthermore, covariate adjustments were insufficient; while some studies adjusted for age, sex, baseline lung function, or smoking history,6,7 others did not account for clinical variables that could affect the outcomes.
Pulmonary rehabilitation, including exercise training, in individuals with ILD is strongly recommended since it is related to improved exercise performance and quality of life.11,12 Beyond these beneficial effects, exercise may be associated with a reduction in overall mortality in individuals with ILD. However, limited information is available regarding whether maintaining exercise is associated with improved survival in individuals with ILD. Therefore, we aimed to evaluate the impact of exercise maintenance on mortality in individuals with ILD using a longitudinal, large population-based database.
Methods
Data source
This retrospective cohort study utilized the Korean National Health Insurance claims-based database, administered by the National Health Insurance Service (NHIS). The NHIS functions as a government-managed universal insurance provider, extending coverage to approximately 97% of the Korean population, encompassing around 50 million individuals.13,14 The NHIS provides health screening exams and reports on sociodemographic data, a self-questionnaire survey, basic clinical laboratory findings, inpatient and outpatient usage, prescription records, and recorded diagnoses based on International Classification of Diseases 10th Revision (ICD-10) codes. Mortality information was also acquired from the NHIS database, supplied by Statistics Korea under the auspices of the Ministry of Strategy and Finance of Korea.15,16 This extensive NHIS database has been widely employed in diverse epidemiological studies aimed at identifying risk factors and prognoses in pulmonary diseases.17,18
Study population
Among 8050 individuals who were newly diagnosed with ILD between January 2010 and December 2016, we included 4042 who underwent two consecutive health examinations within 2 years before and after the ILD diagnosis to measure changes in their exercise habits between the two time points. Among them, those with a follow-up duration of less than 1 year were excluded to ensure enough prognostic observation. In addition, we excluded individuals aged <20 years and those with missing variables on health screening examinations (

Flow chart of the study population.
ILD was defined as claims with any of the ICD-10 diagnosis codes J84. ILD subtypes were categorized as IPF (J84.1 or J84.18) or non-IPF ILD (CTD-related ILD, J84 codes with CTD codes (M05, M06, M30, M31, M32, M33, M34, M35, M36, M45), hypersensitivity pneumonitis (D67), sarcoidosis (D86), and other ILDs that do not classify into the previous subtypes).17,19,20
Study exposure: Exercise
The exposure in our study was the change in exercise habits. The intensity and frequency of regular physical activity were assessed using a self-reported questionnaire during two consecutive health screening examinations. This questionnaire has been used in several high-quality studies and includes the frequency of light, moderate, or vigorous weekly physical activity of varying intensity.21–23
Regular exercise was defined as a moderate-intensity exercise for more than 5 days per week or vigorous-intensity exercise for more than 3 days per week. 24 To evaluate the impact of changes in exercise habits on the risk of mortality, the study population was subdivided into three groups according to changes in exercise performance: non-exercisers were individuals who never performed exercise in either examination, while exercise non-maintainers were those who were non-exercisers in the first examination but exercisers in the second examination or exercisers in the first examination but non-exercisers in the second examination. Exercise maintainers were individuals who exercised in both examinations.
Study outcome: Mortality
The primary outcome was a comparison of the incidence and risk of all-cause mortality. Death-related information was available through December 31, 2022. Thus, for the mortality analysis, we tracked individuals until the date of death or December 31, 2022, whichever came first. We had a 1-year lagging period to evaluate the impact of exercise on the long-term prognosis of ILD.
Covariates
Data for baseline characteristics were collected from the study dataset. Body mass index (BMI) was calculated as body weight divided by the square of height (kg/m2) and was classified into four groups, as recommended for Asians: normal (18.5–22.9 kg/m2), low (<18.5 kg/m2), overweight (23.0–24.9 kg/m2), and obese (⩾25 kg/m2).25,26 Smoking status (never-, ex-, or current smoker) and alcohol status were determined based on a self-reported questionnaire. Income status was divided into the highest 20% (high), the lowest 20% (low), and the rest (middle); individuals supported by the medical aid program were classified as the low-income group. The residential area was classified into metropolitan cities, middle- and small-sized cities, and rural areas. Comorbidities were defined as one or more claims under the following ICD-10 diagnosis codes as major diagnoses within 1 year of the index date27,28: hypertension (I10–I15), diabetes mellitus (E10–E14), dyslipidemia (E78), chronic kidney disease (N18), and airway disease (chronic obstructive pulmonary diseases (J43, J44, except for J430) and asthma (J45, J46)).
Statistical analysis
Descriptive statistics are presented as numbers (percentages) for categorical variables and as the mean (standard deviation) for continuous variables. Group comparisons were conducted using the chi-square (χ2) test for categorical variables and one-way analysis of variance (ANOVA) for continuous variables. The incidence rates of severe exacerbation and death were calculated by dividing the number of incident events by the total follow-up period (per 10,000 person-years (PY)). A cumulative incidence plot was used to compare the incidences of death, and a log-rank test was utilized to evaluate significant differences between the three groups. We used Cox proportional hazards regression analyses to evaluate the risk of death in the non-exercisers versus exercise non-maintainers versus exercise maintainers. The multivariable analysis was adjusted for age, sex, BMI, smoking status, alcohol status, economic status, residential area, and comorbidities (hypertension, diabetes mellitus, dyslipidemia, chronic kidney disease, and airway disease). We also conducted subgroup analyses based on sex, age, BMI, smoking status, economic status, residential area, and alcohol status. A two-sided
The reporting of this study conforms to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement. 29
Results
Baseline characteristics
Of 3850 individuals with ILD, 32% were exercise maintainers (
Baseline characteristics of individuals.
Data are presented as number (percentages) or means (standard deviations).
Income status was divided into the highest 20% (high), the lowest 20% (low), and the rest (middle); individuals supported by the medical aid program were classified as the low-income group.
COPD, chronic obstructive pulmonary disease.
Impact of exercise habits on mortality
During a median follow-up of 7.2 (interquartile range, 5.6–9.2) years, the incidence rate of death was lower in exercisers (341.28 per 10,000 PY) compared to non-exercisers (401.81 per 10,000 PY), yielding an adjusted hazard ratio (aHR) of 0.82 (95% confidence interval (CI): 0.72–0.94; Table 2). In the cumulative incidence plot for death, exercisers have a better survival rate than non-exercisers (Figure 2(a), log-rank
Impact of exercise on risk of mortality in ILD.
Adjusted for age, sex, BMI, smoking status, alcohol status, economic status, residential area, and comorbidities (hypertension, diabetes mellitus, dyslipidemia, chronic kidney disease, and airway disease).
Cox regression analyses could not be performed because the proportional hazards assumption was not satisfied.
BMI, body mass index; CI, confidence interval; HR, hazard ratio; ILD, interstitial lung disease; IPF, idiopathic pulmonary fibrosis; PY, person-years.

Cumulative incidence of death in the ILD cohort, stratified by (a) exercise status and (b) exercise maintenance status.
The incidence rate of death was inversely correlated with the level of exercise adherence (319.56 per 10,000 PY for exercise maintainers vs 360.44 per 10,000 PY for exercise non-maintainers vs 401.81 per 10,000 PY for non-exercisers). The risk of death was lowest among exercise maintainers (aHR: 0.78, 95% CI: 0.66–0.92), followed by exercise non-maintainers (aHR: 0.85, 95% CI: 0.73–0.99), compared to non-exercisers (Table 2). The cumulative incidence plot for death showed similar results (Figure 2(b), log-rank
When analyzed by ILD subtypes, among individuals with non-IPF ILD, exercisers showed a lower rate of death compared to non-exercisers (284.35 vs 360.36 per 10,000 PY,

Cumulative incidence of death in individuals with (a) non-IPF ILD and (b) IPF, stratified by exercise status.
Subgroup analyses
Subgroup analyses regarding the risk of mortality according to exercise among individuals with ILD are shown in Table 3. Age, sex, smoking status, economic status, residential area, and alcohol status exhibited no interaction in the association of exercise with the development of deaths (all
Subgroup analysis of mortality risk in ILD.
Adjusted for age, sex, BMI, smoking status, alcohol status, economic status, residential area, and comorbidities (hypertension, diabetes mellitus, dyslipidemia, chronic kidney disease, and airway disease).
BMI, body mass index, CI, confidence interval; HR, hazard ratio; ILD, interstitial lung disease; PY, person-years.
Subgroup analyses regarding the risk of mortality among individuals with ILD according to exercise maintenance are shown in Supplemental Table 1.
Discussion
Utilizing a nationwide representative dataset, we evaluated the risk of death in nearly 4000 individuals with ILD according to their adherence to exercise. During a median follow-up of 7.2 years, exercisers had a significantly lower risk of death compared to non-exercisers by about 20%. The risk of mortality further decreased as the adherence to exercise maintenance increased in a dose-dependent manner, showing the lowest risk in exercise maintainers.
In individuals with ILD, exercise capacity is reduced, and the mechanisms behind this reduction are multifactorial. The basic physiology involved is known to be as follows. Impaired gas exchange occurs due to the destruction of the pulmonary capillary bed, leading to ventilation-perfusion mismatch and oxygen diffusion limitations. 30 Circulatory limitations arise from pulmonary capillary destruction and pulmonary vasoconstriction, which can result in pulmonary hypertension and cardiac dysfunction in some patients. 31 Ventilatory limitations to exercise may also occur, though they are not considered a major contributor in most patients. 32 Peripheral muscle dysfunction, often due to physical deconditioning, may limit exercise capacity. 33
Previous studies have examined measures of exercise capacity and physical activity and their relationship to prognosis in ILD. One retrospective, international cohort study including 710 fibrotic patients with ILD found that those who improved physical performance due to pulmonary rehabilitation had better survival compared with those who did not improve. 34 In addition, in a single-center study of 169 patients with ILD without pulmonary hypertension, a 6MWT distance of less than 300 m and a 6MWT final oxygen saturation of less than 85% were associated with an increased risk of mortality. 35 However, except for the aforementioned studies, many studies on this issue focused exclusively on IPF. In one review article on the relationship between physical activity, exercise capacity, and mortality risk in ILD, including IPF, both physical activity and exercise capacity were associated with mortality risk. 36 However, the evidence was weak due to the small sample sizes of the included studies, and the follow-up periods were not long enough to estimate the prognosis of ILD. Overcoming these limitations of previous studies, we comprehensively included individuals with various types of ILDs as a whole, with a large sample size (a total of 3850 ILDs in three exercise groups). The follow-up duration was also longer than in previous studies, with a median of 7.2 years.
With these advantages, our analysis adds robust evidence for the relationship between exercise habits and ILD prognosis; beyond merely showing exercise as a predictor of mortality in individuals with ILD, our study findings emphasize that adherence to exercise could provide more benefits in individuals with ILD over a long-term period of the natural disease course. In our study, the risk reduction for death increased as adherence to exercise increased (compared to non-exercisers, 22% reduction in exercise maintainers and 15% in exercise non-maintainers), suggesting that the importance of monitoring exercise capacity or providing rehabilitation should be considered in the long-term management of ILD.
Interestingly, in subgroup analyses by ILD subtypes, the benefit of exercise in the IPF group showed for up to 7 years after follow-up; however, it disappeared thereafter. The reasons for this observation are not fully explainable. However, we carefully suggest that factors affecting survival other than exercise capacity may be more complex in patients with IPF than in those with non-IPF, considering the worst prognosis of IPF among various ILDs. For example, the risk of lung cancer and acute exacerbation in IPF increases as the duration of IPF increases, which may affect mortality regardless of exercise capacity.37,38 Future studies evaluating exercise and mortality in IPF with a long-term follow-up duration should consider this.
Our subgroup analysis further revealed a significant interaction between BMI and economic status in the association between exercise and mortality. Individuals with a lower BMI and higher economic status had a reduced risk of mortality compared to their counterparts. Previous studies have demonstrated that adipose tissue releases pro-inflammatory cytokines, including tumor necrosis factor-alpha, interleukin-6, and leptin, as well as adipokines that contribute to the progression of pulmonary fibrosis.39,40 This may explain why individuals with a low or normal BMI in our study exhibited a lower risk of mortality compared to those who were overweight or obese. Economic status can influence access to exercise resources, health literacy, comorbidity burden, and healthcare utilization,41,42 all of which may impact the interaction between exercise and mortality. Our additional subgroup analyses indicated that individuals in a higher economic status had a lower risk of mortality compared to those in a lower economic status.
Our study has several limitations. First, individuals who underwent health screening examinations were included, which might be associated with healthy user bias. In addition, only subjects with two consecutive health examination results were included in our study, leading to the exclusion of half of the total population and potentially introducing bias. Second, the diagnosis of ILD and other comorbidities relied on ICD-10 codes, potentially leading to overestimations or underestimations of diagnoses. However, the ICD-10 codes used in this study have been validated in many previous studies of ILD.19,20 Third, there was no data on lung function, imaging, treatment, ILD control status, and time since ILD diagnosis. Future studies are needed to determine if the combination of the data types mentioned above can be used to better predict the prognosis of ILD. Fourth, generalizing our findings to other countries or ethnicities may be challenging since our research utilized a dataset from the Korean population. An additional limitation is the absence of a sample size calculation or justification for the selected cohort in the study’s methodology.
Since factors other than exercise may influence clinical outcomes in individuals with ILD, and given the limitations of this study, including the lack of detailed information on ILD severity and treatment, the results should be interpreted with caution.
Conclusion
In conclusion, among individuals with ILD, exercise was associated with reduced risk of mortality, which was particularly beneficial among exercise maintainers. The effect was more pronounced in the non-IPF group than in the IPF group, and individuals with lower BMI and higher economic status had a lower risk of death compared to other groups.
Supplemental Material
sj-docx-1-tar-10.1177_17534666251362380 – Supplemental material for Impact of exercise maintenance on mortality in interstitial lung disease: a population-based retrospective cohort study
Supplemental material, sj-docx-1-tar-10.1177_17534666251362380 for Impact of exercise maintenance on mortality in interstitial lung disease: a population-based retrospective cohort study by Bo-Guen Kim, Min Gu Kang, Sung Jun Chung, Hyung Koo Kang, Jong Seung Kim and Hyun Lee in Therapeutic Advances in Respiratory Disease
Supplemental Material
sj-docx-2-tar-10.1177_17534666251362380 – Supplemental material for Impact of exercise maintenance on mortality in interstitial lung disease: a population-based retrospective cohort study
Supplemental material, sj-docx-2-tar-10.1177_17534666251362380 for Impact of exercise maintenance on mortality in interstitial lung disease: a population-based retrospective cohort study by Bo-Guen Kim, Min Gu Kang, Sung Jun Chung, Hyung Koo Kang, Jong Seung Kim and Hyun Lee in Therapeutic Advances in Respiratory Disease
Footnotes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
