HLA and non-HLA genes in Behçet’s disease: a multicentric study in the Spanish population
- Marco Antonio Montes-Cano1,
- Marta Conde-Jaldón1,
- José Raul García-Lozano1,
- Lourdes Ortiz-Fernández1,
- Norberto Ortego-Centeno2,
- María Jesús Castillo-Palma3,
- Gerard Espinosa4,
- Genaro Graña-Gil5,
- Miguel Angel González-Gay6,
- Ana Celia Barnosi-Marín7,
- Roser Solans8,
- Patricia Fanlo9,
- Teresa Camps10,
- Santos Castañeda11,
- Juan Sánchez-Bursón12,
- Antonio Núñez-Roldán1,
- Javier Martín13 and
- María Francisca González-Escribano1Email author
© Montes-Cano et al.; licensee BioMed Central Ltd. 2013
Received: 4 October 2012
Accepted: 19 September 2013
Published: 4 October 2013
According to genome wide association (GWA) studies as well as candidate gene approaches, Behçet’s disease (BD) is associated with human leukocyte antigen (HLA)-A and HLA-B gene regions. The HLA-B51 has been consistently associated with the disease, but the role of other HLA class I molecules remains controversial. Recently, variants in non-HLA genes have also been associated with BD. The aims of this study were to further investigate the influence of the HLA region in BD and to explore the relationship with non-HLA genes recently described to be associated in other populations.
This study included 304 BD patients and 313 ethnically matched controls. HLA-A and HLA-B low resolution typing was carried out by PCR-SSOP Luminex. Eleven tag single nucleotide polymorphisms (SNPs) located outside of the HLA-region, previously described associated with the disease in GWA studies and having a minor allele frequency in Caucasians greater than 0.15 were genotyped using TaqMan assays. Phenotypic and genotypic frequencies were estimated by direct counting and distributions were compared using the χ2 test.
In addition to HLA-B*51, HLA-B*57 was found as a risk factor in BD, whereas, B*35 was found to be protective. Other HLA-A and B specificities were suggestive of association with the disease as risk (A*02 and A*24) or protective factors (A*03 and B*58). Regarding the non-HLA genes, the three SNPs located in IL23R and one of the SNPs in IL10 were found to be significantly associated with susceptibility to BD in our population.
Different HLA specificities are associated with Behçet’s disease in addition to B*51. Other non-HLA genes, such as IL23R and IL-10, play a role in the susceptibility to the disease.
Behçet’s disease (BD) is a systemic vasculitis characterized mainly by recurrent oral and genital ulceration, although other clinical manifestations, such as skin lesions, ocular, gastrointestinal and neurological disorders, are relatively common . The aetiology of BD remains unclear; nevertheless, imbalances in the innate or adaptive immune response triggered by infectious agents or environmental factors in genetically predisposed individuals could be some of the underlying mechanisms of the disease. Evidence of genetic contributions to the pathogenesis of the disease are based on familial aggregation, predominance in patients with Mediterranean or Asian ancestry and association with the human leukocyte antigen (HLA) B51 (HLA-B51) region in several ethnic groups. Multiple studies investigating associations between HLA and BD have been reported [2, 3]. All of them described an association between the disease and the specificity HLA-B51 independent of the ethnicity of the study population . Nevertheless, associations of this pathology with other HLA-B allele groups are less consistent [5–12]. This lack of replication may be due to the lack or weaker association of BD with other HLA-B allele groups, as well as with the relatively small sample size included in each study, which makes it difficult to find an association with relatively underrepresented allele groups. The typing procedures have evolved to achieve a better resolution. Nevertheless, the basis of the association of HLA with BD remains elusive. Very recently, in a study employing dense genotyping in the HLA region, three peaks of association with the disease were found: one of these single-nucleotide polymorphisms (SNPs) is located in the HLA-A region, another is located in the HLA-B human major histocompatibility complex class I chain-related gene A (HLA-B-MICA) region and another is located in the HLA-C region. After inferring the classical HLA alleles from the SNPs, the authors of that suggested that the robust association between HLA-B*51 and BD is explained by a variant located between the HLA-B and MICA genes . Improved knowledge of the HLA-B molecules associated with BD would help to establish the basis of this association. Contribution of the HLA region has been estimated to represent approximately 20% of the genetic component of this disease . Different studies designed to establish the contribution of other HLA and non-HLA genes have been published. Two different approaches have been used in genetic association studies and in BD. First, the candidate gene strategy has permitted the establishment of the association of BD with HLA-B51 and with other HLA and non-HLA genes . Second, the free hypothesis approach has also been used to assess the genetic component of BD [16–20]. Three of these studies, which were carried out using a large-scale SNP genotyping strategy (genomewide association studies (GWASs)), were recently published [18–20]. In general, these studies confirmed the association with HLA-B, identified a second HLA class I region independently associated with the disease in the HLA-A region and described an association with other non-HLA regions [19, 20]. The two aims of the present study were to further investigate the influence of the HLA region in BD and to explore the relationship with the previously described non-HLA genes associated with BD.
A total of 304 BD patients (134 males and 170 females) with a mean age at onset ± SD of 38.7 ± 13.8 years fulfilled the 1990 International Study Group classification criteria for BD  and they were recruited into the study. A total of 313 bone marrow and blood donors (50% males) were included in the study as normal healthy controls. All the participants were Spanish Europeans recruited from different Spanish hospitals. The study was approved by all local ethical committees of the corresponding hospitals, and all the study participants gave their written informed consent to participate. Clinical features of the patient group were as follows: 100% had oral ulcers; 74.4% had genital ulcers; 49.8% had uveitis; 47.4% had arthritis; and 24.5% had vascular involvement, 18.1% had neurological involvement and 19.1% had gastrointestinal. The distribution of the frequencies of the different markers in the cohorts from different hospitals was not significantly different.
Peripheral blood was obtained from the healthy controls, and peripheral blood or saliva served as the starting material obtained from patients. Samples were obtained from participants after they provided their written informed consent. Genomic DNA was extracted using the QIAamp DNA Mini Kit (Qiagen, Barcelona, Spain) according to the manufacturer’s recommendations and stored at -20°C. The purity of DNA was determined using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). Only DNA samples having a 260 nm to 280 nm absorbance ratio between 1.7 and 2.0 and a final concentration of 10 to 20 ng/μl were considered appropriate. Ten DNA samples from saliva were eliminated because they did not meet these quality criteria.
HLA-A and HLA-B low-resolution typing was carried out by using the polymerase chain reaction sequence-specific oligonucleotide (PCR-SSOP) Luminex method using LABType SSO (One Lambda Inc., Canoga Park, CA, USA) according to the manufacturer’s instructions. Briefly, target DNA was PCR-amplified using HLA-A or HLA-B group-specific primers. The biotinylated PCR products were denatured and hybridized with specific probes bound to color-coded microspheres. Phycoerythrin-conjugated streptavidin was used to label and reveal reactions, and a flow analyser, LABScan 100 (One Lambda Inc), was used to identify fluorescence intensity in each microsphere. The software HLA Fusion 2.0 (One Lambda Inc) was used to assign the HLA typing of each locus. Typing of locus A and locus B was possible for 92% and 93%, respectively, of the DNA samples meeting the quality criteria.
Additionally, an intergene (LOC100129342) and six gene regions (IL23R, IL10, CPVL, TDRD7/C9orf174, UBASH3B and UBAC2) previously reported to be associated with BD in GWASs in candidate gene studies were included [18–20, 22]. The SNPs were selected from among those included in previous publications, taking into account the minor allele frequency (MAF) in the Utah residents with Northern and Western European ancestry from the Centre d’Etude du Polymorphisme Humain collection (CEU) population . Only SNPs with MAF greater than 0.15 in the CEU were included in this study, one each for LOC100129342 (rs11206377), CPVL (rs317711), TDRD7 (rs2061634), UBASH3B (rs4936742) and UBAC2 (rs7999348) and three for each of the two non-HLA loci with the highest evidence of association with BD in different populations, which are IL23R (rs17375018, rs7517847 and rs1343151) and IL10 (rs3024498, rs2222202 and rs1800872). Genotyping of these SNPs was performed using the TaqMan SNP Genotyping Assay (Applied Biosystems, Barcelona, Spain) and the LightCycler 480 System (Roche, Barcelona, Spain).
Phenotypic and genotypic frequencies were estimated by direct counting, and distributions were compared using the χ2 test. Odds ratios (ORs) and 95% confidence intervals (95% CI) were calculated according to Woolf’s method using the StatCalc calculator with the EPI Info software package (Centers for Disease Control and Prevention, Atlanta, GA, USA). In the case of the HLA univariate analysis, P-values ≤0.001 (0.05 divided by 48 which is the number of groups is the number of HLA alleles found in our population) were considered significant, whereas P-values between 0.01 and 0.001 were considered to be suggestive of an association and P-values between 0.05 and 0.01 were considered marginally associated. The adjusted ORs were calculated after constructing a logistic regression model categorizing all the individuals according to the presence or absence of the HLA-A and HLA-B groups with P < 0.05 in the univariate analysis (Epi Info 2002). In the cases of HLA risk factors, the analysis proposed by Svejgaard and Ryder to detect the strongest HLA association was also used. This analysis consisted of different two-tailed Fisher’s exact tests performed in various 2 × 2 tables using the GraphPad Software QuickCalcs tool (http://graphpad.com/quickcalcs/contingency1.cfm; GraphPad Software, La Jolla, CA, USA), starting with the basic data collected in a 2 × 4 table that included the four phenotypic combinations of the two HLA factors under study . The ORs were calculated according to Haldane’s modification of Woolf’s method, and P-values of tests 1 through 4 were corrected by a factor of 6 as these authors suggested.
Gene–gene interactions were evaluated using the nonparametric multifactor dimensionality reduction (MDR) method [25–27]. A tenfold cross-validation was used to estimate the balanced testing accuracy (TA) and cross-validation consistency (CVC) of MDR models. Balanced accuracy is (sensitivity + specificity)/2, where sensitivity is true-positive/(true-positive + false-negative) and specificity is true-negative/(false-positive + true-negative). The MDR model with the highest TA and CVC values was selected as the best model for studying from one to n loci. Statistical significance was evaluated using a 10,000-fold permutation test.
Allelic frequencies of the HLA-A and HLA-B allele groups in Behçet’s disease patients and healthy controls a
OR (95% CI)
2n = 542 (%)
2n = 626 (%)
1.47 (1.14 to 1.91)
0.53 (0.34 to 0.83)
1.70 (1.10 to 2.62)
2n = 548 (%)
2n = 626 (%)
0.58 (0.35 to 0.95)
0.49 (0.32 to 0.76)
0.48 (0.23 to 0.99)
4.11 (2.79 to 6.06)
3.70 (1.78 to 7.83)
0.10 ( 0.0 to 0.71)
Logistic regression analysis a
OR adjusted values (95% CI)
1.18 (0.80 to 1.72)
0.66 (0.38 to 1.11)
1.84 (1.11 to 3.04)
0.59 (0.35 to 1.02)
0.50 (0.30 to 0.84)
0.61 (0.29 to 1.30)
4.96 (3.21 to 7.66)
4.16 (1.98 to 8.76)
0.11 (0.01 to 0.91)
Basic data and 2 × 2 comparisons to detect the strongest human leukocyte antigen association of the risk factors a
One-locus analysis of the non–human leukocyte antigen genes included in this study a
OR (95% CI)
The first finding of the present work is the confirmation of a second HLA-B allele group as a risk factor for BD, at least in European populations. Thus, in addition to B*51, we confirmed association of B*57 with BD in a model of logistic regression analysis, an independent association according to the analysis recommended by Svejgaard and Ryder  and that the model B*51+B*57 performed the best in the two-attribute model in the MDR analysis including HLA and non-HLA factors associated with the disease, Consistent with one previous study describing the same association in a UK cohort, the OR for B*57 was similar to that for B*51 in our study population . Moreover, HLA-B*35 was detected as protective in our cohort. This HLA-B allele group had previously been found to be protective in an Italian cohort with a very limited sample size . Interestingly, our study suggests a protective effect of B*58 because its frequency is significantly decreased among patients before correction and remains associated in the logistic regression analysis. Although this result needs confirmation, it is in agreement with different reports where a slight decrease in the frequency of HLA-B*58 was observed [6, 8, 28]. Both B57 and B58 have lower worldwide frequencies than B51 and B35; therefore, the statistical power to detect associations with the same ORs is worse. Noticeably, the only patient bearing B*58 had also B*51. HLA-B51 and B57 have the Bw4 epitope which is the ligand of the natural killer (NK) cell inhibitory receptor KIR3DL1. Interactions between Bw4 and KIR3DL1 have been proposed as a possible underlying mechanism to explain the relationship between HLA-B and BD. However, a large number of common HLA-B alleles encoding molecules with the Bw4 epitope (for example, B*44, B*49) have never been associated with BD. In this regard, B*35, protective according to our results, encodes molecules without the Bw4 epitope, but B*58 (which our study suggests is protective) encodes HLA molecules with this epitope. It has been described that changes at specific positions of the HLA-B molecule affect the interaction of Bw4 with KIR3DL1 . Among them, position 97 (α2 domain) is very interesting because most of the HLA-B alleles with the Bw4 epitope (including B*58) encode proteins with 97R or 97S. The only exceptions are proteins encoded by alleles of the groups HLA-B*13, HLA-B*27, HLA-B*51, HLA-B*52 and HLA-B*57 . From among these groups, B*51 and B*57 have been associated with BD, and HLA-B*27 was not found to be associated with the disease in our cohort. However, HLA-B*27 has been associated with BD in other studies [8, 10]. The remaining groups (B*13 and B*52) have not been consistently associated with BD. HLA-B*51 and HLA-B*52 differ in only two amino acid positions, one of which is position 67 (α1 domain) (F67S). Changes in position 67 also affect the interaction of Bw4 with KIR3DL1 . The amino acid in position 67 is shared between B*52 and B*13 . Therefore, a possible explanation for the association of HLA-B and BD is that these molecules have, in addition to the Bw4 epitope, amino acid changes in specific positions such as 67 and 97, which modify interactions with KIR3DL1. Thus, determining which of the HLA-B alleles are involved in the disease could contribute to clarification in the etiology of BD.
The HLA-A region has been associated with the disease in GWASs, as well as in candidate gene studies, but it is not clear whether it is the primary associated locus. The distribution of the HLA-A groups was different in our patient and control cohorts. Nevertheless, in the univariate analysis, none of the differences reach the significance level necessary to be considered as associated, although three of them were suggestive of an association, two as risk factors and one as a protective factor. Moreover, HLA-A*24 was independently associated according to the logistic regression analysis and the analysis recommended by Svejgaard and Ryder . Our results, as well as the previously published GWASs, suggest that the weight of the HLA-A region in the susceptibility to BD is weaker than that of the HLA-B region. Relatively few studies have investigated the influence of the HLA-A polymorphism in BD, and the associations described have not been consistent. HLA-A*26 has been described as associated with BD in Asian populations [5, 16], but only one study in a Greek European population with a very limited number of patients reported an association with HLA-A26 . We did not detect differences in the distribution of A*26 between patients and controls. Interestingly, A*24, a suggested risk factor in our work, is one of the few HLA-A specificities with the Bw4 epitope. Taking our data and those reported in other studies to date, we think that, contrary to the situation with HLA-B, a worldwide association of HLA-A with BD is unlikely to be found. Additionally, it is not completely clear whether other genes located near the HLA-A locus, such as HLA-E (ligand of the NK receptor CD94:NKG2A) could be primarily responsible for the association observed. In this way, HLA-E is highly selective in binding peptides derived from HLA class I leader sequences, and this condition regulates its expression on the cell surface . Thus, polymorphisms located in the leader sequence could influence the expression of HLA-E in the cell surface and, as a consequence, the NK activity. Further studies are needed to clarify which HLA-A allele groups are associated in each ethnic group as well as the role of the HLA-A region in BD.
Regarding the bases of the HLA association with BD, a study in which three different peaks with independent effects were identified and replicated in the region was published while our manuscript was under review. That study also suggested that the association with the HLA-B*51 could be due to a noncoding variant located between HLA-B and MICA . The results of the study by Hughes et al. are in agreement with the findings of other previously reported studies . Nevertheless, studies designed to assess the association between MICA and BD have had conflicting results, although most of them support the association of MICA being due to linkage disequilibrium with HLA-B51. Moreover, different GWASs have not found an independent association between MICA and BD [16, 19, 20]. In the study by Hughes et al., HLA typing was deduced by imputation from SNPs. Imputation is a valuable tool for identifying and fine-mapping association between one region and one disease, as well as to extract much more information from existing GWAS data sets. Nevertheless, the accuracy of the imputation is not 100%, depending on the ethnicity of the population and the HLA locus under study . Therefore, it would be interesting to compare these results where the HLA typing was assigned by imputation with those obtained by HLA typing to further clarify this question.
Regarding the non-HLA genes included in the study, IL23R is the most strongly associated in our population because the three SNPs within this region were associated with BD, with ORs ranging from 1.36 to 1.59. One of the SNPs located in IL10 was also significantly associated. According to our sample size, we can detect associations with ORs of about 1.3. The SNPs rs11206377 (intergene), rs3024498 (IL10), rs1800872 (IL10) and rs2061634 (TDR7) had ORs between 1.1 and 1.3 in our cohort; therefore, our study is underpowered in these cases. The rest of the SNPs included do not seem to greatly influence susceptibility to BD in our population, because they have ORs lower than 1.1. Specifically, rs799348, a functional SNP in UBAC2 which has been strongly associated with BD in the Turkish population, has an OR of 1.0 in our population. IL23R and IL10 genes have been found to be associated to BD in GWASs as well as in candidate gene studies, and this association has been replicated in other populations [19, 20, 22, 35]. Polymorphisms in IL23R have been associated with Crohn’s disease, rheumatoid arthritis and ankylosing spondylitis, and variants of IL10 have been associated with ulcerative colitis and systemic lupus erythematosus. Therefore, these genes could be involved in a common pathway of susceptibility to several autoimmune and autoinflammatory diseases [36–38].
Our study demonstrates an association of BD with B*57 independently of B*51, at least in European populations. Regarding HLA-A, taking our data together with those reported to date, a worldwide association of HLA-A with this disease is unlikely to be found. With respect to non-HLA genes, IL-23R and IL10 played a role in this disease in our population.
Utah residents with Northern and Western European ancestry from the Centre d’Etude du Polymorphisme Humain collection
Genomewide association study
Human leukocyte antigen
Minor allele frequency
Multifactor dimensionality reduction
Single nucleotide polymorphism
This work was supported by Fondo de Investigaciones Sanitarias (10/1701), Fondos FEDER, Plan Andaluz de Investigación (CTS-0197 and CTS-180), Red Enfermedades Inflamatorias y Reumáticas RD08/0075/0013 and Consejería de Salud de la Junta de Andalucía (PI0411/2010). LOF is the recipient of a fellowship (FI11/00547). The authors thank Asociación Andaluza de Enfermedades Autoinmunes (AADEA) and all patients and donors enrolled in the present study for their cooperation. Also, the authors thank the following public Spanish institutions for their participation in this study: Servicio Andaluz de Salud (Hospitales Torrecárdenas, Almería; Clínico San Cecilio, Granada; Carlos Haya, Málaga, and Virgen del Rocío and Valme, Sevilla), Servei Catalá de la Salut (Clinic and Vall d’Hebron hospitals, Barcelona), Servizo Galego de Saude (Complejo Hospitalario Universitario, A Coruña), Servicio Cántabro de Salud (Hospital Marqués de Valdecilla, Santander), Servicio Navarro de Salud (Hospital Virgen del Camino, Pamplona), Servicio Madrileño de Salud (Hospital de la Princesa, Madrid) and Consejo Superior de Investigaciones Científicas (IPB López Neyra, Granada).
- Mendes D, Correia M, Barbedo M, Vaio T, Mota M, Gonçalves O, Valente J: Behçet’s disease—a contemporary review. J Autoimmun. 2009, 32: 178-188. 10.1016/j.jaut.2009.02.011.View ArticlePubMedGoogle Scholar
- Koné-Paut I, Geisler I, Wechsler B, Ozen S, Ozdogan H, Rozenbaum M, Touitou I: Familial aggregation in Behçet’s disease: high frequency in siblings and parents of pediatric probands. J Pediatr. 1999, 135: 89-93. 10.1016/S0022-3476(99)70333-1.View ArticlePubMedGoogle Scholar
- Sakane T, Takeno M, Suzuki N, Inaba G: Behçet’s disease. N Engl J Med. 1999, 341: 1284-1291. 10.1056/NEJM199910213411707.View ArticlePubMedGoogle Scholar
- de Menthon M, LaValley MP, Maldini C, Guillevin L, Mahr A: HLA-B51/B5 and the risk of Behçet’s disease: a systematic review and meta-analysis of case–control genetic association studies. Arthritis Rheum. 2009, 61: 1287-1296. 10.1002/art.24642.View ArticlePubMedGoogle Scholar
- Itoh Y, Inoko H, Kulski JK, Sasaki S, Meguro A, Takiyama N, Nishida T, Yuasa T, Ohno S, Mizuki N: Four-digit allele genotyping of the HLA-A and HLA-B genes in Japanese patients with Behçet’s disease by a PCR-SSOP-Luminex method. Tissue Antigens. 2006, 67: 390-394. 10.1111/j.1399-0039.2006.00586.x.View ArticlePubMedGoogle Scholar
- Choukri F, Chakib A, Himmich H, Hüe S, Caillat-Zucman S: HLA-B*51 and B*15 alleles confer predisposition to Behçet’s disease in Moroccan patients. Hum Immunol. 2001, 62: 180-185. 10.1016/S0198-8859(00)00249-4.View ArticlePubMedGoogle Scholar
- Ahmad T, Wallace GR, James T, Neville M, Bunce M, Mulcahy-Hawes K, Armuzzi A, Crawshaw J, Fortune F, Walton R, Stanford MR, Welsh KI, Marshall SE, Jewell DP: Mapping the HLA association in Behçet’s disease: a role for tumor necrosis factor polymorphisms?. Arthritis Rheum. 2003, 48: 807-813. 10.1002/art.10815.View ArticlePubMedGoogle Scholar
- Gül A, Uyar FA, Inanç M, Ocal L, Barrett JH, Aral O, Koniçe M, Saruhan-Direskeneli G: A weak association of HLA-B*2702 with Behçet’s disease. Genes Immun. 2002, 3: 368-372. 10.1038/sj.gene.6363863.View ArticlePubMedGoogle Scholar
- Bettencourt A, Pereira C, Carvalho L, Carvalho C, Patto JV, Bastos M, Silva AM, Barros R, Vasconcelos C, Paiva P, Costa L, Costa PP, Mendonça D, Correia J, Silva BM: New insights of HLA class I association to Behçet’s disease in Portuguese patients. Tissue Antigens. 2008, 72: 379-382. 10.1111/j.1399-0039.2008.01087.x.View ArticlePubMedGoogle Scholar
- Radouane A, Oudghiri M, Chakib A, Naya A, Belhouari A, El Malki A, Bennani S: HLA-B*27 allele associated to Behçet’s disease and to anterior uveitis in Moroccan patients. Ann Biol Clin (Paris). 2011, 69: 419-424.Google Scholar
- Arber N, Klein T, Meiner Z, Pras E, Weinberger A: Close association of HLA-B51 and B52 in Israeli patients with Behçet’s syndrome. Ann Rheum Dis. 1991, 50: 351-353. 10.1136/ard.50.6.351.PubMed CentralView ArticlePubMedGoogle Scholar
- Soto-Vega E, García-Muñoz R, Richaud-Patin Y, Zúñiga-Ramos J, Crispín JC, Díaz-Jouanen E, Flores-Suárez LF, Llorente L, Granados J: Class I and class II MHC polymorphisms in Mexican patients with Behçet’s disease. Immunol Lett. 2004, 93: 211-215. 10.1016/j.imlet.2004.03.017.View ArticlePubMedGoogle Scholar
- Hughes T, Coit P, Adler A, Yilmaz V, Aksu K, Düzgün N, Keser G, Cefle A, Yazici A, Ergen A, Alpsoy E, Salvarani C, Casali B, Kötter I, Gutierrez-Achury J, Wijmenga C, Direskeneli H, Saruhan-Direskeneli G, Sawalha AH: Identification of multiple independent susceptibility loci in the HLA region in Behçet’s disease. Nat Genet. 2013, 45: 319-324. 10.1038/ng.2551.View ArticlePubMedGoogle Scholar
- Yazici H, Fresko I, Yurdakul S: Behçet’s syndrome: disease manifestations, management, and advances in treatment. Nat Clin Pract Rheumatol. 2007, 3: 148-155. 10.1038/ncprheum0436.View ArticlePubMedGoogle Scholar
- Ohno S, Ohguchi M, Hirose S, Matsuda H, Wakisaka A, Aizawa M: Close association of HLA-Bw51 with Behçet’s disease. Arch Ophthalmol. 1982, 100: 1455-1458. 10.1001/archopht.1982.01030040433013.View ArticlePubMedGoogle Scholar
- Meguro A, Inoko H, Ota M, Katsuyama Y, Oka A, Okada E, Yamakawa R, Yuasa T, Fujioka T, Ohno S, Bahram S, Mizuki N: Genetics of Behçet disease inside and outside the MHC. Ann Rheum Dis. 2010, 69: 747-754. 10.1136/ard.2009.108571.View ArticlePubMedGoogle Scholar
- Karasneh J, Gül A, Ollier WE, Silman AJ, Worthington J: Whole-genome screening for susceptibility genes in multicase families with Behçet’s disease. Arthritis Rheum. 2005, 52: 1836-1842. 10.1002/art.21060.View ArticlePubMedGoogle Scholar
- Fei Y, Webb R, Cobb BL, Direskeneli H, Saruhan-Direskeneli G, Sawalha AH: Identification of novel genetic susceptibility loci for Behçet’s disease using a genome-wide association study. Arthritis Res Ther. 2009, 11: R66-10.1186/ar2695.PubMed CentralView ArticlePubMedGoogle Scholar
- Remmers EF, Cosan F, Kirino Y, Ombrello MJ, Abaci N, Satorius C, Le JM, Yang B, Korman BD, Cakiris A, Aglar O, Emrence Z, Azakli H, Ustek D, Tugal-Tutkun I, Akman-Demir G, Chen W, Amos CI, Dizon MB, Kose AA, Azizlerli G, Erer B, Brand OJ, Kaklamani VG, Kaklamanis P, Ben-Chetrit E, Stanford M, Fortune F, Ghabra M, Ollier WE: Genome-wide association study identifies variants in the MHC class I, IL10, and IL23R-IL12RB2 regions associated with Behçet’s disease. Nat Genet. 2010, 42: 698-702. 10.1038/ng.625.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizuki N, Meguro A, Ota M, Ohno S, Shiota T, Kawagoe T, Ito N, Kera J, Okada E, Yatsu K, Song YW, Lee EB, Kitaichi N, Namba K, Horie Y, Takeno M, Sugita S, Mochizuki M, Bahram S, Ishigatsubo Y, Inoko H: Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behçet’s disease susceptibility loci. Nat Genet. 2010, 42: 703-706. 10.1038/ng.624.View ArticlePubMedGoogle Scholar
- International Study Group for Behçet’s Disease: Criteria for diagnosis of Behçet’s disease. Lancet. 1990, 335: 1078-1080.Google Scholar
- Jiang Z, Yang P, Hou S, Du L, Xie L, Zhou H, Kijlstra A: IL-23R gene confers susceptibility to Behçet’s disease in a Chinese Han population. Ann Rheum Dis. 2010, 69: 1325-1328. 10.1136/ard.2009.119420.View ArticlePubMedGoogle Scholar
- International HapMap Project: HapMap DataRel 24 database.http://hapmap.ncbi.nlm.nih.gov/cgi-perl/gbrowse/hapmap24_B36/,
- Svejgaard A, Ryder LP: HLA and disease associations: detecting the strongest association. Tissue Antigens. 1994, 43: 18-27. 10.1111/j.1399-0039.1994.tb02291.x.View ArticlePubMedGoogle Scholar
- Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH: Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001, 69: 138-147. 10.1086/321276.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore JH: Computational analysis of gene-gene interactions using multifactor dimensionality reduction. Expert Rev Mol Diagn. 2004, 4: 795-803. 10.1586/14737184.108.40.2065.View ArticlePubMedGoogle Scholar
- Computational Genetics Laboratory, Institute for Quantitative Biomedical Sciences, Geisel School of Medicine and Dartmouth College: Multifactor Dimensionality Reduction (MDR) version 2 open-source software package.http://www.multifactordimensionalityreduction.org/,
- Kera J, Mizuki N, Ota M, Katsuyama Y, Pivetti-Pezzi P, Ohno S, Inoko H: Significant associations of HLA-B*5101 and B*5108, and lack of association of class II alleles with Behçet’s disease in Italian patients. Tissue Antigens. 1999, 54: 565-571. 10.1034/j.1399-0039.1999.540605.x.View ArticlePubMedGoogle Scholar
- Sanjanwala B, Draghi M, Norman PJ, Guethlein LA, Parham P: Polymorphic sites away from the Bw4 epitope that affect interaction of Bw4+ HLA-B with KIR3DL1. J Immunol. 2008, 181: 6293-6300.PubMed CentralView ArticlePubMedGoogle Scholar
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI): IMGT/HLA Sequence Alignment Tool.http://www.ebi.ac.uk/imgt/hla/align.html,
- Mizuki N, Ohno S, Ando H, Chen L, Palimeris GD, Stavropoulos-Ghiokas E, Ishihara M, Goto K, Nakamura S, Shindo Y, Isobe K, Ito N, Inoko H: A strong association between HLA-B*5101 and Behçet’s disease in Greek patients. Tissue Antigens. 1997, 50: 57-60. 10.1111/j.1399-0039.1997.tb02835.x.View ArticlePubMedGoogle Scholar
- Braud VM, Allan DS, O’Callaghan CA, Söderström K, D’Andrea A, Ogg GS, Lazetic S, Young NT, Bell JI, Phillips JH, Lanier LL, McMichael AJ: HLA-E binds to natural killer cell receptors CD94/NKG2A, B and C. Nature. 1998, 391: 795-799. 10.1038/35869.View ArticlePubMedGoogle Scholar
- Ota M, Mizuki N, Katsuyama Y, Tamiya G, Shiina T, Oka A, Ando H, Kimura M, Goto K, Ohno S, Inoko H: The critical region for Behçet disease in the human major histocompatibility complex is reduced to a 46-kb segment centromeric of HLA-B, by association analysis using refined microsatellite mapping. Am J Hum Genet. 1999, 64: 1406-1410. 10.1086/302364.PubMed CentralView ArticlePubMedGoogle Scholar
- de Bakker PI, Raychaudhuri S: Interrogating the major histocompatibility complex with high-throughput genomics. Hum Mol Genet. 2012, 21 (R1): R29-R36. 10.1093/hmg/dds384.PubMed CentralView ArticlePubMedGoogle Scholar
- Xavier JM, Shahram F, Davatchi F, Rosa A, Crespo J, Abdollahi BS, Nadji A, Jesus G, Barcelos F, Patto JV, Shafiee NM, Ghaderibarim F, Oliveira SA: Association study of IL10 and IL23R-IL12RB2 in Iranian patients with Behçet’s disease. Arthritis Rheum. 2012, 64: 2761-2772. 10.1002/art.34437.View ArticlePubMedGoogle Scholar
- Lees CW, Barrett JC, Parkes M, Satsangi J: New IBD genetics: common pathways with other diseases. Gut. 2011, 60: 1739-1753. 10.1136/gut.2009.199679.View ArticlePubMedGoogle Scholar
- Franke A, Balschun T, Karlsen TH, Sventoraityte J, Nikolaus S, Mayr G, Domingues FS, Albrecht M, Nothnagel M, Ellinghaus D, Sina C, Onnie CM, Weersma RK, Stokkers PC, Wijmenga C, Gazouli M, Strachan D, McArdle WL, Vermeire S, Rutgeerts P, Rosenstiel P, Krawczak M, Vatn MH, Mathew CG, Schreiber S, the IBSEN study group: Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility. Nat Genet. 2008, 40: 1319-1323. 10.1038/ng.221.View ArticlePubMedGoogle Scholar
- Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, Ortmann W, Kosoy R, Ferreira RC, Nordmark G, Gunnarsson I, Svenungsson E, Padyukov L, Sturfelt G, Jönsen A, Bengtsson AA, Rantapää-Dahlqvist S, Baechler EC, Brown EE, Alarcón GS, Edberg JC, Ramsey-Goldman R, McGwin G, Reveille JD, Vilá LM, Kimberly RP, Manzi S, Petri MA, Lee A, Gregersen PK: A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet. 2009, 41: 1228-1233. 10.1038/ng.468.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.