Patient characteristics
We identified 249 patients with uILD diagnosed between 2003 and 2018; 143 patients were included in Aarhus from 2003-2009 and 2012-2018, and 106 patients in Heidelberg from 2012-2018. The majority of the patients were males with a smoking history (Table 1). The median follow-up time was 2.0 years; 57 patients (23%) died during follow-up. At baseline, FVC% predicted was mildly reduced whereas DLCO% predicted was moderately reduced. More patients from Aarhus had a smoking history compared to patients from Heidelberg (p=0.004). In the Aarhus cohort, DLCO%, FVC%, and 6MWD was higher, but only the difference in DLCO% reached statistical significance (p<0.001, p=0.1 and p=0.08, respectively). Baseline characteristics are shown in Table 1.
Table 1. Baseline characteristics of the uILD patients.
|
All patients,
n = 249
|
Aarhus cohort,
n = 143
|
Heidelberg cohort,
n = 106
|
Age, years (IQR)
|
70.0 (60.0-75.0)
|
69.0 (60.0-75.0)
|
70.5 (58.0-77.0)
|
Male gender, %
|
60.2
|
58.7
|
62.3
|
Never smokers, %
Current smokers, %
Former smokers, %
Missing, %
|
34.5
14.9
49.0
1.6
|
27.3
18.9
52.4
1.4
|
44.3
9.4
44.3
1.9
|
Pack years (IQR)
|
25.0 (10.0-40.0)
|
30.0 (15.0-40.0)
|
20.0 (10.0-40.0)
|
Charlson comorbidity index (IQR)
|
0 (0-1)
|
0 (0-1)
|
0 (0-1)
|
FVC, % pred (SD)
|
75.9 (24.6)
|
78.1 (24.9)
|
73.0 (23.9)
|
DLCO, % pred (SD)
|
49.9 (20.3)
|
54.1 (20.8)
|
44.6 (18.3)
|
6MWD, m (SD)
|
357.4 (136.2)
|
372.3 (142.5)
|
338.1 (125.6)
|
Follow up, years (IQR)
|
2.0 (0.8-3.3)
|
1.8 (0.9-2.9)
|
2.5 (0.7-3.6)
|
1-year survival (95% CI)
|
0.89 (0.84-0.92)
|
0.88 (0.81-0.92)
|
0.91 (0.82-0.95)
|
3-year survival (95% CI)
|
0.73 (0.66-0.79)
|
0.67 (0.56-0.76)
|
0.81 (0.70-0.88)
|
Data are presented as frequencies, mean with standard deviation (SD), median with interquartile range (IQR) or survival with 95% confidence intervals (CI). FVC: Forced vital capacity, DLCO: diffusion capacity of the lung for carbon monoxide, 6MWD: distance walked during the 6-minute walk test.
Number of comorbidities
Eighty-five percent of the patients had one or more comorbidities (Figure 1a). The median number of comorbidities was two (IQR 1-3), and the median Charlson comorbidity index was zero (IQR 0-1). The frequency of the registered comorbidities ranged from 1-39%. Arterial hypertension (39%), emphysema (30%), diabetes (19%), gastro-esophageal reflux (18%) and coronary artery disease (17%) were most common. The baseline prevalence of all comorbidities is presented in Figure 1b.
Data are presented as a percentage of all patients. A: Total number of comorbidities per patient. B: Spectrum of comorbidities in the uILD cohort. Multiple comorbidities could be reported.
Survival analysis
The only comorbidity associated with excess mortality was dyslipidemia in both univariate and multivariate analyses (Table 2). None of the five pre-specified comorbidities of special interest (emphysema, ischemic heart disease, pulmonary hypertension, diabetes, and gastro-esophageal reflux disease) were associated with excess mortality. No survival difference was seen in patients with zero comorbidities vs. 1-3 or zero vs. 4 or more comorbidities, nor in analyses of patients stratified by the median number of comorbidities or median Charlson comorbidity index (Table 3).
Table 2. Survival analyses for specific comorbidities
Comorbidities
|
Univariate analysis (95% CI)
|
Multivariate analysis (95% CI)
|
Diabetes (n=47)
|
1.20 (0.63-2.29)
|
0.93 (0.48-1.83)
|
Pulmonary hypertension (n=20)
|
1.66 (0.75-3.68)
|
1.52 (0.67-3.43)
|
Lung cancer (n=2)
|
3.19 (0.44-23.40)
|
2.32 (0.31-17.32)
|
Emphysema (n=74)
|
1.04 (0.56-1.91)
|
1.24 (0.65-2.37)
|
Gastro-esophageal reflux disease (n=46)
|
0.92 (0.46-1.83)
|
0.95 (0.46-1.96)
|
Arterial hypertension (n=98)
|
1.25 (0.74-2.11)
|
1.02 (0.58-1.79)
|
Depression (n=19)
|
1.75 (0.79-3.86)
|
1.75 (0.74-4.13)
|
Valve disease (n=12)
|
1.18 (0.37-3.80)
|
0.85 (0.26-2.77)
|
Atrial fibrillation (n=21)
|
1.66 (0.75-3.69)
|
0.90 (0.38-2.16)
|
Dyslipidemia (n=31)
|
2.16 (1.18-3.96)
|
2.25 (1.19-4.24)
|
Coronary artery disease (n=42)
|
1.35 (0.71-2.56)
|
1.06 (0.53-2.10)
|
Cancer (non-lung) (n=16)
|
1.75 (0.69-4.40)
|
1.16 (0.44-3.06)
|
Osteoporosis (n=21)
|
1.14 (0.49-2.65)
|
1.20 (0.50-2.87)
|
Cerebrovascular disease (n=17)
|
1.85 (0.79-4.32)
|
1.64 (0.68-3.97)
|
Venous thromboembolism (n=17)
|
0.80 (0.25-2.58)
|
0.48 (0.12-2.01)
|
Thyroid disease (n=16)
|
0.64 (0.20-2.05)
|
0.74 (0.17-3.20)
|
Sleep apnea (n=12)
|
*
|
*
|
Chronic heart failure (n=9)
|
1.97 (0.84-4.61)
|
1.41 (0.58-3.42)
|
Chronic renal failure (n=9)
|
0.83 (0.20-3.41)
|
0.58 (0.14-2.42)
|
Liver failure (n=2)
|
*
|
*
|
Data are presented as hazard ratios with 95% confidence intervals (CI). Hazard ratios >1 indicate an association with increased mortality. Multivariate analyses are adjusted for gender, age, and FVC% predicted. *: Number of deaths too low for analysis
Table 3. Survival analyses for number of comorbidities and Charlson comorbidity index
Parameter
|
Univariate analysis (95% CI)
|
Multivariate analysis (95% CI)
|
1-year survival (95% CI)
|
3-year survival (95% CI)
|
Number of comorbidities
|
0
|
Ref.
|
Ref.
|
0.91 (0.74-0.97)
|
0.80 (0.56-0.92)
|
1-3
|
1.07 (0.47-2.41)
|
0.78 (0.34-1.81)
|
0.90 (0.83-0.94)
|
0.74 (0.65-0.82)
|
4 or more
|
1.40 (0.56-3.50)
|
0.90 (0.35-2.31)
|
0.86 (0.71-0.93)
|
0.64 (0.46-0.78)
|
≤ median (0-2)
|
Ref.
|
Ref.
|
0.89 (0.83-0.93)
|
0.79 (0.70-0.86)
|
> median (3 or more)
|
1.46 (0.86-2.49)
|
1.12 (0.64-1.95)
|
0.88 (0.78-0.94)
|
0.63 (0.49-0.74)
|
Charlson comorbidity index
|
≤ median (0)
|
Ref.
|
Ref.
|
0.88 (0.79-0.93)
|
0.77 (0.64-0.85)
|
> median (1 or more)
|
1.17 (0.68-2.03)
|
0.75 (0.42-1.35)
|
0.90 (0.83-0.94)
|
0.71 (0.61-0.79)
|
Data are presented as hazard ratios or percent survivors with 95% confidence intervals (CI). Hazard ratios >1 indicate an association with increased mortality. Multivariate analyses are adjusted for gender, age, and FVC% predicted. Ref.: Reference group
Clusters of comorbidities
Three clusters with different comorbidity profiles were identified (Figure 2, Table 4). Patients in cluster 1 had significantly fewer comorbidities than the entire cohort. Cluster 2 was dominated by patients with a larger total number of comorbidities, primarily cardiovascular and associated diseases, a higher body mass index (BMI) and more severely impaired pulmonary function and exercise capacity based on six-minute walk test Emphysema, cancer and depression were more prevalent in cluster 3, and these patients had a lower BMI and better exercise capacity. No difference in mortality was observed between the clusters (Figure 3, Table 4). A significant annual decline in FVC and DLCO % predicted was observed in cluster 1 and 2, but not in cluster 3 (Table 5).
Cluster borders are indicated by the black lines. Each patient is placed in the same area on all maps. Red colors indicate a high frequency of the specific comorbidity, while blue colors indicate absence of the comorbidity. C1: Cluster 1; C2: Cluster 2; C3: Cluster 3; GERD: Gastro-esophageal reflux disease; CAD: Coronary artery disease; VTE: Venous thromboembolism.
Table 4. Clinical characteristics and prevalence of comorbidities in the three comorbidity clusters.
|
Cluster1:
(Few comorbidities)
n = 77
|
Cluster 2:
(Cardiovascular)
n = 85
|
Cluster 3:
(Emphysema)
n = 87
|
Age, years (IQR)
|
68.0 (57.0-75.0), p=0.22
|
69.0 (58.0-75.0), p=0.99
|
71.0 (61.0-77.0), p=0.23
|
Male, %
|
59.7, p=0.91
|
63.5, p=0.45
|
57.5, p=0.51
|
Never smokers, %
Current smokers, %
Former smokers, %
Missing
|
36.4, p=0.78
13.0, p=0.53
50.6, p=0.86
0.0
|
37.8, p=0.53
12.2, p=0.37
50.0, p=0.96
3.5
|
31.4, p=0.37
19.8, p=0.13
48.8, p=0.83
1.2
|
Pack years (IQR)
|
20.0 (15.0-30.0), p=0.18
|
25.5 (13.0-40.0), p=0.47
|
35.0 (10.0-43.0), p=0.56
|
Body mass index (SD)
|
27.2 (5.3), p=0.14
|
29.8 (4.9), p=0.001†
|
26.8 (4.8), p=0.03*
|
FVC, % predicted (SD)
|
77.4 (21.3), p=0.52
|
71.1 (26.4), p=0.03*
|
79.3 (24.8), p=0.11
|
DLCO, % predicted (SD)
|
53.8 (20.5), p=0.06
|
44.5 (19.7), p=0.006*
|
51.6 (19.8), p=0.36
|
6MWD, m (SD)
|
377.7 (123.6), p=0.19
|
304.2 (139.5), p<0.0001*
|
394.6 (127.2), p=0.005†
|
Total number of comorbidities, n (SD)
|
0.9 (1.1), p<0.0001*
|
3.4 (2.0), p<0.0001†
|
2.0 (1.1), p=0.41
|
Emphysema, %
|
5.7, p<0.0001*
|
25.7, p=0.01*
|
70.3, p<0.0001†
|
Pulmonary hypertension, %
|
3.9, p=0.11
|
14.1, p=0.01†
|
5.7, p=0.33
|
Venous thromboembolism, %
|
3.9, p=0.22
|
11.8, p=0.03†
|
4.6, p=0.31
|
Sleep apnea, %
|
2.6, p=0.28
|
8.2, p=0.07
|
3.4, p=0.46
|
Coronary artery disease, %
|
22.1, p=0.14
|
28.2, p=0.0005†
|
1.1, p<0.0001*
|
Dyslipidemia, %
|
6.5, p=0.06
|
20.0, p=0.009†
|
10.3, p=0.46
|
Arterial hypertension, %
|
0.0, p<0.0001*
|
63.5, p<0.0001†
|
50.6, p=0.008†
|
Atrial fibrillation, %
|
10.4, p=0.46
|
7.1, p=0.58
|
8.0, p=0.87
|
Chronic heart failure, %
|
7.8, p=0.56
|
8.2, p=0.40
|
3.4, p=0.16
|
Heart valve disease, %
|
3.9, p=0.65
|
8.2, p=0.07
|
2.3, p=0.17
|
Cerebrovascular disease, %
|
7.8, p=0.69
|
3.5, p=0.14
|
9.2, p=0.28
|
Chronic renal failure, %
|
1.3, p=0.19
|
8.2, p=0.005†
|
1.1, p=0.13
|
Diabetes, %
|
0.0, p<0.0001*
|
55.3, p<0.0001†
|
0.0, p<0.0001*
|
Osteoporosis, %
|
5.2, p=0.22
|
8.2, p=0.94
|
11.5, p=0.20
|
Gastro-esophageal reflux disease, %
|
0.0, p<0.0001*
|
51.8, p<0.0001†
|
2.3, p<0.0001*
|
Thyroid disease, %
|
7.8, p=0.56
|
9.4, p=0.17
|
2.3, p=0.05
|
Liver failure, %
|
1.3, p=0.56
|
0.0, p=0.31
|
1.1, p=0.66
|
Lung cancer, %
|
0.0, p=0.34
|
2.4, p=0.049†
|
0.0, p=0.30
|
Cancer, %
|
1.3, p=0.03*
|
4.7, p=0.43
|
12.6, p=0.003†
|
Depression, %
|
1.3, p=0.01*
|
7.1, p=0.81
|
13.8, p=0.007†
|
Data are presented as means with standard deviations (SD) or interquartile range (IQR) for continuous variables and frequencies for categorical variables. Significance levels for the comorbidities was based on comparison between the result in one cluster and the rest of the cohort (the two other clusters combined) using the t-test or Wilcoxon Mann-Whitney U test.
*: Significantly lower or less frequent in this cluster compared with the rest of the cohort (the two other clusters combined). †: Significantly higher or more frequent in this cluster compared with the rest of the cohort (the two other clusters combined). FVC: Forced vital capacity, DLCO: diffusion capacity of the lung for carbon monoxide, 6MWD: distance walked during the 6-minute walk test.
Table 5. Survival analyses and changes in pulmonary function during follow-up. Patients were followed for a median time of two years.
Parameter
|
Cluster 1
n = 77
|
Cluster 2
n = 85
|
Cluster 3
n = 87
|
Survival analyses
|
Number of deaths (%)
|
17 (22%)
|
22 (26%)
|
18 (21%)
|
Univariate analysis (95% CI)
|
Ref.
|
1.07 (0.56 to 2.03)
|
0.93 (0.48 to 1.81)
|
Multivariate analysis (95% CI)
|
Ref.
|
0.82 (0.42 to 1.61)
|
0.84 (0.42 to 1.67)
|
Pulmonary function
|
ΔFVC, % predicted
|
-2.59 (-3.44 to -1.74)
|
-2.17 (-3.06 to -1.27)
|
-0.54 (-1.54 to 0.47)
|
ΔDLCO, % predicted
|
-1.62 (-2.60 to -0.64)
|
-1.37 (-2.36 to -0.37)
|
-0.95 (-2.05 to 0.16)
|
Data are presented as frequencies, hazard ratios (Cox regression analyses), or change pr. Year (linear mixed effects models) with 95% confidence intervals. Multivariate analyses are adjusted for gender, age, and FVC% predicted. Δ: Change per year. FVC: Forced vital capacity, DLCO: diffusion capacity of the lung for carbon monoxide. Ref.: Reference group