# Role of serum immunoglobulins for predicting sarcoidosis outcome A cohort study

Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Olivier Decaux,

Mathieu Lederlin, Stéphanie Guillot, Antoinette Perlat, Patrick Jego

► To cite this version:

Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Olivier Decaux, Mathieu Lederlin, et al. Role of serum immunoglobulins for predicting sarcoidosis outcome A cohort study. PLoS ONE, 2018, 13 (4), pp.e0193122. 10.1371/journal.pone.0193122. hal-01771420

HAL Id: hal-01771420

https://univ-rennes.hal.science/hal-01771420

Submitted on 19 Apr 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Role of serum immunoglobulins for predicting sarcoidosis outcome: A cohort study

Nicolas Belhomme<sup>1,2*</sup>, Stéphane Jouneau<sup>2,3,4</sup>, Guillaume Bouzillé<sup>5,6,7,8</sup>, Olivier Decaux<sup>1,2</sup>, Mathieu Lederlin<sup>2,5,9</sup>, Stéphanie Guillot<sup>10</sup>, Antoinette Perlat<sup>1</sup>, Patrick Jégo<sup>1,2,4</sup>

1 Internal Medicine Department, Rennes University Hospital, Rennes, France, 2 University of Rennes 1, Rennes, France, 3 Department of Respiratory Medicine, Rennes University Hospital, Rennes, France. University of Rennes 1, Rennes, France, 4 INSERM-IRSET UMR1085, Rennes, France, 5 INSERM, U1099, Rennes, France, 6 Université de Rennes 1, LTSI, Rennes, France, 7 CHU Rennes, CIC Inserm 1414, Rennes, France, 8 CHU Rennes, Centre de Données Cliniques, Rennes, France, 9 Department of Radiology, Rennes University Hospital, Rennes, France, 10 Department of Respiratory Physiology, Rennes University Hospital, Rennes, France

* nicolas.belhomme@chu-rennes.fr

<figure><img src="image_2.png" /></figure>

## Abstract

### Background

Sarcoidosis is a systemic granulomatous disease which carries variable outcomes. Serum protein electrophoresis is an easily accessible and routinely performed examination at diagnosis, in order to search for hypergammaglobulinemia, which is frequently found, and to rule out other granulomatous diseases such as common variable immunodeficiency. We aimed to assess the impact of baseline immunoglobulin level on the outcome of sarcoidosis.

### Methods

We conducted a retrospective cohort-study, at Rennes University Hospital, in which all newly diagnosed patients for whom a serum protein electrophoresis had been performed at baseline were enrolled, from 2006 to 2014. The main outcome was the need for corticosteroid treatment within 2 years from diagnosis, the secondary outcome was the occurrence of relapse among treated patients.

### Results

Eighty patients were included in the study, and 41.25% of them exhibited an elevated globulins rate. In univariate analysis, an elevated ACE level >70 U/l, Afro-Caribbean origin, and extra-pulmonary involvement, were associated with the need for corticosteroid treatment. In multivariate analysis, only ACE elevation (OR = 1.03, IC95% 1.01–1.05, p = 0.009) and extra-pulmonary involvement (OR = 5.8, IC95% 1.4–24, p = 0.015) were significant. Immunoglobulin level was not associated with the main outcome. Regarding the secondary outcome, none of the studied features were predictive of relapse among the 34 treated patients followed for two years.

<figure><img src="image_3.png" /></figure>

**OPEN ACCESS**

**Citation:** Belhomme N, Jouneau S, Bouzillé G, Decaux O, Lederlin M, Guillot S, et al. (2018) Role of serum immunoglobulins for predicting sarcoidosis outcome: A cohort study. PLoS ONE 13(4): e0193122. https://doi.org/10.1371/journal.pone.0193122

**Editor:** Elliott D. Crouser, Ohio State University Wexner Medical Center, UNITED STATES

**Received:** January 3, 2018

**Accepted:** February 5, 2018

**Published:** April 11, 2018

**Copyright:** © 2018 Belhomme et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The authors received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.
## Conclusions

There was no relation between the immunoglobulin level at diagnosis and the evolution of sarcoidosis. An elevated ACE level and the presence of initial extra-pulmonary involvement were both associated with a more severe course of the disease necessitating a corticosteroid treatment.

## Introduction

Sarcoidosis is a systemic disease of unknown etiology, which is characterized by the development of non-necrotizing epithelioid granulomas in involved organs [1]. The course of the disease is unpredictable: spontaneous remission may occur while 20 to 70% of the patients will need a systemic therapy, for which oral corticosteroids are considered as the first-line option [1-3]. Relapse arises in 37 to 74% of treated patients [4]. Biologically, an elevated ACE is found in 60% of the patients [5]. Other parameters, such as interleukin-2 receptor, neopterin, or chitotriosidase, have been proposed for monitoring the activity of the disease, nevertheless their prognosis value needs to be confirmed, and their accessibility in routine practice is limited [6-8]. Hypergammaglobulinemia is frequently observed among sarcoidosis patients, but its relation to the disease's course has not been established [9,10]. It has been suggested that elevated immunoglobulins or circulating immune-complex may be correlated to the disease's course [11, 12], nonetheless no recent study has aimed to further explore this hypothesis to our knowledge [13]. As immunoglobulin level is highly influenced by many factors such as corticosteroid treatment or intercurrent infections, it is not appropriate for monitoring the disease's activity. However, serum protein electrophoresis (SPE) is frequently performed at diagnosis, and easily accessible in routine practice. Hence, the aim of our study was to investigate the potential role of the immunoglobulin level at the time of diagnosis for predicting the disease's course.

## Methods

### Study design

We conducted a retrospective cohort study in Rennes University Hospital, a tertiary care center referral for rare pulmonary diseases, between January, 1st 2006 and December, 31st 2014. All of the data used in this study existed prior to the study and the request for ethical approval, and was collected as part of standard medical procedure. The authors had access to the identity of included patients.

Ethics approval was received from the Rennes University Hospital ethics committee (number 16.65, 2016, June 6th). The need for consent was waived by the ethics committee for this particular study, and only the non-opposition of participating patients had to be collected. Therefore, all patients received a written information about the study and were given the possibility to decline to participate by completing and returning an opposition form which was specific to this particular study. Minors were not included. All patients newly diagnosed with sarcoidosis during the study period were included. Potential cases were retrieved from the institutional clinical data warehouse, eHOP software [14], and their medical records were individually reviewed. Sarcoidosis was diagnosed in accordance with the American Thoracic Society and World Association of Sarcoidosis and Other Granulomatous disorders (WASOG) criteria [2]. Histological confirmation was not needed for patients presenting with typical
Löfgren's syndrome [1,2]. The time of diagnosis was defined as the date of the first positive histopathological result or the date of the clinical diagnosis.

Inclusion required for a SPE to have been performed in our center within 3 months from diagnosis, before any systemic anti-sarcoidosis medication was started, and for the patients not to have been exposed to corticosteroids within a period of 4 month prior to SPE. Patients with chronic diseases that affect immunoglobulins amount and function (e.g primary immune deficiency, cirrhosis, hemopathy, HIV infection) were excluded. Patients were followed up for 2 years, in accordance with previous studies addressing the disease's prognosis [15,16]. The primary outcome was the need for corticosteroid (CS) treatment. The decision to treat was taken according to guidelines [1,2,17]. The secondary outcome was the association between the baseline gammaglobulin level and the occurrence of relapse. Organ involvement was determined using the 2014 WASOG assessment instrument, and only involvements classified as "highly probable" or "at least probable" were considered [18,19]. Organ involvements were deemed extra-pulmonary if they did not involve lungs or mediastinal lymph nodes. Neither hypercalcemia, nor Löfgren's syndrome or erythema nodosum were considered as extra-pulmonary involvement. Visceral involvement was defined by the involvement of the heart, liver, central nervous system, kidney or spleen. Relapse was defined by the necessity of escalating CS therapy due to progressive disease, or of resuming CS therapy after cessation [20].

Between January 1st 2006 and December 31st 2014, sarcoidosis was diagnosed according to WASOG criteria in 97 patients in our center. SPE had been performed at diagnosis in 82 of them. Among those, one had received corticosteroids within one month prior to diagnosis, and one refused to participate. Consequently, 80 patients were finally included in the analysis. Three patients were excluded during the follow-up period: 2 were lost to follow up, and one developed multiple myeloma requiring chemotherapy. All of these 3 patients were included in the analysis for main outcome, as the severity of their disease brought the need for CS treatment before they were excluded from the analysis. They were not analyzed for secondary outcome as they didn't complete the follow up period (Fig 1).

## Demographics and clinical data collection

We gathered demographic variables (age, ethnicity, gender), clinical variables (histological confirmation, involvement of the organs and manifestations of the disease). Pulmonary function tests (Forced Expiratory Volume in 1 second-VEF₁, Total Lung Capacity-TLC and Diffusing Capacity of the Lung for Carbon Monoxide-DLCO) were analyzed, and roentgenographic stage at baseline was retrieved from the patient's records. Among biological variables, we collected Angiotensin converting Enzyme (ACE) (normal range:20-70 U/l), lymphocytes count (normal range: 1000-4000 cells/µl), gammaglobulin level (normal range: 8-13.5g/l), creatininemia (normal: 0.67-1.18 mg/dl), and calcemia (normal: 2.26-2.62 mmol/l). SPE were all performed in our center using a capillary electrophoresis (Capillarys 2 from Sebia Inc, Norcross, USA, with capillary Protein(e) 6 assay, ref 2003). Normal range was 8-13.5 g/l.

## Statistical analysis

Quantitative variables are given with their mean±SD or with their median and interquartile range (IQR) depending on their distribution. Comparison between quantitative variables were analyzed by conducting Student or Mann Whitney U test. The analysis of predictive factors associated with study outcomes were conducted by logistic regression using Firth method for correction of bias. Results are given with their Odds-Ratio (OR) and 95% confidence interval (CI95%). All variables were studied in univariate analysis. Factors for which the corresponding p-value was <0.10 were included in multivariate analysis and a stepwise selection procedure
<figure><img src="image_4.png" /><figcaption>**Fig 1. Study profile.** CS = Corticosteroids; FU = Follow-Up; SPE = Serum Protein Electrophoresis.</figcaption></figure>

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

was applied to get the final models. Missing data were taken into account with the use of the multiple imputation by chained equations procedure. Five datasets were imputed to perform analyses. We followed Rubin's rules to obtain the final results. The rates of missing data were in every case in the range of 0-6% (variables associated with the highest rates of missing data were TLC and DLCO with a rate of 6%).
<table><caption>**Table 1. Patients characteristics.**</caption><thead><tr><th></th><th></th></tr></thead><tbody><tr><td><strong>Patients characteristics</strong></td><td></td></tr><tr><td>Sex ratio (f/m)</td><td>1.5 (48/32)</td></tr><tr><td>Age (mean±SD)</td><td>39±14.69</td></tr><tr><td><strong>Clinical features</strong></td><td>n(%)</td></tr><tr><td>Ethnicity</td><td></td></tr><tr><td>Caucasian</td><td>72(90)</td></tr><tr><td>Caribbean</td><td>3(3.75)</td></tr><tr><td>African</td><td>5(6.25)</td></tr><tr><td>Histological confirmation</td><td>63(78.75)</td></tr><tr><td>Löfgren' syndrome</td><td>26(32.5)</td></tr><tr><td>Extra-pulmonary involvement</td><td>44(55)</td></tr><tr><td>Periph. adenopathy</td><td>21(26.25)</td></tr><tr><td>Ocular</td><td>10(12.5)</td></tr><tr><td>Arthritis</td><td>9(11.25)</td></tr><tr><td>Cutaneous sarcoid</td><td>9(11.25)</td></tr><tr><td>ENT</td><td>8(10)</td></tr><tr><td>Liver</td><td>7(8.75)</td></tr><tr><td>Spleen</td><td>5(6.25)</td></tr><tr><td>Nervous System</td><td>3(3.75)</td></tr><tr><td>Muscles</td><td>2(2.5)</td></tr><tr><td>Heart</td><td>2(2.5)</td></tr><tr><td>Bones</td><td>1(1.25)</td></tr><tr><td>Respiratory symptoms</td><td></td></tr><tr><td>Dyspnea</td><td>23(28.75)</td></tr><tr><td>Cough</td><td>19(23.75)</td></tr><tr><td><strong>Chest X-ray stage</strong></td><td>n(%)</td></tr><tr><td>Normal</td><td>9(11.25)</td></tr><tr><td>I</td><td>38(47.5)</td></tr><tr><td>II</td><td>29(36.25)</td></tr><tr><td>III</td><td>4(5)</td></tr><tr><td>IV</td><td>0</td></tr><tr><td><strong>Biological data</strong></td><td>Median (Q1-Q3)</td></tr><tr><td>Immunoglobulin level (g/l)</td><td>12.90(10.17-14.57)</td></tr><tr><td>Blood lymphocytes count (/µl)</td><td>1.21(0.87-1.58)</td></tr><tr><td>ACE (U/l) *</td><td>67(44-92.75)</td></tr><tr><td>Calcemia (mmol/l)</td><td>2.39(2.32-2.46)</td></tr><tr><td>Creatininemia (mg/dl)</td><td>0.82(0.68-0.93)</td></tr><tr><td><strong>Pulmonary function tests (% Pred)</strong></td><td>Median (Q1-Q3)</td></tr><tr><td>FEV<sub>1</sub>**</td><td>103.9(91.9-112.9)</td></tr><tr><td>TLC***</td><td>102(92.85-113.95)</td></tr><tr><td>DLCO***</td><td>70.20(62.65-80.90)</td></tr><tr><td><strong>Treatments received</strong></td><td>n(%)</td></tr><tr><td>Corticosteroids</td><td>37(46.25)</td></tr><tr><td>Hydroxychloroquine</td><td>3(3.75)</td></tr><tr><td><strong>Relapse</strong></td><td>10(12.5)</td></tr></tbody></table>

<footnote>* Data available for 78 patients.</footnote>

<footnote>** Data available for 77 patients;</footnote>

<footnote>*** Data available for 75 patients.</footnote>

<footnote>ACE = Angiotensin Converting Enzyme; DLCO = Diffusion Capacity of Lungs for Carbon Monoxide; ENT = Ear, Nose and Throat; FEV₁ = Forced Expiried Volume in one second; TLC = Total Lung Capacity.</footnote>
In order to test specifically the prognostic value of serum immunoglobulin level, this variable was added separately into the logistic regression model, and the area under curve (AUC) of the two predictive models (without immunoglobulin level and after including immunoglobulin level) were compared, and the likelihood ratio was calculated. Analysis were conducted with R version 3.3.3. software.

## Results

### Characteristics of the study population

Mean age of the patients was 39.0±14.69 years, sex ratio (f/m) was 1.5. Regarding ethnicity, 69 (86.25%) patients were Caucasian, 8 (10%) were Afro-Caribbean. Forty-four (55%) patients had an extra-pulmonary involvement at diagnosis. Chest X-ray showed a stage 1 in most of the patients (n = 38, 47.5%). Median immunoglobulin level was 12.90 g/l (10.17–14.57, extremes 8.1–35.9). Thirty-three (41.25%) patients presented a hypergammaglobulinemia >13.5g/l. A CS therapy was needed for 37 (46.25%) patients during the study period. Population characteristics are presented in Table 1.

### Relation of immunoglobulin level to individual features

Median immunoglobulin level was significantly higher in Afro-Caribbean patients (19.45 versus 12.56 g/l, p<0.001) and in patients presenting an extra-pulmonary involvement (13.65 versus 11.85, p = 0.016). Patients presenting a Löfgren's syndrome had a lower immunoglobulin level (10.90 versus 13.55, p = 0.011). We found a relation neither between immunoglobulin level and ACE level (p = 0.11), nor with a decrease in pulmonary function test parameters. Patients who necessitated a corticosteroid treatment during the study period and those who presented a relapse had a significantly higher immunoglobulin level (13.60 versus 12.40 g/l, p = 0.036 and 15.15 versus 12.53g/l, p = 0.014, respectively) (Table 2).

### Results for main outcome

In univariate analysis, neither the immunoglobulin level at diagnosis, nor hypergammaglobulinemia, were associated with CS treatment (OR = 1.1, CI95% 1–1.3, p = 0.059 and OR = 2, CI95% 0.79–5.2, p = 0.137, respectively). Factors that were predictive of the need for treatment were the Afro-Caribbean origin (OR = 24, CI95% 1–554, p = 0.049), an ACE level >70 U/L (OR = 6.2, CI95% 2.2–17, p<0.001), extra-pulmonary involvement (OR = 6.7, CI95% = 2.4–19, p<0.001) including cutaneous localization (OR = 7.7, CI95% = 1.1–55, p = 0.041), ENT (OR = 7.7, IC95% = 1.1–55, p = 0.041) and lymph nodes involvement (OR = 5.6, CI95% 1.7–18, p = 0.004). Löfgren's syndrome was associated with treatment non-necessity (OR = 0.097, CI95% 0.027–0.35, p<0.001). The results of the univariate analysis are shown in Table 3. In multivariate analysis, both extra-pulmonary involvement (OR = 5.8, CI95% 1.4–24, p = 0.0152) and an elevated ACE (OR = 1.03, CI95% = 1.01–1.05, p = 0.009, meaning that each elevation of one unit of ACE increased the probability for treatment of 3%) were predictive of the need for CS (Table 4, Model 1). These findings remained unchanged after including the immunoglobulin level into the logistic regression model (Table 4, Model 2). The likelihood ratio for the two models was 1, reflecting the absence of influence of the immunoglobulin level on the main outcome.

### Results for secondary outcome

Relapse was analyzed in the 34 patients who received CS and were followed for two years. Their characteristics are presented in an additional table (S1 Table). A relapse occurred in 10 of them. The duration of exposure to the secondary outcome was comparable (22 months from
<table><caption>Table 2. Relation between immunoglobulins and individual characteristics.</caption><thead><tr><th>Variables</th><th>Median value of immunoglobulins (Median Q1-Q3)</th><th>p-value</th></tr></thead><tbody><tr><td>Age</td><td></td><td></td></tr><tr><td>&gt;40</td><td>13.0 (10.2-14.5)</td><td></td></tr><tr><td>≤40</td><td>12.7 (10.2-14.6)</td><td>0.93</td></tr><tr><td>Gender</td><td></td><td></td></tr><tr><td>Female n = 48</td><td>12.90 (10.25-14.57)</td><td></td></tr><tr><td>Male n = 32</td><td>13.05 (10.17-14.50)</td><td>0.976</td></tr><tr><td>Ethnicity</td><td></td><td></td></tr><tr><td>Afro-Caribbean n = 8</td><td>19.45 (14.85-25.27)</td><td></td></tr><tr><td>Other n = 72</td><td>12.56 (10.07-14.05)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Chest X-ray stage</td><td></td><td></td></tr><tr><td>0 or 1 n = 47</td><td>12.90 (10.25-14.65)</td><td></td></tr><tr><td>&gt;1 n = 33</td><td>13.00 (9.90-14.20)</td><td>0.957</td></tr><tr><td>Extra-pulmonary involvement</td><td></td><td></td></tr><tr><td>No n = 36</td><td>11.85 (9.83-13.60)</td><td></td></tr><tr><td>Yes n = 44</td><td>13.65 (10.70-15.43)</td><td><strong>0.016</strong></td></tr><tr><td>Löfgren's syndrome</td><td></td><td></td></tr><tr><td>Yes n = 26</td><td>10.90 (9.98-12.90)</td><td></td></tr><tr><td>No n = 54</td><td>13.55 (10.33-15.47)</td><td><strong>0.011</strong></td></tr><tr><td>ACE*</td><td></td><td></td></tr><tr><td>&gt;70 U/l</td><td>13.7 (10.4-15.5)</td><td></td></tr><tr><td>≤70 U/l</td><td>11.6 (9.7-13.5)</td><td>0.11</td></tr><tr><td>Pulmonary function test</td><td></td><td></td></tr><tr><td>FEV<sub>1</sub>**</td><td>≥80%</td><td></td></tr><tr><td></td><td>&lt;80%</td><td>0.83</td></tr><tr><td>TLC***</td><td>≥80%</td><td></td></tr><tr><td></td><td>&lt;80%</td><td>0.11</td></tr><tr><td>DLCO***</td><td>≥75%</td><td></td></tr><tr><td></td><td>&lt;75%</td><td>0.23</td></tr><tr><td>CS therapy</td><td></td><td></td></tr><tr><td>No (n = 43)</td><td>12.40 (9.95-13.90)</td><td></td></tr><tr><td>Yes (n = 37)</td><td>13.60 (10.55-15.45)</td><td><strong>0.036</strong></td></tr><tr><td>Relapse</td><td></td><td></td></tr><tr><td>No (n = 24)</td><td>12.53 (10.15-13.95)</td><td></td></tr><tr><td>Yes (n = 10)</td><td>15.15 (13.82-16.82)</td><td><strong>0.014</strong></td></tr></tbody></table>

ACE = Angiotensin Converting Enzyme; DLCO = Diffusion Capacity of Lungs for Carbon Monoxide; FEV₁ = Forced Expired Volume in one second; TLC = Total Lung Capacity.

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

start of treatment to end of study in patients with relapse versus 20 months for non-relapsing patients, p = 0.24). In univariate analysis, we didn't find any factors that predict relapse. Particularly, neither the immunoglobulin level (OR = 1.1, CI 95% 0.96-1.4, p = 0.122), nor hypergammaglobulinemia (OR = 5.5, CI95% 0.97-32, p = 0.054) were associated with relapse (Table 3). In multivariate analysis, variables which were retained in the model were dyspnea (OR = 5.1, CI95% 0.85-31, p = 0.07) and visceral involvement (OR = 5.7, CI95% 0.87-38,

<footnote>* Data available for 78 patients</footnote>

<footnote>** Data available for 77 patients</footnote>

<footnote>*** Data available for 75 patients.</footnote>
<table><caption>Table 3. Features associated with corticosteroid therapy and relapse in univariate analysis.</caption><thead><tr><th rowspan="2">Variables</th><th colspan="3">Main outcome(n = 80)</th><th colspan="3">Secondary outcome (n = 34)</th></tr><tr><th>OR</th><th>IC95%</th><th>p</th><th>OR</th><th>CI95.</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>1</td><td>0.98-1.04</td><td>0.659</td><td>0.97</td><td>0.9-1</td><td>0.315</td></tr><tr><td>Gender (female)</td><td>1.8</td><td>0.71-4.7</td><td>0.212</td><td>0.69</td><td>0.15-3.2</td><td>0.629</td></tr><tr><td>Ethnicity</td><td>24</td><td>1-554</td><td><strong>0.049</strong></td><td>1</td><td>0.16-6.8</td><td>0.964</td></tr><tr><td>X-ray stage >1</td><td>2.5</td><td>0.98-6.6</td><td>0.055</td><td>0.85</td><td>0.19-3.9</td><td>0.831</td></tr><tr><td>Immunoglobulin level</td><td>1.1</td><td>0.99-1.3</td><td>0.059</td><td>1.1</td><td>0.96-1.4</td><td>0.122</td></tr><tr><td>Immunoglobulins >13.5g/l</td><td>2</td><td>0.79-5.2</td><td>0.137</td><td>5.5</td><td>0.97-32</td><td>0.054</td></tr><tr><td>Lymphocytes count</td><td>0.96</td><td>0.55-1.7</td><td>0.886</td><td>0.67</td><td>0.21-2.2</td><td>0.494</td></tr><tr><td>ACE</td><td>1</td><td>1-1.05</td><td><strong>0.002</strong></td><td>1</td><td>0.99-1.03</td><td>0.143</td></tr><tr><td>ACE >70U/l</td><td>6.2</td><td>2.2-17</td><td><strong>&lt;0.001</strong></td><td>3.3</td><td>0.41-26</td><td>0.256</td></tr><tr><td>Calcemia</td><td>2.9</td><td>0.86-10</td><td>0.084</td><td>0.97</td><td>0.59-1.6</td><td>0.892</td></tr><tr><td>Creatininemia</td><td>5.1</td><td>0.49-52</td><td>0.169</td><td>0.54</td><td>0.022-14</td><td>0.702</td></tr><tr><td>Dyspnea</td><td>1.8</td><td>0.64-4.9</td><td>0.264</td><td>4.1</td><td>0.82-21</td><td>0.083</td></tr><tr><td>Cough</td><td>0.42</td><td>0.13-1.3</td><td>0.141</td><td>1.8</td><td>0.24-14</td><td>0.557</td></tr><tr><td>Löfgren's syndrome</td><td>0.09</td><td>0.03-0.3</td><td><strong>&lt;0.001</strong></td><td>0.29</td><td>0.0076-11</td><td>0.497</td></tr><tr><td>Extra-pulmonary involvement</td><td>6.7</td><td>2.4-19</td><td><strong>&lt;0.001</strong></td><td>9</td><td>0.34-242</td><td>0.183</td></tr><tr><td>Visceral involvement</td><td>6.2</td><td>1.3-30</td><td><strong>0.023</strong></td><td>4.6</td><td>0.84-25</td><td>0.0776</td></tr><tr><td>Arthritis</td><td>0.8</td><td>0.2-3.6</td><td>0.743</td><td>0.29</td><td>0.0076-11</td><td>0.497</td></tr><tr><td>Cutaneous</td><td>7.7</td><td>1.1-55</td><td><strong>0.041</strong></td><td>2.1</td><td>0.35-13</td><td>0.397</td></tr><tr><td>Ocular</td><td>1.6</td><td>0.39-6.8</td><td>0.492</td><td>1.8</td><td>0.24-14</td><td>0.557</td></tr><tr><td>ENT</td><td>7.7</td><td>1.1-55</td><td><strong>0.041</strong></td><td>1</td><td>0.16-6.8</td><td>0.964</td></tr><tr><td>Periph. adenopathy</td><td>5.6</td><td>1.7-18</td><td><strong>0.004</strong></td><td>2.4</td><td>0.5-11</td><td>0.268</td></tr><tr><td>Spleen</td><td>4.2</td><td>0.51-35</td><td>0.181</td><td>7.3</td><td>0.67-80</td><td>0.1</td></tr><tr><td>Liver</td><td>6.5</td><td>0.88-47</td><td>0.0657</td><td>6.2</td><td>0.87-45</td><td>0.068</td></tr><tr><td>Nervous system</td><td>2.2</td><td>0.19-25</td><td>0.524</td><td>2.5</td><td>0.12-49</td><td>0.542</td></tr><tr><td>TLC</td><td>0.97</td><td>0.93-1</td><td>0.044</td><td>0.98</td><td>0.92-1.1</td><td>0.631</td></tr><tr><td>FEV<sub>1</sub></td><td>0.97</td><td>0.94-1</td><td>0.115</td><td>1</td><td>0.96-1.1</td><td>0.667</td></tr><tr><td>DLCO</td><td>0.95</td><td>0.91-1</td><td>0.033</td><td>0.99</td><td>0.93-1.1</td><td>0.834</td></tr></tbody></table>

ACE = Angiotensin Converting Enzyme; DLCO = Diffusion Capacity of Lungs for Carbon Monoxide; FEV₁ = Forced Expired Volume in one second; TLC = Total Lung Capacity.

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

p = 0.06), and none were significant. Once added into the model, the immunoglobulin level was also not significant (OR = 1.1, CI95% 0.91-1.4, p = 0.326). However, it resulted in an improvement of the relapse predicting value of the model. Indeed, the likelihood ratio of the

<table><caption>Table 4. Features associated with corticosteroid therapy in multivariate analysis.</caption><thead><tr><th rowspan="2">Variables</th><th colspan="3">Model 1</th><th colspan="3">Model 2</th></tr><tr><th>OR</th><th>CI95%</th><th>p</th><th>OR</th><th>CI95%</th><th>P</th></tr></thead><tbody><tr><td>Ethnicity</td><td>20</td><td>0.6-659</td><td>0.093</td><td>19</td><td>0.61-618</td><td>0.092</td></tr><tr><td>Extra-pulmonary involvement</td><td>5.8</td><td>1.4-24</td><td><strong>0.015</strong></td><td>5.6</td><td>1.4-22</td><td><strong>0.016</strong></td></tr><tr><td>Chest X-ray stage</td><td>4</td><td>0.99-16</td><td>0.052</td><td>3.8</td><td>0.97-15</td><td>0.056</td></tr><tr><td>ACE</td><td>1.03</td><td>1.01-1.05</td><td><strong>0.009</strong></td><td>1.03</td><td>1.006-1.05</td><td><strong>0.010</strong></td></tr><tr><td>Immunoglobulin level</td><td>-</td><td>-</td><td>-</td><td>0.97</td><td>0.8-1.2</td><td>0.791</td></tr></tbody></table>

ACE = Angiotensin Converting Enzyme.

https://doi.org/10.1371/journal.pone.0193122.t004
two models favored the second one (p = 0.012), and the derivation area under curve was 0.77 (CI95%: 0.56-0.97) for model 2, comprising immunoglobulins, whereas it was 0.63 (CI95%: 0.41-0.85, non-significant) for model 1.

Discussion

At the end, although the patients who required CS therapy and those who experienced a relapse had a higher baseline immunoglobulin level, our results show that such a parameter is not of interest for predicting the need for CS therapy during the first 2 years following diagnosis. Among study parameters, only ACE and extra-pulmonary involvement were found to be associ- ated with the main outcome. In our study, 41.25% of the patients had a hypergammaglobuline- mia at diagnosis, an intermediate value compared to those found in other studies [11,21]. Patients from Afro-Caribbean ethnicity had a higher rate of serum immunoglobulin compared to Caucasians (19.45 versus 12.56 g/l, p<0.001). This finding is in line with the results of Kataria and Stilzbach [11,22]. Group comparison analysis revealed that patients undergoing corticoste- roid therapy had a higher immunoglobulin level than the non-treated patients (13.6 versus 12.4 g/l, p = 0.036). Treated patients who experienced a relapse had a higher rate of immunoglobulin as well (15.15 versus 12.53 g/l, p = 0.014). Multiple logistic regression showed that the serum immunoglobulin level was associated neither with disease severity as assessed by the need for CS therapy, nor with the occurrence of relapse. In multivariate analysis, extra-pulmonary involve- ment at diagnosis was associated with CS therapy (OR = 5.8, CI95% 1.4-24, p = 0.015). We used a stepwise selection procedure to identify the final predictors. This procedure is known to poten- tially fail at identifying the true predictors. However, our findings are consistent with previous studies, which showed that initial extra-pulmonary involvement entails a worsening in the dis- ease's course resulting in additional organs involvements [16,23].

An elevated ACE was also predictive of the need for treatment, and this result is in line with many studies which have discussed the potential role of ACE for monitoring disease activity [24,25]. Nevertheless, whether ACE is a reliable marker of disease activity remains very contro- versial. Moreover, its reliability suffers from significant individual variations even in healthy subjects, mainly due to genetic polymorphism: thus, we should remain very cautious regarding its interpretation [26-28]. Finally, we observed a clear-although not significant-trend on the association between the Afro-Caribbean ethnicity and the need for treatment. As already shown in previous studies, black ethnicity is associated with more severe forms of the disease [29-32], to a more frequent need for systemic therapy [30], and to a higher mortality [33,34]. Although there was also no association between black ethnicity and treatment in a study con- ducted by Baughman [35], the non-significance of this result in the present study may be due to insufficient power as few black patients were included.

There are several limitations to our study. First of all, the number of patients (n = 80) may have impaired the statistical power of our analysis. Indeed, our inclusion criteria were very restrictive, and only the patients whose diagnosis was histologically confirmed or who presented a typical Löfgren's syndrome were considered, in accordance to international WASOG criteria [1,2]. Nonetheless, these criteria may be too restrictive: it has been shown that a typical presen- tation on chest CT-scan, and hyperlymphocytosis with a CD4/CD8 ratio greater than 3.5 on broncho-alveolar lavage, combined with the exclusion of other possible diagnosis, are highly predictive of sarcoidosis [36,37]. Then, the possibility of retaining the diagnosis in such patients should be discussed: we may consider the sarcoidosis as “probable” in cases fulfilling these crite- ria, whereas it would be classified as “certain” once histological proof is obtained. The secondary outcome is limited by the few relapses observed: indeed, only 10 relapses occurred, concerning 29.4% of the 34 treated patients followed-up for 2 years. This is principally due to the definition
employed: we made the choice to consider the relapse only among the patients who received or had received a treatment. As a consensual international definition of relapse is lacking, we aimed to avoid confusion between relapse and disease progression, which could be defined by the need of introducing a treatment due to disease aggravation after a period of stability [17,20]. Including this feature of the disease into the analysis would have biased the main outcome. Furthermore, CS therapy on its own may be associated with a higher risk of relapse [15,17,38], and limiting the assessment of relapse to the treated patients rather than assessing the whole population permitted to prevent this interaction bias. Another explication to the low number of relapses is the duration of follow-up, limited to 2 years. Thereby, only early relapses were studied. However, the risk of relapse lasts for many years over the disease's course, and our results do not reflect the global risk of relapse of sarcoidosis, as it was not our purpose [3,4,17].

Moreover, the risk of relapse is highest in the year following therapy disruption: in our study 20 out of the 34 treated patients were still receiving therapy at the end of follow up, thus carrying a lower potentiality of relapse [4]. On top of that, the risk of relapse is higher in black patients, who were a small number in our cohort. Finally, the relapse rate observed in our study is comparable to the study of Rizzato et al., which was also conducted in an essentially Caucasian population [4]. It is worth noting that the inclusion of the immunoglobulin level variable improved substantially the predicting capacity for relapse of the model, which implies that immunoglobulin level at diagnosis may be of interest for predicting relapse. This hypothesis needs to be further investigated in larger studies.

## Conclusions

In this cohort study including 80 patients followed for two years after sarcoidosis was diagnosed, we didn't find an association between the immunoglobulin level at baseline and the need for CS therapy or the occurrence of relapse. Nevertheless, the patients requiring a treatment and those who experienced a relapse had a higher level of immunoglobulin, as well as Afro-Caribbean patients. The immunoglobulin level statistically improved the capacity of the model to predict relapse, bringing the need for further larger studies to confirm its potential prognosis value. Furthermore, a significant association was found between initial extra-pulmonary involvement or an elevated ACE and the need for CS therapy.

## Supporting information

**S1 Table. Characteristics of the 34 patients who received a corticosteroid therapy and who completed the 2 years follow-up period.**
*Data available for 33/34 patients;
**Data available for 32/34 patients.
(DOCX)

## Author Contributions

**Conceptualization:** Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Olivier Decaux, Mathieu Lederlin, Stéphanie Guillot, Antoinette Perlat, Patrick Jégo.

**Data curation:** Nicolas Belhomme, Mathieu Lederlin, Stéphanie Guillot, Antoinette Perlat.

**Formal analysis:** Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Mathieu Lederlin, Stéphanie Guillot, Antoinette Perlat, Patrick Jégo.

**Investigation:** Stéphane Jouneau, Antoinette Perlat.

**Methodology:** Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Olivier Decaux, Stéphanie Guillot, Antoinette Perlat, Patrick Jégo.
**Supervision:** Nicolas Belhomme, Stéphane Jouneau, Antoinette Perlat.

**Validation:** Nicolas Belhomme, Stéphane Jouneau, Guillaume Bouzillé, Olivier Decaux, Mathieu Lederlin, Antoinette Perlat, Patrick Jégo.

**Writing – original draft:** Nicolas Belhomme, Stéphane Jouneau, Antoinette Perlat.

**Writing – review & editing:** Nicolas Belhomme, Stéphane Jouneau, Antoinette Perlat.

## References

1. Valeyre D, Prasse A, Nunes H, Uzunhan Y, Brillet P-Y, Müller-Quernheim J. Sarcoidosis. Lancet. 2014; 383: 1155–1167. https://doi.org/10.1016/S0140-6736(13)60680-7 PMID: 24090799

2. Statement on sarcoidosis. Joint Statement of the American Thoracic Society (ATS), the European Respiratory Society (ERS) and the World Association of Sarcoidosis and Other Granulomatous Disorders (WASOG) adopted by the ATS Board of Directors and by the ERS Executive Committee, February 1999. Am J Respir Crit Care Med. 1999; 160: 736–755. https://doi.org/10.1164/ajrccm.160.2.ats4-99 PMID: 10430755

3. Gibson GJ, Prescott RJ, Muers MF, Middleton WG, Mitchell DN, Connolly CK, et al. British Thoracic Society Sarcoidosis study: effects of long term corticosteroid treatment. Thorax. 1996; 51: 238–247 PMID: 8779124

4. Rizzato G, Montemurro L, Colombo P. The late follow-up of chronic sarcoid patients previously treated with corticosteroids. Sarcoidosis Vasc Diffuse Lung Dis. 1998; 15: 52–58 PMID: 9572002

5. Iannuzzi MC, Rybicki BA, Teirstein AS. Sarcoidosis. N Engl J Med. 2007; 357: 2153–2165. https://doi.org/10.1056/NEJMra071714 PMID: 18032765

6. Bargagli E, Bennett D, Maggiorelli C, Di Sipio P, Margollicci M, Bianchi N, et al. Human chitotriosidase: a sensitive biomarker of sarcoidosis. J Clin Immunol. 2013; 33: 264–270. https://doi.org/10.1007/s10875-012-9754-4 PMID: 22878841

7. Miyoshi S, Hamada H, Kadowaki T, Hamaguchi N, Ito R, Irifune K, et al. Comparative evaluation of serum markers in pulmonary sarcoidosis. Chest. 2010; 137: 1391–1397. https://doi.org/10.1378/chest.09-1975 PMID: 20081103

8. Harlander M, Salobir B, Zupančič M, Dolenšek M, Bavčar Vodovnik T, Terčelj M. Serial chitotriosidase measurements in sarcoidosis—two to five year follow-up study. Respir Med. 2014; 108: 775–782. https://doi.org/10.1016/j.rmed.2014.02.002 PMID: 24594143

9. Takada K, Ina Y, Noda M, Sato T, Yamamoto M, Morishita M. The clinical course and prognosis of patients with severe, moderate or mild sarcoidosis. J Clin Epidemiol. 1993; 46: 359–366 PMID: 8483000

10. Schoenfeld N, Schmolke B, Schmitt M, Remy N, Ellensohn P, Wahn U, et al. Specification and quantitation of circulating immune complexes in the serum of patients with active pulmonary sarcoidosis. Thorax. 1994; 49: 688–691 PMID: 8066564

11. Kataria YP, Holter JF. Immunology of sarcoidosis. Clin Chest Med. 1997; 18: 719–739 PMID: 9413655

12. Rømer FK, Sølling J. Relationship between circulating immune complexes and angiotensin-converting enzyme in pulmonary sarcoidosis. Acta Med Scand. 1981; 210: 299–303 PMID: 6274169

13. Tzouvelekis A, Kouliatsis G, Anevavis S, Bouros D. Serum biomarkers in interstitial lung diseases. Respir Res. 2005; 6: 78. https://doi.org/10.1186/1465-9921-6-78 PMID: 16042760

14. Bouzillé G, Sylvestre E, Campillo-Gimenez B, Renault E, Ledieu T, Delamarre D, et al. An Integrated Workflow For Secondary Use of Patient Data for Clinical Research. Stud Health Technol Inform. 2015; 216: 913 PMID: 26262215

15. Baughman RP, Nagai S, Balter M, Costabel U, Drent M, du Bois R, et al. Defining the clinical outcome status (COS) in sarcoidosis: results of WASOG Task Force. Sarcoidosis Vasc Diffuse Lung Dis. 2011; 28: 56–64 PMID: 21796892

16. Judson MA, Baughman RP, Thompson BW, Teirstein AS, Terrin ML, Rossman MD, et al. Two year prognosis of sarcoidosis: the ACCESS experience. Sarcoidosis Vasc Diffuse Lung Dis. 2003; 20: 204–211 PMID: 14620163

17. Gottlieb JE, Israel HL, Steiner RM, Triolo J, Patrick H. Outcome in sarcoidosis. The relationship of relapse to corticosteroid therapy. Chest. 1997; 111: 623–631 PMID: 9118698

18. Judson MA. Extrapulmonary sarcoidosis. Semin Respir Crit Care Med. 2007; 28: 83–101. https://doi.org/10.1055/s-2007-970335 PMID: 17330194
19. Judson MA, Costabel U, Drent M, Wells A, Maier L, Koth L, et al. The WASOG Sarcoidosis Organ Assessment Instrument: An update of a previous clinical tool. Sarcoidosis Vasc Diffuse Lung Dis. 2014; 31: 19–27 PMID: 24751450

20. Baughman RP, Judson MA. Relapses of sarcoidosis: what are they and can we predict who will get them? Eur Respir J. 2014; 43: 337–339. https://doi.org/10.1183/09031936.00138913 PMID: 24488991

21. Studdy PR, Bird R, Neville E, James DG. Biochemical findings in sarcoidosis. J Clin Pathol. 1980; 33: 528–533 PMID: 6249857

22. Siltzbach LE, James DG, Neville E, Turiaf J, Battesti JP, Sharma OP, et al. Course and prognosis of sarcoidosis around the world. Am J Med. 1974; 57: 847–852 PMID: 4432869

23. Inoue Y, Inui N, Hashimoto D, Enomoto N, Fujisawa T, Nakamura Y, et al. Cumulative Incidence and Predictors of Progression in Corticosteroid-Naïve Patients with Sarcoidosis. PloS One. 2015; 10: e0143371. https://doi.org/10.1371/journal.pone.0143371 PMID: 26575272

24. Ainslie GM, Benatar SR. Serum angiotensin converting enzyme in sarcoidosis: sensitivity and specific-
ity in diagnosis: correlations with disease activity, duration, extra-thoracic involvement, radiographic
type and therapy. Q J Med. 1985; 55: 253–270 PMID: 2991971

25. Baughman RP, Shipley R, Eisentrout CE. Predictive value of gallium scan, angiotensin-converting enzyme level, and bronchoalveolar lavage in two-year follow-up of pulmonary sarcoidosis. Lung. 1987; 165: 371–377 PMID: 2828772

26. Prasse A, Katic C, Germann M, Buchwald A, Zissel G, Müller-Quernheim J. Phenotyping sarcoidosis from a pulmonary perspective. Am J Respir Crit Care Med. 2008; 177: 330–336. https://doi.org/10.1164/rccm.200705-742OC PMID: 17975200

27. Baughman RP, Culver DA, Judson MA. A concise review of pulmonary sarcoidosis. Am J Respir Crit Care Med. 2011; 183: 573–581. https://doi.org/10.1164/rccm.201006-0865CI PMID: 21037016

28. Biller H, Zissel G, Ruprecht B, Nauck M, Busse Grawitz A, Müller-Quernheim J. Genotype-corrected reference values for serum angiotensin-converting enzyme. Eur Respir J. 2006; 28: 1085–1090. https://doi.org/10.1183/09031936.00050106 PMID: 17138677

29. Israel HL, Karlin P, Menduke H, DeLisser OG. Factors affecting outcome of sarcoidosis. Influence of race, extrathoracic involvement, and initial radiologic lung lesions. Ann N Y Acad Sci. 1986; 465: 609–618 PMID: 3460398

30. Judson MA, Boan AD, Lackland DT. The clinical course of sarcoidosis: presentation, diagnosis, and treatment in a large white and black cohort in the United States. Sarcoidosis Vasc Diffuse Lung Dis. 2012; 29: 119–127 PMID: 23461074

31. Baughman RP, Teirstein AS, Judson MA, Rossman MD, Yeager H, Bresnitz EA, et al. Clinical char-
acteristics of patients in a case control study of sarcoidosis. Am J Respir Crit Care Med. 2001; 164:
1885-1889. https://doi.org/10.1164/ajrccm.164.10.2104046 PMID: 11734441

32. Rabin DL, Thompson B, Brown KM, Judson MA, Huang X, Lackland DT, et al. Sarcoidosis: social predictors of severity at presentation. Eur Respir J. 2004; 24: 601–608. https://doi.org/10.1183/09031936.04.00070503 PMID: 15459139

33. Mirsaeidi M, Machado RF, Schraufnagel D, Sweiss NJ, Baughman RP. Racial difference in sarcoidosis mortality in the United States. Chest. 2015; 147: 438–449. https://doi.org/10.1378/chest.14-1120 PMID: 25188873

34. Tukey MH, Berman JS, Boggs DA, White LF, Rosenberg L, Cozier YC. Mortality among African Ameri-
can women with sarcoidosis: data from the Black Women's Health Study. Sarcoidosis Vasc Diffuse
Lung Dis. 2013; 30: 128–133 PMID: 24071884

35. Baughman RP, Judson MA, Teirstein A, Yeager H, Rossman M, Knatterud GL, et al. Presenting charac-
teristics as predictors of duration of treatment in sarcoidosis. QJM. 2006; 99: 307–315. https://doi.org/
10.1093/qjmed/hcl038 PMID: 16595563

36. Costabel U, Bonella F, Ohshimo S, Guzman J. Diagnostic modalities in sarcoidosis: BAL, EBUS, and PET. Semin Respir Crit Care Med. 2010; 31: 404–408. https://doi.org/10.1055/s-0030-1262207 PMID: 20665390

37. Winterbauer RH, Lammert J, Selland M, Wu R, Corley D, Springmeyer SC. Bronchoalveolar lavage cell populations in the diagnosis of sarcoidosis. Chest. 1993; 104: 352–361 PMID: 8339618

38. Eule H, Weinecke A, Roth I, Wuthe H. The possible influence of corticosteroid therapy on the natural course of pulmonary sarcoidosis. Late results of a continuing clinical study. Ann N Y Acad Sci. 1986; 465: 695–701 PMID: 3524366