Introduction

Breast cancer is a complex multiple process influenced by multiple factors. It is considered that specific gene polymorphisms have effects on gene transcription, mRNA stability and protein activity [1]. Multi-functional cytokines involved in the process of inflammatory and immunological responses are closely associated with the pathogenesis of autoimmune and malignant diseases, making them potential risks for breast cancer [24]. The interleukin 1 gene family located on chromosome 2q14.2 includes three members: IL-1α, IL-1β and IL-1 receptor antagonist (IL-RA) encoded by IL-1RN [5]. Interleukin 1α and IL-1β are potent proinflammatory cytokines, whereas IL-RA is an anti-inflammatory cytokine [6]. Interleukin 1β can be produced by various cells and it modifies the process of host response to microbial invasion, tissue injury and inflammation [7].

The interleukin 1β gene has three potentially functional SNPs: –31 (rs1143627, C>T), –511 (rs16944, C>T) in the promoter region and +3954 (rs1143634, C>T) in exon 5 [8, 9]. So far, many studies have been conducted to assess the relations between the three SNPs (rs16944, rs1143634 and rs1143627) in IL-1β and breast cancer risk [1020]. However, the results remain conflicting. Therefore, we performed this meta-analysis in order to obtain a more precise evaluation of these links.

Material and methods

Literature search

Relevant studies published before June 1st, 2018 were identified through a search in PubMed, Embase, and Web of Science using a combination of the following terms: (“polymorphism” or “SNP” or “single nucleotide polymorphisms”), (“breast cancer” or “breast carcinoma” or “breast tumor”) and (“Interleukin-1” or “IL-1β” or “Interleukin-1 beta”). The references from the eligible articles or textbooks were also manually searched by us to find additional potential sources.

Inclusion and exclusion criteria

The criteria for studies in our meta-analysis are as follows: (a) studies concentrated on relations between IL-1β polymorphisms and breast cancer risk; (b) providing sufficient data for the frequencies of alleles and genotypes; (c) published in English. Studies were excluded when: (a) they were not case-control studies; (b) did not supply complete and essential information; (c) they were meta-analyses, reviews, or editorial articles.

Data extraction

Data were extracted from each publication independently by two authors (Wang and Yuan) based on the inclusion criteria mentioned above. For each study, the data were collected as follows: the first author, year of publication, country of origin, ethnicity, menopausal state, numbers of patients and controls, source of controls, mutation detection methods, genotyping methods, minor allele frequency (MAF), allele and genotype frequencies and the evidence of Hardy-Weinberg equilibrium (HWE) in controls. Disputes were settled by consulting with a third author if disagreements occured.

Statistical analysis

In order to assess the strength of associations between IL-1β gene polymorphisms and breast cancer susceptibility under five genetic models which include the allele model, dominant model, recessive model, homozygous model and the heterozygous model, crude odds ratios (ORs) with their corresponding 95% confidence intervals (CIs) were adopted [2124]. Hardy-Weinberg equilibrium in control groups was estimated using the χ2 test. The pooled OR’s statistical significance was verified using the Z test with a two-tailed p < 0.05 which is regarded as statistically significant. Between-study variations and heterogeneities were evaluated using either Cochran’s Q-statistic or I2 test. When the result is a p-value < 0.05 or I2 > 50%, it indicates the existence of heterogeneity among studies; only in this circumstance was the random effects model (DerSimonian Laird method) used. If not, the fixed effects model (Mantel-Haenszel method) was performed. In order to explore sources of heterogeneity, we performed subgroup analysis with reference to ethnicity, menopausal state, source of controls and genotyping methods. Additionally, to investigate the sensitivity, we removed each study in turn to evaluate the quality and consistency of results. We used Begg’s funnel plot and Egger’s linear regression test to detect publication biases. All analyses were conducted with STATA version 12.0.

Results

Eligible studies

Four studies with controls deviating from HWE were excluded [12, 15, 18, 20]. Ten case-control studies were conducted by us to assess associations between IL-1β polymorphisms and breast cancer risk. As shown in Tables I and II, 6 studies were eligible for rs16944(C>T) including 2454 cases and 2720 controls [11, 12, 14, 16, 18, 19]. Details are as follows: 1) Caucasians, Asians, and Africans were investigated in 4, 2 and 1 studies respectively, 3 of which were associated with menopause. 2) Only 1 study was based on PB regarding the source of controls. 3) Genotype methods included TaqMan, PCR, sequencing and MALDI-TOF. For rs1143634(C>T), it was the same as above, but only two studies were qualified, which contained 996 cases and 620 controls [13, 14]. Finaly, for rs1143627(C>T), just 2413 cases and 2438 controls were involved in five studies [10, 1517, 19]. Genotype distributions in controls of all studies were in accordance with HWE.

Table I

Characteristics of case-control studies included in the meta-analysis

SNPs/First authorYearRacial descentCountryMenopausal stateSource of controlsGenotype methodsCasesControls
rs16944:
 Smith2004CaucasianUKNPBTaqman141261
 Hefler2005CaucasianGermanyNHBPCR269227
 Liu2006AsianChinaNHBPCR365631
 Balasubramanian2006CaucasianUKNHBPCR703489
 Pooja2012AsianIndiaPre-menopausalHBSequence107200
 Pooja2012AsianIndiaPost-menopausalHBSequence93200
 Gong2013CaucasianAmericanPre-menopausalHBMALDL-TOF185163
 Gong2013CaucasianAmericanPost-menopausalHBMALDL-TOF141146
 Gong2013AfricanAmericanPre-menopausalHBMALDL-TOF237195
 Gong2013AfricanAmericanPost-menopausalHBMALDL-TOF213208
 Zuo2018AsianChinaNHBSequenom MassARRAY530628
rs1143634:
 Snoussi2005AfricanTunisiaNHBPCR305200
 Hefler2005CaucasianGermanyNHBPCR269227
 Balasubramanian2006CaucasianUKNHBPCR691420
 Pooja2012AsianIndiaPre-menopausalHBSequence107200
 Pooja2012AsianIndiaPost-menopausalHBSequence93200
 Pooja2012AsianIndiaNHBSequence200200
rs1143627:
 Ito2002AsianJapanNHBPCR-CTPP227185
 Liu2006AsianChinaNHBPCR365631
 Lee2006AsianKoreaPre-menopausalPBPCR-CTPP353290
 Lee2006AsianKoreaPost-menopausalPBPCR-CTPP206215
 Lee2006AsianKoreaNPBPCR-CTPP559505
 Akisik2007AsianTurkeyNNPCR-RFLP126110
 Gong2013CaucasianAmericanPre-menopausalHBMALDL-TOF186166
 Gong2013CaucasianAmericanPost-menopausalHBMALDL-TOF142146
 Gong2013AfricanAmericanPre-menopausalHBMALDL-TOF239195
 Gong2013AfricanAmericanPost-menopausalHBMALDL-TOF216210
 Zuo2018AsianChinaNHBSequenom MassARRAY530628
Table II

Characteristics of case–control studies included in the meta-analysis

SNPs/First authorGenetype distributionHWE test
CasesControls
CTCCCTTTMAFCTCCCTTTMAF
rs16944:
 Smith187956067140.3430921387135390.410.25
 Hefler362176124114310.3328716788111280.370.61
 Liu358372941701010.496995631973051290.450.58
 Balasubramanian972434339294700.31670308232206510.310.61
 Pooja671471243520.311442562594810.360.78
 Pooja701161834410.381442562594810.360.78
 Gong2541168878190.312151117271200.340.71
 Gong1601224962300.22205877163120.290.71
 Gong22125356109720.471782124490610.460.33
 Gong18823834120590.441932234995640.460.24
 Zuo5565041602361340.486286281423441420.50.02
rs1143634:
 Snoussi428182157114340.293069412066140.240.24
 Hefler41512315997130.233371171199990.260.04
 Balasubramanian1062320410242390.23629211231167220.250.24
 Pooja17836752840.173643617612120.09<0.01
 Pooja16026761840.133643617612120.09<0.01
 Pooja340601474670.153643617612120.09<0.01
rs1143627:
 Ito21923558103660.481552152899580.420.18
 Liu379351102175880.485796831333131850.460.98
 Lee344362961521050.4928529570145750.490.99
 Lee20920351107480.5821121943125470.490.02
 Lee5535651472591530.494965141132701220.490.12
 Akisik991531863450.29981222156330.450.75
 Gong1192532079870.321132192171740.340.54
 Gong1221623062500.43872051263710.290.71
 Gong28219684114410.412381527686330.390.31
 Gong26616681104310.382461747792410.410.16
 Zuo5105501362381560.486306261443421420.490.03

Meta-analysis

The critical results of the current meta-analysis are described in Table III. Based on data, three common polymorphisms in IL-1β (rs16944, rs1143634, rs1143627) were not significantly associated with breast cancer risk in all genetic models. To go a step further, the data were stratified into different subgroups according to ethnicity, menopausal state, source of controls and genotyping methods. In terms of analysis by ethnic subgroup, TT genotype of rs16944 polymorphism represented a higher risk of breast cancer compared to CT genotype and CC + TT genotype and this was only significant in Asian people (TT vs. CC + CT: 1.229, 95% CI: 1.063–1.422, p = 0.005, Figure 1; TT vs. CT: 1.211, 95% CI: 1.057–1.388, p = 0.006, Figure 2), not in Caucasian and African populations. Similarly, it is only in Asian people that a significantly decreased breast cancer risk was found in the dominant model (CT + TT vs. CC: OR = 0.944, 95% CI: 0.897–0.994, p = 0.027, Figure 3) for the rs1143627 polymorphism. With patients being stratified according to the menopausal state, we noted that the rs1143627 polymorphism correlated with reduced breast cancer risk among post-menopausal women in three genotype models: the allele model (T vs. C: 0.859, 95% CI: 0.753–0.98, p = 0.024), the recessive model (TT vs. CC + CT: 0.727, 95% CI: 0.576–0.918, p = 0.007, Figure 4), and the homozygous model (TT vs. CC: 0.743, 95% CI: 0.626–0.882, p = 0.001). No significant associations between IL-1β polymorphism and breast cancer risk were found in other stratified analyses based on source of controls and genotyping methods.

Table III

Summary ORs and 95% CI of interleukin-1β polymorphisms and breast cancer risk

VariableAllele modelDominant modelRecessive modelHomozygous modelHeterozygous model
OR(95% CI)P-valueaOR(95% CI)P-valueaOR(95% CI)P-valueaOR(95% CI)P-valueaOR(95% CI)P-valuea
rs16944:
 Overall1.0170.944–1.0960.6531.0010.945–1.060.6091.0710.915–1.2540.3931.0360.898–1.1960.6261.0720.944–1.2180.283
Ethnicity:
 Caucasion0.9860.837–1.160.8620.970.863–1.090.8131.0290.703–1.5070.0820.9880.657–1.4870.9561.0390.847–1.2730.714
 Asian1.0710.98–1.1720.131.010.927–1.1010.4361.2291.063–1.4220.0051.0830.884–1.3270.4431.2171.06–1.3960.005
 African1.0170.927–1.1040.7911.0430.938–1.1610.9790.9370.763–1.150.5331.0410.889–1.2190.6150.9010.743–1.0920.289
rs1143634:
 Overall1.0720.783–1.4660.6651.0360.777–1.3820.8071.270.863–1.870.2261.2770.758–2.1510.3591.2430.86–17960.247
rs1143627:
 Overall0.9590.903–1.0180.1710.9560.912–1.0020.0580.9660.855–1.0920.5810.9180.823–1.0250.1281.0170.92–1.1240.747
Ethnicity:
 Caucasion0.9180.726–1.1590.4720.940.792–1.1150.4750.8790.611–1.2650.4880.8790.617–1.2520.4730.9430.779–1.1430.553
 Asian0.9650.899–1.0340.3120.9440.897–0.9940.0271.0230.881–1.1870.7680.9290.812–1.0640.2891.0840.963–1.220.181
 African0.9880.873–1.1180.8491.0240.924–1.1350.6510.8660.632–1.1860.370.9320.689–1.2590.6450.8450.639–1.1180.239
Menopausal:
 Pre-menopausal1.0250.958–1.0980.4711.0050.952–1.0610.8671.0850.925–1.2720.3261.0370.928-1.1590.5191.080.932–1.2510.308
 Post-menopausal0.8590.753–0.980.0240.9130.787–1.0580.2250.7280.575–0.9220.0080.7430.626-0.8820.0010.8120.653–1.0110.062

a Two-side χ2 test, P < 0.05 was considered statistically significant.

Figure 1

OR of breast cancer in different ethnicities associated with rs16944 in IL-1β gene for the TT genotype compared with the CC + CT genotype

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g001_min.jpg
Figure 2

OR of breast cancer in different ethnicities associated with rs16944 in IL-1β gene for the TT genotype compared with the CT genotype

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g002_min.jpg
Figure 3

OR of breast cancer in different ethnicities associated with rs1143627 in IL-1β gene for the CT + TT genotype compared with the CC genotype

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g003_min.jpg
Figure 4

OR of breast cancer in different menopausal state associated with rs1143627 in IL-1β gene for the TT genotype compared with the CC + CT genotype

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g004_min.jpg

Publication bias

Egger’s test and Begg’s funnel plot were used so that we could evaluate the potential publication bias of the studied literature. No obvious evidence of publication bias was detected in IL-1β (rs16944, rs1143634, rs1143627) (Figure 5). Also, good results were obtained in the sensitivity analysis (Figure 6).

Figure 5

Begg’s funnel plot analysis of publication bias. A – rs16944 in TT vs. CC + CT model. B – rs16944 in TT vs. CT model. C – rs1143627 in CT + TT vs. CC model. D – rs1143627 in TT vs. CC + CT model

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g005_min.jpg
Figure 6

Sensitivity analysis. A – rs16944 in TT vs. CC + CT model. B – rs16944 in TT vs. CT model. C – rs1143627 in CT + TT vs. CC model. D – rs1143627 in TT vs. CC + CT model

https://www.archivesofmedicalscience.com/f/fulltexts/99839/AMS-18-1-99839-g006_min.jpg

Discussion

It is believed that cytokines are strongly connected to cancer pathogenesis accompanying increasing evidence, which suggests that they participate in tumor initiation, growth and metastasis [2]. Numerous studies about cytokine gene polymorphisms have been conducted to investigate their relations with many inflammatory and neoplastic diseases [14]. Emerging studies have reported the associations between IL-1β polymorphisms and breast cancer risk because of IL-1β’s crucial importance in breast cancer development. The size of samples in a single study was relatively small and the controls of some studies deviated from HWE, so the results were controversial. For example, Ito et al. first reported that rs1143627 was significantly associated with breast cancer risk [10] and another case-control study by Liu et al. verified this conclusion in the Chinese population [16]. However, Lee et al. and Akisik and Dalay did not find significant differences [15, 17]. Studies by Snoussi et al. and Pooja et al. revealed highly significant associations between rs1143634 and the aggressive phenotype of breast cancer [13, 18]. On the other hand, Hefler et al. and Balasubramanian et al. did not find significant differences [12, 14]. For rs16944, some studies indicated that rs16944 genotype reduced the risk of breast cancer, while others showed the opposite results.

Liu et al. conducted a meta-analysis to investigate the relations between three polymorphisms in IL-1β gene and the risk of breast cancer [25]. The variant genotype of rs1143627 was found to be associated with a significantly increased breast cancer risk while the polymorphisms rs16944 and rs1143634 did not represent any associations with breast cancer risk, which was inconsistent with our results. As Wang et al. mentioned [26], the data reported by Liu et al. conflicted with the data from some previous publications, so the results provided by Liu et al. were untrustworthy. Given this situation, we conducted an updated meta-analysis to investigate the associations between three polymorphisms in the IL-1β gene and breast cancer risk.

Some potential limitations in the present meta-analysis should be taken into consideration. First of all, compared to the one included in previous studies, this sample size was larger, but it was still relatively small for some SNPs and stratified analyses. Next, we could not conduct haplotype analysis and linkage disequilibrium, because more details about personal information regarding genotypes of the SNPs (rs16944, rs1143634 and rs1143627) in IL-1β gene were unavailable. Furthermore, it is difficult for us to evaluate potential interactions between gene-environment, gene-gene and multiple polymorphic loci in the same gene. Regardless of these limitations, our present meta-analysis includes a much larger number of eligible studies and a stratified analysis.

In conclusion, our current meta-analysis suggests that the rs16944 polymorphism is significantly associated with increased risk of breast cancer among the Asian population, while the rs1143627 polymorphism reduces the breast cancer risk in post-menopausal women. In the future, we need further large-scale and rigorous studies to validate these findings.