This article has Open Peer Review reports available.
Influence of KIR genes and their HLA ligands in the pathogenesis of leprosy in a hyperendemic population of Rondonópolis, Southern Brazil
- Luciana Ribeiro Jarduli1,
- Hugo Vicentin Alves1,
- Fabiana Covolo de Souza-Santana2,
- Elaine Valim Camarinha Marcos2,
- Ana Carla Pereira2,
- Ida Maria Foschiani Dias-Baptista2,
- Vinícius Medeiros Fava3,
- Marcelo Távora Mira3,
- Milton Ozório Moraes4,
- Marcos da Cunha Lopes Virmond2 and
- Jeane Eliete Laguila Visentainer1Email author
© Jarduli et al.; licensee BioMed Central Ltd. 2014
Received: 3 July 2013
Accepted: 4 August 2014
Published: 12 August 2014
The objective of this study was to investigate the association between KIR genes and the immunopathogenesis of leprosy.
The types of KIR and HLA genes were evaluated by PCR-SSOP-Luminex in 408 patients with leprosy and 413 healthy individuals. Statistical analysis was performed using the Chi-square or Fisher’s exact test and stepwise multivariate analysis.
There was a higher frequency of activating KIR genes (KIR2DS1, 2DS2 and 3DS1) together with their HLA ligands in the tuberculoid (TT) group as compared to the lepromatous leprosy (LL) group. KIR2DL2/2DL2-C1 was more frequent in the patient, TT and LL groups than in the control group. Borderline patients presented a higher frequency of inhibitory pairs when compared to the control group, and a higher frequency of activating pairs as compared to the LL group. Multivariate analysis confirmed the associations and demonstrated that being a female is a protective factor against the development of the disease per se and the more severe clinical form.
This study showed that activating and inhibitory KIR genes may influence the development of leprosy – in particular, activating genes may protect against the more aggressive form of the disease – thereby demonstrating the role of NK cells in the immunopathology of the disease.
Leprosy or Hansen’s disease is a chronic, slowly developing infectious disease caused by Mycobacterium leprae, which global prevalence has been decreasing over recent years due to effective multiple drug therapy. However, the detection of new cases remains constant at approximately 250,000 per year, with pockets of high endemicity existing in countries such as Brazil, Angola, Africa, India, Nepal, and the United Republic of Tanzania . In Brazil, leprosy has been diagnosed in all regions of the country, with the Northern and Mid-Western regions being hyperendemic and the Northeastern region having high endemicity .
The most widely used classification model for the disease is that of Ridley-Jopling, which divides the disease into five polar and interpolar forms: Tuberculoid (TT), Borderline-Tuberculoid (BT), Borderline-Borderline (BB), Borderline-lepromatous leprosy (BL) and Lepromatous leprosy (LL) . Initially, the disease is of indeterminate form, which may regress spontaneously or evolve into one of these stages. The clinical manifestations of the disease depend on factors such as the ability of the bacilli to proliferate and the immune response presented by the host .
Over the years, research has been aimed at better understanding the genetics of susceptibility to the different clinical forms of leprosy. The clinical and pathological spectrum of leprosy and the epidemiological, geographical, and ethnic heterogeneity might be explained by genetic differences in the host’s resistance. While some loci affect the intrinsic susceptibility to leprosy itself, others modify the clinical form of the disease . Among the candidate genes for an individual’s susceptibility or resistance to leprosy are the KIR (Killer Immunoglobulin-Like Receptor) genes. These genes span over approximately 150 Kb of the LCR (Leukocyte Receptor Complex) region of chromosome 19q13.4, which encodes KIRs, members of a group of regulatory molecules present on the surface of NK (Natural Killer) cells .
The KIRs are the main functional regulators of NK cells; the activity of these effector cells is determined by the balance between their activation and inhibition that occurs as a result of the binding of KIRs with HLA class I molecules present in all nucleated cells . Most KIRs bind to HLA-C molecules; the KIR2DL2 and KIR2DL3 inhibitory receptors interact with HLA-C of Group 1, while KIR2DL1 recognizes HLA-C of Group 2. In the absence of these inhibitory receptors, NK cells can be activated by activating receptors, such as KIR2DS1 and KIR2DS2 . The genotype may predispose a subject to certain interactions that vary from activation to inhibition depending on the individual’s KIR/HLA profile .
To date, only one study has reported an association between KIR genes and their HLA ligands with leprosy phenotypes, using a population sample of 165 patients from south Brazil . Thus, the aim of this study was to investigate the association between KIR genes and their HLA ligands in the development of leprosy and its clinical forms in patients from the hyperendemic region of Rondonópolis, Mato Grosso, Mid-Western Brazil, ultimately aiming to better understand the immunopathogenic mechanisms of M. leprae.
Population and epidemiological data collection
This study was approved by the Research Ethics Committee of the Instituto Lauro de Souza Lima Bauru, São Paulo, Brazil, and funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). All participants signed an informed consent form agreeing to participate in the study. Data for epidemiological characterization were collected in government healthcare clinics by applying a questionnaire and by revisiting medical records and notification registration forms of patients treated between 2006 and 2009.
This case–control study was conducted using a population sample from Rondonópolis (Mid-Western State of Mato Grosso, Brazil: latitude 16°28′15″ South and longitude 54°38′08″ West), consisting of 408 patients with leprosy (250 men and 158 women) with a median age of 41 years treated in government healthcare clinics of the municipality of Rondonópolis. The patients were classified into distinct groups according to the Ridley-Jopling criteria : LL (5.14%), borderline (BL 15.4%; BB 25.8%; BT 39.4%) and TT (14.7%), with 25 patients (6.12%) presenting the indeterminate form of the disease. The control group consisted of healthy individuals 413 (249 men and 164 women) with a median age of 42 years. The following criteria were used while forming the control group: no family and/or personal history with leprosy and no history of chronic infections, inflammation, or autoimmune diseases. The control group was matched to the case group according to epidemiological characteristics such as ethnicity, gender, age and geographical region.
After collecting 5 mL of peripheral venous blood in EDTA anticoagulant, the extraction and purification of DNA were performed by selective precipitation using the technique described by John et al. .
Typing of the KIR genes and HLA class I alleles
Typing of the KIR genes and HLA class I alleles (HLA-A, -B, -C) was performed by PCR-SSOP (Polymerase Chain Reaction-Sequence Specific Oligonucleotide Probes) using the SSO LabType® kit (One Lambda, San Diego, CA, USA) after evaluating the purity and adjusting the concentration of the DNA (20 ng). The amplified product was hybridized with microspheres bound to probes specific for the KIR gene and HLA class I alleles. The resulting products were analysed using a LUMINEX® flow cytometer and the results were compared by HLA Fusion™ (One Lambda, San Diego, CA, USA).
The frequencies of KIR genes and HLA class I alleles were obtained by direct counting and by comparing these frequencies between patients and controls using the Chi-square test with Yates correction or Fisher’s exact test using a 2 × 2 contingency table with a 95% confidence interval (95% CI). The risk of developing leprosy and its clinical forms was calculated by determining the OR (odds ratio) for P-values < 0.05 according to Wolf , using the statistics program OpenEpi v. 2.3.1 (http://www.openepi.com/OE2.3/Menu/OpenEpiMenu.htm). P values were corrected by the Bonferroni inequality method, by multiplying them by the number of KIR genes analysed.
Subsequently, multivariate logistic regression analysis was performed using stepwise method of data analysis, using age, gender and ethnicity as covariates, employing IBM® SPSS® 20.0 Statistics software (SPSS, IBM Inc., Chicago, IL, USA) with the level of significance set at 5% (P-value ≤ 0.05).
Distribution of the population of patients and healthy individuals according to clinical forms of leprosy, age, gender and ethnic group
N = 408
N = 21
N = 302
N = 60
N = 413
(Mean ± SD)
42.0 ± 16.4
45.7 ± 15.2
42.8 ± 15.9
37.4 ± 16.8
42.0 ± 14.5
11 ( 3.64)
3 ( 5.0)
17 ( 4.11)
There was a predominance of men in the total patients (61.2%), LL (80.9%), Borderline (64.2%) and Control groups (60.2%), while the TT group was composed predominantly of females (58.3%). There was a predominance of Afro-Brazilians followed by Caucasians and Blacks in all groups, with similar proportions of ethnicity in each group.
Distribution of the frequencies of KIR genes in patients with leprosy and their clinical forms and healthy individuals in the population of Rondonopolis, MT
N = 408
N = 21
N = 302
N = 60
N = 413
Frequencies of KIR genes in the presence of their ligands HLA Class I in patients with leprosy and their clinical forms and healthy individuals in the population of Rondonópolis, MT
KIR – HLA Class I
N = 408
N = 21
N = 302
N = 60
N = 413
KIR2DL1 – C2
KIR2DL1 – C2/C2
KIR2DL2 – C1
KIR2DL2 – C1/C1
KIR2DL3 – C1
KIR2DL3 – C1/C1
KIR3DL1 – Bw4
KIR3DL1 – Bw4/Bw4
KIR3DL2 – A*03/*11
KIR2DS1 – C2
KIR2DS1 – C2/C2
KIR2DS2 – C1
KIR2DS2 – C1/C1
KIR3DS1 – Bw4
KIR3DS1 – Bw4/Bw4
KIR2DL2/2DL2 – C1
KIR2DL2/2DL2 – C1/C1
KIR2DL2/2DL3 – C1
KIR2DL2/2DL3 – C1/C1
KIR2DL3/2DL3 – C1
KIR2DL3/2DL3 – C1/C1
Stepwise logistic regression analysis of factors associated with the development of leprosy and its clinical forms
per se vs. controls
KIR2DL1 - C2/C2
KIR2DS2 - C1
LL vs. controls
KIR2DL2/2DL2 - C1
Borderline vs. controls
TT vs. controls
KIR2DL1 - C2/C2
KIR2DL2 - C1
KIR2DL2/2DL2 - C1
TT vs. Borderline
KIR2DL1 - C2/C2
KIR2DL2 - C1
TT vs. LL
Borderline vs. LL
KIR2DL2/2DL2 - C1
When only the HLA ligands of the KIR gene were evaluated, they were in equilibrium of Hardy-Weinberg, however no significant differences in their frequencies were found between groups.
Studying the impact of genes on infectious diseases such as leprosy in endemic regions is extremely important because it provides information about the contributions of host genetic factors and the environment in which the individual lives on the development of the disease. Historically in Brazil, the Mid-Western region has always had high rates of new cases of leprosy and the least favourable evolution of endemics in Brazil . A recent study conducted by Magalhães et al.  assessed the pattern of leprosy in the State of Mato Grosso. The data revealed that the evolution of the disease was related to the process by which the region was settled, with successive migrations causing changes in the epidemiological structure due to the poor living conditions of the migrant population associated with other factors of infection in the disease process.
The current study shows that men are more affected by leprosy than women; the literature presents conflicting results regarding the prevalence of leprosy with respect to gender [16, 17]. Some authors claim that the greater social contact between men and their frequent exposure to risk environments can contribute to the increased number of cases . Others claim that the prevalence of the disease in men is not universal and that, when related to leprosy reactions, women who are pregnancy or breastfeeding are at greater risk due to immunological and genetic characteristics . In the evaluation of ethnic background, the highest frequency was found in Afro-Brazilians followed by Caucasians, which reflects the greater proportion of half-castes in the population of Mato Grosso.
In the analysis of activating KIR genes, the KIR2DS2-C1 pair was significantly more common in the Total Patient and TT Groups than the Control Group. It is important to note the presence of a strong binding imbalance between the KIR2DL2 and KIR2DS2 genes , as was confirmed in this study. The KIR2DL2/2DS2 haplotype frequencies were higher than expected in the Total Patient and Control Groups. Considering the presence of linkage disequilibrium between these two genes, the presence of the association may be attributable to the high frequency of the KIR2DL2 gene for both Groups.
The function of the KIR genes in the immune response is highly dependent on the HLA molecules expressed on the surface of target cells; the KIR-HLA interaction can result in either activation or inhibition of NK cells. The recognition of specific HLA molecules by inhibiting KIRs is well established , with a hierarchy of inhibition of the KIR2DL molecules: 2DL1-C2/C2 has the greatest inhibitory potential, followed by 2DL2-C1 and 2DL3-C1 . In this study, was found a tendency of association for the gene KIR2DL1-C2C2 between the BL and BT forms, which are two opposite clinical forms, being BL nearer to the lepromatous pole. The results suggest that greater inhibition of NK cells may occur in these patients, leading to a more severe manifestation of the disease.
The distribution of the KIR2DL3 and KIR2DL2 inhibitory genes, and whether they were homozygous or heterozygous, was observed in this study, supporting the hypothesis that these genes segregate as alleles from a single locus [21, 22]. Interesting results were found for KIR2DL2/2DL2-C1; its frequency was higher in all patient groups compared to the Control Group, although the difference was only significant for the Total Patient, TT and LL Groups, suggesting that the inhibitory effect of KIR2DL2/2DL2-C1 may contribute to the development of leprosy, but mainly contributes to a worse prognosis in M. leprae infections.
It is important to highlight the results observed in the analysis of the number of pairs of KIR-HLA ligands. The Total Patient Group had a more inhibitory profile due to the higher frequency of only one inhibitory pair than the Control Group, which may favour the development of the disease. A similar result was found for the Borderline Group of patients compared to the Control Group. Interestingly, in the analysis of the activating pairs of ligands, a higher frequency of two pairs was observed in the borderline patients compared to the LL Group. The interpolar forms, termed borderline, comprised the majority of the cases. The analysis of the KIR-HLA ligand pairs in this group showed a more inhibitory profile relative to the Control Group, added to a more activating profile compared to the LL Group. These individuals have a characteristic immune instability against the mycobacteria, which means that there is great variation in their clinical manifestations, whether of the skin, nerves, or systemic involvement. The balance between activation and inhibition of NK cells may explain the ill-defined characteristics observed in these patients.
The results of the multivariate analysis, in addition to confirming the influence of the KIR genes and their HLA ligands on the immunopathology of leprosy and its different forms, demonstrated the influence of gender. Interestingly, the LL Group presented a negative association with respect to gender compared to the Control Group. Thus, there were fewer females in the LL Group compared to the Control Group. On the other hand, there was a positive association with respect to gender on comparing the TT Group: TT vs. Control (β = 0.7; OR = 2.1), TT vs. Borderline (β = 0.9; OR = 2.5) or TT vs. LL (β = 1.7; OR = 5.9). Among the groups, the TT Group was the only one that comprised more women than men. These results suggest that being a woman is a protective factor against the most severe form and spread of the disease, in addition to contributing to its mildest manifestation with susceptibility to the TT form. Specific immune responses to autoimmune and infectious diseases differ between men and women. The female sex hormone, oestrogen, can increase the production of INF-γ and IL-2, thereby influencing the immune response in women. Gonadotropin-releasing hormone (GnRH) is involved in the maturation of the thymus and exerts a potent stimulatory effect, leading to increases in IL-2 and its receptor (IL-2R), INF-γ and activation of CD4+ helper T cells . The male androgens, on the other hand, are related to the inhibition of T and B cell immune response and are found at low levels in men suffering from autoimmune diseases such as systemic lupus erythematosus .
With respect to ethnicity, the multivariate analysis demonstrated a positive association between TT and controls (β 0.3; OR 1.4) and between TT and Borderline (β 0.4; OR 1.5). For both groups, the TT Group had a higher frequency of Afro-Brazilians. It is possible that the migrations that occurred during settlement of the State of Mato Grosso influenced the ethnic composition of the population, whereby individuals classified as of mixed descent result from miscegenation between Afro-Brazilians and Indians.
The significant difference in the percentage of the KIR2DL1 gene in the Total Patient and Borderline Groups compared to the Control Group in the Chi-square analysis was also significant in some multivariate analysis comparisons; however, in both cases the results were not correlated with the development of the disease or its different forms. KIR2DL1-C2/C2, which was not significant in the Chi-square analysis, was significantly different in the comparison between the Total Patient and Control groups; however, this did not result in associations with the disease or its clinical forms. The KIR2DL2/2DL2-C1 inhibitory pair, which was significant in the Chi-square analysis, was also significant in the multivariate analysis (LL vs. Control; Borderline vs. LL), with a tendency being seen in the comparison between the TT vs. Control Groups. The significance of KIR2DL2-C1 in the multivariate analysis may be due to the higher frequency of this pair in the TT Group (45.0%), but without implications for the pathogenesis of the disease. An equivalent result was observed for KIR2DS2-C1, whose significance, found when comparing the Total Patient and Control Groups, is probably due to a linking imbalance with KIR2DL2 as reported for the two groups in this study.
The current study shows that activating and inhibitory KIR genes in the presence of their corresponding HLA ligands may have some effect on the development of leprosy and its clinical forms, seen particularly with the strong presence of activating genes in the TT Group. The balance between activating and inhibitory genes may interfere with the progression of the disease to milder or more aggressive forms, or even cause it to vary between the two poles with indeterminate forms, thus highlighting the role of NK cells in the immunopathology of the disease.
The highlight of this study was identifying the profiles of patients with leprosy from a hyperendemic area of the country and comparing them to a group of extremely well matched healthy individuals.
The authors would like to thank the staff of the Lauro de Souza Lima Institute in Bauru, the Laboratory of Immunogenetics of UEM and the government health clinics of Rondonópolis.
This study was funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Instituto Lauro Souza Lima, Bauru and the Immunogenetics laboratory of the State University of Maringá.
- World Health Organization: Global leprosy situation. Wkly Epidemiol Rec. 2012, 87: 317-328.Google Scholar
- Ministério da Saúde: Situação epidemiológica da hanseníase no Brasil, 2001–2006. http://portal.saude.gov.br/portal/arquivos/pdf/poster_hansen_brasil.pdf,
- Ridley DS, Jopling WH: Classification of leprosy according to immunity. A five-group system. Int J Lepro Other Mycobact Dis. 1966, 34: 255-273.Google Scholar
- Misch EA, Berrigton WR, Vary JC, Hawn TR: Leprosy and the human genome. Microbiol Mol Biol Rev. 2010, 74: 589-620. 10.1128/MMBR.00025-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Schurr E, Alcaïs A, de Léséleuc L, Abel L: Genetic predisposition to leprosy: a major gene reveals novel pathways of immunity to Mycobacterium leprae. Semin Immunol. 2006, 18: 404-410. 10.1016/j.smim.2006.07.005.View ArticlePubMedGoogle Scholar
- Rajagopalan S, Long EO: Understanding how combinations of HLA and KIR genes influence disease. J Exp Med. 2005, 201: 1025-1029. 10.1084/jem.20050499.View ArticlePubMedPubMed CentralGoogle Scholar
- Rajalingam R: Human diversity of killer cell immunoglobulin-like receptors and disease. Korean J Hematol. 2011, 46: 216-228. 10.5045/kjh.2011.46.4.216.View ArticlePubMedPubMed CentralGoogle Scholar
- Winter CC, Gumperz JE, Parham P, Long EO, Wagtmann N: Direct binding and functional transfer of NK cell inhibitory receptors reveal novel patterns of HLA-C allotype recognition. J Immunol. 1998, 161: 571-577.PubMedGoogle Scholar
- Carrington M, Martin MP: The impact of variation at the KIR gene cluster on human disease. Curr Top Microbiol Immunol. 2006, 298: 225-257.PubMedGoogle Scholar
- Franceschi DS, Mazini PS, Rudnick CCC, Sell AM, Tsuneto LT, de Melo FC, Braga MA, Peixoto PRF, Visentainer JEL: Association between killer-cell immunoglobulin-like receptor genotypes and leprosy in Brazil. Tissue Antigens. 2008, 72: 478-482. 10.1111/j.1399-0039.2008.01127.x.View ArticlePubMedGoogle Scholar
- John SW, Weitzner G, Rozen R, Scriver CR: A rapid procedure for extracting genomic DNA from leukocytes. Nucleic Acids Res. 1991, 19: 408-10.1093/nar/19.2.408.View ArticlePubMedPubMed CentralGoogle Scholar
- Woolf B: On estimating the relation between blood group and disease. Ann Hum Genet. 1995, 19: 251-253.View ArticleGoogle Scholar
- Hosmer DW, Lemeshow S: Applied logistic regression. John Wiley & Sons. 1989, 10: 1162-1163.Google Scholar
- Ministério da Saúde: Programa Nacional de Controle da Hanseníase: Informe Epidemiológico. http://www.portal.saude.gov.br/arquivos/pdf/boletim-novembro.pdf,
- Magalhães MCC, Santos ES, Queiroz ML, Lima ML, Borges RCM, Souza MS, Ramos AN: Migração e hanseníase em Mato Grosso. Rev bras epidemiol. 2011, 14: 386-397. 10.1590/S1415-790X2011000300004.View ArticleGoogle Scholar
- Miranzi SSC, Pereira LHM, Nunes AA: Perfil epidemiológico da hanseníase em um município brasileiro, no período de 2000 a 2006. Rev Soc Bras Med Trop. 2010, 43: 62-67. 10.1590/S0037-86822010000100014.View ArticleGoogle Scholar
- Lima HMN, Sauaia N, Costa VRL, da Costa L, Neto GTC, Figueiredo PMS: Perfil epidemiológico dos pacientes com hanseníase atendidos em Centro de Saúde em São Luís, MA. Rev Bras Clin Med. 2010, 8: 323-327.Google Scholar
- Silva AR, Matos WB, Silva CCB, Gonçalves EGR: Hanseníase no município de Buriticupu, estado do Maranhão: busca ativa de casos na população adulta. Rev Soc Bras Med Trop. 2010, 43: 691-694. 10.1590/S0037-86822010000600018.View ArticlePubMedGoogle Scholar
- Lodoen MB, Lanier LL: Natural killer cells as an initial defense against pathogens. Curr Opin Immunol. 2006, 18: 391-398. 10.1016/j.coi.2006.05.002.View ArticlePubMedGoogle Scholar
- Long EO, Rajagopalan S: HLA class I recognition by killer cell Ig-like receptors. Semin Immunol. 2000, 12: 101-108. 10.1006/smim.2000.0212.View ArticlePubMedGoogle Scholar
- Khakoo SI, Carrington M: KIR and disease: a model system or system of models?. Immunol Rev. 2006, 214: 186-201. 10.1111/j.1600-065X.2006.00459.x.View ArticlePubMedGoogle Scholar
- Parham P: MHC class I molecules and KIRs in human history, health and survival. Nat Rev Immunol. 2005, 5: 201-214.View ArticlePubMedGoogle Scholar
- Tanriverdi F, Silveira LF, MacColl GS, Bouloux PM: The hypothalamic-pituitary-gonadal axis: immune function and autoimmunity. J Endocrinol. 2003, 176: 293-304. 10.1677/joe.0.1760293.View ArticlePubMedGoogle Scholar
- Verthelyi D, Petri M, Ylamus M, Klinman DM: Disassociation of sex hormone levels and cytokine production in SLE patients. Lupus. 2001, 10: 352-358. 10.1191/096120301674365881.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/438/prepub
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 credited.