- Research article
- Open Access
Association between the IL-10 and IL-6 polymorphisms and brucellosis susceptibility: a meta-analysis
BMC Medical Genetics volume 21, Article number: 63 (2020)
Brucellosis is a quite normal zoonotic infection, which is caused by immediate contact with animals infected with Brucella or its products. IL-10 (− 1082 G/A, − 819 C/T, − 592C/A) and IL-6 -174 G/C polymorphisms have a great relationship with IL-10 and IL-6 production, which brings about Brucellosis pathogenesis and development. So far, the results of published literatures were controversial. Now, we perform a meta-analysis in different ethnic populations to get a more precise estimate of above polymorphisms with Brucellosis susceptibility.
Both OR and corresponding 95%CI were enrolled to make an assessment of the association strength through extracting genotyping frequency of cases and controls. The χ2-test based Q-statistic and I2 statistics were applied. If there was no evident heterogeneity, the fixed-effects model would be applied. If not, the random-effects model would be used.
The significant associations were only found in Asian population of − 819 loci under three genetic models as follows: (Allele model: OR = 0.60, 95%CI = 0.44–0.82, P = 0.001), (homozygote comparison: OR = 0.24, 95%CI = 0.09–0.62, P = 0.003), (recessive genetic model: OR = 0.22, 95%CI = 0.05–0.91, P = 0.036).
In conclusion, IL-10 − 819 loci polymorphism contributes no risk to Caucasian population but may be associated with decreased risk in Asian population. And IL-10 -1082 G/A, 592 loci and IL-6 -174 G/C polymorphism are not associated with Brucellosis risk.
Brucellosis is a quite normal zoonotic infection, which is caused by immediate contact with animals infected with Brucella or its products. Although the number of Brucellosis patients is relatively small, it remains a severe problem and it is very popular locally in most areas of Asia and Africa [1, 2]. Some patients show onset fever, some show fatigue or sweating. Multiple clinical manifestations can be shown by Brucellosis. What’s worse, nontypical clinical presentation brings difficulties to diagnosis. The exact pathogenesis of Brucellosis remains unknown. It has been reported that cell-mediated immunity is considered to play a crucial role in immunity response to Brucellosis infection . CD4+ and CD8+ T lymphocytes are considered as playing a key role in cellular immunity as they can release IFN-λ and activate the functions in macrophages [4, 5]. Additionally, Interleukin-10 (IL-10) is a crucial cytokine contributing to resist inflammation, which makes various biological effects on multiple types of cell. After the infection of Brucella pathogen, IL-10 can lead to the production drawdown of IFN-λ and inhibition of macrophages function .
IL-6 is another cytokine, not only involving in the process of inflammation response but also regulating multiple biological processes such as metabolism, regeneration and nervous system . It has been demonstrated that polymorphism can influence the expression of cytokine and play a crucial role in infectious diseases [8, 9].
There are quite a few literatures demonstrating that hereditary factors play crucial parts in the development of Brucellosis [10,11,12,13,14]. Apart from IL-6 -174 G/C polymorphism, a large number of documents have focused on IL-10 promoter polymorphisms, including position − 1082(G → A) (rs1800896), − 819(C → T) (rs1800871) and − 592(C → A) (rs1800872). However, inconsistent results were obtained. In the present study, we carry out a meta-analysis to obtain a more precise estimate of the above polymorphisms. In addition, the present study is a multi-ethnic study because ethnic populations may contribute tremendous influence to genetic polymorphisms of IL-10 and IL-6.
The present meta-analysis was performed on account of predefined protocol named Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group . Two databases consisting of PubMed and Embase were put into use by searching keywords as follows: (“IL” or “Interleukin”) and “Brucellosis” and “polymorphism”. The expiration date was the end of February, 2020. We reviewed titles and abstracts of all citations and retrieved studies. We did not limit anything of literatures such as its language, regional culture, publication year and sample size. After these studies were received, their corresponding references were also investigated for other possible studies.
Inclusion and exclusion criteria
The inclusion criteria consisted of three items: (a) a case-control research or short communication; (b) assessment of relationship between IL polymorphisms and Brucellosis susceptibility; (c) providing enough data to judge the final result or the present data can predict the final result. Accordingly, other types of research were excluded such as case report, review article.
The whole information and data were independently extracted by first author (Yueyuan Wu, Shuzhou Yin and Xiaochun Jin) and the final information was reviewed by final referees (Xu Chen and Youtao Zhang). The present study is a multi-ethnic study. Different ethnic descents were categorized as African, Asian, and Caucasian. To get a accurate result, three authors (Shuzhou Yin, Yueyuan Wu and Xiaochun Jin) were asked to check all data and information. If they could not reach an agreement, they would check the above data and information again and have a meeting trying to come to an agreement. If the controversial results still existed, the final referees (Xu Chen and Youtao Zhang) would be asked to give the final decision.
Methodological quality assessment
Three authors (Shuzhou Yin, Yueyuan Wu and Xiaochun Jin) were asked to perform methodological quality evaluation according to predefined assessment criteria (Table 1), which was based on the items of Jiang et al. . The scores with the range of 0(lowest) to 18(highest) were based on several assessment items consisting of credibility of controls, matching criteria, diagnostic criteria, genotyping examination, sample size and Hardy-Weinberg equilibrium. The above assessment items were performed based on popular epidemiological methods and Brucellosis special characteristic. Literatures whose scores less than 12 were considered as “low-quality” literatures. However, the literatures with scores equal to or more than 12 were regarded as “high-quality” studies.
Both the OR and 95%CI were estimated to make the assessment of association power in four different genetic models. For position − 174 G/C of IL-6, there were allele comparison (G versus C), homozygote comparison (GG versus CC), recessive model (GG versus GC/CC), and dominant model (GG/GC versus CC). For position − 1082 G/A of IL-10, there were allele comparison (A versus G), homozygote comparison (AA versus GG), recessive model (AA versus AG/GG), and dominant model (AA/AG versus GG). For position − 819 C/T of IL-10, there were allele comparison (T versus C), homozygote comparison (TT versus CC), recessive model (TT versus TC/CC), and dominant model (TT/TC versus CC). For position − 592 loci of IL-10, there were allele comparison (A versus C), homozygote comparison (AA versus CC), recessive model (AA versus AC/CC), and dominant model (AA/AC versus CC). The χ2-test based Q-statistic and I2 statistics were used. If there was no evident heterogeneity, the fixed-effects model would be applied . If not, the random-effects model would be used . Sensitivity analysis was used to evaluate the stability of the results. Funnel plots and Egger’s test were applied to detect the potential publication bias . All statistics were conducted using Stata software (version 12.0; StataCorp LP, College Station, TX, USA).
A flow chart of the search course is displayed in Fig. 1. Seven literatures were enrolled in our meta-analysis [10, 20,21,22,23,24,25]. The detailed information of seven studies was listed in Table 2.
Quantitative synthesis of data
Detailed results for the relationship between these polymorphisms and Brucellosis susceptibility are displayed in Table 3. P < 0.05 is considered as a significant association. Overall, significant associations were only found in Asian population of − 819 loci under three genetic models as follows: (Allele model: OR = 0.60, 95%CI = 0.44–0.82, P = 0.001) (Fig. 2), (homozygote comparison: OR = 0.24, 95%CI = 0.09–0.62, P = 0.003) (Fig. 3), (recessive genetic model: OR = 0.22, 95%CI = 0.05–0.91, P = 0.036) (Fig. 4).
Sensitivity analysis and publication Bias
Sensitivity Analysis was carried out to indicate single literature’s effect on the final result through continuous removal of individual studies under every genetic model. In the present study, our results could not be affected by any study, suggesting its reliability and robustness. Although slightly asymmetrical funnel plots were discovered (P = 0.806), we could not find any distinct publication bias by Egger’s test under allele model (P = 0.509) or recessive model (P = 0.509).
Previous studies have explored the association of IL-10 and IL-6 polymorphisms with Brucellosis risk. The present meta-analysis consists of 7 studies for IL-10 and IL-6 polymorphisms. To the best of our knowledge, the present research was the first to investigate the relationship between IL-10 and IL-6 polymorphisms and Brucellosis risk. Our meta-analysis shows that IL-10 -819 loci polymorphism contributes no risk to Caucasian population but may be associated with decreased risk in Asian population. And other polymorphisms of IL-10 and IL-6 are not related with Brucellosis susceptibility. For IL-6 -174 G/C polymorphism, Budak et al. reported in 2007 that GG genotype was more common in patients of Brucellosis than healthy controls (40 cases and 50 controls) and concluded that GG genotype was a risk factor in the development of Brucellosis . However, the other four literatures reported no association between IL-6 -174 G/C polymorphism and Brucellosis susceptibility. We consider sample size as an important factor contributing to the above discrepancy. Small-study bias is not an emerging event when performing genetic association studies . The current result of IL-10 -819 loci polymorphism may be attributed to race. It is well-known that different races have discrepant genotype number and allele frequency. So that it is necessary to conduct subgroup analysis by race. However, we made a conservative conclusion. Only two literatures with 250 cases and 141 controls were enrolled in our study. Considering the small sample size of eligible studies, a new updated study needs to be urgently conducted. Regarding the significant association found at IL10–1082 AA genotyping in Caucasian, we feel that it is likely a false positive, which is based on two reasons. Firstly, the P value approximates 0.05 and is a critical value. Secondly, the P values of other genetic models are 0.217, 0.484, and 0.916(P > 0.05).
Extensive variations have been established in the frequencies of cytokine polymorphism among different healthy population, including the − 1082 G/A polymorphism, which has been investigated most widely in healthy populations. The allele frequency show wide difference in different countries and regions. For example, the prevalence rate of IL-10 -1082 G allele was reported to be 42.5, 48.9 3.8 and 13.0% in Iranians, Norwegians, Japanese and Koreans, respectively [27,28,29,30].
HWE has been demonstrated to be an important genetic equilibrium law. And if the genotype distribution of control population is not conforming to HWE law, selection bias would occur. It has been reported that several factors would contribute to HWE deviation including fixed mating, or allele do not achieve equilibrium condition and other possible factors, which may lead to inaccurate results of meta-analysis. In our meta-analysis, the problem contributes little influence on our results. Because genotype distributions of all studies are consistent with HWE, which improves reliability of our results greatly.
Genome-wide association study (GWAS) is undoubtedly an effective method to find the sequence of gene mutations, accordingly screening the specific SNP associated with disease. GWAS opens the door to the study of complex diseases by comparing patients’ SNP sites with control groups in whole genome, which identifies all variant allele mutation frequencies. It is regrettable that there is rare GWAS on Brucellosis. Sankarasubramanian etal pulished a study titled “Identification of genetic variants of Brucella spp. through genome-wide association studies” and they found special SNPs, which were closely related with human’s specificity and virulence. They also thought that origin of these special SNPs was early and might derive from B. abortus evolution . However, we think this is just a new start in exploring species-specific SNPs of Brucellosis. In the future, we should investigate further in the genomes and their roles of Brucellosis.
Despite a lot of work we have done, there are still disadvantages existing, which should be remarkable here. Firstly, the present research is still a small-sample study and we should be careful in the results. Secondly, some miscellaneous factors including age, sex and other environment factors are not considered and calculated for final estimates. So that the results of our meta-analysis are on the account of unadjusted estimates. We should also be careful in the results.
In conclusion, this is the first meta-analysis which investigates the association between IL-10 polymorphism and Brucellosis risk. In conclusion, IL-10 -819 loci polymorphism contributes no risk to Caucasian population but may be associated with decreased risk in Asian population. And IL-10 -1082 G/A, 592 loci and IL-6 -174 G/C polymorphism are not associated with Brucellosis risk.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Restricted Fragment Length Polymorphism
Single nucleotide polymorphism
Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.
Boschiroli ML, Foulongne V, O'Callaghan D. Brucellosis: a worldwide zoonosis. Curr Opin Microbiol. 2001;4(1):58–64.
Golding B, Scott DE, Scharf O, Huang LY, Zaitseva M, Lapham C, Eller N, Golding H. Immunity and protection against Brucella abortus. Microbes Infect. 2001;3(1):43–8.
Ko J, Splitter GA. Molecular host-pathogen interaction in brucellosis: current understanding and future approaches to vaccine development for mice and humans. Clin Microbiol Rev. 2003;16(1):65–78.
Araya LN, Elzer PH, Rowe GE, Enright FM, Winter AJ. Temporal development of protective cell-mediated and humoral immunity in BALB/c mice infected with Brucella abortus. J Immunol. 1989;143(10):3330–7.
Fernandes DM, Baldwin CL. Interleukin-10 downregulates protective immunity to Brucella abortus. Infect Immun. 1995;63(3):1130–3.
Seruga B, Zhang H, Bernstein LJ, Tannock IF. Cytokines and their relationship to the symptoms and outcome of cancer. Nat Rev Cancer. 2008;8(11):887–99.
Bidwell J, Keen L, Gallagher G, Kimberly R, Huizinga T, McDermott MF, Oksenberg J, McNicholl J, Pociot F, Hardt C, et al. Cytokine gene polymorphism in human disease: on-line databases, supplement 1. Genes Immun. 2001;2(2):61–70.
Sepanjnia A, Eskandari-Nasab E, Moghadampour M, Tahmasebi A, Dahmardeh F. TGFbeta1 genetic variants are associated with an increased risk of acute brucellosis. Infect Dis Ther. 2015;47(7):458–64.
Bravo MJ, de Dios CJ, Alonso A, Caballero A. Polymorphisms of the interferon gamma and interleukin 10 genes in human brucellosis. Eur J Immunogenet. 2003;30(6):433–5.
Orozco G, Sanchez E, Lopez-Nevot MA, Caballero A, Bravo MJ, Morata P, de Dios CJ, Alonso A, Martin J. Inducible nitric oxide synthase promoter polymorphism in human brucellosis. Microbes Infect. 2003;5(13):1165–9.
Bravo MJ, Colmenero Jde D, Alonso A, Caballero A. HLA-B*39 allele confers susceptibility to osteoarticular complications in human brucellosis. J Rheumatol. 2003;30(5):1051–3.
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.
Rezazadeh M, Hajilooi M, Haidari M, Rafiei A, Alavi SA, Keramat F. Association of susceptibility to brucellosis and interleukin-4 promoter polymorphism. Scand J Infect Dis. 2006;38(11–12):1045–9.
Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, Rennie D, Moher D, Becker BJ, Sipe TA, Thacker SB. Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis of observational studies in epidemiology (MOOSE) group. Jama. 2000;283(15):2008–12.
Jiang DK, Wang WZ, Ren WH, Yao L, Peng B, Yu L. TP53 Arg72Pro polymorphism and skin cancer risk: a meta-analysis. J Invest Dermatol. 2011;131(1):220–8.
Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. 1959;22(4):719–48.
DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7(3):177–88.
Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50(4):1088–101.
Budak F, Goral G, Heper Y, Yilmaz E, Aymak F, Basturk B, Tore O, Ener B, Oral HB. IL-10 and IL-6 gene polymorphisms as potential host susceptibility factors in brucellosis. Cytokine. 2007;38(1):32–6.
Karaoglan I, Pehlivan S, Namiduru M, Pehlivan M, Kilincarslan C, Balkan Y, Baydar I. 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.
Rasouli M, Kiany S, Behbin M. Interleukin-10 gene polymorphisms and susceptibility to brucellosis in Iranian patients. Iran J Immunol. 2008;5(2):131–5.
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 Investig. 2016;45(2):107–15.
Asaei S, Rasouli M, Moravej A. Interleukin-8 but not interleukin-6 variant may affect susceptibility to brucellosis. Iran J Immunol. 2013;10(3):158–66.
Bravo MJ, Colmenero JD, Queipo-Ortuno MI, Alonso A, Caballero A. TGF-beta1 and IL-6 gene polymorphism in Spanish brucellosis patients. Cytokine. 2008;44(1):18–21.
Chen X, Yan Y, Li P, Yang Z, Qin L, Mo W. Association of GSTP1 -313A/G polymorphisms and endometriosis risk: a meta-analysis of case-control studies. Eur J Obstet Gynecol Reprod Biol. 2013;171(2):362–7.
Myhr KM, Vagnes KS, Maroy TH, Aarseth JH, Nyland HI, Vedeler CA. Interleukin-10 promoter polymorphisms in patients with Guillain-Barre syndrome. J Neuroimmunol. 2003;139(1–2):81–3.
Bagheri M, Abdi-Rad I, Omrani D, Khalkhali HR. Heterogeneity of cytokine single-nucleotide polymorphisms among the Iranian and in the other east-south Asian populations. Transfus Med. 2006;16(3):192–9.
Miyazoe S, Hamasaki K, Nakata K, Kajiya Y, Kitajima K, Nakao K, Daikoku M, Yatsuhashi H, Koga M, Yano M, et al. Influence of interleukin-10 gene promoter polymorphisms on disease progression in patients chronically infected with hepatitis B virus. Am J Gastroenterol. 2002;97(8):2086–92.
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(3–4):143–7.
Sankarasubramanian J, Vishnu US, Gunasekaran P, Rajendhran J. Identification of genetic variants of Brucella spp. through genome-wide association studies. Infect Genet Evol. 2017;56:92–8.
The authors declare that there are no sources of funding to be acknowledged.
Ethics approval and consent to participate
Consent for publication
The authors declare that there are no competing interests associated with the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jin, X., Wu, Y., Yin, S. et al. Association between the IL-10 and IL-6 polymorphisms and brucellosis susceptibility: a meta-analysis. BMC Med Genet 21, 63 (2020). https://doi.org/10.1186/s12881-020-01006-0