Skip to main content

Polymorphisms in IL-1β, vitamin D receptor Fok1, and Toll-like receptor 2 are associated with extrapulmonary tuberculosis



Human genetic variants may affect tuberculosis susceptibility, but the immunologic correlates of the genetic variants identified are often unclear.


We conducted a pilot case-control study to identify genetic variants associated with extrapulmonary tuberculosis in patients with previously characterized immune defects: low CD4+ lymphocytes and low unstimulated cytokine production. Two genetic association approaches were used: 1) variants previously associated with tuberculosis risk; 2) single nucleotide polymorphisms (SNPs) in candidate genes involved in tuberculosis pathogenesis. Single locus association tests and multifactor dimensionality reduction (MDR) assessed main effects and multi-locus interactions.


There were 24 extrapulmonary tuberculosis cases (18 black), 24 pulmonary tuberculosis controls (19 black) and 57 PPD+ controls (49 black). In approach 1, 22 SNPs and 3 microsatellites were assessed. In single locus association tests, interleukin (IL)-1β +3953 C/T was associated with extrapulmonary tuberculosis compared to PPD+ controls (P = 0.049). Among the sub-set of patients who were black, genotype frequencies of the vitamin D receptor (VDR) Fok1 A/G SNP were significantly different in extrapulmonary vs. pulmonary TB patients (P = 0.018). In MDR analysis, the toll-like receptor (TLR) 2 microsatellite had 76% prediction accuracy for extrapulmonary tuberculosis in blacks (P = 0.002). In approach 2, 613 SNPs in 26 genes were assessed. None were associated with extrapulmonary tuberculosis.


In this pilot study among extrapulmonary tuberculosis patients with well-characterized immune defects, genetic variants in IL-1β, VDR Fok1, and TLR2 were associated with an increased risk of extrapulmonary disease. Additional studies of the underlying mechanism of these genetic variants are warranted.

Peer Review reports


A possible genetic predisposition to tuberculosis has been suggested in several studies,[13] but the functional immunologic correlate of the genetic polymorphisms identified is often unclear[4]. We have sought to first identify immunologic defects that may predispose to tuberculosis, then assess for genetic polymorphisms associated with these immunologic defects. Because there is evidence that extrapulmonary tuberculosis is the result of an underlying immune defect, [5] we have focused our search on persons with prior extrapulmonary disease. In previous studies we noted that HIV-seronegative adults with prior extrapulmonary tuberculosis had lower levels of CD4+ lymphocytes, unstimulated cytokine production, and tumor necrosis factor (TNF)-α production in response to lipopolysaccharide (LPS) or LPS + interferon (IFN)-γ than persons with prior pulmonary tuberculosis or latent M. tuberculosis infection[6, 7].

In the current study we have combined two previous study populations in which cytokine responses were well-characterized, [6, 7] to identify genetic polymorphisms associated with extrapulmonary tuberculosis, and by extension, the associated immunologic abnormalities. We make use of a candidate gene approach utilizing single nucleotide polymorphisms (SNPs) and microsatellites previously reported to be associated with tuberculosis, and we also include SNPs in candidate genes that are hypothesized to play a role in tuberculosis pathogenesis.

Because several loci could contribute to the phenotype expressed in complex diseases such as tuberculosis, it is important to search for gene-gene interactions (epistasis), which may give a more accurate prediction of disease risk. Such interactions are difficult to detect using traditional statistical methods such as single locus association tests because those tests were not developed to detect purely interactive effects. Such tests identify genes with main effects and then follow-up analyses assess for interactions between genes that exhibit a main effect. New statistical and computational methods that have better power to detect interactions, including those without main effects, in relatively small sample sizes, are required. Multifactor dimensionality reduction (MDR) analysis is a novel method developed to address this need[8]. In the current study both single locus association tests and MDR analysis were used to detect potential single or multi-locus interactions that predict extrapulmonary tuberculosis.


Clinical and Demographic Characteristics

There were 24 extrapulmonary tuberculosis cases (18 black), 24 pulmonary tuberculosis controls (19 black), and 57 controls (49 black) with latent M. tuberculosis infection. The clinical and demographic characteristics of the study population are in Table 1. There were no significant differences between cases and controls according to age, sex, or race. However, pulmonary controls had lower median body mass index (BMI) at the time of study than either extrapulmonary tuberculosis cases or controls with latent M. tuberculosis infection. Extrapulmonary tuberculosis cases had lower median CD4+ lymphocytes at time of study than controls. The sites of disease among extrapulmonary cases included lymphatic (n = 9), pleural (n = 4), bone/joint (n = 4), pericardial (n = 2), genito-urinary (n = 2), miliary/meningeal (n = 1), gastrointestinal (n = 1), and laryngeal (n = 1).

Table 1 Clinical and Demographic Characteristics of the Study Population

SNPs Associated with Tuberculosis in Previous Studies

Of the 27 SNPs and 3 microsatellites previously reported in the literature as being associated with tuberculosis (Table 2), 25 SNPs and all 3 microsatellites were assessed by PCR-based SNP genotyping. The IL-1β -511 (C→T) and vitamin D receptor (VDR) Bsm1 SNPs were on the GeneChip® Human Mapping 50K Xba Array. Of the 30 polymorphisms tested, 5 SNPs (rs 17235416, rs1800450, rs5743708, rs6265786, rs2393799) were unable to be genotyped after quality control and were therefore excluded. Only markers with < 10% missing data were considered for analysis. No imputation of missing values was performed. For univariate analyses, missing values were removed for analysis; for the MDR analysis, missing observations were ignored, such that attribute construction and accuracy calculations were only performed based on observed values.

Table 2 Polymorphisms Tested that had Previously Been Reported as Associated with Tuberculosis

In single locus association testing of each polymorphism among all study participants, the IL-1β +3953 (C→T) SNP was associated with any tuberculosis (P = 0.044) and extrapulmonary tuberculosis (P = 0.049) compared to PPD+ controls (Table 3). Among black patients, the Fok1 SNP in the VDR gene was associated with extrapulmonary tuberculosis compared to pulmonary tuberculosis controls (P = 0.018).

Table 3 Of the Polymorphisms Tested that had Previously Been Reported as Associated with Tuberculosis (Table 2), the SNPs Associated with Tuberculosis in this Study

In MDR analysis of the 25 polymorphisms among all study participants, the microsatellite in TLR2 (GT) had 61% prediction accuracy for any tuberculosis compared to controls with M. tuberculosis infection (P = 0.038) (Table 3). It also had 63% prediction accuracy for any tuberculosis in blacks (P = 0.047) and 76% prediction accuracy for extrapulmonary tuberculosis in blacks (P = 0.002). Results did not change when the analysis included CD4+ lymphocytes and BMI as potential covariates (results not shown).

SNPs in Candidate Genes

There were 661 SNPs in genes presumed to be associated with tuberculosis pathogenesis (Table 4), of which 613 (93%) had genotyping efficiency > 90% (i.e., < 10% missing data) and were included in the analysis. In the MDR analysis of these 613 SNPs, an intronic SNP in TNF-α (C→T) (rs1811063) predicted tuberculosis with 62% accuracy among all patients (P = 0.05; Table 5). The combination of two SNPs in TLR4 (C→T) (rs1399431) and TNF-α (C→T) (rs7791836) predicted tuberculosis risk with 71% accuracy in blacks (P = 0.02). Results did not change when the analysis included CD4+ lymphocytes and BMI as potential covariates (results not shown). There were no SNPs associated with extrapulmonary tuberculosis compared to controls.

Table 4 Polymorphisms in the Affymetrix GeneChip® Human Mapping 50K Xba Array in Genes Hypothesized to be Associated with Tuberculosis Pathogenesis
Table 5 Of the SNPs Listed in Table 4, the SNPs Associated with Tuberculosis in This Study Population


In this study population we first identified immunologic defects among persons with extrapulmonary tuberculosis, then explored possible associations with candidate genetic variants. In single locus association tests among genetic variants previously associated with tuberculosis susceptibility, IL-1β +3953 was associated with extrapulmonary disease, as well as all forms of tuberculosis, compared to persons with M. tuberculosis infection. Among black patients, the Fok1 SNP in the VDR gene distinguished extrapulmonary from pulmonary disease. The MDR analysis provides important additional results, however, because it is more powerful than traditional logistic regression analysis,[9] adjusts for multiple comparisons via permutation testing, and because the results were more consistent across the groups assessed. MDR also evaluates potential gene-gene and gene-environment interactions, which are important to investigate in complex phenotypes such as tuberculosis pathogenesis. The TLR2 microsatellite predicted extrapulmonary tuberculosis in black patients with 76% accuracy. It also predicted all forms of tuberculosis in blacks, as well as in the full study population. The prediction accuracy estimates the importance of variables in the dataset, with the intent of generalizing the model to the full population.

These results point to the importance of IL-1β +3953, VDR Fok1 A/G, and the TLR2 microsatellite in extrapulmonary tuberculosis risk, with the latter two being particularly significant in blacks, even in such a small sample size. This is important in light of the increased proportion of extrapulmonary tuberculosis in blacks[5]. Further studies are needed to determine if these polymorphisms could account for this epidemiological finding. This study included very few individuals of other racial backgrounds, making comparisons of genetic models between racial groups impossible. Our results are consistent with the previously noted association between IL-1β +3953 and tuberculous pleurisy, and VDR Fok1 and several forms of extrapulmonary tuberculosis--both noted among in Gujarati Asians living in London[10, 11]. They are also consistent with the recently identified link between TLR and the innate immune response to M. tuberculosis: the relationship between TLR signaling, the up-regulation of the VDR, and vitamin D-mediated killing of intracellular M. tuberculosis via the microbial peptide cathelicidin[12]. In that study, blacks had low 25-hydroxyvitamin D levels and low cathelicidin messenger RNA induction. While immediate biological interpretation of our results is not possible, defects in M. tuberculosis recognition and/or subsequent intracellular signaling in the nuclear factor kappa B (NF-κB) pathway related to TLR2 polymorphisms would be consistent with the subtle innate immune defects in the extrapulmonary tuberculosis patients in this study population: low unstimulated cytokine production, decreased TNF-α production in response to LPS or LPS + IFN-γ, and low CD4+ lymphocytes[6, 7, 4].

The results from SNPs in genes hypothesized to play a role in tuberculosis pathogenesis showed that no SNPs were significantly associated with extrapulmonary tuberculosis. There was, however, a SNP in TNF-α that was associated with tuberculosis among all study participants. And among blacks, a combination of SNPs in TNF-α and TLR4 predicted tuberculosis risk with 71% accuracy. These SNPs cannot be interpreted as easily in light of the immunologic findings, which pertain to extrapulmonary tuberculosis. However, they suggest that similar cytokine pathways as those noted above are important in tuberculosis pathogenesis.

It was notable that persons with prior extrapulmonary tuberculosis or latent M. tuberculosis infection had significantly higher body mass index (BMI) than persons with pulmonary tuberculosis. If BMI after completion of treatment is comparable to BMI prior to development of disease, it would suggest that lower BMI may predispose to pulmonary tuberculosis. It should be noted that the median BMI in persons with pulmonary disease was within the normal range, whereas BMI in persons in the other two groups was high. A recent large population-based study in Hong Kong found that as BMI increased, the risk of pulmonary tuberculosis decreased, but the risk of extrapulmonary tuberculosis did not[13]. The protective effect of increased BMI on pulmonary disease in that study persisted even after controlling for confounding variables, including smoking and diabetes mellitus. The above findings suggest that risk factors for extrapulmonary tuberculosis are not ameliorated by increased BMI--which would be consistent with an immunogenetic predisposition to extrapulmonary tuberculosis. Factors related to the M. tuberculosis strain also play a role in extrapulmonary disease[14].

The most important limitation of the current study was the small sample size. This was driven by the intensive immunologic testing performed on specimens from each participant. Statistically significant genetic associations were indeed identified, but the sample was under-powered (with either traditional methods or MDR analysis) to detect other potential associations. Additionally, even for the positive associations detected, reliable estimates of effect size cannot be determined, due to the "winner's curse" related to the small sample size[15]. Table 6 shows the detectable effect sizes with the current sample size for the variants previously associated with tuberculosis. These calculations assume 80% power and an uncorrected alpha of 0.05%. Larger studies are needed in which candidate gene polymorphisms and associated cytokine pathways are fully evaluated in all patients, so that any genetic polymorphisms identified can be interpreted in the context of the cytokine abnormalities noted. Given this, one should consider the results of this study--particularly those from the candidate gene approach and interaction analysis--as hypothesis-generating rather than hypothesis-confirming. While MDR has greater statistical power than traditional methods, its power to detect gene-gene and gene-environment interactions is less than that to detect main effects. In searching for interactive effects, the small sample size did not allow for investigating interactions with more than two loci.

Table 6 Odds Ratios that Could Have Been Detected Given the Sample Size and Minor Allele Frequency in the Study Population

There were other limitations of this study. First, of the SNPs tested that had previously been associated with tuberculosis, most pertained to pulmonary rather than extrapulmonary disease. Due to possible differences in the pathophysiology of these two disease manifestations, [7] one might expect that not all of the SNPs tested would be related to the pathogenesis of extrapulmonary tuberculosis. Second, not all genes that are presumably associated with tuberculosis pathogenesis had SNPs included in the GeneChip® Human Mapping 50K Xba Array, so not all such genes could be assessed. Third, we were unable to directly incorporate cytokine response data into the genetic analyses because the methodology used to quantify cytokine responses was not the same in the two immunologic studies.

Additionally, concerns with multiple testing arise when screening such a large number of genetic variants. The results of the first approach are presented without any corrections for multiple testing. Because these variants have been associated with tuberculosis before, each test represents an independent statistical hypothesis. By using MDR and permutation in the second approach, P values were empirically derived based on the total number of tests at each stage.


The results of this study suggest that an evaluation of the underlying mechanism(s) of the genetic variants in IL-1β +3953, VDR Fok1 A/G, and the TLR2 microsatellite--and their role in the pathogenesis of extrapulmonary tuberculosis--is warranted. Comprehensive immunogenetic studies will contribute to our understanding of tuberculosis pathogenesis, and may allow us to identify persons at highest risk of developing tuberculosis.


Study participants

The study population was pooled from two immunologic studies that have been described previously; the inclusion criteria for both studies were similar and are described in detail elsewhere[6, 7]. Briefly, patients were identified through the Baltimore City Health Department Eastern Chest Clinic and Nashville Metropolitan Health Department Tuberculosis Clinic. In this case-control study, eligibility criteria for case patients included a history of treated culture-confirmed extrapulmonary tuberculosis, age ≥ 18 years old, and HIV-seronegative status. Extrapulmonary disease was defined as any site outside of the pulmonary parenchyma. Exclusion criteria included serum creatinine >2 mg/dL, use of corticosteroids or other immunosuppressive agents at the time of diagnosis or time of study entry, malignancy, or diabetes mellitus. The criteria for pulmonary tuberculosis control patients included HIV-seronegative adults ≥ 18 years old who had completed treatment for culture-confirmed pulmonary tuberculosis, and had no evidence of extrapulmonary tuberculosis. Positive cultures of sputum, bronchoalveolar lavage, or pulmonary parenchyma were required. Controls with latent M. tuberculosis infection were ≥ 18 years old, HIV-seronegative, and had a positive tuberculin skin test (defined as ≥ 10 mm induration after intradermal placement of 5 tuberculin units of PPD) without evidence of active tuberculosis. Participants in this control group were U.S.-born (and therefore not vaccinated with BCG), and were mostly close contacts of tuberculosis cases. Exclusion criteria for both control groups were the same as for the case group. Controls were drawn from the same two clinic populations as cases, and were not related to the cases. Extrapulmonary tuberculosis cases and pulmonary controls completed treatment prior to study entry.

This study was approved by the institutional review boards of the Johns Hopkins Hospital, the Baltimore City Health Department, the National Institutes of Health, Vanderbilt University Medical Center, and the Nashville Metropolitan Health Department. All study participants provided written informed consent.

Sample Collection

CD4+ lymphocytes were quantified by flow cytometry. DNA was extracted from blood samples using the Puregene DNA Isolation kit (Gentra Systems, Minneapolis, MN) following the manufacturer's protocol. Genomic DNA was stored at -70°C until genotyping. Laboratory personnel were blinded to the case-control status of the specimens.

SNP Genotyping

As of November 15, 2006, 27 SNPs were identified in the literature as being associated with tuberculosis (Table 2). SNPs not included in the GeneChip® Human Mapping 50K Xba Array (see below) were genotyped using validated TaqMan SNP genotyping assays from Applied Biosystems. For each SNP, 25 ng DNA was used. The polymerase chain reaction (PCR) was performed in an ABI 9700 thermocycler under the following conditions: 95°C for 10 minutes followed by 50 cycles of 92°C for 15 seconds and 60°C for 1 minute. The 384-well plates were read on an ABI 7900HT sequence detection system according to manufacturer's manual. The probes were labeled with FAM or VIC dye at the 5' end and a minor-groove binder and non-fluorescent quencher at the 3' end.

Microsatellite Genotyping

As of November 15, 2006, 3 microsatellites were identified in the literature as being associated with tuberculosis (Table 2). A 5'(GT)n microsatellite associated with NRAMP1, a VNTR associated with interleukin (IL)-1RA, and a GT repeat polymorphism in the toll-like receptor (TLR) 2 gene were all genotyped using fluorescent fragment analysis [16]. All 3 microsatellites were amplified by PCR. The PCR primers and conditions for the 5'(GT)n and IL-1RA microsatellites were previously reported[17, 10]. The TLR2 microsatellite was amplified using the following conditions: 1 minute at 96°C; 30 cycles of 94°C for 1 minute, 1 minute at 53°C, and 2 minutes at 70°C; and a final elongation period of 10 minutes at 70°C. The primers used to amplify the TLR2 microsatellite were previously published [18]. The forward primers for NRAMP1, IL-1RA and TLR2 microsatellites (5'-ACT CGC ATT AGG CCA ACG AG-3', 5'-CTC AGC AAC ACT CCT AT-3', and 5'-GCA TTG CTG AAT GTA TCA GGG A-3' respectively) were all 5' labeled with FAM purchased from MWG Biotech. Following PCR, the amplicons were electrophoresed on a 3730s DNA analyzer (Applied Biosystems) and analyzed with GeneMapper 4.0 (Applied Biosystems).

Candidate Gene Genotyping

Twenty six candidate genes were identified based on their potential role in tuberculosis immunopathogenesis and the inclusion of SNPs in these genes in the GeneChip® Human Mapping 50K Xba Array (Affymetrix, Inc., Santa Clara, CA). A total of 661 SNPs in these 26 candidate genes were present (Table 4). Genotyping was performed according to the manufacturer's protocol. Briefly, a complexity reduction process was performed where genomic DNA (250 ng) was digested with XbaI, ligated to XbaI adaptor (Affymetrix), and amplified by polymerase chain reaction (PCR) using primers specific to the ligated adaptor. Cycling conditions were an initial denaturation of 94°C for 3 minutes followed by 30 cycles of 94°C for 30 seconds, 60°C for 45 seconds, and 68°C for 60 seconds. A final extension of 68°C for 7 minutes concluded the reactions. PCR products were assayed by gel electrophoresis, purified, fragmented to < 250 bp using dilute DNaseI (Affymetrix), biotin end-labeled with terminal deoxynucleotidyl transferase, and hybridized to the 50K Xba Array at 48°C for 16 hours at 60 rpm. The hybridized arrays were washed and stained on Fluidics Station 450 and scanned with the GeneChip Scanner according to the manufacturer's settings (Affymetrix). The arrays were analyzed with software GDAS version 3.0.2 (Affymetrix), which provides rank scores for the probability of particular genotypes at SNP loci. The scores were AA or BB for homozygous alleles and AB for heterozygous alleles, and confidence scores showed the accuracy of the genotype call. Standard procedures and default analysis parameters for individual DNA samples were employed. An internal control run in parallel did not detect any DNA contamination. All procedures were performed using the same lots of reagents.

Quality Control

Genotype calling was performed using GDAS version 3.0.2 software as noted above. This software uses the Dynamic Models algorithm for genotype calling [19]. All markers with significant (P < 0.05 after Bonferroni correction) departure from Hardy-Weinberg equilibrium were excluded from final analyses. In addition, SNPs in the GeneChip® Human Mapping 50K Xba Array with <90% efficiency (>10% missing data) were excluded from the analysis. Missing data patterns did not deviate from random expectations.

Statistical Analysis

Clinical and demographic characteristics were compared among the three patient groups (extrapulmonary tuberculosis, pulmonary tuberculosis, and PPD+) using the Kruskal-Wallis test for continuous variables and the chi-squared and Fisher exact tests for categorical variables. Missing observations were individually excluded from all stages of analyses (for both the univariate and MDR analyses (no imputation was performed).

Two populations were assessed for genetic factors associated with extrapulmonary tuberculosis: 1) all individuals regardless of race; 2) only black participants. This stratified analysis was performed to minimize spurious associations due to population stratification [20]. There were not enough individuals to perform analyses stratified by other racial groups. The two populations were further divided into three subgroups for association analyses: a) any tuberculosis vs. PPD+; b) extrapulmonary tuberculosis vs. PPD+; and c) extrapulmonary vs. pulmonary tuberculosis.

Genetic analyses were performed in two stages (see below). Single locus chi-squared association tests (genotypic tests, with 2 degrees of freedom) (STATA version 9.0; College Station, TX) were performed on variants previously reported to be associated with tuberculosis susceptibility. MDR analyses were performed both on previously reported variants and on variants in candidate genes. MDR was performed to explore potential single-locus associations and gene-gene interactions between variants[21]. Details of the MDR algorithm and implementation are described in the Additional file 1 and Figure 1.

Figure 1
figure 1

Summary of the general steps to implement the MDR method, adapted from Ritchie and Motsinger 2005 [52]. In step one, the exhaustive list of n combinations are generated from the pool of all independent variables. In step two, for k = 1 to N, the combinations are represented in k-dimensional space, and the number of responders and non-responders are counted in each multifactor cell. In step three, the ratio of responders to non-responders is calculated within each cell. In step four, each multifactor cell in the k-dimensional space is labeled as high-likelihood/high-risk if the ratio of responsive individuals to non-responsive individuals exceeds a threshold and low-likelihood/low-risk if the threshold is not exceeded. In step five the training accuracy is calculated. This is then repeated for each multifactor combination. In step seven, the model with the best training accuracy is selected and evaluated in the test set. In step eight, the testing accuracy of the model is estimated. In step nine a permutation test is conducted to determine the statistical significance of the model(s). Steps 1 through 6 are repeated for each possible cross-validation interval. Bars represent hypothetical distributions of responders (left) and non-responders (right) with each multifactor combination. Dark-shaded cells represent high-likelihood genotype combinations while light-shaded cells represent low-likelihood genotype combinations.

The first stage of the genetic analysis examined variants reported to be associated with tuberculosis susceptibility in previous studies (Table 2). Because each variant represented a distinct statistical hypothesis (since each was being evaluated for replication), no correction for multiple testing was used.

In the second stage of genetic analysis, a nested candidate-gene study was compiled from the Affymetrix Xba genotyping arrays. Only SNPs in genes hypothesized to play a role in tuberculosis pathogenesis were included (see Table 4).

A Linux version of the MDR software was used for data analysis (compiled and benchmarked on a PC with a 600 MHz Pentium-III running Red Hat 2.2.5-15, written in C and compiled with the GNU C compiler). Presently, MDR software is being distributed in a JAVA version with a graphical user interface

For both analysis stages MDR was also performed with genetic factors plus CD4+ lymphocytes and body mass index (BMI). This allowed for exploration of both gene-gene and gene-clinical factor interactions.


  1. Bellamy R, Ruwende C, Corrah T, McAdam KPWJ, Whittle HC, Hill AVS: Variations in the NRAMP1 gene and susceptibility to tuberculosis in West Africans. N Engl J Med. 1998, 338: 640-644. 10.1056/NEJM199803053381002.

    CAS  Article  PubMed  Google Scholar 

  2. Bellamy R, Ruwende C, Corrah T, McAdam KP, Thursz M, Whittle HC, Hill AV: Tuberculosis and chronic hepatitis B virus infection in Africans and variation in the vitamin D receptor gene. J Infect Dis. 1999, 179: 721-724. 10.1086/314614.

    CAS  Article  PubMed  Google Scholar 

  3. Comstock GW: Tuberculosis in twins: a re-analysis of the Prophit survey. Am Rev Resp Dis. 1978, 117: 621-624.

    CAS  PubMed  Google Scholar 

  4. Berrington WR, Hawn TR: Mycobacterium tuberculosis, macrophages, and the innate immune response: does common variation matter?. Immunol Rev. 2007, 219: 167-186. 10.1111/j.1600-065X.2007.00545.x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. Gonzalez OY, Adams G, Teeter LD, Bui TT, Musser JM, Graviss EA: Extra-pulmonary manifestations in a large metropolitan area with a low incidence of tuberculosis. Int J Tuberc Lung Dis. 2003, 7: 1178-1185.

    CAS  PubMed  Google Scholar 

  6. Sterling TR, Dorman SE, Chaisson RE, Ding L, Hackman J, Moore K, Holland SM: HIV-seronegative adults with extrapulmonary tuberculosis have abnormal innate immune responses. Clin Infect Dis. 2001, 33: 976-982. 10.1086/322670.

    CAS  Article  PubMed  Google Scholar 

  7. Antas PRZ, Ding L, Hackman J, Reeves-Hammock L, Shintani AK, Schiffer J, Holland SM, Sterling TR: Decreased CD4+ lymphocytes and innate immune responses in adults with previous extrapulmonary tuberculosis. J Allerg Clin Immunol. 2006, 117: 916-923. 10.1016/j.jaci.2006.01.042.

    CAS  Article  Google Scholar 

  8. Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH: Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001, 69: 138-147. 10.1086/321276.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Williams SM, Ritchie MD, Phillips JA, Dawson E, Prince M, Dzhura E, Willis A, Semenya A, Summar M, White BC, Addy JH, Kpodonu J, Wong LJ, Felder RA, Jose PA, Moore JH: Multilocus analysis of hypertension: a hierarchical approach. Hum Hered. 2004, 57: 28-38. 10.1159/000077387.

    Article  PubMed  Google Scholar 

  10. Wilkinson RJ, Patel P, Llewelyn M, Hirsch CS, Pasvol G, Snounou G, Davidson RN, Toossi Z: Influence of polymorphism in the genes for the interleukin (IL)-1 receptor antagonist and IL-1beta on tuberculosis. J Exp Med. 1999, 189: 1863-1874. 10.1084/jem.189.12.1863.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. Wilkinson RJ, Llewelyn M, Toossi Z, Patel P, Pasvol G, Lalvani A, Wright D, Latif M, Davidson RN: Influence of vitamin D deficiency and vitamin D receptor polymorphisms on tuberculosis among Gujarati Asians in west London: a case-control study. Lancet. 2000, 355: 618-621. 10.1016/S0140-6736(99)02301-6.

    CAS  Article  PubMed  Google Scholar 

  12. Liu PT, Stenger S, Li H, Wenzel L, Tan BH, Krutzik SR, Ochoa MT, Schauber J, Wu K, Meinken C, Kamen DL, Wagner M, Bals R, Steinmeyer A, Zugel U, Gallo RL, Eisenberg D, Hewison M, Hollis BW, Adams JS, Bloom BR, Modlin RL: Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response. Science. 2006, 311: 1770-1773. 10.1126/science.1123933.

    CAS  Article  PubMed  Google Scholar 

  13. Leung CC, Lam TH, Chan WM, Yew WW, Ho KS, Leung G, Law WS, Tam CM, Chan CK, Chang KC: Lower risk of tuberculosis in obesity. Arch Intern Med. 2007, 167: 1297-1304. 10.1001/archinte.167.12.1297.

    Article  PubMed  Google Scholar 

  14. Kong Y, Cave MD, Zhang L, Foxman B, Marrs CF, Bates JH, Yang ZH: Association between Mycobacterium tuberculosis Beijing/W lineage strain infection and extrathoracic tuberculosis: Insights from epidemiologic and clinical characterization of the three principal genetic groups of M. tuberculosis clinical isolates. J Clin Microbiol. 2007, 45: 409-414. 10.1128/JCM.01459-06.

    CAS  Article  PubMed  Google Scholar 

  15. Kraft P: Curses--winner's and otherwise--in genetic epidemiology. Epidemiology. 2008, 19: 649-651. 10.1097/EDE.0b013e318181b865.

    Article  PubMed  Google Scholar 

  16. Weber JL, May PE: Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction. Am J Hum Genet. 1989, 44: 388-396.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Soborg C, Andersen AB, Madsen HO, Kok-Jensen A, Skinhoj P, Garred P: Natural resistance-associated macrophage protein 1 polymorphisms are associated with microscopy-positive tuberculosis. J Infect Dis. 2002, 186: 517-521. 10.1086/341775.

    CAS  Article  PubMed  Google Scholar 

  18. Yim JJ, Lee HW, Lee HS, Kim YW, Han SK, Shim YS, Holland SM: The association between microsatellite polymorphisms in intron II of the human Toll-like receptor 2 gene and tuberculosis among Koreans. Genes Immun. 2006, 7: 150-155. 10.1038/sj.gene.6364274.

    CAS  Article  PubMed  Google Scholar 

  19. Di X, Matsuzaki H, Webster TA, Hubbell E, Liu G, Dong S, Bartell D, Huang J, Chiles R, Yang G, Shen MM, Kulp D, Kennedy GC, Mei R, Jones KW, Cawley S: Dynamic model based algorithms for screening and genotyping over 100 K SNPs on oligonucleotide microarrays. Bioinformatics. 2005, 21: 1958-1963. 10.1093/bioinformatics/bti275.

    CAS  Article  PubMed  Google Scholar 

  20. Reich DE, Goldstein DB: Detecting association in a case-control study while correcting for population stratification. Genet Epidemiol. 2001, 20: 4-16. 10.1002/1098-2272(200101)20:1<4::AID-GEPI2>3.0.CO;2-T.

    CAS  Article  PubMed  Google Scholar 

  21. Motsinger AA, Ritchie MD: The effect of reduction in cross-validation intervals on the performance of multifactor dimensionality reduction. Genet Epidemiol. 2006, 30: 546-555. 10.1002/gepi.20166.

    Article  PubMed  Google Scholar 

  22. Awomoyi AA, Marchant A, Howson JM, McAdam KP, Blackwell JM, Newport MJ: Interleukin-10, polymorphism in SLC11A1 (formerly NRAMP1), and susceptibility to tuberculosis. J Infect Dis. 2002, 186: 1808-1814. 10.1086/345920.

    CAS  Article  PubMed  Google Scholar 

  23. Hoal EG, Lewis LA, Jamieson SE, Tanzer F, Rossouw M, Victor T, Hillerman R, Beyers N, Blackwell JM, van Helden PD: SLC11A1 (NRAMP1) but not SLC11A2 (NRAMP2) polymorphisms are associated with susceptibility to tuberculosis in a high-incidence community in South Africa. Int J Tuberc Lung Dis. 2004, 8: 1464-1471.

    CAS  PubMed  Google Scholar 

  24. Gao PS, Fujishima S, Mao XQ, Remus N, Kanda M, Enomoto T, Dake Y, Bottini N, Tabuchi M, Hasegawa N, Yamaguchi K, Tiemessen C, Hopkin JM, Shirakawa T, Kishi F: Genetic variants of NRAMP1 and active tuberculosis in Japanese populations. International Tuberculosis Genetics Team. Clin Genet. 2000, 58: 74-76. 10.1034/j.1399-0004.2000.580113.x.

    CAS  Article  PubMed  Google Scholar 

  25. Ma X, Dou S, Wright JA, Reich RA, Teeter LD, El Sahly HM, Awe RJ, Musser JM, Graviss EA: 5' dinucleotide repeat polymorphism of NRAMP1 and susceptibility to tuberculosis among Caucasian patients in Houston, Texas. Int J Tuberc Lung Dis. 2002, 6: 818-823.

    CAS  PubMed  Google Scholar 

  26. Li HT, Zhang TT, Zhou YQ, Huang QH, Huang J: SLC11A1 (formerly NRAMP1) gene polymorphisms and tuberculosis susceptibility: a meta-analysis. Int J Tuberc Lung Dis. 2006, 10: 3-12.

    CAS  PubMed  Google Scholar 

  27. Ryu S, Park YK, Bai GH, Kim SJ, Park SN, Kang S: 3'UTR polymorphisms in the NRAMP1 gene are associated with susceptibility to tuberculosis in Koreans. Int J Tuberc Lung Dis. 2000, 4: 577-580.

    CAS  PubMed  Google Scholar 

  28. Liu W, Cao WC, Zhang CY, Tian L, Wu XM, Habbema JD, Zhao QM, Zhang PH, Xin ZT, Li CZ, Yang H: VDR and NRAMP1 gene polymorphisms in susceptibility to pulmonary tuberculosis among the Chinese Han population: a case-control study. Int J Tuberc Lung Dis. 2004, 8: 428-434.

    CAS  PubMed  Google Scholar 

  29. Malik S, Abel L, Tooker H, Poon A, Simkin L, Girard M, Adams GJ, Starke JR, Smith KC, Graviss EA, Musser JM, Schurr E: Alleles of the NRAMP1 gene are risk factors for pediatric tuberculosis disease. Proc Natl Acad Sci USA. 2005, 102: 12183-12188. 10.1073/pnas.0503368102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Tosh K, Campbell SJ, Fielding K, Sillah J, Bah B, Gustafson P, Manneh K, Lisse I, Sirugo G, Bennett S, Aaby P, McAdam KP, Bah-Sow O, Lienhardt C, Kramnik I, Hill AV: Variants in the SP110 gene are associated with genetic susceptibility to tuberculosis in West Africa. Proc Natl Acad Sci USA. 2006, 103: 10364-10368. 10.1073/pnas.0603340103.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Lio D, Marino V, Serauto A, Gioia V, Scola L, Crivello A, Forte GI, Colonna-Romano G, Candore G, Caruso C: Genotype frequencies of the +874T-->A single nucleotide polymorphism in the first intron of the interferon-gamma gene in a sample of Sicilian patients affected by tuberculosis. Eur J Immunogenet. 2002, 29: 371-374. 10.1046/j.1365-2370.2002.00327.x.

    CAS  Article  PubMed  Google Scholar 

  32. Rossouw M, Nel HJ, Cooke GS, van Helden PD, Hoal EG: Association between tuberculosis and a polymorphic NFkappaB binding site in the interferon gamma gene. Lancet. 2003, 361: 1871-1872. 10.1016/S0140-6736(03)13491-5.

    CAS  Article  PubMed  Google Scholar 

  33. Lopez-Maderuelo D, Arnalich F, Serantes R, Gonzalez A, Codoceo R, Madero R, Vazquez JJ, Montiel C: Interferon-gamma and interleukin-10 gene polymorphisms in pulmonary tuberculosis. Am J Respir Crit Care Med. 2003, 167: 970-975. 10.1164/rccm.200205-438BC.

    Article  PubMed  Google Scholar 

  34. Henao MI, Montes C, Paris SC, Garcia LF: Cytokine gene polymorphisms in Colombian patients with different clinical presentations of tuberculosis. Tuberculosis (Edinb). 2006, 86: 11-19. 10.1016/

    CAS  Article  Google Scholar 

  35. Tso HW, Ip WK, Chong WP, Tam CM, Chiang AK, Lau YL: Association of interferon gamma and interleukin 10 genes with tuberculosis in Hong Kong Chinese. Genes Immun. 2005, 6: 358-363. 10.1038/sj.gene.6364189.

    CAS  Article  PubMed  Google Scholar 

  36. Cooke GS, Campbell SJ, Sillah J, Gustafson P, Bah B, Sirugo G, Bennett S, McAdam KP, Sow O, Lienhardt C, Hill AV: Polymorphism within the interferon-gamma/receptor complex is associated with pulmonary tuberculosis. Am J Respir Crit Care Med. 2006, 174: 339-343. 10.1164/rccm.200601-088OC.

    CAS  Article  PubMed  Google Scholar 

  37. Amirzargar AA, Rezaei N, Jabbari H, Danesh AA, Khosravi F, Hajabdolbaghi M, Yalda A, Nikbin B: Cytokine single nucleotide polymorphisms in Iranian patients with pulmonary tuberculosis. Eur Cytokine Netw. 2006, 17: 84-89.

    CAS  PubMed  Google Scholar 

  38. Gomez LM, Camargo JF, Castiblanco J, Ruiz-Narvaez EA, Cadena J, Anaya JM: Analysis of IL1B, TAP1, TAP2 and IKBL polymorphisms on susceptibility to tuberculosis. Tissue Antigens. 2006, 67: 290-296. 10.1111/j.1399-0039.2006.00566.x.

    CAS  Article  PubMed  Google Scholar 

  39. Awomoyi AA, Charurat M, Marchant A, Miller EN, Blackwell JM, McAdam KP, Newport MJ: Polymorphism in IL1B: IL1B-511 association with tuberculosis and decreased lipopolysaccharide-induced IL-1beta in IFN-gamma primed ex-vivo whole blood assay. J Endotoxin Res. 2005, 11: 281-286.

    CAS  PubMed  Google Scholar 

  40. Bellamy R, Ruwende C, Corrah T, McAdam KP, Whittle HC, Hill AV: Assessment of the interleukin 1 gene cluster and other candidate gene polymorphisms in host susceptibility to tuberculosis. Tuber Lung Dis. 1998, 79: 83-89. 10.1054/tuld.1998.0009.

    CAS  Article  PubMed  Google Scholar 

  41. Delgado JC, Baena A, Thim S, Goldfeld AE: Ethnic-specific genetic associations with pulmonary tuberculosis. J Infect Dis. 2002, 186: 1463-1468. 10.1086/344891.

    CAS  Article  PubMed  Google Scholar 

  42. Oral HB, Budak F, Uzaslan EK, Basturk B, Bekar A, Akalin H, Ege E, Ener B, Goral G: Interleukin-10 (IL-10) gene polymorphism as a potential host susceptibility factor in tuberculosis. Cytokine. 2006, 35: 143-147. 10.1016/j.cyto.2006.07.015.

    CAS  Article  PubMed  Google Scholar 

  43. Shin HD, Park BL, Kim YH, Cheong HS, Lee IH, Park SK: Common interleukin 10 polymorphism associated with decreased risk of tuberculosis. Exp Mol Med. 2005, 37: 128-132.

    CAS  Article  PubMed  Google Scholar 

  44. Soborg C, Madsen HO, Andersen AB, Lillebaek T, Kok-Jensen A, Garred P: Mannose-binding lectin polymorphisms in clinical tuberculosis. J Infect Dis. 2003, 188: 777-782. 10.1086/377183.

    CAS  Article  PubMed  Google Scholar 

  45. Bornman L, Campbell SJ, Fielding K, Bah B, Sillah J, Gustafson P, Manneh K, Lisse I, Allen A, Sirugo G, Sylla A, Aaby P, McAdam KP, Bah-Sow O, Bennett S, Lienhardt C, Hill AV: Vitamin D receptor polymorphisms and susceptibility to tuberculosis in West Africa: a case-control and family study. J Infect Dis. 2004, 190: 1631-1641. 10.1086/424462.

    CAS  Article  PubMed  Google Scholar 

  46. Flores-Villanueva PO, Ruiz-Morales JA, Song CH, Flores LM, Jo EK, Montano M, Barnes PF, Selman M, Granados J: A functional promoter polymorphism in monocyte chemoattractant protein-1 is associated with increased susceptibility to pulmonary tuberculosis. J Exp Med. 2005, 202: 1649-1658. 10.1084/jem.20050126.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. Hawn TR, Dunstan SJ, Thwaites GE, Simmons CP, Thuong NT, Lan NT, Quy HT, Chau TT, Hieu NT, Rodrigues S, Janer M, Zhao LP, Hien TT, Farrar JJ, Aderem A: A polymorphism in Toll-interleukin 1 receptor domain containing adaptor protein is associated with susceptibility to meningeal tuberculosis. J Infect Dis. 2006, 194: 1127-1134. 10.1086/507907.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Ogus AC, Yoldas B, Ozdemir T, Uguz A, Olcen S, Keser I, Coskun M, Cilli A, Yegin O: The Arg753GLn polymorphism of the human toll-like receptor 2 gene in tuberculosis disease. Eur Respir J. 2004, 23: 219-223. 10.1183/09031936.03.00061703.

    CAS  Article  PubMed  Google Scholar 

  49. Ben-Ali M, Barbouche MR, Bousnina S, Chabbou A, Dellagi K: Toll-like receptor 2 Arg677Trp polymorphism is associated with susceptibility to tuberculosis in Tunisian patients. Clin Diagn Lab Immunol. 2004, 11: 625-626.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Li CM, Campbell SJ, Kumararatne DS, Bellamy R, Ruwende C, McAdam KP, Hill AV, Lammas DA: Association of a polymorphism in the P2X7 gene with tuberculosis in a Gambian population. J Infect Dis. 2002, 186: 1458-1462. 10.1086/344351.

    CAS  Article  PubMed  Google Scholar 

  51. Fernando SL, Saunders BM, Sluyter R, Skarratt KK, Goldberg H, Marks GB, Wiley JS, Britton WJ: A polymorphism in the P2X7 gene increases susceptibility to extrapulmonary tuberculosis. Am J Respir Crit Care Med. 2007, 175: 360-366. 10.1164/rccm.200607-970OC.

    CAS  Article  PubMed  Google Scholar 

  52. Ritchie MD, Motsinger AA: Multifactor dimensionality reduction for detecting gene-gene and gene-environment interactions in pharmacogenomics studies. Pharmacogenomics. 2005, 6: 823-834. 10.2217/14622416.6.8.823.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

Download references


This work was presented in part at the 2007 Keystone Tuberculosis Symposium, March 20-25, 2007, Vancouver, British Columbia. Abstract 341.

The authors thank the Vanderbilt DNA Sequencing and Genotyping Facility of the Vanderbilt Institute for Integrative Genomics for assistance with performing microsatellite genotyping (Robert Woodhall). They also thank the Vanderbilt Center for Human Genetics Research for polymorphism analysis (DNA Resources Core; Ping Mayo and Cara B. Sutcliffe) and advice regarding study design (Marylyn D. Ritchie, Ph.D.). All array-based genotyping was performed in the Vanderbilt Microarray Shared Resource (John Mote). The Vanderbilt Microarray Shared Resource is supported by the Vanderbilt Ingram Cancer Center (P30 CA68485), the Vanderbilt Digestive Disease Center (P30 DK58404), and the Vanderbilt Vision Center (P30 EY08126).

This work was supported by the National Institutes of Allergy and Infectious Diseases (K23 AI01654, K24 AI065298, and the intramural program), the CAPES Foundation (Brazil), the Johns Hopkins Hospital General Clinical Research Center (M01-RR00052 from the National Center for Research Resources, National Institutes of Health), and the Vanderbilt University Department of Medicine.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Timothy R Sterling.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AAMR lead the analysis and the writing of the manuscript. PRZA contributed to the design of the study, identification of the genetic polymorphisms to be tested, identification and preparation of the samples to be tested for polymorphisms, and wrote the first draft of the manuscript. NOO contributed to the analysis and the writing of the manuscript. SL oversaw the performance of the genotyping and critically reviewed the manuscript. SMH: contributed to the design of the study and critically reviewed the manuscript. TRS conceived and designed the study, and made significant contributions to the writing of the manuscript. All authors read and approved the final manuscript.

Alison A Motsinger-Reif, Paulo RZ Antas contributed equally to this work.

Electronic supplementary material


Additional file 1: Appendix. Description of the Multifactor Dimensionality Reduction (MDR) Algorithm and its Implementation. (DOC 33 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Motsinger-Reif, A.A., Antas, P.R., Oki, N.O. et al. Polymorphisms in IL-1β, vitamin D receptor Fok1, and Toll-like receptor 2 are associated with extrapulmonary tuberculosis. BMC Med Genet 11, 37 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Tuberculosis
  • Pulmonary Tuberculosis
  • Multifactor Dimensionality Reduction
  • Extrapulmonary Tuberculosis
  • Multifactor Dimensionality Reduction Analysis