Background {#Sec1}
==========

Behçet's disease (BD) is a chronic remitting systemic vasculitis of unrevealed etiology \[[@CR1]\]. BD is heterogeneous in onset \[[@CR2]\] that its common clinical manifestations are orogenital ulcers, skin lesions, and uveitis \[[@CR3]\]. BD is frequently observed in the Far East and the Middle East; however, it has distributed along the Old Silk Route due to immigration \[[@CR1], [@CR4]--[@CR6]\]. Although the exact cause of Behçet's disease is mostly indefinite, it is believed that the pivotal pathophysiologic event is an inflammatory reaction to environmental factors in a genetically predisposed host and presence of epigenetic changes \[[@CR7]--[@CR11]\]. Research has evidenced that various cytokine networks play a crucial role in the determination of the disease onset, course of the disease, and outcome \[[@CR12]\]. Interleukin-10 (IL-10) \[Gene ID: 3586\] is a multifunctional cytokine that acts upon other cytokines' signaling pathways \[[@CR13]\]. IL-10 is a Type II cytokine in a family, including IL-19, IL-20, IL-22, IL-26, and IL-29 \[[@CR14]\]. IL-10 is expressed mostly by activated monocytes/macrophages, natural killer cells, dendritic cells, mast cells, T lymphocytes (mainly Th2 subsets), and B lymphocytes \[[@CR15]\].

This immuno-regulatory cytokine was demonstrated to have an imperative cytokine-suppressing role in the autoimmunity and inflammatory responses as a result of its ability to the downregulate antigen presentation and macrophage activation \[[@CR16]--[@CR20]\]. The abundant evidence has revealed that low IL-10 expression is involved in the pathogenesis of inflammatory and autoimmune diseases such as BD \[[@CR21]\]. Furthermore, polymorphisms in multiple immuno-regulatory genes have indicated as a risk predisposition for the developing of BD owing to their effect on the cytokine production \[[@CR22]\]. The human IL-10 gene maps to 1q31--32 and contains five exons and four introns. The IL-10 promoter is highly polymorphic and these single nucleotide polymorphisms (SNPs) appear to be correlated to the expression level of IL-10. Three frequently investigated SNPs of IL-10 gene are − 1082 A to G (rs1800896) -lies within a putative Ets transcription factor binding site-, − 819 T to C (rs1800871) -placed within a putative positive regulatory region-, and − 592 A to C (rs1800872) -is located within a putative STAT-3 binding site and negative regulatory region- polymorphisms \[[@CR23], [@CR24]\]. Rs1518111 (2195 A to G) which is in linkage disequilibrium (LD) with rs1800896 and rs1554286 (2607 T to C) that is in a strong LD with rs1800571 are the intronic variant of second and third exons \[[@CR25], [@CR26]\].

Considering the role of IL-10 in BD and the correlation between the IL-10 gene polymorphisms and IL-10 production, several studies have been carried out to investigate the association of the IL-10 gene polymorphisms with the BD susceptibility. However, contradictory reported results and/or small sample size could lead to low statistical power and false-positive results. Therefore, in order to deal with these ambiguities, and potential biases, such as publication bias and providing increased statistical power, we turned to meta-analysis. In the present study, a meta-analysis of all available published studies was carried out to clarify the association between the IL-10 gene polymorphisms and the BD susceptibility.

Results {#Sec2}
=======

Literature search and study characteristics {#Sec3}
-------------------------------------------

The literature search filtered 2390 English and Persian-language potential articles related to the IL-10 polymorphisms and BD. After the removal of duplicates, other publications were reviewed based on the titles and abstracts. In the end, 24 eligible studies were selected with a full-text review where 5 studies were eliminated (two studies due to non-availability of the data, even no responses could be gotten after having sought the relevant data via email contacts, two were posters and one study was involving family members). A flowchart of the study selection process and its results is presented in Fig. [1](#Fig1){ref-type="fig"}. Fig. 1Flow diagram of studies for inclusion in the systematic review and meta-analysis

We conduct the meta-analysis, involving 10,626 cases and 13,592 controls. Of the 19 articles, studies including different subpopulations were considered as separate studies. Consequently, these groups were independently analyzed. Ten studies reported on -1082G/A, nineteen studies on -819C/T, thirteen studies on −592C/A, thirteen studies on rs1518111G/A and three studies on rs1554286C/T. The baseline characteristics of the eligible studies are described in Table [1](#Tab1){ref-type="table"}. Table 1Characteristics of studies included in the meta-analysisAuthor(year)EthnicitySample size case/controlMean-ageGenotyping methodsStudied polymorphismHWE P~value~\
(control)refCaseControlBaris S(2016)Turkish71/7037.8 ± 1.0--PCR-RFLPrs1800896--27Hu J. stage I(2015)Chinese Han300/35033.60 ± 8.7839.66 ± 11.19PCR-RFLPrs1800896, rs1800872, rs1554286*p* \> 0.0528Hu J. stage II(2015)Chinese Han418/140333.60 ± 8.7839.66 ± 11.19PCR-RFLPrs1800871p \> 0.0528Al-Okaily F(2015)Saudi61/21137.87 ± 12.536 ± 10PCR-ARMSrs1800896, rs1800871, rs1800872--29Talaat R. M(2014)Egyptian87/9734.37 ± 10.36--PCR-SSPrs1800896, rs1800871p \> 0.0530Xavier J. M(2012)Iranian973/63739.1 ± 11.042.0 ± 11.6Sequenom iPlex assayrs1800896, rs1518111, rs1554286*P* \> 0.0126Shahram F(2011)Iranian150/140----PCR-SSPrs1800896, rs1800871, rs1800872p \> 0.0531Ates O(2010)Turkish102/10237.2 ± 8.452 ± 7PCR-ARMSrs1800896, rs1800871, rs1800872--32Dilek K(2009)Turkish97/127----PCR-SSPrs1800896, rs1800871, rs1800872p \> 0.0533Wallace GR(2007)UK63/182----PCR-SSPrs1800896, rs1800871p \> 0.0534Wallace GR(2007)ME115/113----PCR-SSPrs1800896, rs1800871p \> 0.0534Kramer M(2018)Turkish64/68----Sanger techniquers1800871, rs1800872, rs1518111--35Kramer M(2018)Israeli25/20----Sanger techniquers1800871, rs1800872, rs1518111--35Yu H(2017)Chinese Han1206/247534.86 ± 9.935.46 ± 10.3Sequenom Mass Arrayrs1800871, rs1554286--36Carapito R(2015)Iranian552/417----Taqmanrs1800871, rs1518111--37Wu Z(2014)Chinese407/67938.02 ± 12.4438.81 ± 10.45Sequenom Mass Arrayrs1800871, rs1518111p \> 0.0538Khaib Dit Naib O(2013)Western Algeria51/9626 ± 11--Direct sequencingrs1800871, rs1800872--39Kirino Y(2013)Turkish1209/1278----GWASrs1800871, rs1518111*P* \< 0.000140Mizuki N(2010)Japanese611/737----GWASrs1800871, rs1800872P \< 0.0017Mizuki N(2010)Korean124/140----GWASrs1800871, rs1800872P \< 0.0017Mizuki N(2010)Turkish1215/1279----GWASrs1800871, rs1800872P \< 0.0017Afkari B(2018)Iranian47/5838.02 ± 10.2537.4 ± 8.5PCR-RFLPrs1800872p \> 0.0541Montes-Cano MA(2013)Spanish304/31338.7 ± 13.8--Taqmanrs1800872p \> 0.0542Remmers EF(2010)Discovery-Turkish1161/1221----GWASrs1518111*p* \> 0.000018Remmers EF(2010)Rplication-Turkish110/224----GWASrs1518111p \> 0.000018Remmers EF(2010)Replication-Middle Eastern Arab188/163----GWASrs1518111p \> 0.000018Remmers EF(2010)Replication-Greek107/84----GWASrs1518111p \> 0.000018Remmers EF(2010)Replication-UK Caucasian120/119----GWASrs1518111p \> 0.000018Remmers EF(2010)Replication-Korean77/52----GWASrs1518111p \> 0.000018Remmers EF(2010)Replication-Japanese611/737----GWASrs1518111p \> 0.000018

### Association between rs1800896 polymorphism and risk of BD {#Sec4}

Based on ten studies \[[@CR26]--[@CR34]\], including 1970 cases and 1930 controls, the association between IL-10 -1082G/A polymorphism and BD susceptibility was examined. The combined results showed that considering all the studies, the comparisons of allele and genotypes failed to detect any statistical association under the random-effects model (Table [2](#Tab2){ref-type="table"} and Fig. [2](#Fig2){ref-type="fig"}). Table 2Meta-analysis of the association between IL-10 polymorphisms and BD riskpolymorphismNo. of studiesComparisonTest of associationTest of heterogeneityEgger's test (P)Metareg's test (P)OR (95% CI)P-valueP- valueI^2^ (%)Q testrs1800896 A/G10GG vs. AA0.925 (0.684--1.251)0.6130.26619.99.980.7460.689GG + GA vs. AA0.936 (0.733--1.195)0.5960.02056/018.200.8360.866GG vs. GA + AA0.975 (0.660--1.441)0.9000.02554.417.550.6190.640AG vs. AA0.945 (0.717--1.245)0.6870.00762.221.170.8220.870G vs. A0.965 (0.831--1.120)0.6400.08441.015.250.8160.807rs1800871 T/C19CC vs. TT0.466 (0.368--0.589)0.0000.00357.432.880.0810.347CC + CT vs. TT0.692 (0.584--0.820)0.0000.00066.141.330.8420.674CC vs. CT + TT0.557 (0.405--0.767)0.0000.00086.6119.020.0320.123CT vs. TT0.842 (0.572--1.240)0.3840.00093.7220.520.5420.883C vs. T0.691 (0.626--0.762)0.0000.00067.849.730.3630.493rs1800872 A/C13CC vs. AA0.723 (0.440--1.186)0.1990.00072.432.630.1610.360CC + CA vs. AA0.713 (0.535--0.951)0.0210.00462.423.940.1160.276CC vs. CA + AA0.776 (0.583--1.032)0.0820.00460.327.710.2100.367CA vs. AA0.716 (0.546--0.940)0.0160.02453.019.130.1010.313C vs. A0.779 (0.662--0.915)0.0020.00070.233.590.1940.423rs1518111 A/G13GG vs. AA0.570 (0.471--0.690)0.0000.4590.01.560.6220.741GG + GA vs. AA0.697)0.596--0.814)0.0000.3930.01.870.6840.533GG vs. GA + AA0.701 (0.623--0.790)0.0000.5050.03.320.5580.967AG vs. AA0.786 (0.667--0.927)0.0040.7230.00.650.6040.645G vs. A0.738 (0.681--0.800)0.0000.10636.715.790.5000.352rs1554286 T/C3CC vs. TT0.508 (0.372--0.694)0.0000.19838.23.230.7950.868CC + CT vs. TT0.605)0.467--0.785(0.0000.08759.14.890.6260.466CC vs. CT + TT0.665)0.484--0.916(0.0120.04168.76.400.4870.893CT vs. TT0.646 (0.502--0.832)0.0010.11553.84.330.6560.506C vs. T0.582 (0.440--0.770)0.0000.00087.716.230.3840.265 Fig. 2Forest plot of the association of IL-10 polymorphisms with BD. A random-effects model for the OR with 95% CI was used to detect an association between IL-10 polymorphisms with BD. (**a**) rs1800896 (**b**) rs1800871 (**c**) rs1800872 (**d**) rs1518111 (e) rs1554286

### Association between rs1800871 polymorphism and risk of BD {#Sec5}

In total, we identified nineteen studies \[[@CR7], [@CR28]--[@CR40]\] which contained 6854 cases and 9911 controls assessed the effect of the IL-10 -819C/T polymorphism in the occurrence of BD. By using the random-effects model, meta-analysis revealed a significant association between BD and the IL-10 -819C/T polymorphism under the allelic (OR = 0.691, 95% CI: 0.626--0.762, *P* \< 0.001), homozygous (OR = 0.466, 95% CI: 0.368--0.589, *P* \< 0.001), dominant (OR = 0.692, 95% CI: 0.584--0.820, *P* \< 0.001), and recessive (OR = 0.557, 95% CI: 0.405--0.767, *P* \< 0.001) models, while no significant association was found under the heterozygous model (*P* = 0.384) (Table [2](#Tab2){ref-type="table"} and Fig. [2](#Fig2){ref-type="fig"}).

### Association between rs1800872 polymorphism and risk of BD {#Sec6}

For IL-10 -592C/A polymorphism, thirteen studies \[[@CR7], [@CR28], [@CR29], [@CR31]--[@CR33], [@CR35], [@CR39], [@CR41], [@CR42]\] were involved composed of 3132 BD patients and 3638 controls. The results of the combined analysis of the association between -592 polymorphism and BD are summarized in Table [2](#Tab2){ref-type="table"} and Fig. [2](#Fig2){ref-type="fig"}. There were significant associations for C vs. A allele, CC + CA vs. AA and CA vs. AA genotypes; however, other comparisons failed to obtain any significant association.

### Association between rs1518111 polymorphism and risk of BD {#Sec7}

Thirteen studies \[[@CR8], [@CR26], [@CR35], [@CR37], [@CR38], [@CR40]\], including 5519 cases and 5643 controls focused on the relationship between rs1518111 polymorphism and BD risk. The combined results (Table [2](#Tab2){ref-type="table"} and Fig. [2](#Fig2){ref-type="fig"}) showed significant differences in all models between BD cases and healthy controls.

### Association between rs1554286 polymorphism and risk of BD {#Sec8}

In total, three studies \[[@CR26], [@CR28], [@CR36]\] which contained 2421 cases and 3429 controls were performed to assess the importance of rs1554286 polymorphism for BD susceptibility. As shown in Table [2](#Tab2){ref-type="table"} and Fig. [2](#Fig2){ref-type="fig"}, this polymorphism was found to be significantly associated with BD susceptibility under all models.

Heterogeneity and publication bias {#Sec9}
----------------------------------

A funnel plot was used to reveal any publication bias influencing the analysis. The funnel plots did not show a significant sign of asymmetry for IL-10 polymorphisms. For rs1554286, a funnel plot was not performed because it is useless when the number of studies is limited (Fig. [3](#Fig3){ref-type="fig"}). The Metareg's test and Egger's test were also conducted and the findings showed a lack of publication bias among all comparison models except heterozygous (CC vs. CT + TT) model of rs1800871. Fig. 3Begg's funnel plot for publication bias analysis of the IL-10 polymorphisms with BD. Symmetry in Begg's funnel plots demonstrate the absence of publication bias in the studies investigating the association of (**a**) rs1800896 (**b**) rs1800871 (**c**) rs1800872 (**d**) rs1518111

To test heterogeneity among the selected studies, the I^2^ statistic was employed. The heterogeneity was significant under the genetic models of CC vs. CT + TT and CT vs. TT of rs1800871. After excluding the study by Kirino \[[@CR40]\], the heterogeneity was eliminated (I^2^ = 44.3, association *P* \< 0.001) and (I^2^ = 19.6, association *P* \< 0.001), respectively. Significant heterogeneity was detected in C vs. T model of rs1554286 polymorphism. It can be due to a small number of studies included in this analysis. Therefore, the random-effects model was applied to synthesize the data (Table [2](#Tab2){ref-type="table"}).

Discussion {#Sec10}
==========

In order to better understand the possible genetic association of the IL-10 gene polymorphisms and BD, we conducted a systematic review and meta-analysis using pooling the results of the independent usable studies to determine whether the IL-10 polymorphisms have a potential impact on the BD susceptibility.

Behçet's disease is a multi-system inflammatory disorder with complex etiology, therefore genes involved in the immune system and inflammatory responses are potential candidate genes for the BD development. In addition to HLA-B51 as the strongest genetic factor for BD, which explains 20% incidence of BD \[[@CR43], [@CR44]\], GWAS and subsequent genetic studies have identified a number of non-HLA potent loci for the susceptibility and organ damages of this disease, such as cytokines (IL-4 \[[@CR45]\], IL-6 \[[@CR46]\], IL-10 \[[@CR7], [@CR40]\], CCL2 \[[@CR47], [@CR48]\], TNF-α \[[@CR49]\]), ERAP1 \[[@CR40]\], UBASH3B \[[@CR50], [@CR51]\], STAT4 \[[@CR52]\].

Several studies have been conducted on the influence and impression of the interleukin-10 gene and its protein on the probability of the multiple autoimmune and inflammatory diseases, such as BD relying on its chromosomal location and functional relevance.

Research on the effect of IL-10 on various inflammatory models has led to different and sometimes contradictory results even in the same models, indicating that there are intricate influenced agents on the regulation of cellular responses by IL-10. In the experimental study by Rizzo et al. \[[@CR53]\] indicated that IL-10 knockout mice had a susceptibility to develop the Experimental Autoimmune Uveoretinitis (EAU) in comparison to controls, also, IL-10-treatment in induced mice for the EAU showed a resistance to the disease. Other animal studies have also confirmed the protective effect of IL-10 on the ocular inflammation and CNS autoimmune \[[@CR54]--[@CR56]\]. Some studies on BD have observed a reduction in IL-10 protein level in BD patients' peripheral blood and aqueous humor \[[@CR13], [@CR30], [@CR41], [@CR57]\], while in other studies, there has been an increase in the concentration of this protein in plasma and active ulcers of BD \[[@CR58]--[@CR63]\]. On the other hand, there is a study that has not recorded significant changes between patients and controls \[[@CR64]\]. Similar to the protein concentration determination studies, the results of the mRNA expression level were varied among previous studies. Several investigations \[[@CR13], [@CR41], [@CR64]\] have shown low expression of the IL-10 mRNA in B cells and serums of BD patients. In contrast, Ahmed Ben et al. \[[@CR65]\] have reported a high expression of the IL-10 mRNA in oral/genital ulcers, pseudofolliculitis lesions and sites of positive pathergy test. In contrast, Balkan et al. \[[@CR66]\] found no significant difference in the IL-10 mRNA expression.

So far, different factors have been studied in the production and regulatory the role of IL-10. One of these influential factors that have recently been studied on IL-10 is epigenetic agents. The study of Alipour et al. \[[@CR13]\] showed a hyper-methylation on IL-10 gene promoter region. Genetic variants are other factors that account for up to 75% of the IL-10 expression variance in inter-individual human \[[@CR67]\]. Investigation of the polymorphic region in the transcriptional regulation of IL-10 gene may be helpful in predicting inflammatory diseases such as BD. Three common polymorphisms of IL-10 (− 1082G/A, − 819C/T and -592C/A) have been evaluated in detail for their impacts on gene transcription. In vitro studies showed that the highest IL-10 expression has been observed in GCC (− 1082/− 819/− 592) haplotype and lowest expression in AT (− 1082/− 819) haplotype \[[@CR68]--[@CR70]\]. However, studies on BD have not confirmed the effects of the above alleles \[[@CR30]\]. In addition to the promoter variants, the intronic variant of rs1518111 is also likely to influence transcription \[[@CR66]\].

The present meta-analysis was performed to examine the association between the IL-10 polymorphisms (− 1082, − 819, − 592, rs1518111, and rs1554286) and BD susceptibility in different ethnic groups including East Asian, Middle Eastern, European, and African populations. The strongest associations were observed in the comparisons of allele and genotypes of − 819C/T promoter SNP, rs1518111G/A intronic SNP, and rs1554286C/T intronic SNP. The overall results indicated the protective impact of these polymorphisms against the disease. So that, individuals carrying C allele of −819 and rs1554286 have 31 and 42% lower risk of developing BD, respectively. We also achieved significant association in -592C/A promoter SNP, which showed a relatively protective role against BD. However, the comparisons of allele and genotypes of -1082G/A promoter SNP failed to demonstrate any significant association.

However, our results should be interpreted cautiously, with respect to the following limitation reasons. First, limited number of included studies in some studied SNPs; second, heterogeneity and publication bias, for instance, missing articles due to unpublished studies, possibly with negative results, and/or language bias due to writing in other languages (except English and Persian); third, lack of sufficient data on the characteristics and symptoms of the participants in most of the primary studies to perform subgroups analysis, such as based on age, gender and the disease activity; lastly, some genotype distributions of control groups do not follow the HWE.

Conclusions {#Sec11}
===========

Our data demonstrated significant associations of the IL-10 polymorphisms with BD. Accordingly, the polymorphisms of -819, − 592, rs1518111 and rs1554286 play relatively protective roles against Behcet's disease.

Methods {#Sec12}
=======

Search strategy {#Sec13}
---------------

A comprehensive search for the IL-10 polymorphisms and BD was performed in PubMed, Embase, Web of Science, and Scopus databases for articles in English, and SID, Magiran, IranMedex, and IranDoc for articles in Persian up to February 2019 according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) checklist and PICO approach \[[@CR43]\]. The most relevant articles were identified by computerized searches using relevant text words and medical subject headings (MeSH) that included all spellings and combination of "Behcet's syndrome" or "Behcet's disease" and "Interleukin-10" or "IL-10". Besides, we hand-searched reference lists of selected publications and previous meta-analysis study to ensure that all relevant studies and recent reviews were included. The search was limited to case-control studies and Genome-Wide Association Studies (GWAS). No restrictions were placed on race, ethnicity or geographic area publication date.

Inclusion and exclusion criteria {#Sec14}
--------------------------------

A study was included in this meta-analysis if (1) it was published by February 2019, (2) it was a case-control study or GWAS that determined the distributions of the IL-10 -1082G/A, −819C/T, −592C/A, rs1518111G/A, and rs1554286C/T polymorphisms, in patients with BD and controls, (3) it contained original data, and (4) it provided enough data to calculate odds ratios (ORs). Exclusion criteria were as follows: (1) Keywords/search terms/reviews, and abstract scanning criteria, (2) the study was a review or abstract and publications in duplicate, (3) unrelated to each of BD and these polymorphisms, (4) studies involving family members, because their analysis was based on linkage considerations, and (5) papers do not involve human subjects.

Data extraction {#Sec15}
---------------

Two authors (E.S. and N.R.) independently extracted the following data from all studies selected by using a study-specific data extraction form: first author's name, year of publication, ethnicity, mean age of participants, number of cases and controls and method for determination of genotype, studied polymorphisms, Hardy--Weinberg equilibrium (HWE) *P*-values. Any disagreement between authors was resolved by discussion until a consensus was reached. If they failed to reach an agreement, a third author (M.G.) was consulted to resolve the discrepancies. Due to the lack of the mean age and HWE *P*-values for some studies, we were not able to analyze the data based on this information. This study was conducted in the Liver and Gastrointestinal Disease Research Center, Tabriz University of Medical Sciences.

Statistical analyses {#Sec16}
--------------------

All statistical manipulations were performed with STATA version-15.0 software (STATA Corporation, College Station, Texas). The strength of the association between IL-10 polymorphisms and BD was assessed by calculating pooled OR and 95% CI values in the following genetic models: allelic model (G vs. A, C vs. A, C vs. T), dominant model (GG + GA vs. AA, CC + CT vs. TT, CC + CA vs. AA), recessive model (GG vs. GA + AA, CC vs. CT + TT, CC vs. CA + AA), homozygous model (GG vs. AA, CC vs. TT, CC vs. AA) and heterozygous model (AG vs. AA, CT vs. TT, CA vs. AA). Z-test with *P* \< 0.05 was used to authenticate the statistical significance of effect size. Sensitivity analysis was used to investigate the cause of dispersion. Then, outlier studies were removed and analysis was recomputed. In the case of a significant reduction in dispersion (I^2^ index), this study was considered as a dispersion factor. With regard to the heterogeneity of the studies, Cochran's Q test (*p*-value \[phet\] \< 0.10 was considered as statistically significant heterogeneity) and I^2^ statistics (75 ≤ I^2^ \< 100 as extreme heterogeneity, 50 ≤ I^2^ \< 75 as high heterogeneity, 25 ≤ I^2^ \< 50 as moderate heterogeneity, and I^2^ \< 25 as no heterogeneity) were used to assess the degree of heterogeneity \[[@CR71]\]. The random-effects model was used for meta-analysis because it accounts for random variability both within and among studies. Publication bias was investigated by Funnel plot and Egger's test. A *P-value* \< 0.05 was considered statistically significant for all analyses.

BD

:   Behcet's disease

CI

:   Confidence Interval

EAU

:   Experimental Autoimmune Uveoretinitis

GWAS

:   Genome-Wide Association Study

HWE

:   Hardy--Weinberg equilibrium

IL-10

:   Interleukin-10

LD

:   linkage disequilibrium

OR

:   Odds Ratio

SNP

:   Single Nucleotide Polymorphism

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The study protocol was reviewed and approved by the Student Research Committee, Tabriz University of Medical Sciences, Tabriz, Iran (63975).We would like to thank "Student Research Committee of Tabriz University of Medical Sciences". Furthermore, we would like to express our deepest appreciation to Dr. Elaine F. Remmers, Dr.Nobuhisa Mizuki, Dr.Maria Francisca Gonzalez-Escribano and Dr.Joana M Xavier those who provided us with the possibility to complete our data.

Elham Shahriyari participated in conception and design, searching and selecting papers, data extraction, analyzing data, writing and approving final paper. L.V. participated in the selection of title, keywords, and approving final paper. N.R. participated in searching and selecting papers and approving final paper. M.A. participated in methodological and statistical advice, and approving final paper. A.K participated in data extraction and approving final paper. M.G. participated in conception and design, searching and selecting papers, data extraction, analyzing data, writing and approving final paper. All authors read and approved the final manuscript.

Authors' information {#FPar1}
====================

Not applicable.

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

The datasets used and/or analyzed during the current study are available from the corresponding and first authors on reasonable request.

Ethical approval is not required, as our study does not include confidential participant data and interventions. We disseminated our results through previous publications.

The meta-analysis is exempt from informed consent because we collected and analyzed data from previous publications in which informed consent had already been obtained.

The authors declare that they have no competing interests.
