Association between sarcoidosis and HLA polymorphisms in a Czech population from Central Europe: focus on a relationship with clinical outcome and treatment

1. Introduction

Sarcoidosis is a multi-systemic inflammatory disease characterized by non-caseating granulomas of unknown etiology that majorly affect the lung but can affect any organ system. The presentation of sarcoidosis is diverse, and the course of the disease can vary from acute to sub-acute inflammatory Löfgren’s syndrome (LS) to progressive form, which occurs in almost 20% of patients (1) and may eventually include fibrosis.

The prevalence and incidence of sarcoidosis vary substantially by ethnicity, region, sex, and age of onset. The lowest incidence and prevalence are in the Asians; in the Caucasian and African American populations, the incidence is more pronounced. The incidence of sarcoidosis in the Czech population is 3.1 per 100,000 inhabitants; its prevalence is 63.1 per 100,000 (2, 3). In Czech patients, apart from the lungs, another extrapulmonary manifestation has also been reported; extrapulmonary impairment is present in 10–90% of patients, depending on the analysis (4, 5).

The genetic background of individuals has been suggested as an essential factor in sarcoidosis susceptibility. Several genetic association studies have supported the genetic basis of sarcoidosis, its sub-phenotypes, and clinical course (6, 7). Studies also reported different genetic variants based explicitly on the ethnic background of the studied subjects (8). The crucial role of the immunogenetic polymorphic system, HLA (Human Leukocyte Antigen), in susceptibility to sarcoidosis has been shown by some previous investigators enrolling mostly patients of North European descent (7, 9, 10); a previous study in a group of 114 Czech, i.e., Central European patients, was limited to a single class II (HLA-DRB1) locus (11).

To investigate a possible role of HLA polymorphisms across class I and II loci in sarcoidosis in this particular population, we enrolled three hundred one well-defined Czech patients and employed precise NGS (Next Generation Sequencing) genotyping focusing on associations between HLA and sarcoidosis clinical phenotypes.

2. Patients and methods

2.1. Patients and control subjects

A total of 301 sarcoidosis patients enrolled in this study were diagnosed and followed according to the ATS (American Thoracic Society), ERS (European Respiratory Society), and WASOG (World Association of Sarcoidosis and other Granulomatous Disorders) statement on sarcoidosis (12) at the Department of Pneumology and Phtiseology, University Hospital and Faculty of Medicine, Brno, the Czech Republic. The control group consisted of 309 unrelated healthy blood donors from the same geographical region of the Czech Republic, the place of sampling was Olomouc, a city located close (70 km) to Brno. Both patients and control subjects were of Czech ethnicity as assessed by place of birth, surname, and Czech as a native language. For the characteristics of patients and control subjects see Table 1. Study participants provided informed consent with participation in the study, which was approved by the institutional review board–the Ethic committee of the University Hospital and the Faculty Medicine Palacký University in Olomouc, vote on project NV18-05-00134 dated 27.06.2017).

TABLE 1

Table 1. Characteristics of the patients and healthy control subjects enrolled in the study.

Genomic DNA from peripheral blood samples obtained from patients and control subjects was isolated using the Arrow automated extraction system (Isogen Life Science, PW De Meern, Utrecht, Netherlands).

2.2. NGS HLA analysis

To determine HLA alleles, DNA samples obtained from the study subjects were sequenced using the Omixon Holotype HLA 96/11 and 96/7 (Omixon Biocomputing Ltd., Budapest, Hungary) kit on the Illumina MiSeq next-generation sequencing platform (Illumina, San Diego, CA, USA). Long-range PCR amplified the HLA class I and class II loci. All seven amplicons from each sample were pooled into a final 35 μl volume on a new 96-well PCR plate and purified from residual primers and unincorporated nucleotides. Then follows library preparation with the fragmentation of pooled amplicons, fragment end repair, and ligation of sample-specific indexed adaptors. Equal aliquots of indexed sample-specific libraries were combined into one pooled library to carry out magnetic bead-based library clean-up. The concentration of the size selected library was determined on LightCycler 480 II (Roche Diagnostics, Mannheim, Germany) real-time PCR instrument using KAPA Sybr Fast qPCR Master Mix (KAPA Biosystems, Boston, MA, USA). Before sequencing, the library was denatured by NaOH, diluted with hybridization buffer, and a 9 pM library was loaded on a MiSeq flow cell (Illumina, San Diego, CA, USA) sequenced in a single 300 cycle (V2) paired-end sequencing run. Collected reads were exported in fastq format (13).

Allele assignment was performed with the Omixon Twin software v4.1.0.–v4.4.1. and the IPD-IMGT/HLA database Release 3.34−3.45.

The success rate of NGS HLA-typing in all seven loci is presented in Supplementary Table 1.

The HLA-DPB1 locus was excluded from further analyses due to a suboptimal genotyping score in the patients (84%). HLA-DQB1 locus was genotyped in all control samples and in 205 out of 301 patient samples.

2.3. Statistical analysis

Frequencies of the determined HLA alleles were obtained by direct counting (Supplementary Tables 2, 3). Association between the allele frequency and sarcoidosis and its clinical phenotypes was evaluated by two-sided Fisher’s exact probability test, providing odds ratio (OR), 95% confidence interval (CI), and level of significance (p-value < 0.05). Bonferroni’s correction for multiple comparisons was performed for each HLA loci as p_corr = 1-(1-p)ⁿ (n of variants in loci at four digits). The statistical analysis was performed using¹ (9.09.2022), 2 × 2 Contingency Table vassarstats.net (9.09.2022) and excel tables; p < 0.05. For subanalyses between clinical phenotypes, an additional parameter-carriage rate-was employed; this parameter, also known as HLA “antigen frequency,” was determined as a number of individuals carrying a particular variant divided by the number of all individuals tested for the given locus; comparisons using this data are shown in Supplementary Tables 4−7.

3. Results

3.1. Comparison of the distribution of HLA variants between patients and healthy population

We analyzed 6 HLA loci (HLA-A,-B,-C,-DRB1,-DQA1, and -DQB1) in 301 unrelated patients with sarcoidosis diagnosed and followed according to the ATS, ERS, and WASOG statement on sarcoidosis (12) at the Department of Pneumology and Phtiseology, University Hospital and Faculty of Medicine, Brno, Czechia and 309 unrelated Czech healthy control subjects, blood donors from the same geographical region of Czechia. The numbers of the determined variants in the typed loci in patients and control subjects are shown in Supplementary Table 1 and frequencies of determined variants in Supplementary Tables 2, 3.

When comparing the distribution of HLA polymorphisms between sarcoidosis patients and control subjects, we observed 18 HLA variants associated with sarcoidosis on the primary level, 3 of those in HLA class I and 15 in HLA class II, respectively; for numbers that correspond to loci see Supplementary Table 1 (the right column).

Table 2 shows the variants which remained associated with the disease after the correction for multiple comparisons. All were within HLA class II loci. HLA-DQB1*06:02 and HLA-DQB1*06:04 can be considered as susceptibility (risk) factors for sarcoidosis; the frequency of the HLA-DQB1*06:02 in the patients was increased 1.75 times compared with the control subjects, and the HLA-DQB1*06:04 frequency increased four times in the patients. HLA-DRB1*01:01, HLA-DQA1*03:01, and HLA-DQB1*03:02 may confer protection against sarcoidosis: frequency of HLA-DRB1*01:01 decreased in the patients 2.5 times compared with the healthy control subjects, HLA-DQA1*03:01 dropped two times, and HLA-DQB1*03:02 4.5 times.

TABLE 2

Table 2. Summary of the observed associations between HLA alleles and sarcoidosis in the investigated Czech patients (only the alleles associated at least on primary level are shown).

3.2. HLA variants in patients with Löfgren’s syndrome

We performed several subgroup analyses to evaluate the relationship between HLA and sarcoidosis clinical phenotypes. We first compared the patients with Löfgren’s syndrome (LS), a benign sarcoidosis phenotype, and those with non-Löfgren sarcoidosis. Ten variants were associated with LS on the primary level, and six remained significantly associated after correction for multiple comparisons (Table 3). Regarding HLA class I, variants HLA-A*01:01 (p_corr = 3.86 × 10^–04), HLA-B*08:01 (p_corr = 2.50 × 10^–04), and HLA-C*07:01 (p_corr = 2.46 × 10^–04) were overrepresented in the patients with LS (Table 3). When comparing carriage rates in LS and non-LS groups, HLA-B*08:01 (p_corr = 9.34 × 10^–05) and HLA-C*07:01 (p_corr = 0.002) again appeared more frequently (Supplementary Table 4), and can be considered risk factors for LS.

TABLE 3

Table 3. Comparison of the HLA allele frequencies between the patients with Löfgren’s syndrome, LS (n = 55) compared with the patients with non-LS (n = 244).

Regarding HLA class II, 3 variants associated with LS: HLA-DRB1*03:01 (p_corr = 5.16 × 10^–08); HLA-DQA1*05:01 (p_corr = 3.42 × 10^–09) and HLA-DQB1*02:01 (p_corr = 1.67 × 10^–03) (Table 3). Two of these variants were associated with LS when carriage rates were used for comparisons–HLA-DRB1*03:01 (p_corr = 1.93 × 10^–08) and HLA-DQA1*05:01 (p_corr = 0.02) (Supplementary Table 4).

3.3. Comparison of the distribution of HLA variants in the context of different chest radiography stages

Next, we were interested if there is an association between the distinct CXR stage(s) of sarcoidosis and particular HLA variants. Therefore, the observed HLA frequencies were compared between the patient subgroups separated according to the CXR stages. The patients were divided into two groups—in the first group there were patients with CXR stage 1 sarcoidosis, and in the second group, patients with sarcoidosis CXR stage 2−4. We observed nine variants associated with CXR stage 1 on the primary level; two variants remained associated after correction for multiple comparisons (Table 4). Both these variants (HLA-DRB1*03:01 and HLA-DQA1*05:01) occurred three times more often in the CXR 1 patients than in patients with higher CXR stages. In comparisons based on carriage rates, apart from HLA-DRB1*03:01 (p_corr = 6.29 × 10^–05) and HLA-DQA1*05:01 (p_corr = 3.63 × 10^–06) CXR stage 1 associated also with HLA-B*08:01 (p_corr = 0.042) and HLA-C*07:01 (p_corr = 0.011) (Supplementary Table 5).

TABLE 4

Table 4. Comparison of the HLA allele frequencies between the patients with chest X-ray (CXR) stage 1 (n = 95) compared with the patients with stages 2–4 (n = 206); only the alleles associated at least on the primary level are shown.

Of the four variants associated with higher CXR stages on the primary level, two (HLA-DRB1*11:01, HLA-DQA1*05:05) remained associated after correction for multiple comparisons (Table 4). HLA-DRB1*11:01 frequency was four times higher in CXR stages 2−4 than in CXR stage 1 and this variant was associated with higher sarcoidosis CXR stages (p_corr = 0.001) also when comparing carriage rates (Supplementary Table 5).

3.4. HLA variants in the context of clinical disease course

We were also interested if HLA variation can be linked to the clinical course of the disease. We, therefore, compared the patients undergoing remission and the patients in whom the disease persisted; the disease status was assessed 2 years after the diagnosis. Eight variants were associated with remission on the primary level, two variants remain associated after correction for multiple comparisons (Table 5). Patients with variants HLA-DRB1*03:01 or HLA-DQA1*05:01 were three times more likely to develop remission of the disease. Those two variants were significantly associated with remission also in comparison of carriage rates, HLA-DRB1*03:01 (p_corr = 1.65 × 10^–04) and HLA-DQA1*05:01 (p_corr = 6.03 × 10^–05) (Supplementary Table 6). We observed eight variants associated with persistent sarcoidosis on the primary level of significance (Table 5).

TABLE 5

Table 5. Comparison of the HLA allele frequencies between the patients with disease remission (n = 142) compared with the patients with persistent disease course (n = 134); only the alleles associated at least on the primary level are shown.

3.5. HLA variants in the context of corticosteroid treatment

Treatment requirement may be considered as another indicator of sarcoidosis disease course. Our further subanalysis examined whether corticosteroid treatment is connected to a particular HLA variant. We compared the patients who received corticosteroid treatment with those who did not receive it. Nine HLA variants were primarily associated with no need of treatment; three variants (HLA-A*01:01, HLA-DRB1*03:01, and HLA-DQA1*05:01) remained associated after correction for multiple comparisons (Table 6), they occurred twice more often in non-treated patients. The same three variants were associated with non-treated patients also in comparison of carriage rates: HLA-A*01:01 (p_corr = 1.45 × 10^–04), HLA-DRB1*03:01 (p_corr = 0.024), and HLA-DQA1*05:01 (p_corr = 0.007) (Supplementary Table 7). By contrast, one variant associated with the need for treatment: HLA-A*02:01 (p_corr = 0.038) (Supplementary Table 7).

TABLE 6

Table 6. Comparison of the HLA allele frequencies between the patients with treatment (n = 194) compared with the patients without treatment (n = 105); only the alleles associated at least on the primary level are shown.

3.6. HLA variants in patients with an extrapulmonary manifestation of sarcoidosis

Our final subanalysis focused on a possible link between HLA and the extrapulmonary manifestation of sarcoidosis, which was reported in 91 (30%) of our patients. When comparing patients with the extrapulmonary manifestation of sarcoidosis with the patients without it, we observed one variant associated with extrapulmonary manifestation–HLA-DQB1*05:03 (p = 0.033; OR; 95% CI = [2.653; 1.114−6.318]); its frequency in patients with extrapulmonary manifestation was doubled. By contrast, HLA-DRB1*01:01 (p = 0.049; OR; 95% CI = [0.151; 0.02−1.151]) occurred six times more often in patients without extrapulmonary manifestation; both these associations were detected on the primary level.

3.7. Linkage disequilibrium between associated HLA variants

Linkage disequilibrium between the alleles associated with sarcoidosis in Czech patients was calculated. A non-random association was identified between variants associated with the risk of sarcoidosis. The strongest linkage was detected between HLA-DQB1*06:02-HLA-DRB1*15:01. In the context of variants protective from sarcoidosis, the strong association was (among others) between alleles HLA-DQA1*03:01-HLA-DQB1*03:02. Regarding variants associated with higher stages of the disease there was a linkage disequilibrium between variants HLA-DQA1*05:05-HLA-DRB1*11:01. A strong linkage disequilibrium was detected between variants HLA-DQA1*05:01-HLA-DRB1*03:01 which related both to CXR stage 1 and to disease remission. There was also a linkage between variants HLA-A*01:01-HLA-DRB1*03:01 associated with no need for treatment. The variants associated with LS, i.e., HLA-A*01:01, HLA-B*08:01, HLA-C*07:01, HLA-DRB1*03:01, HLA-DQA1*05:01, and HLA-DQB1*02:01 were found in linkage with each other except for HLA-A*01:01- HLA-DQB1*02:01 (Table 7).

TABLE 7

Table 7. Linkage disequilibrium between HLA alleles observed to be associated with sarcoidosis in Czech population among variants being susceptible for sarcoidosis, protective for sarcoidosis, associated with unfavorable prognosis and associated with better prognosis in sarcoidosis.

4. Discussion

This study, performed in the Czech population, investigated the association of HLA polymorphisms with sarcoidosis and focused on associations between HLA and sarcoidosis clinical phenotypes. Observations of an association between sarcoidosis and HLA-DQB1*06:04 as a susceptibility and of HLA-DQA1*03:01 as a protective marker in the Czech population represent novel findings. Further, this study, including altogether 610 subjects, confirmed associations between HLA-DQB1*06:02 (susceptibility) and HLA-DRB1*01:01, HLA-DQB1*03:02 (protection), previously reported in smaller-scale studies in other, namely North-European and also East-Asian populations (13–17). Regarding clinical phenotypes of sarcoidosis, HLA-A*01:01, HLA-B*08:01, HLA-C*07:01, HLA-DRB1*03:01, HLA-DQA1*05:01, and HLA-DQB1*02:01 associated with Löfgren’s syndrome. The variants HLA-DRB1*03:01 and HLA-DQA1*05:01 were associated with disease remission and the presence of CXR stage 1, while the variants HLA-DRB1*11:01 and HLA-DQA1*05:05 associated with CXR stages 2−4. Associations were also observed with other clinically relevant disease features, i.e., non-requirement of corticosteroid treatment and systemic disease: i.e., HLA-A*02:01 was overrepresented in patients needing treatment, and HLA-DQB1*05:03 associated with extrapulmonary sarcoidosis.

4.1. HLA in sarcoidosis patients compared with the healthy population

This first detailed analysis of the HLA variants in 6 loci (HLA-A,-B,-C,-DRB1,-DQA1, and -DQB1) could compare so far the highest number of patients and control subjects from the Czech population.

In line with their role in sarcoidosis immunopathogenesis (1, 8, 18), HLA class II loci emerged as the central location of the disease-associated variants: HLA-DQB1*06:02, here defined as a risk variant; was previously reported to be associated mainly with a worse prognosis in sarcoidosis in Swedish and Dutch patients (14, 19) and linked with variants HLA-DRB1*15:01 and HLA-B*07 (15). In our study, HLA-DRB1*15:01 and HLA-B*07:02 were also observed as risk factors for sarcoidosis at the primary significance level. As these two variants were reported in linkage disequilibrium with the HLA-DQB1*06:02 variant (Table 7), it may be speculated that the HLA-DQB1 locus allele is responsible for the observed association.

In this context, we newly characterized variant HLA-DQB1*06:04 as a risk factor for sarcoidosis development, contrary to a previous report from the Netherlands, where there was no difference in the occurrence of this variant in patients compared to the controls (15). The discrepancy could reflect different cohort sizes and distinct ethnic backgrounds of the investigated patients (8).

Variant HLA-DRB1*01:01 was previously reported as protective against sarcoidosis in the population of the United Kingdom, Netherlands, Finland, and Japan, where it was identified as part of the haplotype DQA1*01:01 ∼ DQB1*05:01 ∼ DRB1*01:01; all these variants we report to be in LD and act as protective against sarcoidosis development (17, 20, 21). Except for HLA-DRB1*01:01, which occurred more frequently in our Czech control population than in our patients, the other two protective variants were HLA-DQA1*03:01 and HLA-DQB1*03:02. HLA-DQB1*03:02 was recently reported as protective from sarcoidosis also in Korean population (13). HLA-DQA1*03:01 has so far never been mentioned in connection with sarcoidosis. A hereby observed association can, however, result from LD between HLA-DQA1*03:01 and HLA-DQB1*03:02, and another protective variant–HLA-DRB1*04:01, which we report to protect on the primary level. These three variants are part of one of the most common haplotypes in Caucasian populations (22). Another common haplotype, HLA-DQB1*06:02 ∼ HLA-DQA1*01:02, relates to the risk of sarcoidosis in this study.

4.2. HLA and clinical phenotypes of sarcoidosis

We determined variant HLA-DRB1*11:01 as a risk factor for more advanced CXR stages of sarcoidosis. This variant was previously reported to relate to sarcoidosis in general across a spectrum of ethnicities (9), but we are the first to connect this variant with higher CXR stages of the disease, specifically for the Czech population.

Regarding Löfgren’s syndrome, we report six variants (HLA-A*01:01∼HLA-B*08:01∼HLA-C*07:01∼HLA-DRB1*03:01 ∼HLA-DQA1*05:01∼HLA-DQB1*02:01), forming 8.1 haplotype, to be associated with LS. The 8.1 ancestral haplotype is most common in Northern Europe and was associated with risk for autoimmune diseases (23), it was also described as associated with sarcoidosis in populations from Northern Europe (18, 24, 25). The linkage of variants HLA-B*08, HLA-DRB1*03:01, and HLA-DQB1*02:01 with good prognosis of sarcoidosis was observed previously in a Croatian population (26). A favorable disease course was also previously linked to HLA-DQB1*02:01 (19, 27).

Some of the variants of the 8.1 haplotype, HLA-DRB1*03:01 and HLA-DQA1*05:01, were associated with a milder disease course, i.e., initial chest X-ray (CXR 1) or remitting disease. Further, together with HLA-A*01:01, these variants occurred more frequently in patients not requiring corticosteroid treatment. Our observations, therefore, extend findings from previous reports from sarcoidosis association studies across populations that showed the link between HLA-DRB1*03:01 and sarcoidosis patients, especially in the group with a better prognosis (11, 27, 28). We can also speculate about the effect of HLA-DRB1*03:01 homozygosity on better prognosis as four of five homozygotes among our Czech sarcoidosis patients had a remitting disease, a similar relationship was previously observed in HLA-DRB1*03:01 homozygous patients with remission of multiple sclerosis (29). The association of variant HLA-B*08:01 with the milder disease and especially with LS has also been previously reported (27, 30). Concerning our observation of HLA-DQA1*05:01 associating with better prognosis, this variant was reported to be overrepresented in Japanese sarcoidosis patients compared with controls (31), and HLA-A*01:01 was reported with the resolving course of disease in Sweden (27).

Regarding systemic disease, previous reports from Swedish, UK, and Turkish populations associated alleles of the HLA-DRB1 locus (HLA- DRB1*04 -DRB1*15) with extrapulmonary sarcoidosis (16, 32, 33). In our study, extrapulmonary involvement was associated with another class II locus, HLA-DQB1: the variant HLA-DQB1*05:03, previously mentioned in connection with sarcoidosis uveitis in UK patients (16), was associated with extrapulmonary disease in our Czech patients. Further, HLA-DQB1*05:03 was associated with persistent disease on the primary level, as well as HLA-DRB1*14:54 and HLA-DQA1*01:04; these three variants occur in LD in the Czech population. In an Indian population, these variants were associated with the risk of sarcoidosis (34).

In context of treatment requirement, there has been scarce information on the relationship between immunogenetic markers and non-requirement for corticosteroid therapy. In a previous study, the presence of the mutant CCR5Delta32 allele carried an increased relative risk of the need for treatment in another cohort of Czech sarcoidosis patients (35). Here we observed that HLA-A*02:01 carriers require treatment two times more often than patients not carrying this variant. As corticosteroid treatment is indicated for patients with significant pulmonary sarcoidosis involvement, who are expected to have a higher risk of future mortality or permanent disability (36), the availability of an additional laboratory marker to be included in the complex decision for treatment initiation may be useful. However, before potential clinical application, this finding must be replicated. On the other hand, no need for treatment was associated with three alleles of the 8.1 haplotype and this observation fits well with the hereby reported link of this haplotype with milder sarcoidosis.

4.3. Study strength, limitations, and conclusions

The main strength of this study is the large number of patients characterized in detail and enrolled in a single medical center, allowing comparisons between adequately sized subgroups according to clinical diagnostic standards. We analyzed data from 301 precisely characterized sarcoidosis patients with representative phenotypes. In this context, a low number of patients in CXR stage 4 was the main limitation, the inability to address the allele primarily responsible for an association in cases of linked allele pairs may also be considered limiting. The methodology employed for HLA-typing represents another substantial advantage over previous studies as the precise NGS technique allowed indeterminate characterization at the allele level. Certain limitation of the clinical utility of hereby newly reported results is their validity in a Czech population, unless further replicated elsewhere. On the other hand, the study results could be applied in the meta-analyses of sarcoidosis immunogenetic data.

In conclusion, while the HLA region has long been consistently associated with sarcoidosis, the majority of the reports focused on populations of Northern/Western Europe origin and pointed to class II loci, namely the HLA-DRB1, in line with its major role in antigen presentation (1, 8, 18). In this context, our study provides further information obtained from patients of the Slavonic population from Central Europe, and re-emphasizes the possible role of another class II locus—HLA-DQB1. Study findings nominate a new susceptibility variant for sarcoidosis, HLA-DQB1*06:04, and extend the previous knowledge on HLA-DQB1*06:02 as a sarcoidosis susceptibility marker to the Czech population. Our study also suggests that the variant HLA-DRB1*11:01, previously reported as a general susceptibility factor, may serve as a marker for the advancing course of the disease. Further, other HLA variants associated with extrapulmonary involvement and the requirement for treatment are reported in our study. Patients with these HLA markers should, therefore, be carefully followed up, screened for extrapulmonary involvement, and considered for treatment.

Finally, a possible contribution of the “autoimmune” 8.1 ancestral haplotype to the better sarcoidosis prognosis, is another important outcome of our study. Association of Löfgren’s syndrome and other milder disease forms with the combination of the following six variants, HLA-A*01:01∼HLA-B*08:01∼HLA-C*07:01∼HLA-DRB1*03:01∼HLA-DQA1*05:01∼ HLA-DQB1*02:01, may be, therefore, considered in the context of the personalized approach to sarcoidosis patients.

Previous presentations

Part of this work has been presented as an abstract at the 2021 ERS (European Respiratory Society) International Congress, September 5.-8., 2021.

Data availability statement

The datasets for this article are not publicly available due to concerns regarding participant/patient anonymity and confidentiality. Requests to access the datasets should be directed to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the University Hospital and the Faculty of Medicine Palacký University in Olomouc, vote on project NV18-05-00134 dated 27.06.2017. The patients/participants provided their written informed consent to participate in this study.

Author contributions

MP and MD: conception and design. MD, MP, JP, and MF-V: provision of study materials. KS, MD, MP, KO, MF-V, JP, LK, and AS: collection and assembly of data, manuscript writing, and final approval of the manuscript. KS, MP, KO, MD, and MF-V: data analysis and interpretation. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Ministry of Health of the Czech Republic (grant number: NV18-05-00134).

Acknowledgments

We thank Jitka Svobodová for her assistance with patients enrollment. We appreciate the suggestions of two reviewers for improvement of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2023.1094843/full#supplementary-material

Footnotes

1. ^ https://www.hiv.lanl.gov/content/immunology/hla/index.html

References

1. Calender A, Weichhart T, Valeyre D, Pacheco Y. Current insights in genetics of sarcoidosis: Functional and clinical impacts. J Clin Med. (2020) 9:2633. doi: 10.3390/jcm9082633

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Arkema EV, Cozier YC. Sarcoidosis epidemiology: Recent estimates of incidence, prevalence and risk factors. Curr Opin Pulm Med. (2020) 26:527–34. doi: 10.1097/MCP.0000000000000715

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kolek V. Epidemiological study on sarcoidosis in Moravia and Silesia. Sarcoidosis. (1994) 11:110–2.

PubMed Abstract | Google Scholar

4. Zurkova M, Kolek V, Tomankova T, Kriegova E. Extrapulmonary involvement in patients with sarcoidosis and comparison of routine laboratory and clinical data to pulmonary involvement. Biomed. Papers. (2014) 158:613–20. doi: 10.5507/bp.2014.026

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Doubková M, Pospíšil Z, Skřičková J, Doubek M. Prognostic markers of sarcoidosis: An analysis of patients from everyday pneumological practice. Clin Respir J. (2015) 9:443–9. doi: 10.1111/crj.12160

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Valeyre D, Bernaudin JF, Uzunhan Y, Kambouchner M, Brillet PY, Soussan M, et al. Clinical presentation of sarcoidosis and diagnostic work-up. Semin Respir Crit Care Med. (2014) 35:336–51. doi: 10.1055/s-0034-1381229

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rivera NV, Ronninger M, Shchetynsky K, Franke A, Nöthen MM, Müller-Quernheim J, et al. High-density genetic mapping identifies new susceptibility variants in sarcoidosis phenotypes and shows genomic-driven phenotypic differences. Am J Respir Crit Care Med. (2016) 193:1008–22. doi: 10.1164/rccm.201507-1372OC

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Van Moorsel CHM, Petrek M, Rivera NV. Unravelling the genetic basis of sarcoidosis. In: F Bonnella, DA Culver, D Israël-Biet editors. Sarcoidosis (ERS Monograph). (Vol. 96), Sheffield: European Respiratory Society (2022). p. 41–56. doi: 10.1183/2312508x.10031320

CrossRef Full Text | Google Scholar

9. Rossman MD, Thompson B, Frederick M, Maliarik M, Iannuzzi MC, Rybicki BA. HLA-DRB1* 1101: A significant risk factor for sarcoidosis in blacks and whites. Am J Hum Genet. (2003) 73:720–35. doi: 10.1086/378097

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Levin AM, Adrianto I, Datta I, Iannuzzi MC, Trudeau S, Li J, et al. Association of HLA-DRB1 with sarcoidosis susceptibility and progression in African Americans. Am J Respir Cell Mol Biol. (2015) 53:206–16. doi: 10.1165/rcmb.2014-0227OC

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Mrazek F, Holla LI, Hutyrová B, Znojil V, Vasku A, Kolek V, et al. Association of tumour necrosis factor−α, lymphotoxin−α and HLA-DRB1 gene polymorphisms with Löfgren’s syndrome in Czech patients with sarcoidosis. Tissue Antigens. (2005) 65:163–71. doi: 10.1111/j.1399-0039.2005.00370.x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. 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 ERS. Am J Respir Crit Care Med. (1999) 160:736–55. doi: 10.1164/ajrccm.160.2.ats4-99

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sikorová K, Moon SJ, Yoon HY, Strnad A, Song JW, Petrek M. HLA class II variants defined by next generation sequencing are associated with sarcoidosis in Korean patients. Sci Rep. (2022) 12:9302. doi: 10.1038/s41598-022-13199-w

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Berlin M, Fogdell-Hahn A, Olerup O, Eklund A, Grunewald J. HLA-DR predicts the prognosis in Scandinavian patients with pulmonary sarcoidosis. Am J Respir Crit Care Med. (1997) 156:1601–5. doi: 10.1164/ajrccm.156.5.9704069

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Voorter CE, Drent M, van den Berg-Loonen EM. Severe pulmonary sarcoidosis is strongly associated with the haplotype HLA-DQB1* 0602–DRB1* 150101. Hum Immunol. (2005) 66:826–35. doi: 10.1016/j.humimm.2005.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

16. 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. doi: 10.1093/hmg/ddq325

PubMed Abstract | CrossRef Full Text | Google Scholar

17. 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. doi: 10.1183/09031936.00035213

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Grunewald J. Role of genetics in susceptibility and outcome of sarcoidosis. Semin Respir Crit Care Med. (2010) 31:380–9. doi: 10.1055/s-0030-1262206

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sato H, Grutters JC, Pantelidis P, Mizzon AN, Ahmad T, Van Houte AJ. 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. doi: 10.1165/rcmb.4782

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Grutters JC, Sato H, Welsh KI, du Bois RM. The importance of sarcoidosis genotype to lung phenotype. Am J Respir Cell Mol Biol. (2003) 29:S59.

Google Scholar

21. Naruse TK, Matsuzawa Y, Ota M, Katsuyama Y, Matsumori A, Hara M, et al. HLA-DQB1* 0601 is primarily associated with the susceptibility to cardiac sarcoidosis. Tissue Antigens. (2000) 56:52–7. doi: 10.1034/j.1399-0039.2000.560107.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kockum I, Sanjeevi CB, Eastman S, Landin-Olsson M, Dahlquist G, Lernmark Å. Complex interaction between HLA DR and DQ in conferring risk for childhood type 1 diabetes. Eur J Immunogenet. (1999) 26:361–72. doi: 10.1046/j.1365-2370.1999.00173.x

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gambino CM, Aiello A, Accardi G, Caruso C, Candore G. Autoimmune diseases and 8.1 ancestral haplotype: An update. HLA. (2018) 92:137–43. doi: 10.1111/tan.13305

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Price P, Witt C, Allock R, Sayer D, Garlepp M, Kok CC, et al. The genetic basis for the association of the 8.1 ancestral haplotype (A1, B8, DR3) with multiple immunopathological diseases. Immunol Rev. (1999) 167:257–74. doi: 10.1111/j.1600-065x.1999.tb01398.x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Candore G, Lio D, Romano GC, Caruso C. Pathogenesis of autoimmune diseases associated with 8.1 ancestral haplotype: Effect of multiple gene interactions. Autoimmun Rev. (2002) 1:29–35. doi: 10.1016/s1568-9972(01)00004-0

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Grubić Z, Žunec R, Peroš-Golubičič T, Tekavec-Trkanjec J, Martinez N, Alilović M, et al. HLA class I and class II frequencies in patients with sarcoidosis from Croatia: Role of HLA-B8,- DRB1* 0301, and- DQB1* 0201 haplotype in clinical variations of the disease. Tissue Antigens. (2007) 70:301–6. doi: 10.1111/j.1399-0039.2007.00904.x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Grunewald J, Eklund A, Olerup O. Human leukocyte antigen class I alleles and the disease course in sarcoidosis patients. Am J Respir Crit Care Med. (2004) 169:696–702. doi: 10.1164/rccm.200303-459OC

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Yanardag H, Tetikkurt C, Bilir M, Yilmaz E. Association of HLA antigens with the clinical course of sarcoidosis and familial disease. Monaldi Arch. Chest Dis. (2017) 87:835. doi: 10.4081/monaldi.2017.835

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wu JS, Qiu W, Castley A, James I, Joseph J, Christiansen FT, et al. Presence of CSF oligoclonal bands (OCB) is associated with the HLA-DRB1 genotype in a West Australian multiple sclerosis cohort. J neurol sci. (2010) 288:63–7. doi: 10.1016/j.jns.2009.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

30. 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. doi: 10.1164/rccm.201503-0418OC

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Ishihara M, Ohno S, Ishida T, Ando H, Naruse T, Nose Y, et al. Molecular genetic studies of HLA class II alleles in sarcoidosis. Tissue Antigens. (1994) 43:238–41. doi: 10.1111/j.1399-0039.1994.tb02331.x

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Darlington P, Gabrielsen A, Sörensson P, Tallstedt L, Padyukov L, Eklund A, et al. HLA-alleles associated with increased risk for extra-pulmonary involvement in sarcoidosis. Tissue Antigens. (2014) 83:267–72. doi: 10.1111/tan.12326

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Ozyilmaz E, Goruroglu Ozturk O, Yunsel D, Deniz A, Hanta I, Kuleci S. Could HLA-DRB1*11 allele be a clue for predicting extra-pulmonary sarcoidosis. Sarcoidosis Vasc Diffuse Lung Dis. (2014) 31:154–62.

PubMed Abstract | Google Scholar

34. Sharma SK, Balamurugan A, Pandey RM, Saha PK, Mehra NK. Human leukocyte antigen-DR alleles influence the clinical course of pulmonary sarcoidosis in Asian Indians. Am J Respir Cell Mol Biol. (2003) 29:225–31. doi: 10.1165/rcmb.2003-0007OC

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Petrek M, Drabek J, Kolek V, Zlamal J, Welsh KI, Bunce M, et al. CC chemokine receptor gene polymorphisms in Czech patients with pulmonary sarcoidosis. Am J Respir Crit Care Med. (2000) 162:1000–3. doi: 10.1164/ajrccm.162.3.2001022

CrossRef Full Text | Google Scholar

36. Baughman RP, Valeyre D, Korsten P, Mathioudakis AG, Wuyts WA, Wells A, et al. ERS clinical practice guidelines on treatment of sarcoidosis. Eur Respir J. (2021) 58:2004079. doi: 10.1183/13993003.04079-2020

CrossRef Full Text | Google Scholar