- Research article
- Open Access
- Open Peer Review
Association study between KIR polymorphisms and rheumatoid arthritis disease: an updated meta-analysis
BMC Medical Geneticsvolume 20, Article number: 24 (2019)
Currently published studies investigating association between the killer cell immunoglobulin-like receptor (KIR) gene polymorphisms and rheumatoid arthritis (RA) reported inconsistent and contradictory results. Hence, we aim to carry out this comprehensive meta-analysis of all eligible studies meeting the inclusion criteria to achieve precise and comprehensive relationships between genetic variations in KIR gene cluster and risk of RA.
Databases of Medline/PubMed and Scopus were searched to investigate case-control studies prior to May 2018. The associations between KIR gene polymorphisms and RA susceptibility were analyzed by computing the odds ratio (OR) and 95% confidence interval (95% CI) for each study.
A total of 11 comparative case-control studies involving 1847 RA patients and 2409 healthy individuals were included in this meta-analysis. Four significant associations of 2DL3 (OR = 0.591, 95% CI = 0.351–0.994; P = 0.047), 2DL5 (OR = 0.716, 95% CI = 0.601–0.853; P < 0.001), 2DS5 (OR = 0.623, 95% CI = 0.393–0.988; P = 0.045), and 3DL3 (OR = 0.324, 95% CI = 0.129–0.814; P = 0.016) genes with decreased RA risk were discovered in this meta-analysis. Although, other KIR receptors including 2DL1, 2DL2, 2DL4, 3DL1, 3DL2, 3DS1, 2DS1-2DS4, and two pseudo gens of 2DP1 and 3DP1 displayed no significant association with predisposition to RA.
These findings provide reliable evidence that 2DL3, 2DL5, 3DL3, and 2DS5 might have a potential protective role for RA.
Rheumatoid arthritis (RA) is a common chronic autoinflammatory disorder, which is characterized by erosive articular lesions and progressive joint destruction, resulting in irreversible bone deformity and movement disabilities [1, 2]. The overall prevalence of RA ranges from 0.5 to 1% in global populations . It predominantly affects women with a female to male ratio of 2:1. It is believed that the combination of genetic and environmental factors as well as epigenomic contributing factors can break the self-tolerance and provoke immune system against itself during autoimmunity [4,5,6,7,8]. Current studies have reported approximately 60% genetic contribution in susceptibility to RA in various ethnic populations [9,10,11]. Human leukocyte antigen (HLA)-DR, one of the histocompatibility genes in HLA class II region, is responsible for about 30% of genetic susceptibility to RA . Thus, several genes outside the HLA region are also pivotal for initiation and development of RA. As shown in previous studies, a number of genes, including PTPN22, PADI4, STAT4, TRAF1, TNFAIP3, CD40, and KIR have been associated with RA and drawn much attention to investigate more genetic elements contributing to RA risk [9, 13,14,15,16,17].
The essential role of natural killer (NK) cells in RA pathogenesis have been widely investigated in many studies. These cells are present in synovial fluid as well as the blood of RA patients and participate in inflammation and, therefore, damage of joints [18, 19]. In normal physiological conditions, NK cells are able to eliminate target cells with little or no expression of HLA class I molecules on them. Beside cytotoxic function, their ability to secrete cytokines underlines the role of NK cells in development of autoimmune settings . The efficient functional activity of NK cells is governed by three distinct receptor families including, C-type lectin-like group, immunoglobulin (Ig)-like superfamily, and natural cytotoxicity receptors (NCR) . The KIRs (CD 158 family) are expressed mainly on NK cells and a few proportions of T lymphocytes. HLA class I molecules are considered to be the major ligands for KIRs . The cluster of KIR genes is mapped on the leukocyte receptor complex (LRC) locus on chromosome 19q13.4 which is extremely polymorphic among individuals and codes for glycoprotein receptors acting in inhibitory and activating manner . In humans, 17 highly homologous genes have been recognized for KIR genes encoding receptors containing 9 inhibitory, 6 activating, and 2 pseudogenes [24, 25]. KIR receptors are categorized based on the number of their extracellular Ig-like domains (2D or 3D) and whether they have short (S) or long (L) cytoplasmic tail. Two pseudogenes (2DP1 and 3DP1) do not encode functional receptors . Inhibitory KIRs contain long intracellular tails and transduce signals through immunoreceptor tyrosine-based inhibitory motifs (ITIM). On the other hand, activating signals are generated through stimulation of a transmembrane adaptor molecule, namely DNAX activation protein of 12 kDa (DAP12), which possess domains for immunoreceptor tyrosine-based activating motifs (ITAM) [27, 28].
The association of KIR gene polymorphisms has already been demonstrated in several autoimmune disorders [29, 30]. Previously, 11 genetic studies have tried to clarify the role of KIR polymorphisms in susceptibility to RA. Nonetheless, findings from these studies seem to be contradictory in some cases. This discrepancy may be due to small sample sizes, insufficient statistical power and ethnic variations among different populations. Meta-analysis is a statistical method that combines data from multiple researches to establish a reliable estimation. Therefore, in this meta-analysis, we intend to provide a comprehensive evaluation of previously published studies about KIR association with RA susceptibility to have a conclusive and consistent understanding of the role of these molecules in RA pathogenesis.
Literature search strategy
In the current meta-analysis, a systematic literature search was performed to identify relevant published articles in which associations between KIR variations and RA susceptibility were examined (up to May 2018). The electronic databases including Medline/PubMed and Scopus were searched using the following search keywords and subject terms: “The killer cell immunoglobulin-like receptors” or “KIR” or “KIRs”, “rheumatoid arthritis,” or “RA” and “polymorphism” or “variation” or “single nucleotide polymorphism” or “SNP”. All searched studies were retrieved, and their references were also reviewed to discover additional pertinent studies. Original data collected in English language and human population studies.
Inclusion and exclusion criteria
A study would be included if met following statements: a) independent case-control studies which assess the association of the KIR polymorphisms and the risk of RA; b) sufficient data for calculating the odds ratios (ORs) and 95% confidence intervals (CIs); c) studies with adequate data representing the genotype and allele frequency of patient and healthy subjects. Accordingly, the exclusion criteria were the following: a) reported duplicate publications and overlapping data b) studies without any information of detailed polymorphism frequencies c) Case reports, letter, review, comment, or abstract that published in journals (Table 1, Fig. 1).
Data extraction and quality assessment
The following items were recorded from all 11 eligible studies: first author’s name, publication year, country of origin, ethnicity of the study population, genotyping method, total number of cases and controls, and frequencies of KIR gene polymorphisms. The results were compared and discrepancies were resolved by consensus. For evaluating the quality of the involved studies the Newcastle-Ottawa Scale (NOS) was used . Total NOS scores ranged from 0 to 9. We regarded a study scored 0–3, 4–6, or 7–9 as a low, moderate, or high-quality study, respectively.
The effect size of association between KIR gene polymorphisms and RA susceptibility was determined by odds ratios (ORs) and their 95% confidence intervals (CI) comparing experimental cases to controls for each study. The phenotypic frequency (pf %) of each KIRs was assessed as the percentage of positive numbers between all samples. Genotypic frequency (gf %) among all samples was computed by means of the formula gf = 1– (1 – pf)1/2. In the current meta-analysis, heterogeneity among studies was measured by Cochran’s Q-statistic (P value< 0.10 was considered statistically significant heterogeneity) and quantified by I-squared (I2) tests . I2 values of 25, 50 and 75% were described as low, moderate, and high heterogeneity, respectively. A P value> 0.10 or I2 value< 50% demonstrate that the pooled OR was calculated using the fixed-effect model; otherwise, a random-effect model was used. Moreover, we implemented funnel plots to survey the existence of potential publication bias. Egger’s linear regression test and Begg’s test were employed to estimate the funnel plot asymmetry, and P < 0.05 was considered to be representative of statistically significant . If heterogeneity presented among studies, sensitivity analysis would be carried out. All statistical analyses were achieved by utilizing STATA (version 14.0; Stata Corporation, College Station, TX) and MedCalc version 13.
Characteristics of eligible studies
Regarding to aforementioned inclusion and exclusion criteria, a total of 11 case-control studies reporting the association between KIR gene and RA were available for this meta-analysis. (Shown in Fig. 1) The 11 published studies contained 1847 RA patients and 2409 healthy controls in various ethnic populations including four European, four Asian and three American countries [12, 33,34,35,36,37,38,39,40,41]. All the recruited papers were published from 2001 to 2016. According to the NOS criteria, all eligible studies had a total score ranging from 7 to 9. A detailed characteristic of the genotype and allele frequencies of various studies involved in this meta-analysis are provided in Table 2.
A summary of meta-analysis findings concerning associations between KIR variations and RA is represented in Table 2. Meta-analysis indicated four markedly significant negative associations of 2DL3 (OR = 0.591, 95% CI = 0.351–0.994; P = 0.047), 2DL5 (OR = 0.716, 95% CI = 0.601–0.853; P < 0.001), 2DS5 (OR = 0.623, 95% CI = 0.393–0.988; P = 0.045), and 3DL3 (OR = 0.324, 95% CI = 0.129–0.814; P = 0.016) with RA susceptibility (Fig. 2). Nevertheless, other KIR genes, including 2DL1, 2DL2, 2DL4, 3DL1, 3DL2, 3DS1, 2DS1-2DS4, and two pseudo gens of 2DP1 and 3DP1 showed no significant correlation with predisposition to RA.
A sensitivity analysis was performed by sequential omission of individual and groups of studies. The pooled ORs did not deviate with sequential omission of any participants or group of studies, indicating that our results were statistically robust (Fig. 3).
Heterogeneity and publication bias
In our meta-analysis, significant inter-study heterogeneity (I2 > 50%; P < 0.10) was found in 2DL1 (I2 = 73.8%; P < 0.001), 2DL2 (I2 = 80.2%; P < 0.001), 2DL3 (I2 = 71.1%; P < 0.001), 2DS1 (I2 = 50.7%; P = 0.047), 2DS2 (I2 = 59.9%; P = 0.01), 2DS5 (I2 = 86.7%; P < 0.001), 3DL1 (I2 = 68.1%; P < 0.001), 3DS1 (I2 = 68.4%; P = 0.002), and 2DP1 (I2 = 78.7%; P = 0.002). Therefore, the random-effects model was used to assess the relationship. On the contrary, other KIR polymorphisms revealed no significant heterogeneity and fixed-effects model were utilized to pool the result. Egger’s linear regression test and Begg’s test were applied to check publication bias. Funnel plots also did not show any apparent asymmetry (Fig. 4, Table 2). No publication bias was established for KIR genes in overall analysis (Egger’s regression test P > 0.1).
We intend to perform this meta-analysis to achieve a comprehensive analysis of the previously published studies about KIR association with RA susceptibility to figure out a conclusive and consistent understanding of the RA pathogenesis underlying KIR gene polymorphisms. The analysis revealed significant association of 2DL3, 2DL5, 2DS5, and 3DL3 genes, which decreased RA risk.
RA is a common autoimmune disorder, which mainly affects the joints and brings about bone erosions. One subject of debate at present is the role of genetic elements, which are involved in the initiation and development of RA. Integration analysis of genome wide association study (GWAS) successfully reported newly natural killer cell mediated cytotoxicity pathway as a contributing genes which are associated with RA .
Human NK cells are CD3−CD56+ lymphocytes of innate immune system and are subdivided into CD56bright and CD56dim phenotypes based on CD56 expression level. While CD56dim NK cells participate in cytotoxic responses, CD56bright subsets are highly involved in cytokine secretion especially interferon (IFN)-γ. According to previous researches, following diminished frequency of NK cells in peripheral blood, CD56bright subsets could migrate through blood vessels to inflamed joints and propagated in the synovial fluid of RA patients [43, 44]. Enhancement of NK cells in synovial fluids resulted in high amount of IFN-γ production which, in turn, contributed to perpetuation of inflammatory milieu in this autoimmune disorder . Of note, IFN-γ per se has ability to activate surrounding macrophages through cell-cell contact to release more inflammatory mediators, mostly tumor necrosis factor (TNF)-α which reinforce inflammatory environment in arthritic joints .
NK cells are also responsible for killing target cells that are not able to express self-MHC molecules . Their function is highly modulated by several activating and inhibitory membrane receptors, particularly KIRs . In 2001, the first case-control research regarding association of KIRs and RA undertaken by Yen et al. that showed the KIR2DS2 involvement in the development of rheumatoid vasculitis . From 2001 up to now, several studies were carried out to delineate the relationship of KIR genes repertoire with susceptibility to RA in diverse ethnic populations, though the findings appear inconsistent and conflicting. Thus, we decided to fulfill this comprehensive meta-analysis by incorporating data of 11 relevant studies in order to reevaluate the interest associations. The overall results divulged that KIR2DL3, KIR2DL5, KIR3DL3, and KIR2DS5 were negatively associated with susceptibility to RA, providing evidence that KIRs might potentially affect RA progression. On the other hand, other KIR genes, including 2DL1, 2DL2, 2DL4, 3DL1, 3DL2, 3DS1, 2DS1-2DS4, and two pseudo gens of 2DP1 and 3DP1 displayed no significant association with RA risk in our study.
KIR genes consist of two haplotypes; The “A” haplotype contains inhibitory genes and one activating gene (KIR2DS4) with protective effect against autoimmunities. “B” haplotype, on the other hand, exhibits more activating gens (KIR2DS1, -2DS2, -2DS3, -2DS4, -2DS5, and -3DS1) and inhibitory genes, KIR2DL2 and KIR2DL5 that increase the risk of autoimmune disorders [24, 26, 48]. It is widely accepted that the proper activation of NK cells depends on the equilibrium of activating and inhibitory KIR receptors. Considering their high degree of polymorphisms in two levels of gene content and allelic variations, they are mostly related to many autoimmune diseases, such as RA, ankylosing spondylitis, psoriasis vulgaris, scleroderma, SLE, diabetes mellitus, Crohn’s disease, Behçet’s disease, and etc. [21, 39, 49,50,51,52,53,54].
In this study, KIR2DL3 showed negative association with the RA risk and in agreement with the results of Parkash et al. and Ramirez-De los Santos et al. studies, it provided protection against RA [12, 38]. In this context, it has been clarified earlier that KIR2DL3 inhibits the production of cytokines like IFN-γ. Of note, NK cells-derived IFN-γ is a fundamental source for Th1 priming in RA . Therefore, lack of IFN-γ secretion by NK cells could reduce T cells autoreactivity .
In addition, our results also suggest that KIR2DL5 and KIR2DS5 were protective against RA. KIR2DL5 and KIR2DS5 both are present at haplotype B, but they function differently in the inhibitory and activating manner, respectively. Meanwhile, their respective ligands have not been discovered successfully .
Interestingly, KIR2DS5 display protective effect in some diseases, despite an activating function. A study by Nowak et al. showed that KIR2DS5 plays a protective role in ankylosing spondylitis, endometriosis, and acute kidney graft rejection with no significant effect in RA. According to their literature, this paradoxical protective effect is likely related to the pronounced role of KIR2DS5 in tolerance induction and, therefore, prevention of autoimmune conditions . It should be noted that this phenomenon might be due to linkage disequilibrium (LD) with some neighboring inhibitory loci in KIR haplotype. Thus, LD can contribute to modulate gene expression and give explanation as a putative cause of protective effect of KIR2DS5.
It has been reported previously that KIR2DL5 significantly had protective effect in SLE patients. Due to inhibitory signals from KIR2DL5, NK cells functions can be suppressed in the onset of SLE. Moreover, the less prevalence of KIR2DL5 associated with the increased predisposition to overall infections in SLE patients . Although their findings could not support the relation of KIR2DL5 in RA individuals, we observed a highly significant association of KIR2DL5 with risk of RA in our population study, which may impress the disease susceptibility in the same manner as does in SLE.
KIR3DL3 is one of the so-called “framework” genes, which is present in all KIR haplotypes. It is located in the centromeric half of haplotype at 5’-end and its ligand is uncertain yet . This inhibitory receptor can prohibit NK cell cytotoxicity via single ITIM motif in its cytoplasmic tail. Unlike previous reports, we revealed that a negative association exist between KIR3DL3 and RA risk.
A meta-analysis undertaken by Li et al. in 2015 to evaluate association between KIR variations and RA risk in East Asians and Caucasians subpopulations. The findings of their survey represented two positive associations, including 2DL1, 2DS1, and 2DL3 as a protective factor for RA risk only in East Asians, but not in Caucasians . In the current study, we incorporated more studies with further ethnic populations to achieve better results. The only consistence results obtained from current meta-analysis and that of Li et al. study, was 2DL3 as protective factor in RA patients . While the previous meta-analysis found two positive associations of KIR genes with RA susceptibility, we found that all significantly associated KIR genes had negative association with RA risk. KIR2DL5, KIR2DS5, and KIR3DL3 conferred a protective role in RA development .
This meta-analysis provides comprehensive information from 11 case-control studies to better understanding of KIR and RA association, but some shortcomings should be considered. First, between-study heterogeneity was observed which might distort analysis. Second, the number of published studies entered in this research was relatively small which restrict the estimation power of the present study. However, considering the increased number of subjects included in this meta-analysis, the obtained results are more reliable than the previous analysis . Additionally, various populations of Europe, Asia, and American nations were included in the current meta-analysis, increasing the general conclusion of the results to populations worldwide.
In consideration of all, present meta-analysis identified that KIR2DL3, KIR2DL5, KIR2DS5, and KIR3DL3 seem to have a potential protective role in RA risk. However, premised on forenamed limitations, further research studies are still required to explore the association between the KIR polymorphisms and RA risk with a larger sample sizes and controls in order to validate these findings.
Genome wide association study
Human leukocyte antigen
Immunoreceptor tyrosine-based activating motifs
Immunoreceptor tyrosine-based inhibitory motifs
Killer cell immunoglobulin-like receptor
Leukocyte receptor complex
Natural cytotoxicity receptors
- NK cells:
Natural killer cells
Tumor necrosis factor
McInnes IB, Schett G. Cytokines in the pathogenesis of rheumatoid arthritis. Nat Rev Immunol. 2007;7(6):429–42.
Mahmoudi M, Aslani S, Fadaei R, ARJIjord J. New insights to the mechanisms underlying atherosclerosis in rheumatoid arthritis. Int J Rheum Dis. 2017;20(3):287–97.
Wells AF, Curtis JR, Betts KA, Douglas K, Du EX, Ganguli A. Systematic literature review and meta-analysis of tumor necrosis factor–alpha experienced rheumatoid arthritis. Clin Ther. 2017;39(8):1680–1694.e1682.
Aslani S, Jafari N, Javan MR, Karami J, Ahmadi M, Jafarnejad M. Epigenetic modifications and therapy in multiple sclerosis. Neuromolecular Med. 2017;19(1):11–23.
Karami J, Mahmoudi M, Amirzargar A, Gharshasbi M, Jamshidi A, Aslani S, Nicknam MJG. Immunity: Promoter hypermethylation of BCL11B gene correlates with downregulation of gene transcription in ankylosing spondylitis patients. Genes Immun. 2017;18(3):170.
Ahmadi M, Gharibi T, Dolati S, Rostamzadeh D, Aslani S, Baradaran B, Younesi V, Yousefi MJB. Pharmacotherapy: Epigenetic modifications and epigenetic based medication implementations of autoimmune diseases. Biomed Pharmacother. 2017;87:596–608.
Foma AM, Aslani S, Karami J, Jamshidi A, Mahmoudi MJIR. Epigenetic involvement in etiopathogenesis and implications in treatment of systemic lupus erythematous. Inflamm Res. 2017;66(12):1057–73.
Rezaei R, Mahmoudi M, Gharibdoost F, Kavosi H, Dashti N, Imeni V, Jamshidi A, Aslani S, Mostafaei S, Vodjgani MJI. IRF 7 gene expression profile and methylation of its promoter region in patients with systemic sclerosis. Int J Rheum Dis. 2017;20(10):1551–61.
Kim K, Bang SY, Lee HS, Bae SC. Update on the genetic architecture of rheumatoid arthritis. Nat Rev Rheumatol. 2017;13(1):13–24.
Mousavi MJ, Jamshidi A, Chopra A, Aslani S, Akhlaghi M, Mahmoudi MJJ. Implications of the noncoding RNAs in rheumatoid arthritis pathogenesis. J Cell Physiol. 2018;234(1):335–47.
Aslani S, Mahmoudi M, Karami J, Jamshidi AR, Malekshahi Z, Nicknam MHJA. Epigenetic alterations underlying autoimmune diseases. Autoimmunity. 2016;49(2):69–83.
Ramirez-De los Santos S, Sanchez-Hernandez PE, Munoz-Valle JF, Palafox-Sanchez CA, Rosales-Rivera LY, Garcia-Iglesias T, Daneri-Navarro A, Ramirez-Duenas MG. Associations of killer cell immunoglobulin- like receptor genes with rheumatoid arthritis. Dis Markers. 2012;33(4):201–6.
Ali AM, Vino S. Genetic markers as therapeutic target in rheumatoid arthritis: a game changer in clinical therapy? Rheumatol Int. 2016;36(11):1601–7.
Messemaker TC, Huizinga TW, Kurreeman F. Immunogenetics of rheumatoid arthritis: understanding functional implications. J Autoimmun. 2015;64(Supplement C):74–81.
Abbasi M, Mousavi MJ, Jamalzehi S, Alimohammadi R, Bezvan MH, Mohammadi H, Aslani SJJocp: strategies toward rheumatoid arthritis therapy; the old and the new. 2018.
Aslani S, Mahmoudi M, Salmaninejad A, Poursani S, Ziaee V, Rezaei NJF. Pathology p: Lack of association between STAT4 single nucleotide polymorphisms and Iranian juvenile rheumatoid arthritis patients. Fetal Pediatr Pathol. 2017;36(3):177–83.
Asadipour M, Hassan-Zadeh V, Aryaeian N, Shahram F, Mahmoudi MJIjord: Histone variants expression in peripheral blood mononuclear cells of patients with rheumatoid arthritis 2017.
Fogel LA, Yokoyama WM, French AR. Natural killer cells in human autoimmune disorders. Arthritis Res Ther. 2013;15(4):216.
Shegarfi H, Naddafi F, Mirshafiey A. Natural killer cells and their role in rheumatoid arthritis: friend or foe? Sci World J. 2012;2012.
Rajalingam R. Human diversity of killer cell immunoglobulin-like receptors and disease. Korean J Hematology. 2011;46(4):216–28.
Kuijpers TW, Vendelbosch S, Berg M, Baeten DL. Killer immunoglobulin receptor genes in spondyloarthritis. Curr Opin Rheumatol. 2016;28(4):368–75.
Salim PH, Jobim M, Jobim LF, Xavier RM. Autoimmune rheumatic diseases and their association with killer immunoglobulin-like receptor genes. Rev Bras Reumatol. 2011;51(4):351–6 362-354.
Hsu KC, Chida S, Geraghty DE, Dupont B. The killer cell immunoglobulin-like receptor (KIR) genomic region: gene-order, haplotypes and allelic polymorphism. Immunol Rev. 2002;190:40–52.
Marsh SG, Parham P, Dupont B, Geraghty DE, Trowsdale J, Middleton D, Vilches C, Carrington M, Witt C, Guethlein LA, et al. Killer-cell immunoglobulin-like receptor (KIR) nomenclature report, 2002. Tissue Antigens. 2003;62(1):79–86.
Tajik N, Shahsavar F, Mousavi T, Radjabzadeh MF. Distribution of KIR genes in the Iranian population. Tissue Antigens. 2009;74(1):22–31.
Middleton D, Gonzelez F. The extensive polymorphism of KIR genes. Immunology. 2010;129(1):8–19.
Williams AP, Bateman AR, Khakoo SI. Hanging in the balance. KIR and their role in disease. Mol Interv. 2005;5(4):226–40.
McVicar DW, Taylor LS, Gosselin P, Willette-Brown J, Mikhael AI, Geahlen RL, Nakamura MC, Linnemeyer P, Seaman WE, Anderson SK, et al. DAP12-mediated signal transduction in natural killer cells. A dominant role for the Syk protein-tyrosine kinase. J Biol Chem. 1998;273(49):32934–42.
Rezaei R, Mostafaei S, Aslani S, Jamshidi A, Mahmoudi MJIjord: Association study between killer immunoglobulin-like receptor polymorphisms and ankylosing spondylitis disease: an updated meta-analysis. 2018.
Fathollahi A, Aslani S, Mostafaei S, Rezaei N, Mahmoudi MJIR: The role of killer-cell immunoglobulin-like receptor (KIR) genes in susceptibility to inflammatory bowel disease: systematic review and meta-analysis 2018:1–10.
Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol. 2010;25(9):603–5.
Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ (Clinical research ed). 1997;315(7109):629–34.
Yen JH, Moore BE, Nakajima T, Scholl D, Schaid DJ, Weyand CM, Goronzy JJ. Major histocompatibility complex class I-recognizing receptors are disease risk genes in rheumatoid arthritis. J Exp Med. 2001;193(10):1159–67.
Yen JH, Lin CH, Tsai WC, Wu CC, Ou TT, Hu CJ, Liu HW. Killer cell immunoglobulin-like receptor gene's repertoire in rheumatoid arthritis. Scand J Rheumatol. 2006;35(2):124–7.
Majorczyk E, Pawlik A, Luszczek W, Nowak I, Wisniewski A, Jasek M, Kusnierczyk P. Associations of killer cell immunoglobulin-like receptor genes with complications of rheumatoid arthritis. Genes Immun. 2007;8(8):678–83.
Middleton D, Meenagh A, Wright GD. No association in frequency of KIR receptors in patients with rheumatoid arthritis from Northern Ireland. Tissue Antigens. 2007;69(6):577–82.
Nowak I, Majorczyk E, Wisniewski A, Pawlik A, Magott-Procelewska M, Passowicz-Muszynska E, Malejczyk J, Ploski R, Giebel S, Barcz E, et al. Does the KIR2DS5 gene protect from some human diseases? PLoS One. 2010;5(8):e12381.
Prakash S, Alam S, Bharadwaj U, Aggarwal A, Mishra RN, Agrawal S. Associations of killer cell immunoglobulin like receptors with rheumatoid arthritis among north Indian population. Hum Immunol. 2014;75(8):802–7.
Nazari M, Mahmoudi M, Rahmani F, Akhlaghi M, Beigy M, Azarian M, Shamsian E, Akhtari M, Mansouri R. Association of Killer Cell Immunoglobulin- like Receptor Genes in Iranian patients with rheumatoid arthritis. PLoS One. 2015;10(12):e0143757.
Nishimura WE, Sachetto Z, Costallat LT, Yazbek MA, Londe AC, Guariento EG, Marques SB, Bertolo MB. The role of KIR2DL3/HLA-C*0802 in Brazilian patients with rheumatoid vasculitis. Clinics (Sao Paulo, Brazil). 2015;70(6):408–12.
Velarde-de la Cruz EE, Sanchez-Hernandez PE, Munoz-Valle JF, Palafox-Sanchez CA, Ramirez-de Los Santos S, Graciano-Machuca O, Garcia-Iglesias T, Montoya-Buelna M, Ramirez-Duenas MG. KIR2DL2 and KIR2DS2 as genetic markers to the methotrexate response in rheumatoid arthritis patients. Immunopharmacol Immunotoxicol. 2016;38(4):303–9.
Han J, Geng Y, Deng X, Zhang Z. Subclinical synovitis assessed by ultrasound predicts flare and progressive bone erosion in rheumatoid arthritis patients with clinical remission: a systematic review and metaanalysis. J Rheumatol. 2016;43(11):2010–8.
Dalbeth N, Callan MFC. A subset of natural killer cells is greatly expanded within inflamed joints. Arthritis Rheumatism. 2002;46(7):1763–72.
Pridgeon C, Lennon GP, Pazmany L, Thompson RN, Christmas SE, Moots RJ. Natural killer cells in the synovial fluid of rheumatoid arthritis patients exhibit a CD56bright,CD94bright,CD158negative phenotype. Rheumatology. 2003;42(7):870–8.
Gianchecchi E, Delfino DV, Fierabracci A. NK cells in autoimmune diseases: linking innate and adaptive immune responses. Autoimmun Rev. 2018;17(2):142–54.
Ahern DJ, Brennan FM. The role of natural killer cells in the pathogenesis of rheumatoid arthritis: major contributors or essential homeostatic modulators? Immunol Lett. 2011;136(2):115–21.
Vivier E, Raulet DH, Moretta A, Caligiuri MA, Zitvogel L, Lanier LL, Yokoyama WM, Ugolini SJs. Innate or adaptive immunity? The example of natural killer cells. 2011;331(6013):44–49.
Kulkarni S, Martin MP, Carrington M. The yin and Yang of HLA and KIR in human disease. Semin Immunol. 2008;20(6):343–52.
Mahmoudi M, Jamshidi AR, Karami J, Mohseni A, Amirzargar AA, Farhadi E, Ahmadzadeh N, Nicknam MH. Analysis of killer cell immunoglobulin-like receptor genes and their HLA ligands in Iranian patients with ankylosing spondylitis. Iran J Allergy Asthma Immunol. 2016;15(1):27–38.
Łuszczek W, Mańczak M, Cisło M, Nockowski P, Wiśniewski A, Jasek M, Kuśnierczyk P. Gene for the activating natural killer cell receptor, KIR2DS1, is associated with susceptibility to psoriasis vulgaris. Hum Immunol. 2004;65(7):758–66.
Mahmoudi M, Fallahian F, Sobhani S, Ghoroghi S, Jamshidi A, Poursani S, Dolati M, Hosseinpour Z, Gharibdoost F. Analysis of killer cell immunoglobulin-like receptors (KIRs) and their HLA ligand genes polymorphisms in Iranian patients with systemic sclerosis. Clin Rheumatol. 2017;36(4):853–62.
Akhtari M, Farazmand A, Mahmoudi M, Akbarian M, Ahmadzadeh N, Mirkazemi Z, Mostafaei S, Jamshidi AR. Analysis of killer cell immunoglobulin-like receptors and their human leukocyte antigen-ligands gene polymorphisms in Iranian patients with systemic lupus erythematosus. Lupus. 2016;25(11):1244–53.
Liang HL, Ma SJ, Tan HZ. Association between killer cell immunoglobulin-like receptor (KIR) polymorphisms and systemic lupus erythematosus (SLE) in populations: a PRISMA-compliant meta-analysis. Medicine (Baltimore). 2017;96(10):e6166.
Diaz-Pena R, Vidal-Castineira JR, Moro-Garcia MA, Alonso-Arias R, Castro-Santos P. Significant association of the KIR2DL3/HLA-C1 genotype with susceptibility to Crohn's disease. Hum Immunol. 2016;77(1):104–9.
Flodstrom-Tullberg M, Bryceson YT, Shi FD, Hoglund P, Ljunggren HG. Natural killer cells in human autoimmunity. Curr Opin Immunol. 2009;21(6):634–40.
Morandi B, Bougras G, Muller WA, Ferlazzo G, Munz C. NK cells of human secondary lymphoid tissues enhance T cell polarization via IFN-gamma secretion. Eur J Immunol. 2006;36(9):2394–400.
Campbell KS, Purdy AK. Structure/function of human killer cell immunoglobulin-like receptors: lessons from polymorphisms, evolution, crystal structures and mutations. Immunology. 2011;132(3):315–25.
Kimoto Y, Horiuchi T, Tsukamoto H, Kiyohara C, Mitoma H, Uchino A, Furugo I, Yoshizawa S, Ueda A, Harashima S, et al. Association of killer cell immunoglobulin-like receptor 2DL5 with systemic lupus erythematosus and accompanying infections. Rheumatology (Oxford, England). 2010;49(7):1346–53.
Martin MP, Single RM, Wilson MJ, Trowsdale J, Carrington M. KIR haplotypes defined by segregation analysis in 59 Centre d’Etude Polymorphisme Humain (CEPH) families. Immunogenetics. 2008;60(12):767–74.
Li X, Xia Q, Fan D, Cai G, Yang X, Wang L, Xin L, Ding N, Hu Y, Liu L, et al. Association between KIR gene polymorphisms and rheumatoid arthritis susceptibility: a meta-analysis. Hum Immunol. 2015;76(8):565–70.
This work was supported by Deputy of Research, Tehran University of Medical Sciences (Grant No. 97-03-41-40378). The funders were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.