Patient characteristics
A total of 150 patients with sarcoidosis were screened, of whom 87 were excluded because they lacked assessable serum samples at diagnosis (n = 83) or had insufficient available data (n = 4). As a result, 63 patients with sarcoidosis were included in the study (Fig. 1). The patient characteristics are shown in Table 1. The median age was 58 years (interquartile range [IQR], 45–68 years) and 38 (60.3%) patients were female. The most frequently affected organs were the lungs (98.4%), followed by the eyes (52.4%), skin (15.9%), and heart (6.3%). Forty-five (71.4%) patients had two or more affected organs. On chest X-rays, 21 (33.3%), 40 (63.5%), and 36 (57.1%) patients had bilateral hilar lymphadenopathy (BHL), lung parenchymal involvement, or both, respectively. Among 61 patients who underwent pulmonary function tests, 55 (90.2%) had normal pulmonary function (forced vital capacity, % predicted [%FVC] ≥ 80.0). The median observation time was 5.5 years (IQR, 4.0–10.3 years). The 38 healthy controls had a median age of 34 years (IQR, 27–43 years), a median BMI of 22.2 kg/m2 (IQR, 19.6–23.7 kg/m2), a female proportion of 39.5%, and no smoking history (Supplementary Table 3). The healthy controls were significantly younger and had a lower proportion of smoking history than the patients with sarcoidosis (both p < 0.001)
Table 1
| All patients n = 63 | Single organ involvement, n = 18 | Multiple organ involvement n = 45 | P-value* |
Age, years | 58 (45–68) | 59 (44–75) | 57 (45–65) | 0.330 |
Sex, female | 38 (60.3) | 9 (50.0) | 29 (64.4) | 0.394 |
Body mass index, kg/m2 | 21.4 (20.1–23.8) | 21.9 (20.4–23.1) | 21.4 (19.9–23.8) | 0.897 |
Smoking, ever-smoker | 31 (49.2) | 9 (50.0) | 22 (48.9) | 1.000 |
Number of affected organs, 1/2/3/4/5 | 18/30/13/1/1 | 18/0/0/0/0 | 0/30/13/1/1 | < 0.001 |
Affected organs | | | | |
Lungsa | 62 (98.4) | 17 (94.4) | 45 (100.0) | 0.286 |
Eyes | 33 (52.4) | 1 (5.6) | 32 (71.1) | < 0.001 |
Skin | 10 (15.9) | 0 | 10 (22.2) | 0.051 |
Heart | 4 (6.3) | 0 | 4 (8.9) | 0.317 |
Othersb | 13 (20.6) | 0 | 13 (28.9) | 0.013 |
Serum ACE, IU/L | 19.3 (14.5–25.3) | 16.6 (14.2–21.5) | 21.4 (15.1–25.4) | 0.122 |
Radiographic stage, 0/Ⅰ/Ⅱ/Ⅲ/Ⅳ | 1/21/36/4/1 | 1/6/9/1/1 | 0/15/27/3/0 | 0.802 |
Pulmonary function tests (n = 61) | | | | |
FVC, L | 2.76 (2.34–3.40) | 3.00 (2.43–3.49) | 2.75 (2.33–3.30) | 0.389 |
FVC, % predicted | 94.9 (88.1-100.7) | 96.4 (90.8-104.9) | 94.1 (87.5-100.3) | 0.187 |
FEV1, L | 2.10 (1.82–2.79) | 2.55 (1.92–3.36) | 2.10 (1.81–2.66) | 0.385 |
FEV1, % predicted | 90.3 (77.4-100.2) | 89.5 (78.2–98.8) | 91.1 (76.1-100.5) | 0.915 |
FEV1/FVC, % | 77.9 (74.3–83.8) | 78.2 (74.8–82.8) | 77.9 (72.9–85.1) | 0.812 |
Bronchoalveolar lavage (n = 62) | | | | |
Total cells, ×105/ml | 0.96 (0.61–1.48) | 1.00 (0.51–1.51) | 0.93 (0.68–1.47) | 0.561 |
Lymphocytes, % | 11.9 (8.2–22.9) | 19.0 (5.9–26.5) | 11.5 (8.3–19.9) | 0.852 |
CD4/CD8 ratio | 5.02 (3.30–7.75) | 4.58 (3.30–6.49) | 5.16 (3.34–7.91) | 0.841 |
Data are presented as median (interquartile range) or number (%). |
ACE, angiotensin-converting enzyme; FVC, forced vital capacity; FEV1, forced expiratory volume in 1 second. |
Radiographic stages 0, Ⅰ, Ⅱ, Ⅲ, and Ⅳ represent normal appearance, bilateral hilar lymphadenopathy (BHL) alone, BHL and lung parenchymal involvement, lung parenchymal involvement without BHL, and pulmonary fibrosis, respectively. |
aLungs include pulmonary hilar and mediastinal lymphadenopathy. |
bOthers include muscle, liver, thyroid gland, spleen, nerve, and extramediastinal lymphadenopathy. |
*Comparison between patients with single and multiple organ involvement. |
Clinical Outcomes Of Patients With Sarcoidosis
Eight (12.7%) patients received treatment for sarcoidosis from the time of diagnosis, comprising 7 with corticosteroid alone and 1 with a combination of corticosteroid plus azathioprine. After the treatment, 2 had remission, 2 had stable disease, and 4 had disease progression (Fig. 1). The remaining 55 patients were initially followed without treatment, of whom 20 (36.4%) had spontaneous remission, 25 (45.5%) had stable disease, 6 (10.9%) had disease progression, and 4 (7.3%) were lost to follow-up during the observation period (Fig. 1). Among 6 patients with disease progression after the initial follow-up, 2 (33.3%), 3 (50.0%), and 1 (16.7%) had worsening of lung involvement, development of non-pulmonary lesions, and both, respectively, with 5 receiving corticosteroid alone and 1 receiving no treatment. Among the 5 patients who received corticosteroid, 4 had remission, but 1 had disease progression. One patient who did not receive treatment after disease progression did not deteriorate further or require any treatment.
Associations Of Lcfa Levels With Demographic Characteristics
Among the 63 patients with sarcoidosis, there were weak or moderate correlations between age and n-3 PUFAs (r = 0.35) or n-3/n-6 ratio (r = 0.49), and between BMI and MUFAs (r = 0.25) (Supplementary Fig. 1). The LCFA levels were not correlated with sex (Supplementary Table 4). Among the LCFAs, there were strong correlations between n-3 PUFAs and n-6 PUFAs (r = 0.80), and between SFAs and MUFAs (r = 0.92). The correlations among single LCFAs are presented in Supplementary Fig. 2.
Comparisons Of Serum Lcfa Levels Between Patients With Sarcoidosis And Healthy Controls
The patients with sarcoidosis had significantly lower levels of n-3 PUFAs (p < 0.001) and n-6 PUFAs (p < 0.001) than the healthy controls (Fig. 2A, B). When the LCFAs were evaluated separately, the levels of all n-3 PUFAs (linolenic acid, eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid) were significantly lower in the patients with sarcoidosis (p < 0.001 for all) (Supplementary Fig. 3). Likewise, the levels of all n-6 PUFAs other than eicosadienoic acid (linoleic acid, γ-linolenic acid, dihomo-γ-linolenic acid, arachidonic acid, docosatetraenoic acid) were significantly lower in the patients with sarcoidosis (p < 0.001, p = 0.016, p = 0.007, p < 0.001, and p < 0.001, respectively) (Supplementary Fig. 3). However, there were no significant differences in the levels of n-3/n-6 ratio, SFAs, and MUFAs between the two groups (Fig. 2C–E). Furthermore, the levels of myristic acid (p = 0.004), myristoleic acid (p < 0.001), and palmitoleic acid (p < 0.001) were significantly higher in the patients with sarcoidosis, while the levels of other SFAs and MUFAs did not differ between the two groups (Supplementary Fig. 3).
On univariate logistic analysis, higher levels of MUFAs and lower levels of n-3 PUFAs, n-6 PUFAs, and n-3/n-6 ratio were predictive of sarcoidosis, as were age and sex (Table 2). On multivariate logistic analysis, lower levels of n-3 PUFAs, n-6 PUFAs, and n-3/n-6 ratio were predictive of sarcoidosis, as was age (Table 2). When the LCFAs were evaluated separately, higher levels of myristic acid, myristoleic acid, palmitoleic acid, and eicosenoic acid, and lower levels of arachidic acid, behenic acid, linolenic acid, eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid, linoleic acid, γ-linolenic acid, dihomo-γ-linolenic acid, arachidonic acid, and docosatetraenoic acid were predictive of sarcoidosis on univariate logistic analysis (Supplementary Table 5). After adjustment by age and sex, higher levels of myristoleic acid and palmitoleic acid, and lower levels of stearic acid, arachidic acid, behenic acid, lignoceric acid, linolenic acid, eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid, linoleic acid, γ-linolenic acid, dihomo-γ-linolenic acid, arachidonic acid, and docosatetraenoic acid were independently predictive of sarcoidosis (Supplementary Table 5).
Table 2
Logistic regression analyses for the diagnosis of sarcoidosis
| Univariate analysis | Multivariate analysis |
Set A | Set B | Set C |
Variables | OR (95%CI) | P-value | OR (95%CI) | P-value | OR (95%CI) | P-value | OR (95%CI) | P-value |
Agea | 1.11 (1.06–1.16) | < 0.001 | 1.20 (1.10–1.31) | < 0.001 | 1.11 (1.05–1.18) | < 0.001 | 1.11 (1.06–1.16) | < 0.001 |
Sex, male | 0.43 (0.19–0.98) | 0.044 | 2.66 (0.44–16.30) | 0.289 | 1.64 (0.38–7.01) | 0.505 | 0.81 (0.26–2.59) | 0.727 |
BMIb | 0.97 (0.83–1.13) | 0.671 | | | | | | |
n-3 PUFAs, high | 0.05 (0.01–0.18) | < 0.001 | 0.01 (0.00-0.05) | < 0.001 | | | | |
n-6 PUFAs, high | 0.06 (0.02–0.19) | < 0.001 | | | 0.03 (0.01–0.14) | < 0.001 | | |
n-3/n-6 ratio, high | 0.34 (0.13–0.90) | 0.030 | | | | | 0.25 (0.07–0.83) | 0.023 |
SFAs, high | 2.18 (0.85–5.61) | 0.105 | | | | | | |
MUFAs, high | 3.26 (1.14–9.34) | 0.028 | 1.13 (0.15–8.71) | 0.906 | 4.27 (0.87-21.00) | 0.074 | 1.59 (0.38–6.58) | 0.526 |
The cutoff value for each long-chain fatty acid was determined by the Youden Index in receiver operating characteristic curve analysis (Supplementary Table 1). The Akaike Information criteria for Sets A, B, and C were 58.6, 76.1, and 96.4, respectively. OR, odds ratio; CI, confidence interval; BMI, body mass index; SFAs, saturated fatty acids; PUFAs, polyunsaturated fatty acids; MUFAs, monounsaturated fatty acids. |
aPer 1-year increase. |
bPer 1-kg/m2 increase. |
Associations Of Serum Lcfa Levels With Multiple Organ Involvement In Patients With Sarcoidosis
Among the patients with sarcoidosis, those with multiple organ involvement had significantly lower levels of n-3 PUFAs and n-3/n-6 ratio than those with single organ involvement (Fig. 2A, C). When the LCFAs were evaluated separately, the levels of all n-3 PUFAs except linolenic acid (eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid) were significantly lower in the patients with multiple organ involvement (p = 0.032, p = 0.021, and p = 0.033, respectively) (Supplementary Fig. 3). Meanwhile, there were no significant differences in the levels of n-6 PUFAs, SFAs, and MUFAs between the patients with multiple and single organ involvement (Fig. 2B, D, E).
On univariate logistic analysis, lower levels of SFAs, n-3 PUFAs, n-6 PUFAs, and n-3/n-6 ratio were predictive of multiple organ involvement (Table 3). On multivariate logistic analysis, lower levels of SFAs and n-3/n-6 ratio were predictive of multiple organ involvement (Table 3). When the LCFAs were evaluated separately, lower levels of palmitic acid, stearic acid, palmitoleic acid, eicosatrienoic acid, eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid, linoleic acid, γ-linolenic acid, and dihomo-γ-linolenic acid were predictive of multiple organ involvement on univariate logistic analysis, (Supplementary Table 6). After adjustment by age, sex, and serum angiotensin converting enzyme level, lower levels of palmitoleic acid, eicosatrienoic acid, eicosapentaenoic acid, docosapentaenoic acid, docosahexaenoic acid, γ-linolenic acid, and dihomo-γ-linolenic acid were independently predictive of multiple organ involvement (Supplementary Table 6).
Table 3
Logistic regression analyses for multiple organ involvement in sarcoidosis
| Univariate analysis | Multivariate analysis |
| Set 1 | Set 2 | Set 3 |
Variables | OR (95%CI) | P-value | OR (95%CI) | P-value | OR (95%CI) | P-value | OR (95%CI) | P-value |
Agea | 0.98 (0.95–1.02) | 0.339 | 0.99 (0.94–1.04) | 0.716 | 0.98 (0.93–1.02) | 0.300 | 1.01 (0.95–1.06) | 0.815 |
Sex, male | 0.55 (0.18–1.67) | 0.293 | 0.64 (0.15–2.67) | 0.537 | 0.53 (0.13–2.19) | 0.383 | 0.72 (0.16–3.25) | 0.664 |
BMIb | 1.01 (0.82–1.25) | 0.930 | | | | | | |
Serum ACE, ≥ 21.4 IU/L | 1.08 (0.99–1.18) | 0.096 | 2.41 (0.64–9.14)) | 0.196 | 2.40 (0.64–9.03) | 0.195 | 2.53 (0.62–10.30) | 0.195 |
Radiographic stage, ≥ Ⅱ | 1.27 (0.41–3.95) | 0.676 | | | | | | |
Pulmonary function test | | | | | | | | |
FVC % predictedc | 0.97 (0.93–1.02) | 0.240 | | | | | | |
FEV1/FVCc | 1.00 (0.98–1.01) | 0.692 | | | | | | |
Bronchoalveolar lavage | | | | | | | | |
Lymphocytesc | 1.01 (0.97–1.04) | 0.767 | | | | | | |
CD4/CD8 ratio, ≥ 3.5 | 0.99 (0.29–3.39) | 0.992 | | | | | | |
n-3 PUFAs, high | 0.23 (0.07–0.74) | 0.013 | 0.33 (0.09–1.24) | 0.101 | | | | |
n-6 PUFAs, high | 0.30 (0.09–1.05) | 0.059 | | | 0.43 (0.11–1.67) | 0.225 | | |
n-3/n-6 ratio, high | 0.16 (0.05–0.53) | 0.003 | | | | | 0.14 (0.03–0.61) | 0.009 |
SFAs, high | 0.16 (0.03–0.76) | 0.022 | 0.25 (0.05–1.35) | 0.107 | 0.25 (0.05–1.34) | 0.105 | 0.16 (0.03–0.96) | 0.045 |
MUFAs, high | 1.00 (0.18–5.69) | 1.000 | | | | | | |
The cutoff value for each long-chain fatty acid was determined by the Youden Index in receiver operating characteristic curve analysis (Supplementary Table 2). The Akaike Information criteria for Sets 1, 2, and 3 were 74.1, 75.4, and 69.1, respectively. OR, odds ratio; CI, confidence interval; BMI, body mass index; ACE, angiotensin-converting enzyme; FVC, forced vital capacity; FEV1, forced expiratory volume in 1 second; SFAs, saturated fatty acids; PUFAs, polyunsaturated fatty acids; MUFAs, monounsaturated fatty acids. |
aPer 1-year increase. |
bPer 1-kg/m2 increase. |
cPer 1% increase. |
Associations Of Serum Lcfa Levels With Other Outcomes In Patients With Sarcoidosis
There were no significant associations between lung parenchymal involvement (radiographic stage ≥ 2) and n-3 PUFAs, n-6 PUFAs, n-3/n-6 ratio, SFAs, or MUFAs (Supplementary Fig. 4). Among 51 patients who were followed up without treatment, there were no significant associations between disease progression and n-3 PUFAs, n-6 PUFAs, n-3/n-6 ratio, SFAs, or MUFAs (Supplementary Fig. 5).