Neither the BD-associated genetic risk locus within the HLA-B/MICA region nor being on immunosuppressive medications explained the differences between patients and controls.
We found that the HLA-B*51 allele and the rs76546355/rs116799036 MHC SNP are independent genetic risk factors for BD in Iranian, and that positivity for the rs76546355/rs116799036 risk allele, but not for B*51, does correlate with specific demographic characteristics or clinical manifestations in BD patients.
The alterations in the nature and affinity of HLA-B*51·peptide complexes probably affect T-cell and natural killer cell recognition, providing a sound basis for the joint association of ERAP1 and HLA-B*51 with BD.
The repeated RPE analysis after separating HLA-B alleles carrying B51-KIR binding sequence as distinct alleles within a broad-type allele group revealed B*2702 allele as the only allele showing an association with BD after the deletion of B*51.
In the current study of Behçet disease (BD), nonsynonymous variants (NSVs) identified by deep exonic resequencing of 10 genes found by GWAS (IL10, IL23R, CCR1, STAT4, KLRK1, KLRC1, KLRC2, KLRC3, KLRC4, and ERAP1) and 11 genes selected for their role in innate immunity (IL1B, IL1R1, IL1RN, NLRP3, MEFV, TNFRSF1A, PSTPIP1, CASP1, PYCARD, NOD2, and TLR4) were evaluated for BD association.
The significant increase of HLA-A*2602 and B*3901 in the patient group without HLA-B*51 suggests that these two alleles might also have some secondary influence on the onset of BD.
In this study, HLA-B*5101 gene from 37 individuals including Japanese, Turkish, Jordanian and Iranian patients and healthy controls were fully sequenced to further clarify the B*5101 gene in association with BD.
Along with human leukocyte antigen gene encoding B*51 (HLA-B*51) and areas including the major histocompatibility complex class I, genome-wide association studies have recognized numerous other BD susceptibility genes including those encoding interleukin (IL)-10, IL-12 receptor β 2 (IL-12RB2), IL-23 receptor (IL-23R), C-C chemokine receptor 1 gene, signal transducer and activator of transcription 4 (STAT4), endoplasmic reticulum aminopeptidase (ERAP1), and genes encoding killer cell lectin-like receptor family members (KLRC4-KLRK1).
Since the identification of HLA-B*5101 (and more recently of MICA) as a susceptibility locus for BD, the identification of additional genetic locus/loci, whether inside, or perhaps more importantly outside the MHC has clearly stalled.
The role of the HLA-B*51 gene has been confirmed in recent years, although its contribution to the overall genetic susceptibility to Behçet's disease was estimated to be only 19%.
The evaluation of the IL-2 gene polymorphism (p=0.0065) and IL-10 gene polymorphism (p=0.0483) distributions with respect to age of BD onset revealed a statistically significant distribution.
The alterations in the nature and affinity of HLA-B*51·peptide complexes probably affect T-cell and natural killer cell recognition, providing a sound basis for the joint association of ERAP1 and HLA-B*51 with BD.
Our findings suggest that gene-gene interactions between HLA-B*51 and ERAP1 variants is important for BD development, however, ERAP1 variants which interact with HLA-B*51 may differ among disease phenotypes or populations.
Therefore, although they were not statistically significant, our data were consistent with an association between ERAP1 and BD as well as with an epistatic interaction between ERAP1 and HLA-B in the Spanish population.
The prevalent distribution in a specific geographical area spanning the Mediterranean basin and Asia, the close association with human leukocyte antigen B*51 in different ethnic groups, and the familial clustering of BD are hallmarks accounting for the strong contribution of a genetic background.