# Familial vs. sporadic sarcoidosis: BTNL2 polymorphisms, clinical presentations, and outcomes in a French cohort

Yves Pacheco, Alain Calender, Dominique Israël-Biet, Pascal Roy, Serge Lebecque, Vincent Cottin, Diane Bouvry, Hilario Nunes, Pascal Sève, Laurent Pérard, et al.

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

## To cite this version:

Yves Pacheco, Alain Calender, Dominique Israël-Biet, Pascal Roy, Serge Lebecque, et al. Familial vs. sporadic sarcoidosis: BTNL2 polymorphisms, clinical presentations, and outcomes in a French cohort. Orphanet Journal of Rare Diseases, 2016, 11 (1), 10.1186/s13023-016-0546-4. doi: 10.1595465

HAL Id: hal-01595465

https://hal.science/hal-01595465

Submitted on 26 Sep 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
RESEARCH

Open Access

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

CrossMark

# Familial vs. sporadic sarcoidosis: BTNL2 polymorphisms, clinical presentations, and outcomes in a French cohort

Yves Pacheco<sup>1,11*</sup>, Alain Calender<sup>2</sup>, Dominique Israël-Biet<sup>3</sup>, Pascal Roy<sup>4</sup>, Serge Lebecque<sup>5</sup>, Vincent Cottin<sup>6</sup>, Diane Bouvry<sup>7</sup>, Hilario Nunes<sup>7</sup>, Pascal Sève<sup>8</sup>, Laurent Pérard<sup>9</sup>, Gilles Devouassoux<sup>8</sup>, Nathalie Freymond<sup>1</sup>, Chahira Khouatra<sup>6</sup>, Benoît Wallaert<sup>10</sup>, Raphaelle Lamy<sup>2</sup>, Mad-Hélénie Elsensohn<sup>4</sup>, Claire Bardel<sup>4</sup>, Dominique Valeyre<sup>7</sup> and GSF group

**Abstract**

**Background:** The occurrence of familial forms of sarcoidosis (OMIM 181100) suggests a genetic predisposition. The involvement of butyrophilin-like 2 (BTNL2) gene (rs2076530 variant) has to be investigated.

**Results:** The study performed independent analyses of BTNL2 polymorphism, clinical phenotypes, and outcomes in familial vs. sporadic presentations in 256 sporadic and 207 familial cases from 140 families. The logistic multivariate model showed that a young age at diagnosis and the combination of lung and skin involvement at diagnosis may distinguish sporadic from familial sarcoidosis (p = 0.016 and p = 0.041). We observed also that Sarcoid Clinical Activity Classification (SCAC) profiles were significantly different between familial and sporadic cases (p = 0.0497). Variant rs2076530 was more frequent in patients than in controls (OR = 2.02; 95% CI: [1.32–3.09]) but showed no difference between sporadic and familial cases and no difference according to the clinical phenotype or the outcome.

**Conclusion:** Despite a significant difference in BTNL2 polymorphism between sarcoid patients and controls, there was no such difference between familial and sporadic sarcoidosis cases and no correlation between BTNL2 polymorphism and disease severity or outcome. Thus, BTNL2 difference cannot be considered as a key marker for disease classification or patient management.

**Keywords:** Sarcoidosis, Medical genetics, Human BTNL2 protein, Candidate gene association study, Classification

**Background**

Sarcoidosis is a rare multisystemic granulomatous disorder of still unknown origin. It has various presentations, severities, treatments, and outcomes [1]. Genetic studies have searched for polymorphisms associated with the risk of developing the disease and for genotype-phenotype correlations. Among several candidate genes, BTNL2 (butyrophilin-like 2) has been intensively studied. A recent meta-analysis has confirmed that BTNL2 rs2076530 polymorphism contributes to the risk of sarcoidosis [2]. Meanwhile, it has been difficult to

investigate the clinical phenotypes and outcomes of sarcoidosis in search for genotype-phenotype correlations. Nevertheless, interesting suggestions about classical radiographic staging [3], factor analysis [4, 5], SCAC (Sarcoid Clinical Activity Classification) [6], and WASOG (World Association for Sarcoidosis and Other Granulomatous Disorders) clinical outcome status classification [7] are available.

We investigated BTNL2 gene in familial forms of sarcoidosis to assess the role of this gene as a major indicator of hereditary predisposition to the disease [8–10] and find out whether it can be a useful genetic marker for clinical management and prognosis.

<footnote>* Correspondence: yves.pacheco@univ-lyon1.fr</footnote>

<footnote n="1">¹Hospices Civils de Lyon, Centre Hospitalier Lyon Sud, Université Claude Bernard - Lyon 1, EA-7426, Lyon, France</footnote>

<footnote n="11">¹¹Université Claude Bernard Lyon 1 - EA-7426, 165 Chemin du Grand Revoyet, F-69495 Pierre Benite, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
# Methods

## Participant inclusion

SARCFAM (a national project on familial sarcoidosis) is a prospective observational cohort study started in 2008. It involves 28 French university departments of internal medicine or pulmonology that care for sarcoidosis patients.

Three inclusion criteria were: i) clinical and paraclinical features consistent with sarcoidosis; ii) histopathological evidence of non-caseating granuloma, except Löfgren syndrome; and, iii) exclusion of any other chronic disease.

The evidence for non-caseating granulomatous lesions was obtained by mediastinoscopy (26.4%), or bronchial/pulmonary (34.9%), accessory salivary gland (12.3%), peripheral lymph node (10.4%), skin (9.8%), conjunctiva (2.5%), kidney (1.2%) or other (2.5%) biopsies.

All 463 cases included were diagnosed according to the Joint Statement of the American Thoracic Society (ATS), the European Respiratory Society (ERS), and the (WASOG) [11]. The disease type (familial vs. sporadic) was determined by structured interview.

The control population for the genetic investigation consisted of 430 DNA samples from a single genetic database of healthy human reference population. The mean age (41.5 ± 17.5 years) and geographic origins [European: *n* = 324 (75.3%) and Sub-Saharan African plus Caribbean: *n* = 106 (24.7%)] were quite similar to those of the study sarcoid patients.

## Clinical assessment

All patients with sporadic sarcoidosis or index cases with familial sarcoidosis were diagnosed and followed-up in the 28 clinical centers. A specific case report form was filled in with demographic data, medical history, disease onset (symptomatic or asymptomatic), detailed organ involvement, treatments, and outcome. A detailed familial history of sarcoidosis was also examined to assess genetic predisposition. Familial cases were offered a DNA sampling and the related forms were filled in the center where the index case was diagnosed.

Clinical, biological, and imaging data were also collected at diagnosis and each follow-up visit until December 2012. These data included: i) chest X-ray staging (0: no involvement, 1: isolated hilar lymphadenopathy, 2: hilar lymphadenopathy with pulmonary infiltration, 3: isolated pulmonary infiltration without fibrosis, 4: pulmonary fibrosis); ii) pulmonary function testing: forced expiratory volume in one second (FEV1), forced vital capacity (FVC), FEV1/FVC ratio, and total lung capacity (TLC); iii) biological variables: serum calcium and creatinine, bronchoalveolar lavage cell count; and, iv) information about treatment: systemic treatment with corticosteroids or non-steroid immunosuppressive agents.

The disease outcome was examined using SCAC with six progression patterns [6] and, whenever possible, a classification of outcome in four categories: 1) recovery within 3 years; 2) recovery between 3 and 5 years; 3) no recovery at 5 years; 4) death.

## Genetic study

Blood DNAs were obtained from nearly all patients and tested for BTNL2 gene rs2076530 polymorphism. Primer oligonucleotides for polymerase chain reaction and DNA sequencing were located within the reference sequence surrounding rs2076530, the major single-nucleotide polymorphisms (SNPs) that induce truncation of the BTNL2 protein [8]. A 490-bp amplicon was produced and sequenced using primers 5'-AATGCACAGAGC ATGGAGGTGAG-3' and 5'-GAAGATACTGGAAAAG ATACAAG-3'.

The quality control of PCR products was performed on a LabChipGX system (PERKIN ELMER™) and the products purified by NucleoFast® 96 PCR Clean-up Kit from MACHEREY-NAGEL™. The sequencing was performed by Big Dye terminator v1.1 after purification with BigDye® XTerminator™ Purification Kit. Sequence delineation and base calling used an automated fluorescent DNA sequencer (Applied Biosystems™, model 3130xl).

## Statistical analysis

Familial and sporadic cases were compared using a generalized linear mixed model. To allow for family links, the model included correlations between patients' random effects [12].

In a first application of the model, univariate submodels were constructed and compared with a null model (model with intercept only) using a likelihood ratio test to test the effects of the main variables (Table 1). The clinical characteristics of BTNL2 G/G vs. BTNL2 A/A plus BTNL2 A/G genotypes were compared with the same method (Table 2). Because the study was mainly exploratory, only raw *p*-values are shown; however, a Bonferroni correction for multiple testing was applied and when raw *p*-values were significant, adjusted *p*-values were also calculated [13].

In another application, a multivariate model was built with the most important variables collected at diagnosis (age, sex, organs involved, chest X-ray stage, and BTNL2 polymorphism) to test their association with familial cases.

# Results

## The study population

Among the 463 patients, 256 presented with a sporadic disease and 207 belonged to 140 families (one patient from each of 85 families, two patients from each of 48 families, three patients from each of five families, four
<table><caption>**Table 1** Epidemiological, clinical, and biological characteristics of familial and sporadic sarcoid patients</caption><thead><tr><th>Data</th><th>Number</th><th>Familial cases</th><th>Sporadic cases</th><th>Univariate model raw p value*</th><th>Multivariate model p value</th></tr></thead><tbody><tr><td>Epidemiological data</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td>463</td><td>207 (44.7%)</td><td>256 (55.3%)</td><td></td><td></td></tr><tr><td>Sex</td><td>463</td><td></td><td></td><td>0.49</td><td></td></tr><tr><td>Male</td><td>201</td><td>86 (41.6%)</td><td>115 (44.9%)</td><td></td><td>0.20</td></tr><tr><td>Female</td><td>262</td><td>121 (58.4%)</td><td>141 (55.1%)</td><td></td><td>Ref.</td></tr><tr><td>Age at diagnosis</td><td>424</td><td></td><td></td><td>0.001 [0.021]</td><td>0.016</td></tr><tr><td>Mean ± SD</td><td></td><td>38.2 ± 13.1</td><td>43.0 ± 15.5</td><td></td><td></td></tr><tr><td>Median [q1; q3]</td><td></td><td>35.5 [29.3; 46.2]</td><td>40.3 [30.6; 54.0]</td><td></td><td></td></tr><tr><td>Ethnicity</td><td>303</td><td></td><td></td><td>0.85</td><td></td></tr><tr><td>European</td><td>235</td><td>112 (77.2%)</td><td>123 (77.9%)</td><td></td><td></td></tr><tr><td>Caribbean or Sub-Saharan African</td><td>68</td><td>33 (22.8%)</td><td>35 (22.1%)</td><td></td><td></td></tr><tr><td>Smokers</td><td>381</td><td></td><td></td><td>0.43</td><td></td></tr><tr><td>Yes</td><td>156</td><td>56 (38.4%)</td><td>100 (42.6%)</td><td></td><td></td></tr><tr><td>No</td><td>225</td><td>90 (61.6%)</td><td>135 (57.4%)</td><td></td><td></td></tr><tr><td>Clinical data at diagnosis</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Disease onset</td><td>401</td><td></td><td></td><td>0.015 [0.28]</td><td></td></tr><tr><td>Asymptomatic</td><td>137</td><td>45 (27.3%)</td><td>92 (39.0%)</td><td></td><td></td></tr><tr><td>Symptomatic</td><td>264</td><td>120 (72.7%)</td><td>144 (61.0%)</td><td></td><td></td></tr><tr><td>Löfgren syndrome</td><td>378</td><td></td><td></td><td>0.26</td><td></td></tr><tr><td>Yes</td><td>58</td><td>29 (17.8%)</td><td>29 (13.5%)</td><td></td><td></td></tr><tr><td>No</td><td>320</td><td>134 (82.2%)</td><td>186 (86.5%)</td><td></td><td></td></tr><tr><td>Number of organs involved</td><td>414</td><td></td><td></td><td>0.144</td><td></td></tr><tr><td>1</td><td>158</td><td>69 (39.9%)</td><td>89 (36.9%)</td><td></td><td></td></tr><tr><td>2</td><td>117</td><td>36 (20.8%)</td><td>81 (33.6%)</td><td></td><td></td></tr><tr><td>3</td><td>70</td><td>29 (16.8%)</td><td>41 (17.0%)</td><td></td><td></td></tr><tr><td>≥ 4</td><td>69</td><td>39 (22.5%)</td><td>30 (12.5%)</td><td></td><td></td></tr><tr><td>Organs involved</td><td>414</td><td></td><td></td><td>&lt;0.001 [0.003]</td><td></td></tr><tr><td>Lungs only</td><td>148</td><td>68 (39.31%)</td><td>80 (33.2%)</td><td></td><td>0.65</td></tr><tr><td>Lungs and Eyes</td><td>53</td><td>10 (5.78%)</td><td>43 (17.84%)</td><td></td><td>0.076</td></tr><tr><td>Lungs and Skin</td><td>21</td><td>15 (8.67%)</td><td>6 (2.49%)</td><td></td><td>0.041</td></tr><tr><td>Other organs</td><td>192</td><td>80 (46.24%)</td><td>112 (46.47%)</td><td></td><td>Ref</td></tr><tr><td>X-ray staging</td><td>456</td><td></td><td></td><td>0.20</td><td>0.78</td></tr><tr><td>0</td><td>68</td><td>40 (19.8%)</td><td>28 (11.0%)</td><td></td><td></td></tr><tr><td>1</td><td>162</td><td>66 (32.7%)</td><td>96 (37.8%)</td><td></td><td></td></tr><tr><td>2</td><td>157</td><td>69 (34.2%)</td><td>88 (34.7%)</td><td></td><td></td></tr><tr><td>3</td><td>48</td><td>15 (7.4%)</td><td>33 (13.0%)</td><td></td><td></td></tr><tr><td>4</td><td>21</td><td>12 (5.9%)</td><td>9 (3.5%)</td><td></td><td></td></tr><tr><td>PFT at first visit</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>TLC (%)</td><td>215</td><td></td><td></td><td>0.93</td><td></td></tr><tr><td>&lt; 80%</td><td>47</td><td>18 (22.2%)</td><td>29 (21.6%)</td><td></td><td></td></tr><tr><td>≥ 80%</td><td>168</td><td>63 (77.8%)</td><td>105 (78.4%)</td><td></td><td></td></tr><tr><td>FEV1/FVC</td><td>242</td><td></td><td></td><td>0.43</td><td></td></tr><tr><td>&lt; 70%</td><td>55</td><td>18 (20.0%)</td><td>37 (24.3%)</td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 1** Epidemiological, clinical, and biological characteristics of familial and sporadic sarcoid patients (Continued)</caption><tbody><tr><td>≥ 70%</td><td>187</td><td>72 (80.0%)</td><td>115 (75.7%)</td><td></td></tr><tr><td>Biology at first visit</td><td></td><td></td><td></td><td></td></tr><tr><td>Bronchoalveolar lavage lymphocytes (%)</td><td>95</td><td></td><td></td><td>0.171</td></tr><tr><td>&lt; 10</td><td>23</td><td>3 (11.6%)</td><td>20 (29.0%)</td><td></td></tr><tr><td>[10; 20]</td><td>17</td><td>5 (19.2%)</td><td>12 (17.4%)</td><td></td></tr><tr><td>&gt; 20</td><td>55</td><td>18 (69.2%)</td><td>37 (53.6%)</td><td></td></tr><tr><td>Calcemia</td><td>256</td><td></td><td></td><td>0.80</td></tr><tr><td>&lt; 2.25</td><td>52</td><td>21 (22.6%)</td><td>31 (19.0%)</td><td></td></tr><tr><td>[2.25 - 2.5]</td><td>182</td><td>64 (68.8%)</td><td>118 (72.4%)</td><td></td></tr><tr><td>&gt; 2.5</td><td>22</td><td>8 (8.6%)</td><td>14 (8.6%)</td><td></td></tr><tr><td>Treatment</td><td></td><td></td><td></td><td></td></tr><tr><td>Treatment</td><td>463</td><td></td><td></td><td>&lt;0.001 [0.008]</td></tr><tr><td>None</td><td>136</td><td>79 (38.2%)</td><td>57 (22.3%)</td><td></td></tr><tr><td>At least one</td><td>327</td><td>128 (61.8%)</td><td>199 (77.7%)</td><td></td></tr><tr><td>Type of treatment</td><td>327</td><td></td><td></td><td>0.188</td></tr><tr><td>Corticosteroids only</td><td>186</td><td>67 (52.3%)</td><td>119 (59.8%)</td><td></td></tr><tr><td>Other treatment</td><td>141</td><td>61 (47.7%)</td><td>80 (40.2%)</td><td></td></tr><tr><td>Pattern of disease course</td><td></td><td></td><td></td><td></td></tr><tr><td>SCAC</td><td>414</td><td></td><td></td><td>0.003 [0.0497]</td></tr><tr><td>1- symptomatic - no treatment</td><td>48</td><td>30 (18.1%)</td><td>18 (7.3%)</td><td></td></tr><tr><td>2- symptomatic - treatment ≤ 12 mo.</td><td>47</td><td>18 (10.8%)</td><td>29 (11.7%)</td><td></td></tr><tr><td>3- symptomatic - treatment > 12 mo.</td><td>172</td><td>72 (43.4%)</td><td>100 (40.3%)</td><td></td></tr><tr><td>4- asymptomatic - no treatment</td><td>51</td><td>19 (11.5%)</td><td>32 (12.9%)</td><td></td></tr><tr><td>5- asymptomatic - treatment ≤ 12 mo.</td><td>16</td><td>7 (4.2%)</td><td>9 (3.6%)</td><td></td></tr><tr><td>6- asymptomatic - treatment > 12 mo.</td><td>80</td><td>20 (12.0%)</td><td>60 (24.2%)</td><td></td></tr><tr><td>Classification of outcome</td><td>138</td><td></td><td></td><td>0.165**</td></tr><tr><td>1- recovery within 3 years</td><td>72</td><td>24 (57.14%)</td><td>48 (50%)</td><td></td></tr><tr><td>2- recovery between 3 and 5 years</td><td>32</td><td>6 (14.29%)</td><td>26 (27.08%)</td><td></td></tr><tr><td>3- no recovery at 5 years</td><td>30</td><td>12 (28.57%)</td><td>18 (18.75%)</td><td></td></tr><tr><td>4- death</td><td>4</td><td>0 (0%)</td><td>4 (4.17%)</td><td></td></tr><tr><td>Genetic data</td><td></td><td></td><td></td><td></td></tr><tr><td>BTNL2</td><td>436</td><td></td><td></td><td>0.17***</td></tr><tr><td>BTNL2 A/A</td><td>201</td><td>89 (47.8%)</td><td>112 (44.8%)</td><td>0.39***</td></tr><tr><td>BTNL2 A/G</td><td>194</td><td>84 (45.2%)</td><td>110 (44.0%)</td><td></td></tr><tr><td>BTNL2 G/G</td><td>41</td><td>13 (7.0%)</td><td>28 (11.2%)</td><td></td></tr></tbody></table>

patients from one family, and seven patients from another one). Most of the 303 patients with known ancestors' birth countries were European (78%); the others were Caribbean or Sub-Saharan African (22%). This distribution did not differ from that of the controls. Table 1 shows a comparison of patient epidemiological, clinical, and biological characteristics between sporadic and familial cases.

### Impact of disease type

Univariate comparisons between sporadic and familial cases revealed no significant differences in terms of sex ratio, geographical origin, number of involved organs at diagnosis, distribution of chest X-ray stages at diagnosis, respiratory function tests, biological data, or smoking habits. Also, the incidence of Löfgren syndrome was not different between familial and sporadic cases (17.8% vs.

<footnote>*Significant raw p-values are followed by corrected p-values (Bonferroni correction) in brackets</footnote>

<footnote>**Small "Death" category excluded</footnote>

<footnote>***G/G vs. G/A + A/A</footnote>
<table><caption>**Table 2** Comparison of clinical characteristics according to the BTLN2 genotype (G/G vs. A/A plus G/A)</caption><thead><tr><th rowspan="2">Data</th><th rowspan="2">Number</th><th colspan="3">BTNL2</th><th rowspan="2">Univariate model raw p value *,**</th></tr><tr><th>A/A</th><th>G/A</th><th>G/G</th></tr></thead><tbody><tr><td>Epidemiological data</td><td></td><td></td><td></td><td></td><td>0.59</td></tr><tr><td>Ethnicity</td><td>292</td><td></td><td></td><td></td><td></td></tr><tr><td>European</td><td>226</td><td>103 (75.7%)</td><td>104 (80%)</td><td>19 (73.1%)</td><td></td></tr><tr><td>Caribbean or Sub-Saharan African</td><td>66</td><td>33 (24.3%)</td><td>26 (20%)</td><td>7 (26.9%)</td><td></td></tr><tr><td>Clinical data at diagnosis</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Disease onset</td><td>389</td><td></td><td></td><td></td><td>0.91</td></tr><tr><td>Asymptomatic</td><td>130</td><td>60 (34.29%)</td><td>59 (33.15%)</td><td>11 (30.56%)</td><td></td></tr><tr><td>Symptomatic</td><td>259</td><td>115 (65.71%)</td><td>119 (66.85%)</td><td>25 (69.44%)</td><td></td></tr><tr><td>Löfgren syndrome</td><td>368</td><td></td><td></td><td></td><td>0.042 [0.42]</td></tr><tr><td>Yes</td><td>58</td><td>25 (14.71%)</td><td>32 (19.39%)</td><td>1 (3.03%)</td><td></td></tr><tr><td>No</td><td>310</td><td>145 (85.29%)</td><td>133 (80.61%)</td><td>32 (96.97%)</td><td></td></tr><tr><td>Number of organs involved</td><td>397</td><td></td><td></td><td></td><td>0.074</td></tr><tr><td>1</td><td>152</td><td>70 (38.46%)</td><td>73 (41.01%)</td><td>9 (24.32%)</td><td></td></tr><tr><td>2</td><td>115</td><td>55 (30.22%)</td><td>49 (27.53%)</td><td>11 (29.73%)</td><td></td></tr><tr><td>3</td><td>66</td><td>29 (15.93%)</td><td>33 (18.54%)</td><td>4 (10.81%)</td><td></td></tr><tr><td>≥ 4</td><td>64</td><td>28 (15.39%)</td><td>23 (12.92%)</td><td>13 (35.14%)</td><td></td></tr><tr><td>Organs involved</td><td>397</td><td></td><td></td><td></td><td>0.082</td></tr><tr><td>Lungs only</td><td>143</td><td>68 (37.36%)</td><td>67 (37.64%)</td><td>8 (21.62%)</td><td></td></tr><tr><td>Lungs and Eyes</td><td>52</td><td>23 (12.64%)</td><td>20 (11.24%)</td><td>9 (24.33%)</td><td></td></tr><tr><td>Lungs and Skin</td><td>20</td><td>11 (6.04%)</td><td>8 (4.49%)</td><td>1 (2.7%)</td><td></td></tr><tr><td>Other sites</td><td>182</td><td>80 (43.96%)</td><td>83 (46.63%)</td><td>19 (51.35%)</td><td></td></tr><tr><td>X-ray staging</td><td>432</td><td></td><td></td><td></td><td>0.26</td></tr><tr><td>0</td><td>57</td><td>25 (12.5%)</td><td>26 (13.47%)</td><td>6 (15.38%)</td><td></td></tr><tr><td>1</td><td>158</td><td>77 (38.5%)</td><td>65 (33.68%)</td><td>16 (41.03%)</td><td></td></tr><tr><td>2</td><td>149</td><td>67 (33.5%)</td><td>68 (35.23%)</td><td>14 (35.9%)</td><td></td></tr><tr><td>3</td><td>47</td><td>24 (12%)</td><td>22 (11.4%)</td><td>1 (2.56%)</td><td></td></tr><tr><td>4</td><td>21</td><td>7 (3.5%)</td><td>12 (6.22%)</td><td>2 (5.13%)</td><td></td></tr><tr><td>PFT at first visit</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>TLC (%)</td><td></td><td></td><td></td><td></td><td>0.098</td></tr><tr><td>&lt; 80%</td><td>47</td><td>23 (24.47%)</td><td>22 (23.91%)</td><td>2 (9.52%)</td><td></td></tr><tr><td>≥ 80%</td><td>160</td><td>71 (75.53%)</td><td>70 (76.09%)</td><td>19 (90.48%)</td><td></td></tr><tr><td>Pattern of disease course</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SCAC</td><td>398</td><td></td><td></td><td></td><td>0.83</td></tr><tr><td>1- symptomatic - no treatment</td><td>46</td><td>21 (11.6%)</td><td>22 (12.1%)</td><td>3 (8.3%)</td><td></td></tr><tr><td>2- symptomatic - treatment ≤ 12 mo.</td><td>47</td><td>19 (10.5%)</td><td>22 (12.1%)</td><td>6 (16.7%)</td><td></td></tr><tr><td>3- symptomatic - treatment > 12 mo.</td><td>166</td><td>75 (41.4%)</td><td>74 (40.9%)</td><td>17 (47.2%)</td><td></td></tr><tr><td>4- asymptomatic - no treatment</td><td>48</td><td>19 (10.5%)</td><td>26 (14.4%)</td><td>3 (8.3%)</td><td></td></tr><tr><td>5- asymptomatic - treatment ≤ 12 mo.</td><td>14</td><td>6 (3.3%)</td><td>7 (3.9%)</td><td>1 (2.8%)</td><td></td></tr><tr><td>6- asymptomatic - treatment > 12 mo.</td><td>77</td><td>41 (22.7%)</td><td>30 (16.6%)</td><td>6 (16.7%)</td><td></td></tr><tr><td>Classification of outcome</td><td>137</td><td></td><td></td><td></td><td>0.015 [0.14]***</td></tr><tr><td>1- recovery within 3 years</td><td>72</td><td>29 (52.73%)</td><td>36 (52.17%)</td><td>7 (53.85%)</td><td></td></tr><tr><td>2- recovery between 3 and 5 years</td><td>32</td><td>9 (16.36%)</td><td>17 (24.64%)</td><td>6 (46.15%)</td><td></td></tr><tr><td>3- no recovery at 5 years</td><td>29</td><td>15 (27.27%)</td><td>14 (20.29%)</td><td>0 (0%)</td><td></td></tr><tr><td>4- death</td><td>4</td><td>2 (3.64%)</td><td>2 (2.9%)</td><td>0 (0%)</td><td></td></tr></tbody></table>

<footnote>*Significant raw p-values are followed by corrected p-values (Bonferroni correction) in brackets</footnote>

<footnote>**Small "Death" category excluded</footnote>

<footnote>***G/G vs. G/A + A/A</footnote>
13.5%, respectively). However, the age at diagnosis was lower in familial cases (p = 0.021 after Bonferroni correction) and the combinations of organs involved at diagnosis were significantly different (p = 0.003 after Bonferroni correction).

Sporadic cases needed treatment more often than familial cases (p = 0.008 after Bonferroni correction). Among treated patients, 52.3% of familial and 59.8% of sporadic cases received exclusively corticosteroids. Other treatments (methotrexate, azathioprine, aminoquinolines, or anti-TNF alpha) were given to 47.7% of familial and 40.2% of sporadic cases. The SCAC was significantly different between familial and sporadic cases (p = 0.0497 after Bonferroni correction) but the classification of outcome did not show a significant difference.

A multivariate model was then constructed with the following variables assessed at diagnosis: age, sex, chest X-ray stage, organs involved, and BTNL2 polymorphism. The model was run on the 390 patients without missing data (233 sporadic and 157 familial cases, the former being the reference). The OR estimates and their 95% CIs were 0.81 [0.69; 0.96] for age (per ten-year increment) and 0.73 [0.45; 1.18] for sex (males vs. females). Considering an additive effect for ordinal variables and a dominant genetic model for BTNL2 polymorphism (GG being the wild-type homozygous genotype), the OR estimates were: 1.04 [0.80; 1.34] for chest X-ray stage (per stage increment), 1.46 [0.62; 3.45] for BTNL2 (G/G vs. G/A plus A/A), 1.13 [0.67; 1.88] for organs involved “lungs only”, 0.47 [0.21; 1.08] for “lungs and eyes”, 3.16 [1.05; 9.54] for “lungs and skin” (“other organs” being the reference). In agreement with the univariate analyses, familial cases were significantly younger at diagnosis than sporadic ones (p = 0.016) and the organs involved in familial cases were more often “lungs and skin” (p = 0.041).

## Impact of BTNL2

In the control population (n = 430), the genotypes observed for SNP rs2076530 (normal allele vs. mutated variant A) were respectively 75 (17.4%) G/G, 143 (33.2%) G/A, and 212 (49.3%) A/A. The relative distribution of genotype G/G vs. G/A plus A/A was significantly different between sarcoidosis cases and the control population (odds ratio, OR = 2.02; 95% confidence interval, CI: [1.32–3.09]). The same comparison performed between familial and sporadic cases showed no statistically significant difference.

Table 2 shows that, in the whole population of patients, there is no correlation between BTNL2 polymorphism and the various clinical phenotypes, the SCAC, or the classification of outcome.

## Discussion

The present study was based on clinical and genetic investigations in a nation-wide cohort of both familial and

sporadic sarcoidosis cases. Comparisons between familial and sporadic cases showed some clinical differences but, mostly, the absence of difference in BTNL2 rs2076530 variant frequency. Moreover, the study found no genotype-phenotype correlation between BTNL2 and clinical phenotype, SCAC, or the classification of outcome.

Familial aggregation of sarcoidosis cases is already well known [14]; it raised the issues of predisposing genes and genotypes involved in a Mendelian inheritance of the clinical phenotype. Along with ACCESS [15], we assessed the clinical phenotypes and the outcomes in a large cohort of sarcoid patients comparing sporadic and familial cases. In agreement with the literature [16–19], we observed no difference in sex ratio or geographic origin between familial and sporadic cases. However, in contrast with other reports [16, 18, 19], age at diagnosis was significantly lower in the present study familial cases. This observation suggests that diagnosing sarcoidosis in individuals closely related to index cases in predisposed families may be made earlier than in sporadic cases. A previous work [17] compared the clinical profiles of patients whose diagnosis was made incidentally with symptomatic cases but could not conclude because of small sample sizes. In agreement with Sharma et al. [6], we found a difference between familial and sporadic cases regarding the disease onset (symptomatic vs. asymptomatic); this difference was statistically significant as per the raw *p* value but not significant after Bonferroni correction.

In addition, in agreement with other studies [16, 17], we found no differences in terms of chest X-ray staging or number of organs involved at diagnosis between familial and sporadic cases but the combinations of organs involved at diagnosis was significantly different between the two groups in univariate analyses. Also, the multivariate analysis showed a significant difference between familial and sporadic cases regarding lungs and skin combination. According to Delaveri et al. [20] such an observation may be in relation with the genetic predisposing background. Finally, as also reported by Sharma et al. [16], familial and sporadic cases did not differ in terms of drugs used (corticosteroids alone vs. immunosuppressive drugs).

Several investigators have designed comprehensive scoring systems to assess the disease impact and outcome [3–8]. Two of them have integrated the need for more than a 12-month treatment period [4, 6]. In the present study, we found a significant difference between familial and sporadic cases in terms of disease progression patterns (as per the SCAC) [6] but only a trend toward a difference in the classification of outcome classification according to the duration of the disease. This may reflect a lack of power regarding the classification of outcome. Indeed, only 138 patients had a sufficient follow-up to be correctly classified, whereas 344 patients
could be classified with SCAC. A retrospective study of patients with an isolated lung involvement and 6–144 months of follow-up [5] has identified four disease phenotypes. We cannot easily refer to this series because of the predominance, in our series, of a multi-organ involvement.

The BTNL2-related OR we found here is similar to previously reported ones [2, 8–10]. However, our study adds the original findings that although BTNL2 polymorphism was strongly associated with the disease, its presence was not predictive of the development of a rather familial or rather sporadic case and not associated with the clinical phenotypes, the SCAC, or the outcome. This result differs from others that suggested that BTNL2 might be a predisposing factor for persistent or progressive sarcoidosis [10, 21, 22]. Our data may be interpreted in the light of the low OR of BTNL2 rs2076530 (=2) we found; that is, the BTNL2 gene is not an explanatory factor of disease inheritance. The role of the rs2076530 variant is probably part of a set of genetic and environmental factors associated with the risk of acquiring the disease.

Combined effects of HLA-DRB1/DBP1 haplotypes and the BTNL2 variant have been previously reported [23, 24]. Here, we cannot exclude that the BTNL2-sarcoidosis association we observed is the consequence of gene proximity in the 6p21 area with a synergistic effect of some HLA class II haplotypes and the truncated form of BTNL2 on T-cell activation and proliferation process [22, 25]. Actually, one limitation of the present study is the lack of extensive screening for HLA-DRB1/DBP1 in the patients. In familial cases, we have checked for SNPs within the HLA-DRB1/DBP1 genes in 9 unrelated index cases of familial sarcoidosis and found no common genotypes; this suggests that these index cases do not have a common ancestor and that there is no founder effect. We have also screened the full sequence of BTNL2 in a subset of 40 patients from 40 families and found no other variants than those described in the last version of the NCBI dbSNP database. Actually, we analysed the BTNL2 gene for other variants and found that the rs2076530 A/A variant genotype cosegregates with the rs2076520 synonymous V313V, rs28362679 missense S334L, and rs41441651 missense D336N polymorphisms and, to date, we cannot exclude that these polymorphisms may have a complementary role with that of rs2076530 in the malfunction of BTNL2. Taken together these data demonstrate clearly that, with an OR = 2.0, the rs2076530 splicing variant of BTNL2 should be considered as a genetic risk factor for sarcoidosis but cannot be considered as a major gene explaining a Mendelian inheritance in the familial form of the disease. As others, we suggest strongly that BTNL2 screening belongs to the panel of biomarkers in the diagnosis of sarcoidosis together with HLA class II haplotyping and further coming genes. Despite conflicting opinions, the 6p21.3 region carrying HLA and BTNL2 genes is still

considered to be associated with an increased risk of TH1/TH17 diseases (tuberculosis, leprosy, beryllium disease, ulcerative colitis, or Crohn disease) [26–30]. Here, the BTNL2 rs2076530 splice variant could not distinguish sporadic from familial forms and we still have to understand pathways through which the truncated form of BTNL2 may be involved in the pathogenesis of sarcoidosis.

Sarcoidosis seems to be a multifactorial disease that involves various genetic and environmental pathogenic factors. To date, strong roles are suggested for different SNPs in various genes mostly identified by association studies: HLA-DP subregion, BTNL2, Annexin A11, and, more recently Toll-like receptors, coiled-coil domain-containing protein 88B, Ataxin-2/SH2B adapter protein 3, interleukin IL12B, and Beta-mannosidase/Nuclear factor NF-kappa-B p105 subunit [8, 31–36]. More than ten genes might then be involved in the predisposition to sarcoidosis and next generation sequencing studies will probably provide new candidates for the molecular pathways of granuloma generation. Finally, sarcoidosis might not be a single disease entity but rather a reaction to various triggering events.

Because the present study data were collected from both prevalent and incident cases (to increase the number of familial cases and recoverable DNAs), two biases may be considered. Firstly, some data were missing because of their retrospective collection and others because the 28 centres might not have given the same number of visits or the same delays between visits to all patients. Secondly, ancestry data relied on the ancestors' birth countries and several patients had ancestors from distinct continents. However, BTNL2 polymorphism was shown linked to sarcoidosis regardless of ethnicity in Caucasian, Afro-American, and Japanese patients [37]. Lastly, we did not use the clinical outcome status (COS) of the WASOG Task Force [7] because the COS became available four years after the start of the present study.

## Conclusions

This study is an original comparison of sarcoidosis phenotypes and genotypes between familial and sporadic cases in a very large cohort. Despite a significant difference in BTNL2 polymorphism between sarcoid patients and controls, the lack of correlations between polymorphism, inheritance, clinical phenotypes, and outcome argues against the consideration of this single genetic difference as a practical marker for patient classification or optimized individual patient management.

## Abbreviations

BTNL2: Butyrophilin-like 2 (gene); CI: Confidence interval; FEV1: Forced expiratory volume in one second; FVC: Forced vital capacity; OR: Odds ratio; SARCFAM: Sarcoidose Familiale (study); SCAC: Sarcoid Clinical Activity Classification; SNP: Single-nucleotide polymorphism; TLC: Total lung capacity; WASOG: World Association for Sarcoidosis and Other Granulomatous Disorders
## Acknowledgements

The authors thank Juthier L, Ouaglal H, and Castaldi MH for their technical assistance and Jean Iwaz (Hospices Civils de Lyon) for the revision of the final versions of this manuscript.

The authors also thank the members of the GSF group: Gouita A (CHU Marseille Nord), Chanez P (CHU Marseille), Peyrat-Detis E and Prevot G (CHU Toulouse), Diab S (CHU Grenoble), Prieur F (CHU ST Etienne), Locatelli-Sanchez M, Gaillot-Drevon M, Kiakouama L, Broussolle C, Varron L, Barron A, Kuntz J, Ledoux K, Le Scanff MC, and Bouaziz (CHU Lyon Sud), Pavic M (CH Desgenettes Lyon), Ninet J and Cordier JF (CHU Lyon Est), Vincent M (CH ST Joseph-St Luc, Lyon), Coudurier M (CH Chamberby), GAGNADOUX F (CHU Angers), Jounieaux F (CHU Lille), Terce G (CHU Lille), Diot P (CHU Tours), Marchand-Adams S (CHU Tours), Morel H (CH ST Malo), Jouveau S (CHU Rennes), Olivier-Faivre L (CHU Dijon), Lacronique J (APHP, Hôpital Cochin), Juvin K (APHP, Hôpital Européen Georges Pompidou), Crestani B, Aubier M, Borie R, Roussel R and Chapelon Abric C (APHP, Hôpital Bichat), Hervier B, Piette JC, Janssens T, and Saadoun D (APHP, CH Pitié Salpêtrière), Nathan N and Clément A (APHP, CH Trousseau), Ducroix JP (CHU Amiens), Guillaumat C and Watkin E (CH Sud Francilien), Freynet M, Naccache JM and Uzunhan Y (APHP, Hôpital Avicenne), and Limal N (APHP, Hôpital Henri Mondor).

## Funding

The initial project ("Étude Clinique et Génétique des Formes Familiales de Sarcoïdose") was supported by the French Ministère des Affaires Sociales et de la Santé" (Projet Hospitalier de Recherche Clinique National: PHRC national 2008 N° SP/AO 12/D50604; PHRC national 2012 N°12-027-0309) and by the Comité Français des Maladies Respiratoires (N° DG2013-15KM). The project promoter is Hospices Civils de Lyon.

The funding bodies did not intervene in the study design, data collection and analysis, result interpretation, or drafting the manuscript.

## Availability of data and materials

The data are available on request from the corresponding author.

## Authors' contributions

Study design and monitoring: YP, AC, DI, PR, SL, and DV. Data collection: YP, AC, DI, VC, DB, HN, PS, LP, GD, NF, CK, BW, and DV. Genetic laboratory analyses: AC, SL, RL. Data analysis: PR, ME, and CB. Result interpretation and manuscript drafting: YP, AC, DI, PR, SL, VC, DB, HN, PS, LP, GD, NF, CK, B W, ME, CB, and DV. Approval of final manuscript: All authors.

## Competing interests

The authors declare that they have no competing interests.

## Consent for publication

Not applicable.

## Ethics approval and consent to participate

SARCFAM study obtained Agreement N° 914377 from the Commission Nationale de l'Informatique et des Libertés to collect clinical and genetic information. All included patients gave written informed consents for information collection and use. The genetic study on BTNL2 was performed in accordance with the recommendations of the French bioethic laws.

## Author details

¹Hospices Civils de Lyon, Centre Hospitalier Lyon Sud, Université Claude Bernard - Lyon 1, EA7426, Lyon, France. ²Hospices Civils de Lyon, Hôpital Edouard Herriot, Plateforme de génétique moléculaire, CR-21076, Université Claude Bernard - Lyon 1, INSERM U1052, Lyon, France. ³AP-HP, Hôpital Européen Georges Pompidou, Service de Pneumologie, Centre de Compétence des maladies pulmonaires rares, Université René Descartes - Paris 5, Paris, France. ⁴Hospices Civils de Lyon, Service de Biostatistique, Université Claude Bernard - Lyon 1, CNRS UMR 5558, Lyon, France. ⁵Hospices Civils de Lyon, Hôpital Lyon-Sud, Université Claude Bernard - Lyon 1, INSERM U1052 - CNRS UMR5286, Lyon, France. ⁶Hospices Civils de Lyon, Hôpital Louis Pradel, Université Claude Bernard - Lyon 1, UMR 754, Lyon, France.

⁷AP-HP, Hôpital Avicenne, Université Paris13, EA2363, COMUE Sorbonne Paris Cité, Paris, France. ⁸Hospices Civils de Lyon, Hôpital Croix-Rousse, Université Claude Bernard - Lyon 1, Lyon, France. ⁹Hospices Civils de Lyon, Hôpital Edouard Herriot, Université Claude Bernard - Lyon 1, Lyon, France. ¹⁰Centre hospitalo-universitaire de Lille, Service de Pneumologie et

Immunoallergologie, Centre de Compétence Maladies Pulmonaires Rares, Université Lille 2, Lille, France. ¹¹Université Claude Bernard Lyon 1 - EA-7426, 165 Chemin du Grand Revoyet, F-69495 Pierre Benite, France.

Received: 1 June 2016 Accepted: 21 November 2016

Published online: 03 December 2016

## References

1. Valeyre D, Prasse A, Nunes H, Uzunhan Y, Brillet PY, Müller-Quernheim J. Sarcoidosis. Lancet. 2014;383:1155-67.

2. Lin Y, Wei J, Fan L, Cheng D. BTNL2 gene polymorphism and sarcoidosis susceptibility: a meta-analysis. PLoS ONE. 2015;10:e0122639.

3. Scadding JG. Prognosis of intrathoracic sarcoidosis in England. A review of 136 cases after five years' observation. BMJ. 1961;2:1165-72.

4. Wasfi YS, Rose CS, Murphy JR, et al. A new tool to assess sarcoidosis severity. Chest. 2006;129:1234-45.

5. Rodrigues SC, Rocha NA, Lima MS, Arakaki JS, Coletta EN, Ferreira RG, et al. Factor analysis of sarcoidosis phenotypes at two referral centers in Brazil. Sarcoidosis Vasc Diffuse Lung Dis. 2011;28:34-43.

6. 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-6.

7. Baughman RP, Nagai S, Balter M, et al. Defining the clinical outcome status (COS) in sarcoidosis: results of WASOG Task Force. Sarcoidosis Vasc Diffuse Lung Dis. 2011;28:56-64.

8. Valentonyte R, Hampe J, Huse K, Rosenstiel P, Albrecht M, Stenzel A, et al. Sarcoidosis is associated with a truncating splice site mutation in BTNL2. Nat Genet. 2005;37:57-64.

9. Rybicki BA, Walewski JL, Maliarik MJ, Kian H, Iannuzzi MC, ACCESS Research Group. The BTNL2 gene and sarcoidosis susceptibility in African and Whites. Am J Human Genet. 2005;77:491-9.

10. Li Y, Wollnik B, Pabst S, Lennarz M, Rohmann E, Gillissen A, et al. BTNL2 gene variant and sarcoidosis. Thorax. 2006;61:273-4.

11. American Thoracic Society, European Respiratory Society, World association of Sarcoidosis and Other Granulomatous Disorders. Statement on sarcoidosis: joint statement of the American Thoracic Society (ATS), European Respiratory Society (ERS), 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-55.

12. Vazquez AJ, Bates DM, Rosa GJ, Gianola D, Weigel KA. Technical note: An R package for fitting generalized linear mixed models in animal breeding. J Anim Sci. 2010;88:497-504.

13. Shaffer JP. Multiple Hypothesis Testing. Ann Rev Psychol. 1995;46:561-84.

14. Rybicki BA, Iannuzzi MC, Frederick MM, Thompson BW, Rossman MD, Bresnitz EA, et al. ACCESS Research Group. Familial aggregation of sarcoidosis. A case-control etiologic study of sarcoidosis (ACCESS). Am J Respir Crit Care Med. 2001;164:2085-91.

15. Rossman MD, Kreider ME. Lesson learned from ACCESS. Proc Am Thorac Soc. 2007;4:453-6.

16. Sharma OP, Neville E, Walker AN, James DG. Familial sarcoidosis: a possible genetic influence. Ann N Y Acad Sci. 1976;278:386-400.

17. Neville E, Walker AN, James DG. Prognostic factors predicting the outcome of sarcoidosis: an analysis of 818 patients. Q J Med. 1983;52:525-33.

18. Brennan NJ, Crean P, Long JP, Fitzgerald MX. High prevalence of familial sarcoidosis in an Irish population. Thorax. 1984;39:14-8.

19. Nassif X, Valeyre D, Loiseau A, Battesti JP. Familial sarcoidosis. Apropos of 22 families. Ann Med Interne (Paris). 1985;136:611-4.

20. Delaveri A, Rapti A, Poulou M, Fylaktou E, Tsipi M, Roussos C, et al. BTNL2 gene SNPs as a contributing factor to sarcoidosis pathogenesis in a cohort of Greek patients. Meta Gene. 2014;2:619-30.

21. Wijnen PA, Voorter CE, Nelemans PJ, Verschakelen JA, Bekers O, Brent M. Butyrophilin-like 2 in pulmonary sarcoidosis: a factor for susceptibility and progression. Hum Immunol. 2011;72:342-7.

22. Wennerström A, Pietinalho A, Lasota J, Salli K, Surakka I, Seppänen M, et al. Major histocompatibility complex class II and BTNL2 associations in sarcoidosis. Eur Respir J. 2013;42:550-3.

23. Grunewald J, Brynedal B, Darlington P, Nisell M, Cederlund K, Hillert J, et al. Different HLA-DRB1 allele distributions in distinct clinical subgroups of sarcoidosis patients. Respir Res. 2010;11:25.
24. Sato H, Grutters JC, Pantelidis P, Mizzon AN, Ahmad T, Van Houte AJ, et al. HLA-DQB1 0201. A marker for good prognosis in British and Dutch patients with sarcoidosis. Am J Respir Cell Mol Biol. 2002;27:406–12.

25. Sato H, Woodhead FA, Ahmad T, Grutters JC, Spagnolo P, van den Bosch JM, et al. Sarcoidosis HLA class II genotyping distinguishes differences of clinical phenotype across ethnic groups. Hum Mol Genet. 2010;19:4100–11.

26. Johnson CM, Traherne JA, Jamieson SE, Tremelling M, Bingham S, Parkes M, et al. Analysis of the BTNL2 truncating splice siter mutation in tuberculosis, leprosy and Crohn's disease. Tissue Antigens. 2007;69:236–41.

27. Mochida A, Kinouchi Y, Negoro K, Takahashi S, Takagi S, Nomura E, et al. Butyrophilin-like 2 gene is associated with ulcerative colitis in the Japanese under strong linkage disequilibrium with HLA-DRB1 1502. Tissue Antigens. 2007;70:128–35.

28. Sato H, Spagnolo P, Silveira L, Welsh KJ, du Bois RM, Newman LS, et al. BTNL2 allele associations with chronic beryllium disease in HLA-DPB1 GLU69-negative individuals. Tissue Antigens. 2007;70:480–6.

29. Ali S, Chopra R, Aggarwal S, Srivastava AK, Kalaiarasan P, Malhotra D, et al. Association of variants in, BAT1-LTA-TNF-BTNL2 genes within 6p21-3 region show graded risk to leprosy in unrelated cohorts of Indian population. Human Genet. 2012;131:703–16.

30. Möller M, Kwiatkowski R, Nebel A, van Helden PD, Hoal EG, Schreiber S. Allelic variation in BTNL2 and susceptibility to tuberculosis in a South African population. Microbes Infect. 2007;9:522–8.

31. Schürmann M, Bein G, Kirsten D, Schlaak M, Müller-Quernheim J, Schwinger E. HLA-DQB1 and HLA-DPB1 genotypes in familial sarcoidosis. Respir Med. 1998;92:649–52.

32. Hofmann S, Franke A, Fisher A, Jacobs G, Nothnagel M, Gaede KJ, et al. Genome-wide association study identifies ANXA11 as a new susceptibility locus for sarcoidosis. Nat Genet. 2008;40:1103–6.

33. Veltkamp M, van Moorsel CH, Rijkers GT, Ruven HJ, Grutters JC. Genetic variation in the Toll-like receptor gene cluster (TLR10-TLR1-TLR6) influences disease course in sarcoidosis. Tissue Antigens. 2012;79:25–32.

34. Fischer A, Schmid B, Ellinghaus D, Nothnagel M, Gaede KJ, Schürmann M, et al. A novel sarcoidosis risk locus for Europeans on chromosome 11q13.1. Am J Respir Crit Care Med. 2012;186:877–85.

35. Fischer A, Ellinghaus D, Nutsua M, Hofmann S, Montgomery CG, Iannuzzi MC, et al. Identification of immune-relevant factors conferring sarcoidosis genetic risk. Am J Respir Crit Care Med. 2015;192:727–36.

36. Smith G, Brownell I, Sanchez M, Prystowsky S. Advances in the genetics of sarcoidosis. Clin Genet. 2008;73:401–12.

37. Suzuki H, Ota M, Meguro A, Katsuyama Y, Kawagoe T, Ishihara M, et al. Genetic characterization and susceptibility for sarcoidosis in Japanese patients: risk factors of BTNL2 gene polymorphisms and HLA class II alleles. Invest Ophtalmo Vis Sci. 2012;53:7109–15.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central