ORIGINAL RESEARCH article
Sec. Autoimmune and Autoinflammatory Disorders
Association of Functional Polymorphisms of KIR3DL1/DS1 With Behçet's Disease
- 1Department of Immunology, Hospital Universitario Virgen del Rocío, Instituto de Biomedicina de Sevilla (IBiS), CSIC, Universidad de Sevilla, Seville, Spain
- 2Department of Internal Medicine, Hospital Clínico San Cecilio, Granada, Spain
- 3Department of Internal Medicine, Hospital Universitario Virgen del Rocío, Seville, Spain
- 4Department Autoimmune Diseases, Hospital Universitari Clínic, Barcelona, Spain
- 5Department of Rheumatology, Complejo Hospitalario Universitario A Coruña, A Coruña, Spain
- 6Department of Rheumatology, Hospital Universitario de Valme, Seville, Spain
- 7Department of Immunology, Hospital Universitari Son Espases, Palma de Mallorca, Spain
- 8Department of Internal Medicine, Autoimmune Systemic Diseases Unit, Hospital Vall d'Hebron, Universidad Autonoma de Barcelona, Barcelona, Spain
- 9Department of Rheumatology, Hospital Universitario Marqués de Valdecilla, Santander, Spain
- 10Department of Internal Medicine, Complejo Hospitalario Torrecárdenas, Almería, Spain
- 11Department of Internal Medicine, Hospital Universitario Central de Asturias, Asturias, Spain
- 12Department of Internal Medicine, Hospital Virgen del Camino, Pamplona, Spain
- 13Deparment of Internal Medicine, Hospital Universitari Mútua Terrassa, Terrassa, Spain
- 14Department of Rheumatology, Hospital Clínico San Carlos, Madrid, Spain
- 15Department of Internal Medicine, Hospital Regional Universitario de Málaga, Málaga, Spain
- 16Department of Rheumatology, Hospital de la Princesa, IIS-Princesa, Madrid, Spain
- 17Department of Rheumatology, Hospital Universitario Doctor Peset, Valencia, Spain
- 18Instituto de Parasitología y Biomedicina “López-Neyra”, CSIC, PTS Granada, Granada, Spain
Behçet's disease (BD) is an immune-mediated vasculitis related to imbalances between the innate and adaptive immune response. Infectious agents or environmental factors may trigger the disease in genetically predisposed individuals. HLA-B51 is the genetic factor stronger associated with the disease, although the bases of this association remain elusive. NK cells have also been implicated in the etiopathogenesis of BD. A family of NK receptors, Killer-cell Immunoglobulin-like Receptor (KIR), with a very complex organization, is very important in the education and control of the NK cells by the union to their ligands, most of them, HLA class I molecules. This study aimed to investigate the contribution of certain KIR functional polymorphisms to the susceptibility to BD. A total of 466 BD patients and 444 healthy individuals were genotyped in HLA class I (A, B, and C). The set of KIR genes and the functional variants of KIR3DL1/DS1 and KIR2DS4 were also determined. Frequency of KIR3DL1*004 was lower in patients than in controls (0.15 vs. 0.20, P = 0.005, Pc = 0.015; OR = 0.70; 95% CI 0.54–0.90) in both B51 positive and negative individuals. KIR3DL1*004, which encodes a misfolded protein, is included in a common telomeric haplotype with only one functional KIR gene, KIR3DL2. Both, KIR3DL1 and KIR3DL2 sense pathogen-associated molecular patterns but they have different capacities to eliminate them. The education of the NK cells depending on the HLA, the balance of KIR3DL1/KIR3DL2 licensed NK cells and the different capacities of these receptors to eliminate pathogens could be involved in the etiopathogenesis of BD.
Behçet's disease (BD) [OMIM #109650] is a rare, chronic and systemic vasculitis characterized by recurrent oral and genital ulcers, although other clinical manifestations, such as skin lesions, ocular, gastrointestinal, and neurological disorders are relatively common. This chronic disorder is an immune-mediated disease in which imbalances between the innate and adaptive immune response triggered by infectious agents or environmental factors in genetically predisposed individuals have been suggested as the underlying mechanisms of the disease (1). Among the genetic factors, HLA-B51 has been associated with the disease in many populations (2); more recently, other genes related to the immune system such as IL23R, IL10, STAT4, CCR1, CCR3, KLRC4, ERAP1, TNFAIP3, and FUT2 have been associated with this disease (3).
Natural killer (NK) cells are lymphocytes of the innate immune system with several activators and inhibitor receptors on their surface. Among them, the Killer-cell Immunoglobulin-like Receptor (KIR) gene family consisted of a set of 15 genes: KIR2DL1, KIR2DL2, KIR2DL3, KIR2DL4, KIR2DL5A KIR2DL5B, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DL1, KIR3DL2, KIR3DL3, and KIR3DS1 and 2 pseudogenes: KIR2DP1 and KIR3DP1 located near to the leukocyte receptor complex in the chromosome 19. The KIR gene content is variable among individuals being a major contributor to the KIR diversity in the human population and the NK cell repertoires (4). At present, 625 genotypes (different sets of KIR genes) named by numbers assigned sequentially have been found in the human population1. Based on their gene content, two kinds of KIR haplotype groups, A and B (that give rise to the genotypes AA and Bx), have been described. Although there is no single specific criterion to distinguish them, from a practical point of view, one individual is considered as Bx when its set of KIR genes included at least one of the following genes: 2DL2, 2DL5, 3DS1, 2DS1, 2DS2, 2DS3, 2DS5, and it is assigned as AA when all these genes are absent. The KIR genes encoded activator or inhibitor molecules, in general, the molecules with long cytoplasmic tails (L) contains two immune tyrosine-based inhibitory motifs (ITIM) which transduce inhibitory signals to the NK cell, whereas the molecules with short cytoplasmic tails (S) possess a positively charged amino acid residue in their transmembrane region which allows them to associate with a DAP12 molecule generating an activation signal (4). The A Haplotypes have fewer activator genes than haplotypes B, in fact, the only two activator genes present in AA individuals are KIR2DS4 and KIR2DL4, but only KIR2DS4 binds HLA classical molecules. Besides, each KIR gene has a variable number of alleles, some of them with functional significance. Indeed, a different expression in the cell membrane of diverse KIR3DL1 alleles, ranging from practically null (in KIR3DL1*004, which encodes a misfolded receptor mostly retained inside the cell) to high expression, or production of soluble forms in the case of KIR2DS4 (5, 6).
The KIR-ligands are HLA class I molecules and the interaction between the pair ligand-receptor is fundamental for the regulation of NK cell activity. All the classical HLA-class I molecules (HLA-A, B and C) present peptides to the CD8+ T cells, however, not all of them are ligands of the KIR molecules. In this sense, all the HLA-C molecules are ligands of KIR and they can be grouped as C1 (Asn80) or C2 (Lys80) which are ligands of KIR2DL2 and KIR2DL1, respectively. Nevertheless, only some HLA-A and B molecules are KIR ligands, specifically, those bearing the Bw4 epitope are KIR3DL1 ligands and A3 and A11 have been described to be KIR3DL2 ligands (4). HLA-B51 and other HLA molecules associated with BD have the Bw4 epitope and interactions between Bw4 and KIR3DL1 have been proposed as one of the possible underlying mechanisms to explain the relationship between HLA class I and BD (7). However, a large number of common HLA-B alleles encoding molecules with the Bw4 epitope (e.g., B*44, B*49, etc.) have never been associated with this disease, even those with isoleucine at position 80 (Bw4-80I), which are the strongly interacting ligands. The Bw4 epitope spans the residues 77–83 but it has been described that changes at specific positions in the HLA-B molecule outside this epitope, in particular, position 97, affect the interaction of Bw4 with KIR3DL1 and, interestingly, the position 97 of HLA-B has been found strongly associated with susceptibility to BD (8–10).
With this picture, great complexity and a relative lack of the knowledge of the polymorphism of the KIR system and its relationship with its ligands, the objective of this study was to investigate the contribution of certain KIR functional polymorphisms to the susceptibility to BD.
Materials and Methods
The study included a total of 466 BD-unrelated patients (44.2% males) who fulfilled the 1990 International Study Group classification criteria for BD (11), and 444 unrelated healthy individuals (50% males) included as the control group. All the subjects were Spanish European recruited from 17 Spanish hospitals across the country. 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 the following: 100% had oral ulcers, 59.4% genital ulcers, 53.9% uveitis, 42% arthritis, and 21% vascular, 18.2% neurological, 16.3% positive pathergy test, and 15.4% gastrointestinal involvement.
Peripheral blood collected in EDTA tubes was obtained from the healthy controls, and peripheral blood or saliva served as starting material from patients. 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/280 absorbance ratio of 1.7–2.0 and a final concentration of 10–20 ng/μl were considered appropriate. A total of 14 DNA samples from saliva were eliminated because they did not meet these quality criteria.
The HLA class I (A, B, and C) genotyping was carried out using a PCR-SSOP Luminex method, LABType SSO (One Lambda Inc., Canoga Park, CA), according to the manufacturer's instructions. Briefly, target DNAs were PCR-amplified using group-specific primers (HLA-A, -B, or -C) and the biotinylated-PCR products were denatured and hybridized with specific probes bound to colored-coded microspheres. Phycoerythrin conjugated to Streptoavidin was used to label and reveal reactions which were read in a flow analyzer, LABScanTM100, to identify fluorescent intensity on each microsphere. The software HLA fusion 2.0 (One Lambda Inc.) was used to assign the HLA typing of each locus. This method allows a medium resolution genotyping of each class I HLA gene. After class I genotyping, the samples were classified as A3/A11, Bw4, Bw4-80I, C1, and C2.
The set of KIR genes presents in each sample was determined also using a PCR-SSOP Luminex method, KIR SSO Genotyping Test (One Lambda Inc., Canoga Park, CA), according to the manufacturer's instructions. The principle of the method is the same as that employed for HLA-class I genotyping with the corresponding group-specific primers and bound-probes. Also, to distinguish the presence/absence of KIR genes, this method allows discriminating between some groups of alleles, specifically, between the group of KIR2DS4 alleles with the full sequence (KIR2DS4*001/011/014/015) which encode functional proteins and the group with a deletion of 22 base pairs (KIR2DS4*003/004/006/007/008/009/010/012/013) which encode non-functional variants (6).
KIR3DL1/S1 individuals were genotyped in the rs149123986 by real-time PCR using TaqMan SNP Genotyping Assays (Applied Biosystems, Barcelona, Spain) in a LightCycler 480 (Roche, Barcelona, Spain). This single nucleotide polymorphism (SNP, A/G) allows to distinguish two groups of KIR3DL1 alleles with different expression patterns: the KIR3DL1*004 (rs149123986G) with a very low or null membrane expression (3DL1Null) and the rest of alleles (rs149123986A), which have a normal although variable membrane expression (3DL1Exp)2.
Phenotypic and genotypic frequencies were estimated by direct counting and distributions were compared using the χ2 test. This test was applied to check the association of (a) different ligands: A3/11, B51, Bw4, Bw4-80I, C1, and C2; (b) the KIR genes, haplotypes and genotypes; (c) the KIR3DL1/DS1 variants in the full group and the subgroups B51 Bw4, Bw4-80I; and (d) the KIR2DS4 variants in the total groups and among individuals AA and Bx and stratified according to the presence/absence of the combination ligand and receptor. The P-values were corrected by Bonferroni's adjustment (Pc) considering the number of tests in each case. Those comparisons with Pc-values <0.05 were considered associated whereas those with P < 0.05 but with Pc > 0.05 were considered suggestive of association. The odds ratios (ORs) and 95% confidence intervals (95% CI) were calculated using the web software OpenEpi (Open Source Epidemiologic Statistics for Public Health, Versión 3.01. www.OpenEpi.com).
HLA Molecules KIR-Ligands
A total of 435 BD patients (96.2%) and 439 controls (98.9%) were fulfilled genotyped and included in the statistical analysis. The distribution of different HLA molecules KIR-ligands is displayed in Table 1. The strongest associations were: HLA-B51, conferring risk and A3/11, conferring protection. The frequency of A3/11 was decreased in both patient groups, B51 positive and negative (18 and 21%, respectively, vs. 26 and 33% in their corresponding controls), although statistical significance was reached only in the negative group. Bw4 and Bw4-80I confer risk but with a lower OR than B51. Among the group B51 negative (250 patients and 370 controls), the distribution of Bw4 and Bw4-80I was not significantly different in patients and controls (Bw4 76.8% in patients vs. 76.2 in controls; Bw4-80I 56.0% in patients vs. 56.2% in controls, P > 0.05 in both cases). No differences in the distribution of C1 and C2 were detected in any case.
The distribution of the KIR genes in our cohort of patients and controls is summarized in Table 2. All the KIR genes and pseudogenes were detected in this study and four of them, corresponding to the centromeric (3DL3 and 3DP1) and telomeric (2DL4 and 3DL2) framework genes, were found in all the samples. The distribution of genes was similar in patients and controls and statistically significant differences were not reached in any case. Regarding the distribution of the KIR haplotypes, Bx was the most common in both groups and no significant differences between patients and controls were observed. Concerning the KIR genotypes, a total of 51 genotypes were found in our cohort, but only 12 had a frequency upper 2% in at least one of the groups. The genotype 2 was slightly less frequently represented in patients than in controls (8.7 vs. 14.1%, respectively, P = 0.01, Pc > 0.05; OR 0.58; 95% CI 0.38–0.89).
Functional KIR Polymorphyms
Data regarding the polymorphism of KIR3DL1/S1 are displayed in Table 3. The distribution of individuals KIR3DL1/S1 was not significantly different in patients and controls (KIR3DL1+DS1+ 35.6 vs. 39.9%, KIR3DL1+DS1– 58.3 vs. 56.3%, KIR3DL1-DS1+ 6.0 vs. 3.9%, respectively). Nevertheless, the distribution of the genotypes of rs149123986 was significantly different in both groups (p = 0.01 in 2 × 3 contigency Table) with an over-representation of AA in patients (73.5% vs. 64.3 in controls, Pc = 0.006; OR = 1.54; 95% CI 1.16–2.07) and a down-representation of AG (23.0 vs. 31.7%, Pc = 0.015, OR = 0.65; 95% CI 0.48–0.88). The distribution of the individuals in the six possible groups taking into account both, the KIR3DL1/S1 and the rs149123986, was different in patients and controls (p = 0.002 in 2 × 6 Table) with down-representation of 3DL1Null/3DL1 (6.2 vs. 11.4%, p = 0.007, Pc = 0.042; OR = 0.51; 95% CI 0.32–0.84) and a trend to an over-representation of 3DL1Exp/3DL1Exp (38.1 vs. 31.9%, p = 0.052) individuals among patients. The distribution of the three allele groups: 3DL1Exp, 3DL1Nul, and 3DS1 was different in patients and control (p = 0.01, in 2 × 3 contigency Table) having the group of patients a decreased frequency of 3DL1Null (0.15 vs. 0.20, P = 0.005, Pc = 0.015; OR = 0.70; 95% CI 0.54–0.90) and a slight increased frequency of 3DL1Exp (0.61 vs. 0.56, P = 0.03, Pc > 0.05). Results were similar among B51 positive and negative individuals with a suggestive down-representation of 3DL1Null/3DLS1 in both groups (B51positive: 5.9% in patients vs. 15.9% in controls P = 0.009 Pc > 0.05; OR = 0.34; 95% CI 0.13–0.83 and B51 negative: 6.4% in patients vs. 10.5% in controls P = 0.04 Pc > 0.05; OR = 0.58; 95% CI 0.3–1.05) (Table 4).
Table 4. Frequency of the KIR3DL1 functional polymorphism in BD patients and controls stratified according to the HLA-B51.
Regarding the deletion of 22 bp of 2DS4, differences in the distribution of individuals with at least one copy of the 2DS4Full between patients and controls were not detected (42.7 vs. 37.5%, P = 0.13). Since 2DS4 is the only activator binding classical HLA molecules in AA individuals, patients and controls were stratified according to their haplotype group in AA or Bx but no differences were detected (Data not showed).
The main finding of this study is the association of the KIR3DL1*004 with susceptibility to BD. This allele is protective in the development of the disease and, according to our results, its effect is independent of HLA-B51.
Previous studies analyzed the relationship between the distribution of the KIR genes and the disease. Similarly to our results, no differences between patients and controls have been reported in these studies (12–15). Only one of these studies, performed in the Iranian population, investigated KIR genotypes. That study reported a different distribution of various KIR genotypes, including genotype 2, although the direction of the association was contrary to that found in the present study (15).
KIR3DL1/DS1 has focused on interest because of its function as the receptor of the Bw4 epitope. The distribution of KIR3DL1 and 3DS1 in BD was evaluated in a large cohort of Turkish. Similarly to our results, no association was found in the whole cohort neither in subgroups with B51, Bw4, or Bw4-80I (16) but this study did not investigate variants on this gene. Variants in KIR3DL1 have been reported as associated with the disease in two high throughput studies (10, 17), although none of them was designed to investigate the possible involvement of functional variants in this disease. Very recently, during the submission of this paper, a study that addresses this question in a smaller cohort composed of ethnically-mixed patients and controls was published (18). The main conclusion of that study is the same as ours, an association of the functional alleles of KIR3DL1, independent of HLA-B51 with the disease. Nevertheless, they are some discrepancies in the associated genotypes between these two studies. When the allele frequencies of the study by Petrushkin et al. are deduced of their data, results are similar to ours regarding the allele KIR3DL1*004, which is underrepresented in patients (with similar OR), although the statistical significance is only reached in our cohort. Nevertheless, we have opposite results concerning KIR3DL1Exp and KIR3DS1. In our cohort, the alleles that encode KIR3DL1Exp are significantly overrepresented among patients whereas, in its cohort, it is KIR3DS1 which tends to be overrepresented among patients. In our opinion, the distribution of the allelic frequencies reinforces the idea of association with the allele KIR3DL1*004. Discrepancies could be due to different reasons, among them, methods used in KIR3DL1 genotyping, approaches, stratification of cohorts with a relatively small number of individuals and ethnical differences. Besides, there are previous studies in other class-I pathologies such as psoriatic disease (PD) and ankylosing spondylitis (AS), in which a protective association of the KIR3DL1*004 allele has been reported (19–21).
KIR3DL1*004 encodes a misfolded receptor which is retained in the cell, although it has been described that a small quantity of folded protein can be exported to the membrane (22). One possible explanation to the protective effect of KIR3DL1*004 in these pathologies is that this misfolded molecule retains intracellularly the corresponding class I-risk molecules. In this case, differences in the distribution of KIR3DL1*004 among patients with and without the risk molecule would be expected, nevertheless, no differences were found in the present study neither in other pathologies (20, 21, 23) and no evidence supporting the intracellular interaction between KIR3DL1 and MHC class I molecules have been found (5). Interestingly, KIR3DL1*004 is within a common telomeric haplotype Tel-A1-V which, also, bears two alleles encoding secreted receptors, KIR2DL4*008 and KIR2DS4*006 and another that encodes a membrane receptor KIR3DL2*003/005/011/020 (24)3. Therefore, this Tel-A1-V haplotype, protective in BD, carries non-functional alleles with the only exception of KIR3DL2. Consequently, KIR3DL1*004 could be a marker of the Tel-A1-V haplotype and the protective association found in our study could be explained by the fact that KIR3DL2 is the only functional gene in this haplotype.
KIR3DL2 is an inhibitor framework gene, found in 100% of the population and it has many allelic forms encoding many different proteins (161 and 111, respectively, in the last database), the relative expression levels of these different proteins are unknown. The KIR3DL2 HLA-ligands are A3 and A11 and there is some evidence that KIR3DL2 can bind these HLA molecules in vivo (25). The education of the NK cells is based on the interaction between the inhibitory receptors and the self-MHC-I molecules, only those individuals bearing the HLA-ligand molecules have licensed NK cells that express the corresponding receptor and therefore, only the individuals A3 and A11 have mature NK cells expressing KIR3DL2. Noticeably, similarly to the findings in the present cohort, different studies in several populations have found a protective association of HLA-A3 with BD, independent of HLA-B51 (9, 10, 26). This protective effect of HLA-A3, not yet explained, could be based on the presence of licensed KIR3DL2 NK cells in these individuals.
KIR molecules having 3 extracellular domains (KIR3DL1 and KIR3DL2) are pattern recognition receptors (PRR), they sense pathogen-associated molecular patterns (PAMPs), specifically single-stranded microbial DNA molecules containing CpG motifs (CpG-ODN). KIR3DL2 binds CpG-ODN by the D0 domain and the complexes are internalized and recognized by TLR9 located in the endosomal compartment what results in cell activation. Nevertheless, in spite of KIR3DL1 binds CpG-ODN, the complexes are not internalized (27). Thus, one possible explanation to the protective effect of the Tel-A1-V haplotype is that the lack of expression of KIR3DL1 in these individuals permits a better efficiency to eliminate pathogens by KIR3DL2 NK cells. In other words, NK cells of individuals with the KIR3DL1Exp would be lesser efficient in eliminating pathogens resulting in a perpetuation of the release of inflammatory cytokines. In this sense, a protective effect of KIR3DL1*004 against the progression of the VIH infection to AIDS has been described (28). This mechanism would explain why the association is independent of B51 and Bw4 in our study but also psoriasis. Moreover, some HLA-Bw4 molecules, such as B27 (associated with ankylosis spondylitis) and B51 could be more efficient to license KIR3DL1 NK cells and this could be related to its association with inflammatory diseases. In this sense, the activity of ERAP1, an enzyme that has been related in an epistatic way with the HLA molecules involved in these pathologies and that can conditioner the repertoire of peptides binding to them, could modify the role of the diverse HLA-Bw4 molecules in these diseases (29).
Lastly, the Tel-A1-V haplotype includes KIR2DS4*006, which has the deletion of 22 base pairs, but according to our results, this deletion does not influence the disease susceptibility.
In conclusion, our results suggest a protective role of a telomeric KIR haplotype with only one functional gene in BD. Therefore, the education of the NK cells depending on the HLA, the balance of licensed KIR3DL1/KIR3DL2 NK cells and the different capacities of these receptors to eliminate pathogens could be involved in the etiopathogenesis of BD.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
The studies involving human participants were reviewed and approved by CEI de los Hospitales Universitarios Vírgen Macarena-Virgen del Rocío. The patients/participants provided their written informed consent to participate in this study.
M-FG-E: experimental design. ÁC-N: experimental procedures performance. M-AM-C and J-RG-L: experimental support. ÁC-N, M-AM-C and M-FG-E: data analysis and preparation of the manuscript. NO-C, F-JG-H, GE, GG-G, JS-B, M-RJ, RS, RB, A-CB-M, RG, PF, MR-C, LR-R, TC, SC, J-JA-S, and JM: samples and clinical data providers.
This work was supported by Fondo de Investigaciones Sanitarias, Instituto de Salud Carlos III (ISCIII, 13/01118 and 16/01373), Fondos FEDER and Plan Andaluz de Investigación (CTS-0197).
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.
2. 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–96. doi: 10.1002/art.24642
5. Taner SB, Pando MJ, Roberts A, Schellekens J, Marsh SG, Malmberg KJ, et al. Interactions of NK cell receptor KIR3DL1*004 with chaperones and conformation-specific antibody reveal a functional folded state as well as predominant intracellular retention. J Immunol. (2011) 186:62–72. doi: 10.4049/jimmunol.0903657
6. Middleton D, Gonzalez A, Gilmore PM. Studies on the expression of the deleted KIR2DS4*003 gene product and distribution of KIR2DS4 deleted and nondeleted versions in different populations. Hum Immunol. (2007) 68:128–34. doi: 10.1016/j.humimm.2006.12.007
8. 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–300. doi: 10.4049/jimmunol.181.9.6293
9. Ombrello MJ, Kirino Y, de Bakker PI, Gul A, Kastner DL, Remmers EF. Behcet disease-associated MHC class I residues implicate antigen binding and regulation of cell-mediated cytotoxicity. Proc Natl Acad Sci USA. (2014) 111:8867–72. doi: 10.1073/pnas.1406575111
10. Ortiz-Fernández L, Carmona FD, Montes-Cano MA, García-Lozano JR, Conde-Jaldón M, Ortego-Centeno N, et al. Genetic analysis with the immunochip platform in Behçet disease. Identification of residues associated in the HLA class I region and new susceptibility Loci. PLoS ONE. (2016) 11:e0161305. doi: 10.1371/journal.pone.0161305
13. Arayssi TK, El Hajj N, Shamseddine W, Ibrahim G, Nasr J, Sabbagh AS, et al. Killer cell immunoglobulin-like receptor genotypes in Behçet's disease patients: any role for the 3DP1*001/002 pseudogene? Genet Test Mol Biomarkers. (2009) 13:319–24. doi: 10.1089/gtmb.2008.0108
14. Moon SJ, Oh EJ, Kim Y, Kim KS, Kwok SK, Ju JH, et al. Diversity of killer cell immunoglobulin-like receptor genes in uveitis associated with autoimmune diseases: ankylosing spondylitis and Behçet disease. Ocul Immunol Inflamm. (2013) 21:135–43. doi: 10.3109/09273948.2012.754905
15. Mohammad-Ebrahim H, Kamali-Sarvestani E, Mahmoudi M, Beigy M, Karami J, Ahmadzadeh N, et al. Association of killer cell immunoglobulin-like receptor (KIR) genes and their HLA ligands with susceptibility to Behçet's? disease. Scand J Rheumatol. (2018) 47:155–63. doi: 10.1080/03009742.2017.1340510
17. Kim SJ, Lee S, Park C, Seo JS, Kim JI, Yu HG. Targeted resequencing of candidate genes reveals novel variants associated with severe Behçet's uveitis. Exp Mol Med. (2013) 45:e49. doi: 10.1038/emm.2013.101
18. Petrushkin H, Norman PJ, Lougee E, Parham P, Wallace GR, Stanford MR, et al. KIR3DL1/S1 allotypes contribute differentially to the development of Behçet disease. J Immunol. (2019) 203:1629–35. doi: 10.4049/jimmunol.1801178
20. Berinstein J, Pollock R, Pellett F, Thavaneswaran A, Chandran V, Gladman DD. Association of variably expressed KIR3dl1 alleles with psoriatic disease. Clin Rheumatol. (2017) 36:2261–6. doi: 10.1007/s10067-017-3784-5
21. Díaz-Peña R, Vidal-Castiñeira JR, Alonso-Arias R, Suarez-Alvarez B, Vicario JL, Solana R, et al. Association of the KIR3DS1*013 and KIR3DL1*004 alleles with susceptibility to ankylosing spondylitis. Arthritis Rheum. (2010) 62:1000–6. doi: 10.1002/art.27332
22. Pando MJ, Gardiner CM, Gleimer M, McQueen KL, Parham P. The protein made from a common allele of KIR3DL1 (3DL1*004) is poorly expressed at cell surfaces due to substitution at positions 86 in Ig domain 0 and 182 in Ig domain 1. J Immunol. (2003) 171:6640–9. doi: 10.4049/jimmunol.171.12.6640
23. Zvyagin IV, Mamedov IZ, Britanova OV, Staroverov DB, Nasonov EL, Bochkova AG, et al. Contribution of functional KIR3DL1 to ankylosing spondylitis. Cell Mol Immunol. (2010) 7:471–6. doi: 10.1038/cmi.2010.42
25. Augusto DG, O'Connor GM, Lobo-Alves SC, Bass S, Martin MP, Carrington M, et al. Pemphigus is associated with KIR3DL2 expression levels and provides evidence that KIR3DL2 may bind HLA-A3 and A11 in vivo. Eur J Immunol. (2015) 45:2052–60 doi: 10.1002/eji.201445324
26. Takeuchi M, Mizuki N, Meguro A, Ombrello MJ, Kirino Y, Satorius C, et al. Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behçet's disease susceptibility. Nat Genet. (2017) 49:438–43. doi: 10.1038/ng.3786
27. Sivori S, Falco M, Carlomagno S, Romeo E, Soldani C, Bensussan A, et al. A novel KIR-associated function: evidence that CpG DNA uptake and shuttling to early endosomes is mediated by KIR3DL2. Blood. (2010) 116:1637–47. doi: 10.1182/blood-2009-12-256586
29. Guasp P, Barnea E, González-Escribano MF, Jiménez-Reinoso A, Regueiro JR, Admon A, et al. The Behçet's disease-associated variant of the aminopeptidase ERAP1 shapes a low-affinity HLA-B*51 peptidome by differential subpeptidome processing. J Biol Chem. (2017) 292:9680–9. doi: 10.1074/jbc.M117.789180
Keywords: Behçet's disease, HLA, KIR, NK cells, functional polymorphisms
Citation: Castaño-Núñez Á, Montes-Cano M-A, García-Lozano J-R, Ortego-Centeno N, García-Hernández F-J, Espinosa G, Graña-Gil G, Sánchez-Bursón J, Juliá M-R, Solans R, Blanco R, Barnosi-Marín A-C, Gómez de la Torre R, Fanlo P, Rodríguez-Carballeira M, Rodríguez-Rodríguez L, Camps T, Castañeda S, Alegre-Sancho J-J, Martín J and González-Escribano M-F (2019) Association of Functional Polymorphisms of KIR3DL1/DS1 With Behçet's Disease. Front. Immunol. 10:2755. doi: 10.3389/fimmu.2019.02755
Received: 19 August 2019; Accepted: 11 November 2019;
Published: 29 November 2019.
Edited by:Kutty Selva Nandakumar, Southern Medical University, China
Reviewed by:Piotr Kusnierczyk, Ludwik Hirszfeld Institute of Immunology and Experimental Therapy (PAN), Poland
Mirko Spiroski, Scientific Foundation SPIROSKI, North Macedonia
Copyright © 2019 Castaño-Núñez, Montes-Cano, García-Lozano, Ortego-Centeno, García-Hernández, Espinosa, Graña-Gil, Sánchez-Bursón, Juliá, Solans, Blanco, Barnosi-Marín, Gómez de la Torre, Fanlo, Rodríguez-Carballeira, Rodríguez-Rodríguez, Camps, Castañeda, Alegre-Sancho, Martín and González-Escribano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: María-Francisca González-Escribano, email@example.com