Association Between Cytokine Gene Polymorphisms and Human Susceptibility to Brucellosis

authors:

avatar Faezeh Kabiri 1 , avatar Ramazan Rajabnia ORCID 2 , avatar Abazar Pournajaf ORCID 3 , avatar Mehrdad Gholami 3 , avatar Meysam Hasannejad Bibalan 3 , avatar Aziz Kassani 4 , avatar Sajad Yaghoubi 5 , avatar Abdollah Ardebili ORCID 6 , 7 , *

Faculty of Dentistry, Qom University of Medical Sciences, Qom, IR Iran
Department of Microbiology, Faculty of Medicine, Babol University of Medical Sciences, Babol, IR Iran
Department of Microbiology, Faculty of Medicine, Iran University of Medical Sciences, Tehran, IR Iran
Department of Community Medicine, School of Medicine, Dezful University of Medical Sciences, Dezful, IR Iran
Division of Microbiology, Department of Pathobiology, School of Public Health, Tehran University of Medical Sciences, Tehran, IR Iran
Laboratory Sciences Research Center, Golestan University of Medical Sciences, Gorgan, IR Iran
Department of Microbiology, Faculty of Medicine, Golestan University of Medical Sciences, Gorgan, IR Iran

how to cite: Kabiri F, Rajabnia R, Pournajaf A, Gholami M, Hasannejad Bibalan M, et al. Association Between Cytokine Gene Polymorphisms and Human Susceptibility to Brucellosis. Arch Pediatr Infect Dis. 2017;5(3):e57274. https://doi.org/10.5812/pedinfect.57274.

Abstract

Background:

Allelic single nucleotide polymorphisms (SNPs) in cytokine-encoding genes can affect the degree of cytokine production, and may be related to the tendency to infectious illnesses as well as various clinical consequences.

Objectives:

The aim of this work was to evaluate the possible role of SNPs in the regions of the IL-10 (-592), IL-15 (-367), IL-18 (-656), IL-12 (+1188), IFN-γ (+874), TNF-α (-308), and TNF-β (+252) genes in susceptibility or resistance to brucellosis and its crucial complications.

Methods:

In a period of one year, 125 patients with acute brucellosis referring to 3 large public teaching hospitals were enrolled in this study. We studied the SNPs of IL-10, IL-15, IL-18, IL-12, IFN-γ, and TNF-α/β genes using the allele specific polymerase chain reaction (AS-PCR) with sequence-specific primers.

Results:

Frequency of GG genotype in the TNF-α and TNF-β-encoding genes increased significantly by 52% and 31.2% in patient and control groups, respectively. For IFN-γ, TA genotype was found highly enhanced in patients (60%), while the frequency of AA and TT genotypes were higher in controls (23.2% and 26.4%, respectively). The AA and CC SNPs in IL-12 were dominant in both patient (78.4%) and control (14.4%) groups. In the patient group, the GG and TT genotypes had a higher frequency for genes encoding IL-15 (33.6%) and IL-18 (89.6%).

Conclusions:

Based on the present study, some SNPs within the several cytokine genes, including TNF-α/β (-308/+252), IFN-γ (+874), IL-15 (-367), IL-18 (-656), and IL-12 (+1188) are related to the susceptibility or resistance to brucellosis. In order to approve the biological consequence of our results, additional investigations should be carried out in larger population groups.

1. Background

The Brucella species are small facultative intracellular gram‐negative aerobic bacteria with ability to infect both humans and animals (1). Brucellosis is usually a common zoonotic infection that causes worldwide disease afflicting more than half a million individuals annually (2). B. melitensis and B. abortus are the most common agents of human brucellosis in many countries (3). This illness is endemic in many Asian areas, such as Iran, Turkey, the Arabian Peninsula and Mediterranean countries, Indian subcontinent and Central and South America (4, 5). The actual prevalence of brucellosis in the world is indefinite because of poor reporting systems and lack of access to reliable diagnostic tests in many developing countries (4). More than 500,000 cases with brucellosis are reported to the world health organization (WHO) globally every year, most of whom are from developing countries (6). According to the center for disease control (CDC), the incidence rate of brucellosis in Iran was 39 per 100,000 and 30 per 100,000 populations in 2005 and 2007, respectively. Also, the incidence of brucellosis has increased to 130 per 100,000 populations in western Iran in recent years (6).

The Brucella spp. enter macrophage-monocyte lineage cells, survive, multiply within them, and spread in mononuclear phagocytes (7). Consequently, acquired cell-dependent immune, determined by the T-helper1 (Th1) lymphocyte activation and subsequent activation of macrophages, plays a crucial role in the protection against this infectious disease (7). Production of cytokine index can be considered as T-helper- immune cell, interleukin-2 (IL-2), interferon gamma (IFN-γ), or T-helper-cell-class-2 (Th-2) responses inducing humoral type of immunity (IL-4/5/6/10) and Th-3 type determined by TGF-β (8).

Apart from environmental elements and pathogen strain differences, host genetic determinants are also major causes of susceptibility to or consequence of infectious illnesses. Expression and release of cytokines are dependent, at least in part, on genetic variation within the promoter region or other regulatory sequences of the cytokine genes (8, 9). It has been demonstrated that some single nucleotide polymorphisms (SNPs) in the genes encoding certain cytokines could not only rise susceptibility to some infectious diseases, but also change the course and prediction of the disease (9). For example, it is possible that TNF-α gene SNPs have an impact on the rate of TNF-α synthesis, which in turn might have an effect on inflammatory responses. The substitution of Guanine (G) or Adenine (A) at position +252 of TNF-β and -308 (A/G) genotype of TNF-α is associated with the susceptibility to autoimmune disease and microbial infections (10).

2. Objectives

There are limited data about the effects of polymorphisms on the outcome of brucellosis in the infected patients in Iran. In the present study, we aimed to determine the possible association between cytokine gene SNPs and either susceptibility to or development of focal complications of brucellosis, and also demonstrate that the systemic immune response to Brucella antigens might be enhanced by mutation in the cytokine encoded genes.

3. Methods

3.1. Study Population

This cross-sectional study was performed during a one-year period of time in 2015 - 2016. The study population comprised two patient and control groups. 125 patients with confirmed acute brucellosis were admitted to three teaching therapeutic centers, including Imam Khomeini, Shariati, and Rasoul-e-Akram hospitals. The mean age of the study population was 49 ± 1.5 years, with a range of 17 to 76 years. The isolates were collected from the patients in various age groups: 12 - 22 years (n = 9), 23 - 32 years (n = 21), 33 - 43 years (n = 39), 44 - 54 years (n = 32), 55 - 65 years (n = 17), and 66 - 76 years (n = 7). 72 (57.6%) patients were male and 53 (42.4%) were female. 5 - 10 mL whole blood specimens were obtained from each participant and kept in EDTA‐containing tubes for DNA extraction. Brucellosis was identified based on the clinical manifestation (e.g. Night sweats, fever, debility, arthralgia, weight loss, lymphadenopathy, splenomegaly, malaise, and myalgia), and/or positive blood cultures and serological methods. A single high titer (≥ 1:160) of standard agglutination test (SAT) was considered as a prognostic standard for positive serological results that was approved by a high titer (≥ 1:160) of 2-mercaptoethanol test (2ME) at the time of infection.

The control group comprised healthy blood donors with no history of brucellosis or genetic disorders. They were matched with the patients for sex, age, and geographic area; they had a similar history and were at the same risk of exposure to brucellosis. The study was approved by the institutional review board and the research ethics committee of Islamic Azad University, Qom branch.

3.2. Determination of Cytokine Gene Polymorphism

SNPs for TNF-α/β, IFN-γ, and IL-10/12/15/18 were analyzed by allele specific polymerase chain reaction (AS-PCR). Cellular DNA was obtained from EDTA-treated peripheral venous blood by salting out method (11). The concentration and the quality of the extracted DNA were assessed using a Nanodrop spectrophotometer (ND-1000; Thermo Scientific, Wilmington, DE, USA). Primer sequences used for gene amplification are shown in Table 1.

Table 1.

The Nucleotide Sequences of the Primers Used in This Study

Target GenesPrimer Sequence (5'→3')References
TNF-αF, TCTCGGTTTCTTCTCCATCG(12)
R, ATAGGTTTTGAGGGGCATGG
TNF-βF, AACCTGCTGCTCACCTTGTT(13)
R, CAGTGCAAAGGCTCCAAAGA
IFN-γF, TGAACGCTACACACTGCATC(14)
R, CACTCGGATGAACTCATTGA
IL-10F, CCTAGGTCACAGTGACGTGG(15)
R, GGTGAGCACTACCTGACTAGC
IL-12F, TTTGGAGGAAAAGTGGAAGA(16)
R, AACATTCCATACATCCTGGC
IL-15F, TCT TCA ATA CTT AAG GAT TTAC(17)
R, CAA AAG AGT GGG ATA AGT GA
IL-18F, AGGTCAGTCTTTGCTATCATTCCAGG(18)
R, CTGCAACAGAAAGTAAGCTTGCGGAGAGG

AS-PCR was conducted for the amplification of target genes; a volume of 1.0 µL of total extracted DNA was added to a whole volume of 25 µL PCR reaction mixture including 2.0 µL of 10×PCR buffer, 1.3 µL MgCl2 (50 mM), 0.5 µL dNTPs (10 mM), 1.0 µL of each primer, 0.5 µL of Taq DNA polymerase (5 U/µL) (Amplicon Co., Denmark), and 17.7 µL double-distilled water (ddH2O). PCR was performed in a thermal gradient cycler (Eppendorf Co., Germany) according to the following procedure: initial denaturation at 95°C for 5 minutes, then 33 cycles with denaturation at 94°C for 30 seconds, annealing at 61°C for 55 seconds, extension at 72°C for 32 seconds, and final extension at 72°C for 40 seconds. The amplified products were visualized using a UV transilluminator, following electrophoresis on 1% agarose gel stained with Gel RedTM.

3.3. Statistical Analysis

Statistical analysis was performed using the Statistical Package for the Social Sciences (SPSS) software, ver. 23 (IBM). Chi-square test and logistic regression model were applied for analyzing data. Also, odds ratio (OR) and 95% confidence interval (CI 95%) were calculated. In addition, odds ratios were adjusted for age and sex variables in logistic regression model. Analyses were interpreted as statistically significant at P < 0.05.

4. Results

As shown in Table 2, all SNPs were calculated in the Hardy Weinberg equilibrium in the patient and control groups using chi-square test. Polymorphism in the TNF-α gene at position -308 showed that the frequency of GG genotype increased significantly in patients compared to the control group (52% vs. 27%; P = 0.016, OR = 0.789, 95% CI = 0.345 - 1.890). The patient group in comparison with the control group showed a higher frequency of the AA genotype at position +252 in TNF-β (34.4% vs. 27.2%; P = 0.038, OR = 1.455, 95% CI=0.091 - 1.734) and GG genotype of TNF-β (+252) (27.2% vs. 31.2%; P = 0.047, OR = 1.234, 95% CI = 0.764 - 2.986).

Table 2.

Distribution of the Cytokine Gene Polymorphisms in Patients with Brucellosis

Cytokine GeneGenotypeNo. (%) of PatientsNo. (%) of Healthy ControlOdds Ratio95% CIP Value
TNF-α (-308)GG65 (52)34 (27.2)0.7890.345 - 1.8900.016a
AG35 (28)36 (28.8)1.70.678 - 5.8970.449
AA25 (20)55 (44)1.90.238 - 29.3450.659
TNF-β (+252)AG57 (45.6)43 (34.4)0.8710.677 - 1.9840.578
GG34 (27.2)39 (31.2)1.2340.764 - 2.9860.047a
AA34 (27.2)43 (34.4)1.4550.091 - 1.7340.038a
IFN-γ (+874)TA75 (60)63 (50.4)3.210.452 - 3.1670.012a
AA25 (20)29 (23.2)0.9870.178 - 1.7430.073
TT25 (20)33 (26.4)0.7840.356 - 1.9370.023a
IL-10 (-592)AA69 (55.2)80 (64)2.3450.823 - 8.6780.419
AC12 (9.6)29 (23.2)0.7941.456 - 25.1210.549
CC44 (35.2)16 (12.8)1.9231.876 - 42.2310.713
IL-12 (+1188)AA98 (78.4)89 (71.2)0.7830.178 - 0.9790.016a
AC16 (12.8)18 (14.4)2.7651.345 - 12.6540.457
CC11 (8.8)18 (14.4)1.3511.981 - 13.4510.017a
IL-15 (-367)AA43 (34.4)39 (31.2)1.2310.567 - 1.9860.246
GG42 (33.6)37 (29.6)1.2780.987 - 2.6420.039a
GA40(32)49(39.2)0.7690.327 - 1.5320.065
IL-18 (-656)TT112 (89.6)104 (83.2)0.9450.231 - 1.7630.017a
TG10 (8)14 (11.2)1.9340.654 - 9.8080.091
GG3 (2.4)7 (5.6)2.2350.712 - 26.4200.044a

Polymorphism of IFN-γ gene (at position +874) indicated an increased prevalence of the TA genotype in the patient group (60%, P = 0.012, OR = 3.21, 95% CI = 0.452-3.167). However, the frequency of AA and TT genotypes was higher in the control than the patient group. SNPs in the gene encoding IL-10 (at position -592) showed no significant difference between patient and control groups for all three genotypes studied.

The AA genotype at position +1188 of gene encoding IL-12 was more common in the patients than control group, while the CC genotype was more prevalent in the control group. The most common genotypes in the gene encoding IL-15 (at position -367) were GG (in the patient group) and GA (in the control group). At position -656 in the gene encoding IL-18, TT genotype was more frequent in patients, while both GG and TG genotypes were more common in controls. According to our data, other cytokine genotypes studied showed no significant between-group differences.

5. Discussion

Research to detect SNPs within the genes encoding cytokines, and determining and comparing the cytokine genotype(s) in patients infected with Brocella spp. or healthy individuals can improve our knowledge regarding brucellosis immunology and facilitate the development of new therapeutic and vaccination strategies. Recently, several reports have shown that SNPs in cytokine genes could be an important factor in resistance and/or susceptibility to brucellosis (9-11). To better understanding of brucellosis immunology, we investigated the association of gene polymorphisms of cytokines with this infectious disease.

It has been reported a noteworthy relationship between AA genotype in the promoter region of IFN-γ at the position +874 and various human diseases, such as brucellosis, hepatitis B, and tuberculosis (19). According to study by Bravo et al., the IFN-γ AA genotype was substantially higher in brucellosis patients than healthy people (20). In another study conducted by Budak et al., no association was demonstrated between IFN-γ +874 polymorphism and the risk of acquiring brucellosis in humans (21). In agreement with Budak et al., our results indicated that IFN-γ polymorphisms (at positions +874) were notably prevalent in patients than the control group.

There are several genetic polymorphic sites which are known to be associated with the production of TNF, including TNF-α (-308 G/A, -238 G/A) and TNF-β (+252 A/G) (22, 23). In some previous studies, TNF-α polymorphisms at positions -308 and -238 have been associated with susceptibility to brucellosis (24). Rasouli et al., showed that the frequency of AA genotypes of TNF-β and distribution of A allele were significantly higher in patients than controls (25). It appears that both SNPs reported by Rasouli and coworkers are associated with low production of TNF-α, suggesting the increased susceptibility to infection. Reza et al., demonstrated that TNF-α -308 (A/A) genotype had a higher frequency in the population of patients in comparison with controls (26). Their results indicated that, although the frequency of allele in the two groups was not statistically remarkable, TNF-α polymorphism at nucleotide -308 (A/A) could be involved in the susceptibility to brucellosis (26). In addition, Caballero et al., found that there is no link between the TNF-α -308 (A/A) genotypes and brucellosis (27). Our results indicated that the GG genotypes of TNF-α and -β were significantly higher in patients than healthy people. It seems that individuals who inherit A allele as homozygous (AA) possibly produce lower levels of TNF-α, which could cause the lack of proper immune response at the early stages of Brucella infection, leading the incidence of a full-blown disease. Also, in persons with AG and GG haplotypes, due to more TNF-α production, the bacteria are controlled and the disease cannot develop.

IL-12 is a heterodimeric cytokine that has a key role in the promotion of type 1 immune response. Moreover, it has a protective effect against the Brucella infection (28). In a study conducted by Kamali-Sarvestani et al., IL-12 A allele, which is associated with higher production of IL-12, was significantly more frequent in the controls than the patients. Moreover, AA genotype was significantly more frequent in the controls than the patients (24). Based on these results, it can be concluded that individuals who acquire the AA genotype may produce higher levels of IL-12 which can lead to the initiation of CMI response to brucellosis. On the contrary, the results of the current study indicated that the AA genotypes of IL-12 (+1188) polymorphisms were significantly more frequent in patients with brucellosis compared to the healthy group.

IL-10 is a cytokine with anti-inflammatory action which can strongly prevent the release of cytokines, such as IL-2 and IFN-γ (27). Additionally, IL-10 has an obvious comitogenic effect on proliferation of B and T cells and enhances B cell antibody development (29). Karaoglan et al., showed that the rates of CT and CC genotypes of IL-10 (-819) were more likely in patient and control groups, respectively (30). In the study by Bravo et al., the association between IL-10 gene polymorphisms and susceptibility to brucellosis was not observed (20). Similar to the latter study, our results indicated no significant differences in the genotype distribution of IL-10 (-592) polymorphisms between the patients and controls.

IL-15 in synergy with IL-12 enhances IFN-γ production via NK and T cells and stimulates Th1 responses against the intracellular pathogens (31). In our study, three genotypes of IL-15 (-367), AA/GG and GA had nearly the same distribution in both groups of patient and control, while the GG genotype of IL-15 (-592) was significantly higher in the patient group compared to the control group. In the study by Kalani et al., no significant difference was observed between the frequency of alleles and genotype polymorphisms of IL-15 (-367) in the controls and patients, suggesting that there is not a significant relationship between IL-15 (-367) gene polymorphisms and susceptibility and/or resistance to brucellosis (22).

Some studies showed that SNPs of the IL-18 gene can lead to the increased rates of cytokine expression (23). Similar to findings by the Rasouli et al., our study showed that TT genotypes of IL-18 (-656) have a significantly increased frequency in the patients compared to the controls (32). These findings suggest that TT genotype of IL-18 (-656) could be involved as a potential risk factor in the susceptibility to brucellosis.

Our results also indicated that some genotypes of cytokine polymorphism, including IFN-γ (+874) TA, TNF-α (-308) GG, TNF-β (+252) GG, IL-12 (+1188) AA, IL-15 (-367) GG, and IL-18 (-656) TT are criteria for susceptibility to brucellosis. On the other hand, cytokine genotypes, such as IFN-γ (+874) AA and TT, IL-12 (+1188) CC, TNF-β (+252) AA, and IL-18 (-656) GG may play a protective role against this infection. It is important to evaluate the dependency of these SNPs in cytokine genes on expanding brucellosis in the Iranian population.

Acknowledgements

References

  • 1.

    Guerra H. The brucellae and their success as pathogens. Crit Rev Microbiol. 2007;33(4):325-31. [PubMed ID: 18033597]. https://doi.org/10.1080/10408410701647644.

  • 2.

    Dean AS, Crump L, Greter H, Schelling E, Zinsstag J. Global burden of human brucellosis: a systematic review of disease frequency. PLoS Negl Trop Dis. 2012;6(10). e1865. [PubMed ID: 23145195]. https://doi.org/10.1371/journal.pntd.0001865.

  • 3.

    Diaz Aparicio E. Epidemiology of brucellosis in domestic animals caused by Brucella melitensis, Brucella suis and Brucella abortus. Rev Sci Tech. 2013;32(1):43-51. 53-60. [PubMed ID: 23837364]. https://doi.org/10.20506/rst.32.1.2187.

  • 4.

    Pappas G, Memish ZA. Brucellosis in the middle East: a persistent medical, socioeconomic and political issue. J Chemother. 2007;19(3):243-8. [PubMed ID: 17594917]. https://doi.org/10.1179/joc.2007.19.3.243.

  • 5.

    Gwida M, Al Dahouk S, Melzer F, Rosler U, Neubauer H, Tomaso H. Brucellosis - regionally emerging zoonotic disease? Croat Med J. 2010;51(4):289-95. [PubMed ID: 20718081]. https://doi.org/10.3325/cmj.2010.51.289.

  • 6.

    Eini P, Keramat F, Hasanzadehhoseinabadi M. Epidemiologic, clinical and laboratory findings of patients with brucellosis in Hamadan, west of Iran. J Res Health Sci. 2012;12(2):105-8. [PubMed ID: 23241521].

  • 7.

    von Bargen K, Gorvel JP, Salcedo SP. Internal affairs: investigating the Brucella intracellular lifestyle. FEMS Microbiol Rev. 2012;36(3):533-62. [PubMed ID: 22373010]. https://doi.org/10.1111/j.1574-6976.2012.00334.x.

  • 8.

    Rodriguez-Zapata M, Matias MJ, Prieto A, Jonde MA, Monserrat J, Sanchez L, et al. Human brucellosis is characterized by an intense Th1 profile associated with a defective monocyte function. Infect Immun. 2010;78(7):3272-9. [PubMed ID: 20404074]. https://doi.org/10.1128/IAI.01385-09.

  • 9.

    Davoudi S, Amirzargar AA, Hajiabdolbaghi M, Rasoolinejad M, Soodbakhsh A, Jafari S, et al. Th-1 cytokines gene polymorphism in human brucellosis. Int J Immunogenet. 2006;33(5):355-9. [PubMed ID: 16984280]. https://doi.org/10.1111/j.1744-313X.2006.00626.x.

  • 10.

    Kazemi S, Saidijam M, Hashemi SH, Karami M, Vaisi-Raygani A, Alikhani MY. Analysis of IL-10 and IL-6 Gene Polymorphisms and Their Serum Levels in Patients with Brucellosis: A Case Control Study. Immunol Invest. 2016;45(2):107-15. [PubMed ID: 26849072]. https://doi.org/10.3109/08820139.2015.1096285.

  • 11.

    Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215. [PubMed ID: 3344216]. https://doi.org/10.1093/nar/16.3.1215.

  • 12.

    Lakhanpal M, Singh LC, Rahman T, Sharma J, Singh MM, Kataki AC, et al. Study of single nucleotide polymorphisms of tumour necrosis factors and HSP genes in nasopharyngeal carcinoma in North East India. Tumour Biol. 2016;37(1):271-81. [PubMed ID: 26198046]. https://doi.org/10.1007/s13277-015-3767-6.

  • 13.

    Gonen-Korkmaz C, Sevin G, Gokce G, Arun MZ, Yetik-Anacak G, Yildirim G, et al. Different responses of RT-PCR amplifications after tumor necrosis alpha (TNFa) induction at nuclear factor kappa B (NFkB) gene silenced LNCaP cells. FASEB J. 2013;27(1 Supplement):1105.7.

  • 14.

    Furini AA, Capobianco MP, Storti-Melo LM, Cunha MG, Cassiano GC, Machado RL. Cytokine gene polymorphisms are not associated with anti-PvDBP, anti-PvAMA-1 or anti-PvMSP-119 IgG antibody levels in a malaria-endemic area of the Brazilian Amazon. Malar J. 2016;15(1):374. [PubMed ID: 27435973]. https://doi.org/10.1186/s12936-016-1414-3.

  • 15.

    Ding Q, Shi Y, Fan B, Fan Z, Ding L, Li F, et al. The interleukin-10 promoter polymorphism rs1800872 (-592C>A), contributes to cancer susceptibility: meta-analysis of 16,785 cases and 19,713 controls. PLoS One. 2013;8(2). e57246. [PubMed ID: 23460834]. https://doi.org/10.1371/journal.pone.0057246.

  • 16.

    Sortica VA, Cunha MG, Ohnishi MD, Souza JM, Ribeiro-Dos-Santos AK, Santos NP, et al. IL1B, IL4R, IL12RB1 and TNF gene polymorphisms are associated with Plasmodium vivax malaria in Brazil. Malar J. 2012;11:409. [PubMed ID: 23217179]. https://doi.org/10.1186/1475-2875-11-409.

  • 17.

    Pavkova Goldbergova M, Nemec P, Lipkova J, Jarkovsky J, Gatterova J, Ambrozkova D, et al. Relation of IL-6, IL-13 and IL-15 gene polymorphisms to the rheumatoid factors, anti-CCP and other measures of rheumatoid arthritis activity. Int J Immunogenet. 2014;41(1):34-40. [PubMed ID: 23773307]. https://doi.org/10.1111/iji.12065.

  • 18.

    Wang M, Zhu XY, Wang L, Lin Y. The -607C/A polymorphisms in interleukin-18 gene promoter contributes to cancer risk: evidence from a meta-analysis of 22 case-control studies. PLoS One. 2013;8(10). e76915. [PubMed ID: 24130810]. https://doi.org/10.1371/journal.pone.0076915.

  • 19.

    Ben-Ari Z, Mor E, Papo O, Kfir B, Sulkes J, Tambur AR, et al. Cytokine gene polymorphisms in patients infected with hepatitis B virus. Am J Gastroenterol. 2003;98(1):144-50. [PubMed ID: 12526950]. https://doi.org/10.1111/j.1572-0241.2003.07179.x.

  • 20.

    Bravo MJ, de Dios Colmenero J, Alonso A, Caballero A. Polymorphisms of the interferon gamma and interleukin 10 genes in human brucellosis. Eur J Immunogenet. 2003;30(6):433-5. [PubMed ID: 14675398]. https://doi.org/10.1111/j.1365-2370.2003.00419.x.

  • 21.

    Budak F, Goral G, Heper Y, Yilmaz E, Aymak F, Basturk B, et al. IL-10 and IL-6 gene polymorphisms as potential host susceptibility factors in Brucellosis. Cytokine. 2007;38(1):32-6. [PubMed ID: 17544674]. https://doi.org/10.1016/j.cyto.2007.04.008.

  • 22.

    Kalani M, Rasouli M, Moravej A, Kiany S, Rahimi HR. Association of interleukin-15 single nucleotide polymorphisms with resistance to brucellosis among Iranian patients. Tissue Antigens. 2011;78(5):352-8. [PubMed ID: 21988722]. https://doi.org/10.1111/j.1399-0039.2011.01775.x.

  • 23.

    Folwaczny M, Glas J, Torok HP, Tonenchi L, Paschos E, Bauer B, et al. Polymorphisms of the interleukin-18 gene in periodontitis patients. J Clin Periodontol. 2005;32(5):530-4. [PubMed ID: 15842270]. https://doi.org/10.1111/j.1600-051X.2005.00711.x.

  • 24.

    Kamali-Sarvestani E, Rasouli M, Mortazavi H, Gharesi-Fard B. Cytokine gene polymorphisms and susceptibility to cutaneous leishmaniasis in Iranian patients. Cytokine. 2006;35(3-4):159-65. [PubMed ID: 16950634]. https://doi.org/10.1016/j.cyto.2006.07.016.

  • 25.

    Rasouli M, Kiany S, Moravej A, Kalani M. Interleukin-12 and tumor necrosis factor-β gene polymorphisms as genetic susceptibility factors for brucellosis in Iranian patients. Iran Red Crescent Med J. 2010;2010(3):266-71.

  • 26.

    Reza JS, Alireza R, Mostafa A, Alireza K, Mehrdad H. Association of tumor necrosis factor-alpha-308 (G-->A) polymorphism and susceptibility to brucellosis. J Microbiol Immunol Infect. 2009;42(1):22-6. [PubMed ID: 19424555].

  • 27.

    Caballero A, Bravo MJ, Nieto A, Colmenero JD, Alonso A, Martin J. TNFA promoter polymorphism and susceptibility to brucellosis. Clin Exp Immunol. 2000;121(3):480-3. [PubMed ID: 10971514]. https://doi.org/10.1046/j.1365-2249.2000.01331.x.

  • 28.

    Sathiyaseelan J, Goenka R, Parent M, Benson RM, Murphy EA, Fernandes DM, et al. Treatment of Brucella-susceptible mice with IL-12 increases primary and secondary immunity. Cell Immunol. 2006;243(1):1-9. [PubMed ID: 17184756]. https://doi.org/10.1016/j.cellimm.2006.10.003.

  • 29.

    Jubier-Maurin V, Boigegrain RA, Cloeckaert A, Gross A, Alvarez-Martinez MT, Terraza A, et al. Major outer membrane protein Omp25 of Brucella suis is involved in inhibition of tumor necrosis factor alpha production during infection of human macrophages. Infect Immun. 2001;69(8):4823-30. [PubMed ID: 11447156]. https://doi.org/10.1128/IAI.69.8.4823-4830.2001.

  • 30.

    Karaoglan I, Pehlivan S, Namiduru M, Pehlivan M, Kilincarslan C, Balkan Y, et al. TNF-alpha, TGF-beta, IL-10, IL-6 and IFN-gamma gene polymorphisms as risk factors for brucellosis. New Microbiol. 2009;32(2):173-8. [PubMed ID: 19579695].

  • 31.

    Yoshikai Y, Nishimura H. The role of interleukin 15 in mounting an immune response against microbial infections. Microbes Infect. 2000;2(4):381-9. [PubMed ID: 10817640]. https://doi.org/10.1016/S1286-4579(00)00329-4.

  • 32.

    Rasouli M, Kalani M, Moravej A, Kiany S. Interleukin-18 single nucleotide polymorphisms contribute to the susceptibility to brucellosis in Iranian patients. Cytokine. 2011;54(3):272-6. [PubMed ID: 21393015]. https://doi.org/10.1016/j.cyto.2011.02.011.