- Research article
- Open Access
Evaluation of host genetics on outcome of tuberculosis infection due to differences in killer immunoglobulin-like receptor gene frequencies and haplotypes
BMC Genetics volume 16, Article number: 63 (2015)
Outcome of Mycobacterium tuberculosis (Mtb) infection is affected by virulence of the infecting strain of Mtb, host environment, co-morbidities, and the genetic composition of the host, specifically the presence or absence of genes involved in immune responses/regulation. It is hypothesized that specific killer immunoglobulin-like receptor (KIR) genes may be associated with Mtb infection and clinical outcome. This cross-sectional study examined the KIR gene frequencies, profiles, and haplotypes of individuals with active tuberculosis, latent tuberculosis infection, compared to TB and HIV negative healthy controls.
Analysis of KIR gene frequencies revealed differences among disease status groups, suggesting that enrichment or depletion of specific KIR genes may direct the disease outcome. Mtb infected individuals were more likely to have a centromeric-AA haplotype compared to controls.
The differences in KIR gene frequencies and haplotypes may result in differential cytokine expression, contributing to different disease outcomes, and suggest a genetic influence on Mtb susceptibility and pathogenesis.
Tuberculosis (TB) incidence in the world and in Canada overall are declining toward goals set by the World Health Organization , however, certain populations and/or regions within Canada continue to have rates of tuberculosis exceeding the national average. In 2012, the global burden of TB was estimated at 8.6 million incident cases (122/100,000 population) . Canada reported 1686 new active cases of TB (4.8/100,000) in 2012 . The incidence of TB in Manitoba was more than twice the national rate at 138 cases (10.9/100,000).
The immune response to an intracellular pathogen such as Mycobacterium tuberculosis (Mtb) involves natural killer (NK) cells to bridge the innate and adaptive immune response to infection . NK cells are important in early Mtb infection, as they are capable of activating phagocytic cells at the site of infection , and are producers of INF-γ, which functions to activate macrophages . The activity of NK cells is controlled by a balance of inhibitory and stimulatory signals generated when human leukocyte antigen (HLA) class I ligands bind to killer immunoglobulin-like receptors (KIRs) on the NK cell surface . This highly specific recognition system is controlled by the integration of signals generated by a multitude of inhibitory and activating KIRs, which inhibit or activate, respectively, cytotoxicity and secretion of cytokines ultimately leading to death of the targeted cell . Both inhibition and activation involve a number of signalling molecules, as previously described [8, 9].
There is extensive genomic diversity in KIR genes in humans. Currently, a database and online repository for immune gene frequencies in worldwide populations reports 517 different KIR genotypes [10, 11]. It is believed that this variation may affect resistance or susceptibility to a number of pathogens through ligand-receptor interactions and the downstream signalling and/or cytokine release that follows [12, 13]. Genetic susceptibility or resistance to infectious diseases, in conjunction with environmental and host risk factors, is thought to determine disease progression [14–16].
Present literature indicates that the outcome of Mtb infection is affected not only by virulence of the infecting strain of M. tuberculosis , but also by host environment, disease co-morbidities, and the genetic composition of the host, specifically the presence or absence of genes that regulate the immune system [14, 16, 18–20]. Following Mtb infection, approximately 10 % of individuals will develop active TB (ATB) during their lifetime, while the majority of individuals will exhibit latent TB infection (LTBI) [21, 22]. LTBI refers to the condition in which Mtb remains viable in the macrophage but retains a small amount of metabolic activity . It is not currently known which genes and/or immune components regulate an individual’s disease outcome following exposure (ATB, LTBI, or exposed uninfected). Present literature captures only those studies focusing on genetic profiles among active tuberculosis vs. uninfected individuals. In the majority of these studies, the control group contains both individuals with LTBI infection as identified by a positive Tuberculin skin test (TST), and those with uninfected status.
The novel aspect of this study is to identify unique profiles among the LTBI population, diagnosed using the Interferon gamma release assay (IGRA)' as there is identify twice in this sentence. Differences in KIR profiles and haplotypes may be associated with Mtb infection status [24–26] and play a role in altered TB disease progression and disease outcomes. In this cross-sectional study, we examined the enrichment or depletion of KIR genes in individuals from Manitoba with ATB infection, LTBI and controls, and further explored the association between Mtb infection status and KIR profiles and haplotypes.
The 209 samples consisted of whole blood from individuals living in Manitoba. The sampling was performed at hospital and community TB clinics in Winnipeg, Manitoba, Canada between November 3, 2009 and March 29, 2011 and was cross-sectional in nature. The study was approved by the Health Research Ethics Board at the University of Manitoba (H2008:301). All study participants provided written informed consent following consultation with a study nurse. ATB infection (n = 59) was confirmed by mycobacterial culture. LTBI (n = 46) was identified using the interferon-gamma release assay (IGRA) (QuantiFERON®-TB-Gold, Qiagen). Healthy IGRA negative HIV negative individuals (n = 104) were used as controls and consisted largely of individuals participating in routine occupational health screening, and immigration screening. All individuals within the specified time period who consented to genetic testing were included in this study. Exclusion criteria included those individuals with HIV co-infection, and anyone who exhibited an indeterminate IGRA response. Participant demographics for ethnicity, age, and gender can be seen in Table 1.
DNA extraction and replication
Genomic DNA was extracted using Qiagen DNA Mini Kit as per manufacturer’s instructions (Qiagen, Louisville, KY). The samples were subjected to whole genome replication using the Qiagen Repli-G mini kit as per manufacturer’s instructions to increase DNA concentration of the testing sample.
The concentration of DNA was normalized to 100 μg/mL at 260 nm using the SmartSpec Plus spectrophotometer (Bio-Rad, Mississauga, ON). KIR genotyping was performed by sequence-specific primer polymerase chain reaction using the Miltenyi Biotec KIR Typing Kit (Auburn, CA) as previously described . The amplicons were visualized with UV light (Bio-Rad Gel Doc EZ Imager, Mississauga, ON) following gel electrophoresis at 13 V/cm on a 2 % agarose gel containing ethidium bromide. The KIR typing kit allows for detection of all known human KIR genes and alleles [27, 28]. KIR2DL5A and KIR2DL5B are collectively referred to as KIR2DL5 for this paper.
Data for each individual was entered into BioNumerics software version 5.0 (Applied Maths, Belgium) as binary character data. All KIR genes were combined into a single KIR profile for each individual and clustered to identify prevalent profiles among specified groups using the categorical co-efficient and unweighted pair group method with arithmetic mean (UPGMA) . KIR gene frequencies were tabulated by direct counts from the clustered profiles to determine frequency within a defined group. Differences between Mtb infection status groups were estimated using the two-tailed Fisher’s exact test (GraphPad Software, La Jolla, CA). A P-value ≤0.05 was considered statistically significant. Haplotype designation was determined as previously described .
KIR gene frequencies
In order to determine the differences in KIR gene frequencies between different disease status groups (ATB, LTBI, and controls), the KIR gene frequency data obtained was analyzed and compared. All 209 samples consistently contained the framework genes KIR2DL4, KIR3DL2, KIR3DL3, and the pseudogenes KIR2DP1 and KIR3DP1.
Five KIR genes (KIR2DL2, KIR2DL5, KIR2DL5B, KIR2DS2, and KIR2DS3) differed significantly (P ≤ 0.05) in frequency between disease status groups (Table 2). Two genes differed between individuals with Mtb infection (LTBI and ATB) vs. controls, KIR2DL2 (33.33 % vs. 55.77 %, P = 0.0014) and KIR2DS2 (34.29 % vs. 54.81 %, P = 0.0035). However, the underlying differences can be exposed when analyzing LTBI and ATB separately. KIR2DL5 and KIR2DL5B (both 73.91 % vs. 51.92 %, P = 0.0125) were present in higher frequency in individuals within the LTBI group as compared to controls. KIR2DL2 (27.12 % vs. 55.77 %, P = 0.0005), KIR2DS2 (27.12 % vs. 54.82 %, P = 0.0010), and KIR2DS3 (8.47 % vs. 30.77 %, P = 0.0009) were present in a lower frequency in individuals with ATB compared to controls. Lastly, gene frequencies of KIR2DL5 (73.91 % vs. 49.15 %, P = 0.0156), KIR2DL5B (73.91 % vs. 49.15 %, P = 0.0156), and KIR2DS3 (39.13 % vs. 8.47 %, P = 0.0002) differed significantly between latently and actively infected individuals, respectively.
KIR gene profiles
Forty-three KIR profiles (genotypes) were identified in this study (Fig. 1). These profiles ranged in their frequency of distribution from as high as 24.88 % (52/209) to as low as 0.48 % (1/209). Twenty-two of the 43 profiles identified were unique to a single individual. The most prominent genotypes were #7, (18/209, 8.2 %), #8 (52/209, 24.9 %), #12 (11/209, 5.3 %), #18 (12/209, 5.7 %), and #36 (25/209, 12.0 %). Eight KIR genotypes were shared between all three disease status groups (ATB, LTBI, controls; genotypes # 7, 8, 12, 17, 18, 30, 36, 39) (Fig. 2). Excluding those shared with individuals from the control group, two genotypes were shared between the LTBI and ATB groups (genotypes #11, 22). Nineteen genotypes were exclusive to the control group. Five genotypes were exclusive to those individuals with ATB (genotypes #1, 4 9, 23, 37), and represented 6/59 (10.2 %) active cases. Three genotypes were exclusive to those individuals with LTBI (genotypes # 13, 21, 28), and represented 3/46 (6.5 %) latent cases.
Those individuals with Mtb infection (LTBI or ATB) were found within 25 of the 43 genotypes, most predominantly in #8 (25/105, 23.8 %) and #36 (20/105, 19.0 %). Over 42 % of Mtb infected individuals were found in these two genotypes.
In addition to gene frequency variation, there is haplotypic variation due to the different number and kinds of KIR genes . Both LTBI (34.8 %, P = 0.0004) and ATB (32.3 %, P = 0.0005) infected individuals were significantly more likely to have AA-AB haplotypes than controls (9.6 %; Table 3). Additionally, those individuals with ATB were less likely to have an AB-AB haplotype compared to controls (8.5 % vs. 21.2 %, P = 0.0476). The majority of the Mtb infected cases are represented in the AA-AA (LTBI – 21.7 %, ATB – 35.6 %) and AA-AB (LTBI – 34.8 %, ATB – 32.3 %) haplotypes. Overall, 66.67 % of Mtb infected individuals had a centromeric-AA haplotype (LTBI – 58.7 %, ATB – 72.9 %), compared with only 44.2 % of controls (P = 0.0014).
This study was designed to investigate the KIR gene frequencies in individuals from Manitoba with ATB, LTBI, and a control group, as described in methods. Additionally, KIR profiles and haplotypes were analyzed. KIR genes may influence disease outcome (latent vs. active) which is controlled in part by an organized immune response.
When determining KIR gene frequencies, framework genes and pseudogenes were present in 100 % of the samples, as expected . KIR2DL2, KIR2DL5, KIR2DL5B, KIR2DS2, and KIR2DS3 differed significantly between Mtb status groups.
Mahfouz et al. and Mendez et al. both found KIR2DL3 to be the only statistically significant KIR gene frequency to differ between ATB patients and controls (higher in ATB patients; P = 0.03 and P = 0.02, respectively) [24, 25]. In our study, KIR2DL3 occurred in only a slightly higher frequency of individuals with ATB (98.31 %) compared to controls (91.35 %, P = 0.0957; NS). Pydi et al. found gene frequencies of KIR2DS1, KIR2DS5, KIR3DL1, and KIR2DL3 to be higher in ATB patients compared to controls (P < 0.0001) . Although KIR2DS5 and KIR2DL3 gene frequencies differed slightly between ATB and control groups in our study, none of the observed differences were statistically significant. Lastly, Lu et al. found gene frequencies of KIR2DS1, KIR2DS3, and KIR3DS1 to be significantly higher in ATB patients compared to controls (P < 0.05) . In our study, frequencies of KIR2DS1 and KIR3DS1 were not statistically different between those with ATB compared to controls. However, when looking at KIR2DS3, our study showed a decreased frequency in individuals with ATB compared controls (8.47 % vs. 30.77 %, P = 0.0009). Observed differences in KIR frequencies between TB status groups may be due in large part to ethnicity, as in each of the above mentioned studies, a different country of origin was involved, and KIR are known to differ among ethnic groups . Given that the Manitoba population is very heterogeneous consisting of many foreign born immigrant individuals, as well as indigenous and Canadian born populations, future studies matching TB status groups by ethnicity will help to more clearly define the role of KIR genes in TB pathogenicity. Additionally, group definitions may have played a role; the control groups in the above published studies contained either TST positive individuals or had no data on TB reactivity. These control group then reflect a mixture of both uninfected and LTBI infected participants. In contrast to this, with the use of the IGRA testing, we were able to clearly distinguish our LTBI and control groups.
With the exception of KIR2DL2, our study found an increased presence of inhibitory KIR (KIR2DL5, KIR2DL5B) in LTBI individuals and a decreased presence of activating KIR (KIR2DS2, KIR2DS3) in ATB infected individuals. This may suggest that the enrichment or depletion of specific KIR genes predisposes an individual to progressing to ATB disease by means of an inadequate cytotoxic response to the pathogen. LTBI refers to the condition in which Mtb remains viable in the macrophage but only retains a small amount of metabolic activity . The inability of the immune system to maintain the infection in a latent state results in ATB infection.
Forty-three different gene profiles were identified in the 209 samples, of which, 25 profiles contained Mtb cases. The profiles containing the most Mtb cases were also prevalent in the control group, suggesting an unlikely correlation between profile/genotype and TB status. Many profiles were unique to individuals with Mtb, however there is little to be concluded from those profiles containing only a few individuals. Extrapolation of these findings via continued sampling is warranted to determine the importance of KIR profiles.
Those individuals with LTBI and ATB were more likely to have an AA-AB haplotype than controls. This was the haplotype that contained the most Mtb cases (35/105, 33.3 %). Two-thirds (66.7 %) of individuals with TB had a centromeric-AA halpotype, compared to only 44.2 % of controls (P = 0.0014). A centromeric-AA haplotype represents the haplotype with the fewest number of activating genes. It is hypothesized that this lack of activating genes may prevent the appropriate release of M. tuberculosis killing cytokines .
A limitation to this study is the lack of longitudinal data among our LTBI status group. As we do not know when an individual became LTBI (IGRA positive) it is possible that some of these individuals went on to develop primary or secondary TB, however we do not have access to this data. It is known that these were healthy individuals with low risk for the development of ATB. Another limitation is the unknown TB exposure of our control group (IGRA negative), we can make no conclusions in regards to the role KIR plays in TB susceptibility. This group is used as a reference comparison group for our LTBI and ATB populations.
In summary, major differences can be seen in KIR gene frequencies across Mtb disease status groups. KIR haplotype frequencies differ between these groups as well. The differences in KIR gene frequencies and/or haplotypes may result in differential cytokine expression, contributing to different disease outcomes, and suggest a genetic influence on Mtb susceptibility and pathogenesis. The skewed distribution of A-containing centromeric haplotypes (containing fewer activating genes), along with the increased presence of TB disease in these haplotypes, suggests a correlation. Further investigation is needed to characterize the subtleties of these differences by way of sequencing of specific KIR genes and/or KIR-HLA association studies taking into account different ethnic populations.
Availability of supporting data
The data set supporting the results of this article is included within the article.
Killer immunoglobulin-like receptor
Human leukocyte antigen
Latent tuberculosis infection
Tuberculin skin test
Unweighted pair group method with arithmetic mean
World Health Organization. Global Tuberculosis Report 2013. Geneva, Switzerland: WHO Press; 2013.
Tuberculosis in Canada 2012 - Pre-Release – Public Health Agency of Canada [http://www.phac-aspc.gc.ca/tbpc-latb/pubs/tbcan12pre/tab-eng.php#tab1]
Martin MP, Carrington M. KIR locus polymorphisms: genotyping and disease association analysis. Methods Mol Biol. 2008;415:49–64.
Raja A. Immunology of tuberculosis. Indian J Med Res. 2004;120:213–32.
Flynn JL, Chan J. Immunology of tuberculosis. Annu Rev Immunol. 2001;19:93–129.
Vilches C, Gardiner CM, Parham P. Gene structure and promoter variation of expressed and nonexpressed variants of the KIR2DL5 gene. J Immunol. 2000;165:6416–21.
Rajalingam R. Overview of the killer cell immunoglobulin-like receptor system. Methods Mol Biol. 2012;882:391–414.
Kelley J, Walter L, Trowsdale J. Comparative genomics of natural killer cell receptor gene clusters. PLoS Genet. 2005;1:129–39.
Braun K, Larcombe L, Orr P, Nickerson P, Wolfe J, Sharma M. Killer Immunoglobulin-Like Receptor (KIR) Centromeric-AA Haplotype Is Associated with Ethnicity and Tuberculosis Disease in a Canadian First Nations Cohort. PLoS ONE. 2013;8:e67842.
The Allele Frequency Net Database [KIR Genotype Reference] [http://allelefrequencies.net/kir6001a.asp]
Gonzalez-Galarza FF, Christmas S, Middleton D, Jones AR. Allele frequency net: a database and online repository for immune gene frequencies in worldwide populations. Nucleic Acids Res. 2011;39(Database issue):D913–9.
Bashirova AA, Martin MP, McVicar DW, Carrington M. The killer immunoglobulin-like receptor gene cluster: tuning the genome for defense. Annu Rev Genomics Hum Genet. 2006;7:277–300.
Carrington M, Norman P. The KIR Gene Cluster. Bethesda, MD: National Centre for Biotechnology Information; 2003.
Delgado JC, Baena A, Thim S, Goldfeld AE. Ethnic‐Specific Genetic Associations with Pulmonary Tuberculosis. J Infect Dis. 2002;186:1463–8.
Middleton D, Gonzelez F. The extensive polymorphism of KIR genes. Immunology. 2010;129:8–19.
Morris GAJ, Edwards DRV, Hill PC, Wejse C, Bisseye C, Olesen R, et al. Interleukin 12B (IL12B) genetic variation and pulmonary tuberculosis: a study of cohorts from The Gambia, Guinea-Bissau, United States and Argentina. PLoS ONE. 2011;6:e16656.
Petrelli D, Kaushal Sharma M, Wolfe J, Al-Azem A, Hershfield E, Kabani A. Strain-related virulence of the dominant Mycobacterium tuberculosis strain in the Canadian province of Manitoba. Tuberculosis (Edinb). 2004;84:317–26.
Awomoyi AA, Marchant A, Howson JMM, McAdam KPWJ, Blackwell JM, Newport MJ. Interleukin-10, polymorphism in SLC11A1 (formerly NRAMP1), and susceptibility to tuberculosis. J Infect Dis. 2002;186:1808–14.
Lykouras D, Sampsonas F, Kaparianos A, Karkoulias K, Tsoukalas G, Spiropoulos K. Human genes in TB infection: their role in immune response. Monaldi Arch Chest Dis. 2008;69:24–31.
Uciechowski P, Imhoff H, Lange C, Meyer CG, Browne EN, Kirsten DK, et al. Susceptibility to tuberculosis is associated with TLR1 polymorphisms resulting in a lack of TLR1 cell surface expression. J Leukoc Biol. 2011;90:377–88.
Public Health Agency of Canada. Canadian Lung Assocation/Canadian Thoracic Society: Canadian Tuberculosis Standards, 6th edition. Canada: Minister of Health; 2007.
Versalovic J: Manual of Clinical Microbiology. 10th edition. Washington, DC: ASM Press; 2011.
Vernon A. Treatment of latent tuberculosis infection. Semin Respir Crit Care Med. 2013;34:67–86.
Mahfouz R, Halas H, Hoteit R, Saadeh M, Shamseddeen W, Charafeddine K, et al. Study of KIR genes in Lebanese patients with tuberculosis. Int J Tuberc Lung Dis. 2011;15:1688–91.
Méndez A, Granda H, Meenagh A, Contreras S, Zavaleta R, Mendoza MF, et al. Study of KIR genes in tuberculosis patients. Tissue Antigens. 2006;68:386–9.
Shahsavar F, Mousavi T, Azargon A, Entezami K. Association of KIR3DS1 + HLA-B Bw4Ile80 Combination with Susceptibility to Tuberculosis in Lur Population of Iran. Iran J Immunol. 2012;9:39–47.
KIRtyping_23June08.indd - IM0001493.ashx.
KIR Typing Kits - Miltenyi Biotec [https://www.miltenyibiotec.com/en/products-and-services/macsmolecular/nucleic-acid-research/dna-research/kir-typing-kits.aspx]
Clifford HT: An Introduction to Numerical Classification. New York, NY: Academic Press; 1975.
Robinson J, Mistry K, McWilliam H, Lopez R, Marsh SGE. IPD--the Immuno Polymorphism Database. Nucleic Acids Res. 2010;38(Database):D863–9.
Hsu KC, Liu X-R, Selvakumar A, Mickelson E, O’Reilly RJ, Dupont B. Killer Ig-like receptor haplotype analysis by gene content: evidence for genomic diversity with a minimum of six basic framework haplotypes, each with multiple subsets. J Immunol. 2002;169:5118–29.
Pydi SS, Sunder SR, Venkatasubramanian S, Kovvali S, Jonnalagada S, Valluri VL: Killer Cell Immunoglobulin like Receptor Gene Association with Tuberculosis. Hum Immunol 2012;74:85–92.
Lu C, Bai XL, Shen YJ, Deng YF, Wang CY, Fan G, et al. Potential implication of activating killer cell immunoglobulin-like receptor and HLA in onset of pulmonary tuberculosis. Scand J Immunol. 2012;76:491–6.
The authors wish to acknowledge the study nurse Ms. Carmen Lopez and research technician Ms. Christine Mesa. The authors would also like to thank Drs. T. Blake Ball, Linda Larcombe, and Pamela Orr for research support.
The authors declare that they have no competing interests.
MS and JW participated in the design and coordination of the study. KB performed the laboratory and statistical analysis, and drafted the manuscript. KB performed this work as a graduate student under mentorship of MS. SK was responsible for obtaining study ethics. All authors provided input into the final manuscript.