Association of IL-4 and IL-10 Polymorphisms With Preterm Birth Susceptibility: A Systematic Review and Meta-Analysis

Objective: Preterm birth (PTB) is a typical inflammatory disease with unclear pathogenesis. The studies investigating the relationship between anti-inflammatory factors IL-4 and IL-10 gene polymorphisms and PTB produced conflicting results. This systematic review and meta-analysis aimed to summarize the effects of IL-4 and IL-10 gene polymorphisms and clarify their possible association with PTB.

Methods: A systematic literature review was conducted using PubMed, Web of Science, and Cochrane library (up to 02 April 2022). The MeSH terms, related entry terms, and other names in “Gene” database were used to find relevant articles. A fixed- or random-effects model was used to calculate the significance of IL-4 and IL-10 gene polymorphisms, depending on study heterogeneity. The odds ratios (OR) and 95% confidence intervals (CIs) were calculated in the allele, recessive, dominant, co-dominant, and over-dominant models. The Eggers publication bias plot was used to graphically represent the publication bias.

Results: Polymorphisms in two interleukins (IL-4-590C/T (rs2243250) = 5 and IL-10-592A/C (rs1800872), -819T/C (rs1800871) and -1082A/G (rs1800896) = 16) were found in 21 articles. Overall, only the over-dominant gene model AA + GG vs. AG revealed significant association between IL-10-1082A/G (rs1800896) and PTB (OR [95% CI] = 0.87 [0.76, 0.99], p = 0.04). However, in the allele model, recessive model, dominant model, co-dominant model, and over-dominant model, the polymorphisms for IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872), and IL-10-819T/C (rs1800871) were not found to be associated with the risk of PTB. In gene models, no statistically significant association was found between IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872), IL-10-819T/C (rs1800871), and IL-10-1082A/G (rs1800896) polymorphisms and PTB in subgroup analyses by racial or control group Hardy-Weinberg Equilibrium (HWE) p-value. Eggers’s publication bias plot and heterogeneity test (I²<50%, p = 0.05) of IL-10-1082A/G (rs1800896) suggested that the funnel asymmetry could be due to publication bias rather than heterogeneity.

Conclusion: The current study suggests that the over-dominant gene model AA + GG vs. AG of IL-10-1082A/G (rs1800896) polymorphism may be associated with genetic susceptibility to PTB and may have a protective function against PTB risk. There was unclear association found between IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872) and IL-10-819T/C (rs1800871) polymorphisms and PTB. Due to the limitations of included studies and the risk of publication bias, additional research is required to confirm our findings.

Systematic Review Registration: https://inplasy.com/inplasy-2022-4-0044, identifier INPLASY202240044.

Introduction

Preterm birth (PTB) is defined by the World Health Organization (WHO) as babies born alive before 37 weeks of pregnancy (1). According to the new estimates, PTB prevalence in 2014 ranged from 8.7% to 13.4% of all live births, with approximately 15 million preterm babies born each year (2). PTB is the leading cause of death in children under five years globally (3). Furthermore, premature babies are at a higher risk of short- and long-term complications caused by the immaturity of multiple organ systems, such as cerebral palsy, intellectual disabilities, vision and hearing impairments, and impaired cognitive development (4, 5). PTB has become a global public health issue, with mounting evidence indicating a syndrome attributed to various pathological processes (5). Many studies have revealed that genetic variations in pro-inflammatory cytokines such as tumor necrosis factor-α (TNF-α) and interleukin-1α (IL-1α) are associated with an increased risk of PTB (6). However, the relationship between genetic polymorphisms in anti-inflammatory cytokines and PTB risk remains controversial. The transition from a quiescent to a pro-inflammatory environment during pregnancy initiates labor, confirmed by the increasing production of pro-inflammatory cytokines, such as interleukins and TNF-α during labor (7). Cytokines are required during labor to initiate and regulate uterine contractions, cervical ripening, and fetal membrane rupture (7).

In contrast to the up-regulation of PTB-related pro-inflammatory cytokines, anti-inflammatory cytokines such as IL-10 and IL-4 can provide a compensatory protective role during pregnancy by limiting macrophage production of pro-inflammatory cytokines (8, 9). PTB risk is associated with a shift in the cytokine profile toward pro-inflammatory and decreased anti-inflammatory cytokines. However, in women without identifiable infectious agents, high levels of pro-inflammatory cytokines and low levels of anti-inflammatory cytokines may be caused by other factors (e.g., a genetically determined predisposition to up-regulated or down-regulated synthesis of pro-inflammatory and anti-inflammatory cytokines due to gene polymorphisms) (10). Numerous studies have linked cytokines and diseases, particularly some key single nucleotide polymorphisms (SNPs) in these cytokine genes that may interfere with gene expression, potentially impacting disease pathogenesis.

Various cells secrete IL-4 and IL-10. They play an important role in maintaining tissue homeostasis during infection and inflammation by limiting excessive inflammatory responses, enhancing innate immunity, and promoting tissue repair mechanisms (11–13). Therefore, IL-4 and IL-10 may play an important role in PTB, characterized by infection and inflammation (14, 15). In vitro experiments revealed that IL-4 and IL-10 could down-regulate lipopolysaccharide (LPS)-mediated inflammatory response in human pregnancy-related tissues (14). Animal studies also found that IL-10 could effectively prevent LPS-induced PTB (15). Therefore, understanding the roles of IL-4 and IL-10 in PTB is beneficial for clinical transformation research and may provide potential targets for PTB treatment.

Many studies have found gene polymorphisms in cytokines like IL-4 and IL-10 may play important roles in PTB. However, the findings are contradictory. As a result of our extensive literature search and preliminary screening, we chose the two most studied anti-inflammatory factors, IL-4 and 10, for further investigation. This systematic review and meta-analysis aimed to determine the relationship between SNPs of anti-inflammatory cytokine gene IL-4 and IL-10 and PTB.

Methods

Search Strategy

A systematic literature search was conducted using PubMed, Web of Science, and Cochrane Library, with articles published until 02 April 2022. To ensure a comprehensive literature search, the reference lists of included studies and related reviews were manually searched and screened. There was no language restriction. Two authors (XL and XY) independently searched all databases using the following search terms: (“Polymorphism, Genetic” OR “Polymorphisms, Genetic” OR “Genetic Polymorphism” OR “Genetic Polymorphisms” OR “Gene Polymorphism” OR “Gene Polymorphisms” OR “Polymorphism, Gene” OR “Polymorphisms, Gene” OR “Polymorphism (Genetics)” OR “Polymorphisms (Genetics)” OR “Genetic Variation” OR “Genetic Variations” OR “Variations, Genetic” OR “Variation, Genetic” OR “Diversity, Genetic” OR “Diversities, Genetic” OR “Genetic Diversities” OR “Genetic Diversity”) AND (“Birth, Premature” OR “Births, Premature” OR “Premature Births” OR “Preterm Birth” OR “Birth, Preterm” OR “Births, Preterm” OR “Preterm Births” OR “Premature Birth” OR “Preterm Premature Rupture of the Membranes” OR “PPROM” OR “spontaneous preterm birth” OR “premature labor” OR “Obstetric Labor, Premature” OR “Labor, Premature Obstetric” OR “Preterm Labor” OR “Labor, Preterm” OR “Premature Obstetric Labor” OR “Labor, Premature” OR “Premature Labor” OR “Premature delivery” OR “Preterm delivery”) AND ((“CSIF” OR “TGIF” OR “GVHDS” OR “IL-10” OR “IL10A” OR “Interleukin 10” OR “IL10” OR “CSIF-10” OR “Cytokine Synthesis Inhibitory Factor”) OR (“BSF1” OR “BCGF1” OR “BSF-1” OR “Interleukin 4” OR “Interleukin-4” OR “B-Cell Growth Factor-1” OR “B Cell Growth Factor 1” OR “B-Cell Growth Factor-I” OR “B Cell Growth Factor I” OR “B-Cell Proliferating Factor” OR “B Cell Proliferating Factor” OR “B-Cell Stimulating Factor-1” OR “B Cell Stimulating Factor 1” OR “B-Cell Stimulatory Factor 1” OR “B-Cell Stimulatory Factor-1” OR “BCGF-1” OR “Binetrakin” OR “IL-4” OR “IL4” OR “Mast Cell Growth Factor-2” OR “Mast Cell Growth Factor 2” OR “MCGF-2” OR “B Cell Stimulatory Factor-1” OR “B Cell Stimulatory Factor 1”).

Inclusion and Exclusion Criteria

The following were the inclusion criteria for our systematic review and meta-analysis: 1) case-control or cohort studies; 2) the case group of patients diagnosed with PTB or preterm premature rupture of membranes (PPROM), and the control group of term birth healthy individuals; 3) investigating the association between IL-4 and IL-10 polymorphisms and PTB or PPROM; 4) sufficient data supporting the genotype distribution were provided for the calculation of odds ratios (ORs) and corresponding 95% confidence intervals (CIs), and 5) full-text articles can be obtained. The following were the exclusion criteria: 1) studies that did not meet any of the inclusion criteria; and 2) studies that contained repeated data.

Study Selection and Data Extraction

Using Endnote (EndNote 2020), duplicate studies were excluded during data selection. The titles and abstracts of included studies were further screened, and the full text was reviewed using the inclusion and exclusion criteria. The following characteristics were extracted from each eligible literature by two reviewers (XL and XY) independently following the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidance (16): first author’s surname, publication year, country of the study, age of subjects, the phenotype of case group, single nucleotide polymorphism (SNP) genotyping method, numbers of cases and controls, PTB time of case group, Newcastle-Ottawa Scale (NOS) quality score of each study, Hardy-Weinberg Equilibrium (HWE) p-value of controls and the genotype frequencies of IL-4 and IL-10 gene polymorphisms in cases and controls. Discrepancies were resolved through discussion until consensus was reached, and if inconsistencies emerged, an expert (CM) was referred.

Quality Assessment

The two reviewers (XL and XY) independently assessed the quality of the methods for the included literature using NOS to determine the quality of non-randomized studies in meta-analyses (17). NOS assesses quality utilizing a star rating system; studies with scores ranging from 0 to 9 stars and ≥7 stars were considered high quality. Disputes between the two reviewers were resolved in the same manner described above.

Assessment of Heterogeneity and Publication Bias

According to Cochrane Q test, statistical heterogeneity was assessed using the standard χ² test (α = 0.1) and I² test. A fixed-effects model (Mantel–Haenszel method) was used to confirm collective effectiveness If p ≥ 0.05 and I² ≤ 50% were met. Random-effects models was used if p < 0.05 or I² > 50%. Eggers’s publication bias plot was used to represent the publication bias graphically. Trim-and-fill method was used to adjust the publication bias. Besides, publication bias (p < 0.05) was carefully discussed in our discussion section and listed as one of the limitations of our meta-analysis.

Statistical Analysis

The meta-analysis was performed using Stata (version 17.0) software. ORs and 95% CIs were calculated to assess the association between IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872) (referred as T/G in “dpSNP” database), IL-10-819T/C (rs1800871) (referred as A/G in “dpSNP” database), and IL-10-1082A/G (rs1800896) (referred as T/C in “dpSNP” database) polymorphisms and PTB in the allele model T-allele vs. C-allele, C-allele vs. A-allele, C-allele vs. T-allele and G-allele vs. A-allele, recessive model TT vs. CT + CC, CC vs. AA + AC, CC vs. TT + TC and GG vs. AA + AG, dominant model CT+TT vs. CC, AC + CC vs. AA, TC + CC vs. TT and AG + GG vs. AA, co-dominant model TT vs. CC and CT vs. CC, CC vs. AA and AC vs. AA, CC vs. TT and TC vs. TT, AG vs. AA and GG vs. AA and over-dominant model CC + TT vs. CT, AA + CC vs. AC, TT + CC vs. TC and AA + GG vs. AG respectively. According to HWE, in the case of a pair of alleles, the relationship between gene p (dominant) and gene q (recessive) gene frequency is as follows: p + q = 1, p2 + 2pq + q2. If only genotype or allele frequency information was provided in the included study, we used this formula to calculate the number of specific gene models. Subgroup analysis was used to investigate the potential sources of heterogeneity by racial or by HWE p-value of controls. The goodness-of-fit χ2 test assessed the control group’s HWE deviation. We performed sensitivity analyses to determine the robustness of our findings. For sensitivity analysis, a leave-one-out meta-analysis was used. A p-value < 0.05 was deemed statistically significant.

Trial Sequential Analysis

Trial sequential analysis (TSA) was performed by using the TSA-Trial Sequential Analysis Viewer (version 0.9.5.10 β, Copenhagen Trial Unit, Copenhagen, Denmark) (18). The required information size and TSA monitoring boundaries were then generated automatically.

Expression Quantitative Trait Loci Assessment for SNPs

We searched for expression quantitative trait Loci (eQTL) for the SNPs, using one publicly available database (QTLbase: http://www.mulinlab.org/qtlbase/index.html). Results were presented as nominal p values for each SNP.

Results

Study Selection

According to PRISMA flow diagram, a systematic literature search of PubMed, Web of Science, and Cochrane library yielded 191 studies, including 26 duplicate records (Figure 1, and Table S1). We excluded 124 studies after reviewing the titles and abstracts of 165 non-duplicate articles, for the following reasons: 1) studies were review articles (n = 16); 2) studies were meeting reports (n = 10); 3) studies were about PTB complications (n = 9); and 4) other unrelated studies (n = 92). Of the remaining 38 studies, 15 were excluded because they lacked eligible data, and one study was excluded because its data was duplicated (19). Finally, 21 eligible studies (IL-4 = 5, IL-10 = 16, contains four repetitive studies) with a total sample size ranging from 29 to 559 cases were included in this study. These studies included PTB or PPROM cases as well as healthy controls. Using these studies, we investigated the relationship between alleles and genotypes of the IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872), IL-10-819T/C (rs1800871), and IL-10-1082A/G (rs1800896) genetic polymorphisms and the risk of PTB.

FIGURE 1

Figure 1 The PRISMA flow diagram. PRISMA, Preferred Reporting Items for Systematic Reviews and Meta-Analyses.

Study Characteristics

Table 1 summarizes the key characteristics of included studies in our systematic review meta-analysis. From 2004 to 2020, these studies included 538 cases and 1123 controls for IL-4-590C/T (rs2243250), 787 cases and 1756 controls for IL-10-592A/C (rs1800872), 1225 cases and 1769 controls for IL-10-819T/C (rs1800871) and 1612 cases and 3386 controls for IL-10-1082A/G (rs1800896), respectively. Because not all genotypes were identified, the number of genotypes was somewhat inconsistent with the number of cases and controls. Most case groups were PTB patients, with some studies including PPROM patients. Most studies had NOS scores greater than or equal to 7, indicating high quality.

TABLE 1

Table 1 Main characteristics of included studies.

Results of Meta-Analysis

Five studies in our review investigated the association between IL-4-590C/T (rs2243250) polymorphism and PTB. No significant association was identified between the overall risk of PTB and SNP under different genetic models (Table 2 and Figure S1). A racial subgroup analysis revealed no significant association between PTB and IL-4-590C/T (rs2243250) polymorphism under different genetic models (Figure S2).

TABLE 2

Table 2 Main results of the total and subgroup analysis.

Seven studies evaluated the relationship between IL-10-592A/C (rs1800872) polymorphism and PTB. No significant association was found between the overall risk of PTB and rs1800872 under different genetic factors models (Table 2 and Figure S3). Subgroup analysis by racial and HWE p-value of controls also indicated no significant association between PTB and the IL-10-592A/C (rs1800872) polymorphism under different genetic models (Figures S4 and S5).

Six studies were included to investigate the relationship between IL-10-819T/C (rs1800871) polymorphism and PTB, and no significant association was found between the overall risk of PTB and rs1800871 under different genetic models (Table 2, and Figure S6). Subgroup analysis by HWE p-value of controls also revealed no significant association between PTB risk and the IL-10-819T/C (rs1800871) polymorphism under different genetic models (Figure S7).

Eleven studies were included to investigate the link between the IL-10-1082A/G (rs1800896) polymorphism and PTB risk. The over-dominant gene model AA + GG vs. AG revealed a significant association between IL-10-1082A/G (rs1800896) and PTB (OR [95% CI] = 0.87[0.76, 0.99], p = 0.04) (Table 2 and Figure 2). Under the over-dominant gene model, this polymorphism had a protective effect against PTB susceptibility. In the Austrian, American, Brazilian, Brazilian, and German subgroups, there was no significant association between PTB risk and rs1800896 under different genetic models (Figure S8). However, a potential significant association between IL-10-1082A/G (rs1800896) and PTB was found from the allele model G-allele vs. A-allele, recessive model GG vs. AA + AG, dominant model AG + GG vs. AA, and over-dominant model AA + GG vs. AG (Figure S8) in Indian. However, it should be noted that the sample size is small, and only one Indian study is included.

FIGURE 2

Figure 2 Forest plots of IL-10-1082A/G(rs1800896) polymorphism and PTB. The allele model (G vs A), recessive model (GG vs. AA+AG), dominant model (AG+GG vs AA), co-dominant model (AG vs. AA) and over-dominant model (AA+GG vs. AG). The count for genotypes, weight, OR, 95% confidence interval for each study. The fixed effect and random effect models were respectively utilized according to heterogeneity.

In terms of IL-10-1082A/G (rs1800896) publication bias, there was a relatively obvious asymmetry of the funnel plot (Figure 3), indicating some degree of publication bias. It is unlikely to be the publication bias caused by heterogeneity because of low heterogeneity, as displayed in Figure 2. Funnel plot assessment of IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872), and IL-10-819T/C (rs1800871) by visual inspection was symmetrical (Figure S9), indicating no significant publication bias.

FIGURE 3

Figure 3 Funnel plots of publication bias of IL-10-1082A/G(rs1800896) polymorphism.

A sensitivity analysis of the included studies was performed using leave-one-out forest plots according to a random- or fixed-effects model to investigate the sensitivity of our results. The sequential exclusion of studies had no significant impact on all of the above findings, indicating their robustness (Figures 4 and S10).

FIGURE 4

Figure 4 Sensitivity analysis: Leave-one-out meta-analysis of IL-10-1082A/G(rs1800896) polymorphism and PTB (AA+GG vs. AG).

TSA Analyses

To improve the robustness of our findings, we performed TSA on the over-dominant gene model (AA + GG vs. AG) of IL-10-1082A/G (rs1800896), the allele model T-allele vs. C-allele of IL-4-590C/T (rs2243250), C-allele vs. A-allele of IL-10-592A/C(rs1800872) and C-allele vs. T-allele of IL-10-819T/C (rs1800871) polymorphisms (Figures 5 and S11). TSA was performed in all gene models of each polymorphism, and similar results were found, so the above results were used as an example. The findings indicate that the cumulative Z-curve for the over-dominant gene model (AA + GG vs. AG) of IL-10-1082A/G (rs1800896) polymorphism crossed the trial sequential monitoring boundary (type I error 5%, Z score = 1.96) before reaching the required information size (TSA = 4376) (Figure 5). Similar patterns were observed in the allele model of IL-10-592A/C (rs1800872) and IL-10-819T/C (rs1800871) polymorphisms (TSA = 5142 and 34107, respectively) (Figure S11). However, the cumulative Z-curve for the allele model of IL-4-590C/T (rs2243250) polymorphism did not cross the trial sequential monitoring boundary (type I error 5%, Z score = 1.96) and did not reach the required information size (TSA = 10223) (Figure S11). As a result, the cumulative evidence was sufficient to support the conclusions for IL-10-1082A/G (rs1800896), IL-10-592A/C (rs1800872), and IL-10-819T/C (rs1800871) polymorphisms and PTB in the corresponding gene model. Furthermore, more research is required to confirm the true relationship between IL-4-590C/T (rs2243250) polymorphism and PTB in different gene models.

FIGURE 5

Figure 5 TSA between IL-10-1082A/G(rs1800896) polymorphism and PTB risk (AA+GG vs. AG).

eQTLs for SNPs

To determine whether the above SNPs show evidence as eQTLs in blood, we examined the “QTLbase” for eQTLs. The effective allele and p value for eQTLs of both four SNPs were statistically significant (p < 10⁻⁸) (Table S2).

Discussion

In this systematic review and meta-analysis, the over-dominant gene model AA+GG vs. AG revealed a significant association between IL-10-1082A/G (rs1800896) and PTB. Furthermore, due to an insufficient number of studies for meta-analysis, many other IL-4 and IL-10 gene polymorphisms were excluded from this study.

Human parturition is an inflammatory process, and when the inflammatory balance is disrupted, the formation of an inflammation cascade can lead to the occurrence of PTB (20). Anti-inflammatory cytokines are essential for pregnancy maintenance (21). Moreover, the anti-inflammatory cytokines IL-4 and IL-10 inhibit the production of pro-inflammatory cytokines, and decreased production of IL-4 and IL-10 may increase susceptibility to infections and lead to PTB (22). When comparing preterm and term placentas, the researchers discovered that term placentas had higher IL-4 and IL-10 (23). Furthermore, studies have indicated that the inflammatory response of PTB is genetically controlled (24). As a result, SNPs in IL-4 and IL-10 may result in lower anti-inflammatory cytokines, shifting the cytokine balance toward inflammation. Therefore, SNPs in IL-4 and IL-10 had been extensively studied as potential candidates for predicting PTB risk (10, 25).

In multifetal pregnancies, IL-4 C > T substitution at position -590 had been associated with increased production of IL-4 and susceptibility to PTB (26). Similarly, Heinzmann et al. and Belousova et al. found an association between IL-4-590C/T (rs2243250) gene polymorphism and the risk of PTB (8, 27). Lyubomirskaya et al., on the other hand, found no significant associations between IL-4 SNPs and risk or protective effect for PPROM in the population of Zaporizhzhia region of Ukraine (28). This meta-analysis found no significant association between IL-4-590C/T (rs2243250) gene polymorphism and PTB. However, IL-4 590 T allele frequency was disproportionately high in certain ethnic groups (29). Our race-based subgroup analysis also found no significant link between the IL-4-590C/T (rs2243250) gene polymorphism and PTB. This could be due to the small number of included studies and the considerable heterogeneity among studies, which requires to be confirmed by further research.

Another well-known anti-inflammatory cytokine that can suppress pro-inflammatory cytokines is IL-10 (10). Genetic variations in IL-10 may affect gene expression and thus its levels, resulting in an inflammatory imbalance in pregnant women and PTB (30). Polymorphisms in IL-10-592A/C (rs1800872), IL-10-819T/C (rs1800871) and IL-10-1082A/G (rs1800896) are common SNPs in PTB studies (9). However, studies on IL-10 gene polymorphisms and PTB have yielded contradictory results (9, 31–34). Menon et al. found that the difference in the interaction of IL-10 SNPs and PTB may have a genetic basis (35). Simultaneously, Moura et al., Stonek et al., and Mattar et al. suggested that IL-10-1082A/G (rs1800896) polymorphism was not a genetic marker for identifying women at increased risk of PTB (25, 36, 37). However, our findings demonstrate that IL-10-1082A/G (rs1800896) over-dominant gene model AA+GG vs. AG has a significant protective effect against PTB. Similarly, Kerk et al. revealed that IL-10-1082A/G (rs1800896) might play a role in severe PTB that lasts less than 29 weeks (38). No significant association was found between other gene models of IL-10-592A/C (rs1800872) and IL-10-819T/C (rs1800871) and PTB. Furthermore, no racial differences were found in the subgroup analysis.

Eskdale et al. found that the ability to secrete IL-10 varies in human according to the genetic composition of the IL-10 locus (39). Similarly, Suárez et al. found that IL-10 (mRNA and serum protein) levels showed remarkable inter-individual variations, which are genetically controlled by polymorphic variants at the gene promoter (40). There was a trend of higher IL-10 levels among individuals homozygous for the 1082 G genotype, but no difference in the expression of IL-10 between AG and AA genotype groups (40). The study also found that the concentration of IL-10 transcripts was quite dispersed (interquartile range = 13.85 and 18.87) in individuals with -1082 GG and GA genotypes, indicating that other genetic or environmental factors were also involved (40). In addition, in vitro and animal studies demonstrated that IL-10 could improve the outcome of LPS-induced PTB model (14, 15). Clinical studies revealed that lower IL-10 levels were significantly associated with PTB and had the potential to be used as a biomarker for predicting PTB (41). The eQTL assessment showed that the four SNPs (rs2243250, rs1800872, rs1800871 and rs1800896) in this study were significantly correlated with the expression of the corresponding genes (IL-4 and IL-10), which also proved that the protective effect of the over-dominant gene model AA+GG vs AG of IL-10-1082A/G (rs1800896) polymorphism on PTB may be due to the increase of IL-10 expression in blood. Overall, these findings suggested that the protective effect of AA+GG vs. AG of IL-10-1082A/G (rs1800896) genotype on PTB may be due to an increase in IL-10 levels. Furthermore, the difference between -1082 AG and AA genotype could be due to other unknown genetic or environmental factors.

Our study has some limitations as well. First, the number of studies involving IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872) and IL-10-819T/C (rs1800871) is small. Second, we only performed subgroup analysis on different countries but did not analyze among different races in countries, which may have obscured the correlation caused by race differences. Third, there was a risk of publication bias for IL-10-1082A/G (rs1800896), and data for some genotypes were obtained by calculating, which may differ from the original data. Finally, the absence of publication bias indicated a significant heterogeneity in the included studies of IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872), and IL-10-819T/C (rs1800871).

In conclusion, our findings suggest that the over-dominant gene model AA + GG vs. AG of IL-10-1082A/G (rs1800896) polymorphism may be associated with genetic susceptibility to PTB and may have a protective role. On the other hand, IL-4-590C/T (rs2243250), IL-10-592A/C (rs1800872) and IL-10-819T/C (rs1800871) polymorphisms are unrelated to PTB. Due to the limitations of included studies and the risk of publication bias, future studies in larger populations are required to confirm our findings.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author Contributions

X-LC conceived the idea for this meta-analysis. All authors (X-LC, X-YZ, N-XX, S-CC, C-MX) developed the methodology for the meta-analysis. The manuscript was drafted by X-LC, and revised by X-YZ, N-XX, S-CC, and C-MX. All authors contributed to the article and approved the submitted version.

Funding

This project is supported by the National Natural Science Foundation of China (NO.82171677).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors want to thank those authors who shared the full text of their articles. Thank “Home for Researchers (www.home-for-researchers.com)”.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.917383/full#supplementary-material

References

1. WHO: Recommended Definitions, Terminology and Format for Statistical Tables Related to the Perinatal Period and Use of a New Certificate for Cause of Perinatal Deaths. Modifications Recommended by FIGO as Amended October 14, 1976. Acta Obstet Gynecol Scand (1977) 56(3):247–53.

PubMed Abstract | Google Scholar

2. Chawanpaiboon S, Vogel JP, Moller AB, Lumbiganon P, Petzold M, Hogan D, et al. Global, Regional, and National Estimates of Levels of Preterm Birth in 2014: A Systematic Review and Modelling Analysis. Lancet Glob Health (2019) 7(1):e37–46. doi: 10.1016/s2214-109x(18)30451-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Liu L, Oza S, Hogan D, Chu Y, Perin J, Zhu J, et al. Global, Regional, and National Causes of Under-5 Mortality in 2000-15: An Updated Systematic Analysis With Implications for the Sustainable Development Goals. Lancet (2016) 388(10063):3027–35. doi: 10.1016/s0140-6736(16)31593-8

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Paules C, Pueyo V, Martí E, de Vilchez S, Burd I, Calvo P, et al. Threatened Preterm Labor is a Risk Factor for Impaired Cognitive Development in Early Childhood. Am J Obstet Gynecol (2017) 216(2):157.e1–7. doi: 10.1016/j.ajog.2016.10.022

CrossRef Full Text | Google Scholar

5. Romero R, Dey SK, Fisher SJ. Preterm Labor: One Syndrome, Many Causes. Science (2014) 345(6198):760–5. doi: 10.1126/science.1251816

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Couceiro J, Matos I, Mendes JJ, Baptista PV, Fernandes AR, Quintas A. Inflammatory Factors, Genetic Variants, and Predisposition for Preterm Birth. Clin Genet (2021) 100(4):357–67. doi: 10.1111/cge.14001

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Strizhakov AN, Belousova VS, Svitich OA. Clinical Significance of Toll-Like Receptors in the Pathogenesis of Preterm Labour. (2016). doi: 10.20953/1726-1678-2016-1-35-40

CrossRef Full Text | Google Scholar

8. Belousova VS, Svitich OA, Timokhina EV, Strizhakov AN, Bogomazova IM. Polymorphism of the IL-1 Beta, TNF, IL-1RA and IL-4 Cytokine Genes Significantly Increases the Risk of Preterm Birth. Biochem-Moscow (2019) 84(9):1040–6. doi: 10.1134/s0006297919090062

CrossRef Full Text | Google Scholar

9. Pandey M, Awasthi S, Singh U, Mahdi AA. Association of IL-10 Gene Polymorphism (-819c > T, -592c > A and -1082G > A) With Preterm Birth. Indian J Pediatr (2018) 85(2):93–101. doi: 10.1007/s12098-017-2496-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhu Q, Sun J, Chen Y. Preterm Birth and Single Nucleotide Polymorphisms in Cytokine Genes. Transl Pediatr (2014) 3(2):120–34. doi: 10.3978/j.issn.2224-4336.2014.03.02

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ho IC, Miaw SC. Regulation of IL-4 Expression in Immunity and Diseases. Adv Exp Med Biol (2016) 941:31–77. doi: 10.1007/978-94-024-0921-5_3

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bosurgi L, Cao YG, Cabeza-Cabrerizo M, Tucci A, Hughes LD, Kong Y, et al. Macrophage Function in Tissue Repair and Remodeling Requires IL-4 or IL-13 With Apoptotic Cells. Science (2017) 356(6342):1072–6. doi: 10.1126/science.aai8132

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ouyang W, Rutz S, Crellin NK, Valdez PA, Hymowitz SG. Regulation and Functions of the IL-10 Family of Cytokines in Inflammation and Disease. Annu Rev Immunol (2011) 29:71–109. doi: 10.1146/annurev-immunol-031210-101312

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Bryant AH, Spencer-Harty S, Owens SE, Jones RH, Thornton CA. Interleukin 4 and Interleukin 13 Downregulate the Lipopolysaccharide-Mediated Inflammatory Response by Human Gestation-Associated Tissues. Biol Reprod (2017) 96(3):576–86. doi: 10.1095/biolreprod.116.145680

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Terrone DA, Rinehart BK, Granger JP, Barrilleaux PS, Martin JN Jr., Bennett WA. Interleukin-10 Administration and Bacterial Endotoxin-Induced Preterm Birth in a Rat Model. Obstet Gynecol (2001) 98(3):476–80. doi: 10.1016/s0029-7844(01)01424-7

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JP, et al. The PRISMA Statement for Reporting Systematic Reviews and Meta-Analyses of Studies That Evaluate Healthcare Interventions: Explanation and Elaboration. Bmj (2009) 339:b2700. doi: 10.1136/bmj.b2700

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wells GA, Shea B, O'Connell J. The Newcastle-Ottawa Scale (NOS) for Assessing The Quality of Nonrandomised Studies in Meta-Analyses. (2014).

Google Scholar

18. . Trial Sequential Analysis (TSA) [Computer Program]. Version 0.9.5.10 Beta. The Copenhagen Trial Unit, Centre for Clinical Intervention Research (2021), The Capital Region, Copenhagen University Hospital–Rigshospitalet.

Google Scholar

19. Liubomyrska KS, Kamyshnyi OM, Krut YY. Association Between Single Nucleotide Polymorphism of Immunoregulatory Genes and Preterm Premature Rupture of Membranes in Preterm Labour. Pathologia (2018) 2):187–93. doi: 10.14739/2310-1237.2018.2.141400

CrossRef Full Text | Google Scholar

20. Boyle AK, Rinaldi SF, Norman JE, Stock SJ. Preterm Birth: Inflammation, Fetal Injury and Treatment Strategies. J Reprod Immunol (2017) 119:62–6. doi: 10.1016/j.jri.2016.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Engel SAM, Olshan AF, Savitz DA, Thorp J, Erichsen HC, Chanock SJ. Risk of Small-for-Gestational Age is Associated With Common Anti-Inflammatory Cytokine Polymorphisms. Epidemiology (2005) 16(4):478–86. doi: 10.1097/01.ede.0000164535.36412.6b

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Escobar-Arregoces F, Latorre-Uriza C, Velosa-Porras J, Roa-Molina N, Ruiz AJ, Silva J, et al. Inflamatory Response in Pregnant Women With High Risk of Preterm Delivery and its Relationship With Periodontal Disease: A Pilot Study. Acta Odontol Latinoam (2018) 31(1):53–7. doi: 10.1155/2018/7027683

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tranchot-Diallo J, Gras G, Parnet-Mathieu F, Benveniste O, Marcé D, Roques P, et al. Modulations of Cytokine Expression in Pregnant Women. Am J Reprod Immunol (1997) 37(3):215–26. doi: 10.1111/j.1600-0897.1997.tb00218.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lyubomirskaya ES, Kamyshnyi AM, Krut YY, Smiianov VA, Fedoniuk LY, Romanyuk LB, et al. SNPs and Transcriptional Activity of Genes of Innate and Adaptive Immunity at the Maternal-Fetal Interface in Woman With Preterm Labour, Associated With Preterm Premature Rupture of Membranes. Wiad Lek (2020) 73(1):25–30.

PubMed Abstract | Google Scholar

25. Mattar R, de Souza E, Daher S. Preterm Delivery and Cytokine Gene Polymorphisms. J Reprod Med (2006) 51(4):317–20.

PubMed Abstract | Google Scholar

26. Kalish RB, Vardhana S, Gupta M, Perni SC, Witkin SS. Interleukin-4 and -10 Gene Polymorphisms and Spontaneous Preterm Birth in Multifetal Gestations. Am J Obstet Gynecol (2004) 190(3):702–6. doi: 10.1016/j.ajog.2003.09.066

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Heinzmann A, Mailaparambil B, Mingirulli N, Krueger M. Association of Interleukin-13/-4 and Toll-Like Receptor 10 With Preterm Births. Neonatology (2009) 96(3):175–81. doi: 10.1159/000210091

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lyubomirskaya K, Krut Y, Sergeyeva L, Khmil S, Lototska O, Petrenko N, et al. Preterm Premature Rupture of Membranes: Prediction of Risks in Women of Zaporizhzhia Region of Ukraine. Pol Merkur Lekarski (2020) 48(288):399–405.

PubMed Abstract | Google Scholar

29. Nguyen DP, Genc M, Vardhana S, Babula O, Onderdonk A, Witkin SS. Ethnic Differences of Polymorphisms in Cytokine and Innate Immune System Genes in Pregnant Women. Obstet Gynecol (2004) 104(2):293–300. doi: 10.1097/01.AOG.0000133486.85400.5e

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Speer EM, Gentile DA, Zeevi A, Pillage G, Huo DZ, Skoner DP. Role of Single Nucleotide Polymorphisms of Cytokine Genes in Spontaneous Preterm Delivery. Hum Immunol (2006) 67(11):915–23. doi: 10.1016/j.humimm.2006.08.291

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Han SH, Lee NR, Kim HJ, Kang YD, Kim JS, Park JW, et al. Association Between the IL-6, IL-10, and Tnfα Gene Polymorphisms and Preterm-Birth in Korean Women. Genes Genomics (2020) 42(7):743–50. doi: 10.1007/s13258-020-00946-4

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ramos BRD, Mendes ND, Tanikawa AA, Amador MAT, dos Santos NPC, dos Santos SEB, et al. Ancestry Informative Markers and Selected Single Nucleotide Polymorphisms in Immunoregulatory Genes on Preterm Labor and Preterm Premature Rupture of Membranes: A Case Control Study [Article]. BMC Pregnancy Childbirth (2016) 16:11. doi: 10.1186/s12884-016-0823-1

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Nuk M, Orendi K, Rosenberger S, Griesbacher A, Holzapfel-Bauer M, Resch B, et al. Genetic Variations in Fetal and Maternal Tumor Necrosis Factor-α and Interleukin 10: Is There an Association With Preterm Birth or Periventricular Leucomalacia? J Perinatol (2012) 32(1):27–32. doi: 10.1038/jp.2011.37

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Annells MF, Hart PH, Mullighan CG, Heatley SL, Robinson JS, Bardy P, et al. Interleukins-1, -4, -6, -10, Tumor Necrosis Factor, Transforming Growth Factor-Beta, FAS, and Mannose-Binding Protein C Gene Polymorphisms in Australian Women: Risk of Preterm Birth. Am J Obstet Gynecol (2004) 191(6):2056–67. doi: 10.1016/j.ajog.2004.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Menon R, Fortunato SJ, Edwards DRV, Williams SM. Association of Genetic Variants, Ethnicity and Preterm Birth With Amniotic Fluid Cytokine Concentrations. Ann Hum Genet (2010) 74:165–83. doi: 10.1111/j.1469-1809.2010.00562.x

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Moura E, Mattar R, de Souza E, Torloni MR, Goncalves-Primo A, Daher S. Inflammatory Cytokine Gene Polymorphisms and Spontaneous Preterm Birth [Article]. J Reprod Immunol (2009) 80(1-2):115–21. doi: 10.1016/j.jri.2008.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Stonek F, Metzenbauer M, Hafner E, Philipp K, Tempfer C. Interleukin-10-1082 G/A Promoter Polymorphism and Pregnancy Complications: Results of a Prospective Cohort Study in 1,616 Pregnant Women. Acta Obstet Gynecol Scand (2008) 87(4):430–3. doi: 10.1080/00016340801995657

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kerk J, Dördelmann M, Bartels DB, Brinkhaus MJ, Dammann CE, Dörk T, et al. Multiplex Measurement of Cytokine/Receptor Gene Polymorphisms and Interaction Between Interleukin-10 (-1082) Genotype and Chorioamnionitis in Extreme Preterm Delivery. J Soc Gynecol Investig (2006) 13(5):350–6. doi: 10.1016/j.jsgi.2006.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Eskdale J, Gallagher G, Verweij CL, Keijsers V, Westendorp RG, Huizinga TW. Interleukin 10 Secretion in Relation to Human IL-10 Locus Haplotypes. Proc Natl Acad Sci USA (1998) 95(16):9465–70. doi: 10.1073/pnas.95.16.9465

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Suarez A, Castro P, Alonso R, Mozo L, Gutierrez C. Interindividual Variations in Constitutive Interleukin-10 Messenger RNA and Protein Levels and Their Association With Genetic Polymorphisms. Transplantation (2003) 75(5):711–7. doi: 10.1097/01.TP.0000055216.19866.9A

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Wommack JC, Ruiz RJ, Marti CN, Stowe RP, Brown CE, Murphey C. Interleukin-10 Predicts Preterm Birth in Acculturated Hispanics. Biol Res Nurs (2013) 15(1):78–85. doi: 10.1177/1099800411416225

PubMed Abstract | CrossRef Full Text | Google Scholar