Skip to main content

The interactions between interleukin-1 family genes: IL1A, IL1B, IL1RN, and obesity parameters



Obesity has been recognized as a worldwide growing problem, producing many pathologies including the promotion of “proinflammatory state.” The etiology of human obesity is still only partially understood; however, the genetic background has been proved. Its nature is complex, and currently, it appears that the combined effects of the interactions among multiple genes should receive more attention. Due to the fact that obesity promotes proinflammatory conditions, in this study, we investigated the genetic polymorphism of IL-1 family genes in healthy people with normal and elevated body mass index (BMI) and fat %.


The single-nucleotide polymorphisms (SNPs) within the IL1A -889C > T (rs1800587), IL1B + 3954 T > C (rs1143634), and IL1RN -87G > A (rs2234677) genes alone were associated neither with BMI nor fat % values in tested group. The associations between SNP–SNP interaction and BMI for the IL1B × IL1RN interactions were significant for dominant model (p = 0.02) and codominant model (p = 0.03). The same SNP-SNP interaction (IL1B × IL1RN) was associated also with fat % for codominant (p = 0.01) and recessive (p = 0.002) models.


This study further confirmed that IL-1 family genes are involved in genetic background of obesity. It has been shown that interaction IL1B × IL1RN was associated with both BMI and fat % with rare T allele protecting form higher values. Thus, even if certain polymorphisms in single genes of IL-1 family cannot be defined as related to obesity in examined population, the genetic interrelationships should be analyzed.

Peer Review reports


In 20th century, obesity has been understood as a health problem and has reached epidemic proportions globally. According to the World Health Organization (WHO) data, it contributes to at least 2.8 million deaths each year. Once associated with high-income countries, obesity is now also prevalent in low- and middle-income countries, among people of any age [1]. From 1975 to 2016, an alarming, over fourfold growth (from 4 to 18%) in the prevalence of overweight or obesity was reported also in children and adolescents up to 19 years old. In 2017, over 4 million overweight or obese people died [2], and for the first time in human history, it was shown that the number of obese people worldwide exceeded underweight ones [3]. Thus, obesity represents a serious burden for public health.

Currently, the understanding of obesity is wider than only overweight and involves more aspects that make obesity a complex disease with multiple etiology and background. According to Rajender et al. (2014), human obesity arises from the combined effects of the interactions among multiple genes, environmental factors, and behavior. Thus, the management and prevention of this problem is especially challenging [4].

The etiology of obesity is complex and still not completely understood, but the role of environmental and genetic factors, as well as individual habits are widely accepted. The genetic background of obesity was proposed years ago, and it has been shown that the heritability of body mass index (BMI) ranges from 40 to 70% [5], and multiple genes and associations are involved. BMI is the most frequently used proxy for body composition; however, it is known that this value underestimates adiposity (body fat percentage). Since 1980, the American Army has implemented an Army Weight Control Program (AWCP), which originally was based on the soldier’s BMI but now body fat percentage (fat standard) is considered more important. Both values are tracked for soldiers and limits are implemented for military services [6].

Obesity is understood not only as abnormally high body weight (and so that high BMI) but also as a disease resulting in abnormal metabolic and immune functions. It has been shown that in obese individuals, the proinflammatory state, expressed as chronic, low-grade inflammation, is created, and may result in numerous inflammatory conditions [7,8,9]. It has been proved that the adipose tissue acts as a potent endocrine organ and a source of cytokines [7, 10,11,12], as indicated by increased concentrations of IL-6 and IL-1 in obese patients [13, 14].

The IL-1 is a cytokine family that has been shown to play a key role in multiple innate and acquired immune processes, including inflammation, fibrosis, autoimmune responses, vascular diseases, neoplasia, and others [7, 15, 16]. The IL-1 family consists of 11 members, including 7 with proinflammatory activity (IL-1α, IL-1β, IL-18, IL-33, IL-36α, IL-36β, and IL-36γ) and 4 with antagonistic (IL-1Ra, IL-36Ra, and IL-38) or anti-inflammatory (IL-37) properties. The IL-1 receptor family consists of 10 members: IL-1R1, IL-1R2, IL-1R accessory protein (IL-1RAcP), IL-18Rα, IL-18Rβ, ST2 (IL-33R), IL-36R (previously IL-1Rrp2), single-Ig IL-1R-related molecule (SIGIRR or TIR8), three-Ig-domain-containing IL-1R related 2 (TIGIRR-2 or IL-1RAPLI), and TIGIRR-1 (IL-1RAPL2) [7, 15, 17]. The balance between agonists and antagonists has been reported as being crucial for numerous inflammatory diseases, and the role of IL-1α, IL-1β, and IL-1Ra has been discussed the most frequently [7, 16]. The IL-1 family (except IL-18 and IL-33) is encoded as a cluster located in chromosome 2 (including IL1A, IL1B, and IL1RN) in locus 2q13-21 with three polymorphic genes [7, 16, 18]. Genetic polymorphisms of these genes have been indicated as potentially important in inflammatory diseases, including inflammatory bowel diseases, alopecia areata, psoriasis, lupus erythematosus, Alzheimer’s disease, joint degeneration, and metabolic syndrome, as well as obesity [7, 19,20,21,22,23].

Taking into account the role of IL-1 cytokines and receptors in the inflammatory response and the fact that proinflammatory/anti-inflammatory balance is clearly dysregulated in obese people, we hypothesized that variations and associations among IL-1 family genes may be important as promoting obesity. Thus, the aim of this study was to investigate genetic polymorphisms and interrelationships among genes of IL-1 family in healthy people (without inflammatory disorders) in the context of commonly used obesity measurements: BMI and the percentage of fat mass.



The study population has been selected to maintain maximal homogeneity and to eliminate as many environmental factors as possible. Therefore we have selected 101 military professionals, which lived together, ate the same food, and had very similar daily schedule. For each analysis, the participants were divided into 2 groups differentiated by either BMI or fat % (Table 1). The groups based on BMI involved 76 control (CONBMI) individuals with BMI values between 20.0 and 25.0 and 25 overweight (OVERBMI) ones, with BMI > 25.0. The groups based on fat % were composed of 87 control (CONFat) subjects with fat % below 20.0% and 14 overweight (OVERFat) ones with fat % over 20.0%.

Table 1 Body parameters (expressed as mean ± SD) in the analyzed groups

Genetic analyses for the groups

The measured frequencies of the IL1A rs1800587, IL1B rs1143634, and IL1RN rs2234677 genotypes did not differ significantly from the Hardy-Weinberg equilibrium expectations in the OVERBMI and OVERFat groups (p values 0.55–1 and 0.35–1, respectively). For the CONBMI and CONFat groups, genotype IL1B rs1143634 and IL1RN rs2234677 frequencies were not significantly different from the Hardy-Weinberg equilibrium expectations (p values ranged 0.28–0.56 and 0.19–0.59, respectively), but for the IL1A rs1800587 genotype, where the p value was 0.01, expectations were fulfilled neither for CONBMI nor CONFat group. For whole group, consisting of either CONBMI and OVERBMI, or CONFat and OVERFat, frequencies of the IL1B rs1143634 and IL1RN rs2234677 genotypes did not differ significantly from the Hardy-Weinberg equilibrium expectations (p values ranged from 0.32 to 0.81), and for the IL1A rs1800587 genotype, the expectations were not fulfilled (p = 0.003) (Table 2).

Table 2 Frequencies for ILA, ILB, and ILRN genotypes in OVER, CON, and ALL groups, probabilities that the genotype frequencies do not differ from Hardy-Weinberg expectations (HWE) and minor allele frequency (MAF)

Table 3 summarizes the results of the association analysis of SNPs in the IL1A, IL1B, and IL1RN genes and BMI values. IL1A polymorphism rs1800587 was analyzed only in codominant model due to the lack of AA genotype. The IL1A rs1800587, IL1RN rs2234677, and IL1B rs1143634 genotypes were not related with BMI value (Table 3).

Table 3 Analysis of IL1A (rs1800587), IL1B (rs1143634), IL1RN (rs2234677) polymorphisms in relation to BMI. Due to the lack of AA genotype, only codominant model was analyzed (rs1800587). OR—odds ratio; CI—confidence interval; AIC—the Akaike information criterion

Gene–Gene interactions IL1B rs1143634 × IL1A rs1800587, IL1RN rs2234677 × IL1A rs1800587 and IL1B rs1143634 × IL1RN rs2234677 (only pairwise interactions were considered) were investigated for the same genetic models as for single-gene analyses (except the overdominant model). The frequency of the remaining genotypes did not differ significantly between groups. No association was found between SNP–SNP interaction and BMI for the IL1B rs1143634 × IL1A rs1800587 and IL1RN rs2234677 × IL1A rs1800587 interactions in the dominant and the codominant models, and the recessive model was not applied for these interactions.

A significant association was found between SNP–SNP interaction and BMI for the IL1B rs1143634 × IL1RN rs2234677 interaction in the dominant model, p = 0.02 and the codominant model, p = 0.03. For both models, the odds of being OVERBMI for TC x GG was about 16 times lower than for CCxGG (Fisher test, p = 0.005). No association was found between SNP–SNP interaction and BMI for the IL1B rs1143634 × IL1RN rs2234677 interaction in the recessive model (Table 4). However, the odds of being OVERBMI for TC x GG-AG was 3.6 times lower than for CCxGG-AG (Fisher test, p = 0.025).

Table 4 Analysis of IL1B (rs1143634) × IL1RN (rs2234677) interaction in relation to BMI (dominant model). OR—odds ratio; CI—confidence interval, NA – not applicable

The IL1A rs1800587, IL1RN rs2234677, and IL1B rs1143634 genotypes were not related with fat %. Pairwise gene-gene interactions with fat % were investigated for the same genetic models as for BMI in single gene analyses (except overdominant model). A significant association was found between the SNP-SNP interaction and OVERFat for the IL1B rs1143634 x IL1RN rs2234677 interaction in the codominant model, p = 0.01 and the recessive model, p = 0.002 (Table 5).

Table 5 Analysis of the IL1B(rs1143634) x IL1RN (rs2234677) interaction in relation to fat % (codominant model. OR—odds ratio, CI—confidence intervals, NA—not applicable

The odds of being OVERFat for T/C x G/G (ILB x IL RN) was over 7 times lower than for reference combination C/C x G/G (ILB x ILRN). Similarly, the odd of being OVERFat for C/C x A/G (ILB x ILRN) was 1.5 times lower than for C/C x G/G (ILB x ILRN). All individuals representming T/C x A/G and C/C x A/A (ILB x ILRN) were in CONFat group (Fisher test, p = 0.09). The odds of being OVERFat for T/C x G/G-A/G (ILB x ILRN) was 10 times lower than for reference combination C/C x G/G - A/G (ILB x ILRN). All participants representing C/C x A/A (ILB x ILRN) were in CONFat group (Table 5, Fisher test, p = 0.01).

The most pronounced differences in the frequencies of genotype occurring in either the OVERBMI and CONBMI or OVERFat and CONFat groups were observed for the CC/GG genotype (IL1B rs1143634/ IL1RN rs2234677), which occurred in 40% of the OVERBMI group and 42.86% of the OVERFat group and in 19,74% of the CONBMI group and 21.84% of the CONFat group. In contrast, the CC/AG genotype (IL1B rs1143634/ IL1RN rs2234677) occurred only in 4% of the OVERBMI group and 7.14% of the OVERFat group and 31.59% of the CONBMI group and 27.59% of the CONFat group (Figs. 1 and 2).

Fig. 1
figure 1

Respective genotype percentages in the examined BMI groups

Fig. 2
figure 2

Respective genotype percentages in the examined groups fat %


Obesity is characterized by an abnormal increase in the total amount of triglyceride stored in adipose tissue, resulting from a chronic positive imbalance between energy intake and expenditure [24, 25]. In clinical practice, it is the most often measured as body mass index (BMI), which gives a surrogate measure of overall obesity. According to World Health Organization (WHO) recommendations, body mass index (BMI) higher than 25 is interpreted as indicating overweight, and BMI over 30 indicates obesity [1]; however, it still can underestimate adiposity.

Our study involved the unique, homogenous groups. All participants were students of Military University, so the environmental factors, including the diet, daily schedules, and activities have been the same for several months before the study was conducted. The only factors that distinguished the students for groups were overweighting parameters: BMI (although not the most accurate measurement, but still widely used) and fat %, representing better measure of adiposity. This design of the study allows to avoid the environmental influence and seems the best for the investigation of genetic background of overweighting. The OVERBM group was bigger than OVERFat which might have been interpreted as suggestion that higher BMI resulted from higher muscle mass. Such phenomenon can be observed in the athletes practicing certain disciplines, but it also results in a very low value of fat%. None of the soldiers examined in our study had high BMI and very low fat%, so both parameters indicated overweighting due to fat content, BMI can understood as more general and fat% as more precise measure of adiposity.

Family studies have previously estimated parent–offspring and sibling correlations in agreement with heritabilities with BMI of 20 to 45% [26, 27]. Then, genome-wide association studies (GWAS) allowed to identify to date approximately 140 obesity susceptibility genes, proved to be associated with measures of adiposity (BMI, fat percentage, and/or waist circumference) [28]. The recent GWAS meta-analysis identified 97 BMI-associated loci (including 56 novel ones) in a study involving 339 224 European adults, accounting for 2.7% of BMI variation [29].

About 30 years ago, obesity was introduced in relation with the expression of proinflammatory cytokines. Hotamisligil et al. [30] made the seminal observation that the expression of mRNA for proinflammatory cytokine tumor necrosis factor alpha (TNFα) was high in adipose tissue in several rodent models of obesity, and when TNFα was neutralized, insulin effect was enhanced. Over the years, a close relationship between metabolic pathways and inflammation was identified, including the role of macrophages and adipocytes which secrete many adipocytokines, acting as messenger molecules. Adipocytokines include inflammatory cytokines, such as IL-1 family, IL-6, IL-8, IL-12, interferons (IFNs), TNFα, and transforming growth factor beta (TGFβ), as well as leukemia inhibitory factor (LIF), monocyte chemoattractant protein (MCP-1), macrophage inflammatory protein (MIP-1), leptin, and resistin [31]. IL-1β has been indicated as the prominent cytokine affecting adipocyte metabolic genes and promoting inflammation and metabolic dysfunction in human adipose tissue [32]. The IL-1 family (except IL-18 and IL-33) is encoded as a cluster and polymorphisms of IL1A and IL1B genes have been described as related to obesity, however, conflicting results have been reported [22, 23, 33, 34].

In the studies regarding IL1A, polymorphism -889 C > T (rs1800587) has been investigated the most frequently. Um et al. (2011) revealed that two IL1A polymorphisms -889 C > T (rs1800587) and + 4845 G > T (rs17561) were associated with an increase in BMI in Korean women and lower values were noted in rare allele T carriers. Surprisingly, although fat mass, increased proportionally to BMI, the percentage of body fat was not associated with investigated polymorphisms [23]. The effect of rs1800587 on the transcriptional activity of IL1A was additionally tested in vitro using pre-adipocyte 3 T3-L1 cells and in vivo in obese C57BL/6 J mice and the ones injected with IL-1α. The results indicated that IL1A -889C > T (rs1800587) was a functional polymorphism of IL1A associated with obesity [23, 35]. In Mexican population of adolescents, the same polymorphisms in the IL1A gene, constituting haplotype rs17561T-rs1800587T were associated with slight increase in BMI values in males but not females [33]. However, other study, involving males with ischemic heart disease in Western Australia revealed no interaction between IL1A rs1800587 polymorphism and BMI, although waist circumference, posing a measure of visceral fat mass was larger in TT carriers [22].

Some studies have reported also IL1B polymorphisms as related to obesity. The analysis of + 3953 C/T polymorphism revealed that the T-allele frequency was significantly lower in the overweight (BMI 25–29.9) Korean women and Caucasian men in Sweden [34, 36]. In another Swedish study that involved the population of older men it has been found that IL1B − 31 T > C SNP but not IL1B + 3953 C > T was associated with total fat [37]. In contrast, in Mexican study, none of the examined IL1B polymorphisms: rs1143634 (+3954C > T), rs1143627 (–31C > T), and rs16944 (–511 T > C) were related to BMI [33]. Another pattern has been reported in the population with ischemic heart disease in Western Australia, where IL1B + 3954 T > C (rs1143634) SNPs were studied and TT homozygotes had larger waist circumference and the largest value was noted in individuals with two copies of the IL1A:IL1B T:T haplotype [22].

Polymorphisms of IL-1 receptor family have been reported as either related [11, 34] or not related (rs419598 and rs2234663) to BMI [33, 37].

In our study, the analyzed polymorphisms: IL1A -889C > T (rs1800587), IL1B + 3954 T > C (rs1143634), and IL1RN -87G > A (rs2234677) alone were related neither to BMI value nor fat %. It might have been surprising taking into account the results of the studies presented above, particularly in regard to IL1B + 3953 (rs1143634) that has been suspected to be responsible for body composition in various populations [34, 36]. Although Strandberg et al. (2006) postulated that the association between IL-1 and body fat regulation in humans is robust and not substantially affected by ethnicity, gender, or age [34], it seems more complicated and the results of available studies are at least partially dependent on population parameters. When investigating the relation between genetic background and obesity measurements (including BMI and fat mass), it must be taken into consideration that polymorphisms pose the constant traits and do not change during life, but BMI and fat % may change and depend on many factors including hormonal changes, environmental factors, and individuals’ habits. Thus, if group selection is based on obesity measurements, homogenicity of population is crucial.

Our study revealed the association that according to the authors’ knowledge has not been tested previously. It is currently widely accepted that the interaction between genes, known as epistasis, is the crucial aspect of genetic architecture of complex multicausal phenotypes, both quantitative and qualitative [38,39,40,41]. It can explain the phenotype distributions which could not result from the effect of dominant and recessive gene activity,

Taking the above into consideration, apart from the analysis of polymorphism variants for each of the three genes, in our study, the analysis of the interaction of nonallelic genes was also conducted, limited to the two-gene system (IL1B × IL1RN, IL1B × IL1A, IL1RN × IL1A) for codominant, dominant, and recessive models. In dominant and codominant models the chance of being overweight (for the combination TC × GG) was much (even 16 times) lower than for the reference combination CC × GG. In addition, the chance of being overweight for genotype combination TC × GG-AG (recessive model) was over 3 times lower in comparison to the reference combination CC × GG-AG, independently of IL1RN. These results suggest that the occurrence of rare wild allele T as a heterozygote TC of gene IL1B (rs1143634) protected the young men from high BMI. Additionally, association between the SNP-SNP interaction for the IL1B rs1143634 x IL1RN rs2234677 interaction in the codominant model T/C × G/G (IL1B × IL1RN) and C/C × A/G (IL1B × IL1RN) and the recessive model T/C × G/G-A/G were related also to the fat % in these young males.

In the codominant model, the odds of being OVERFat for T/C x G/G (ILB x IL1RN) was over 7 times lower than for reference combination C/C x G/G (ILB x ILRN); additionally, the odds of being OVERFat for C/C x A/G (ILB x ILRN) was 1.5 times lower than for C/C x G/G (ILB x ILRN). In the recessive model, the odds of being OVERFat for T/C x G/G-A/G (ILB x IL1RN) was 10 times lower than for reference combination C/C x G/G-A/G (IL1B x IL1RN).

The mechanisms by which IL-1 family cytokines affect BMI or fat % are not completely understood yet. It is widely accepted that obesity is characterized by a variety of cytokine-mediated inflammatory responses of chronic, low-grade nature, which result in metabolic and immune derangement [7, 11, 12]. Among IL-1 family, the role of IL-1β has been found pivotal [42]. Over 15 years ago, it was stated that IL-1β was released by human adipocytes and was regulated by TNFα in obesity. Increased concentrations of IL-1β and TNFα synergistically affect lipid metabolism by the regulation of leptin production and release in adipose tissue [36]. Later hypothesis indicated also that IL-1β mediate regulation of leptin effects at hypothalamic level [34].

Currently, it is widely believed that IL-1β production is strictly related to the transcription and subsequent storage of inactive pro–IL-1β, which is converted to an active form by caspase-1 [7, 12, 42, 43]. Active IL-1β downregulates energy generation and expenditure in adipocyte mitochondria leading to fat accumulation and weight gain [44]. The mechanism of this process involves affecting mitochondrial oxidative phosphorylation by suppressing respiratory chain supercomplex formation, without disturbances in other mitochondrial functions. Supercomplexes formation allows to maintain proton gradient necessary for energy production and so that oxygen consumption and fatty acid oxidation are diminished in IL-1β stimulated cells. The pathway involves activation of IL-1R by IL-1β and further recruitment of MyD88. IL-1 signals are propagated downstream of IL-1R through Myddosome, comprising MyD88 and the IL-1R-associated kinases (IRAK). Recently, the pivotal role of IRAK 2, due to its location within mitochondria has been indicated [45], posing the critical link in the modulation of mitochondrial energy metabolism by IL-1β.

The role of IL-1 family polymorphisms in obesity remains poorly understood and only scanty discussed, mostly in the context of increased cytokine production by the carriers of certain alleles [4, 23, 36, 37]. However, it frequently remains unclear which promoters are more potent. For example IL1B − 31 T > C polymorphism has been found important but it was not clear if T or C allele was more potent promoter of IL-1β production and complex interactions with other polymorphisms such as the IL1B − 511 C > T SNP have been suggested [37]. Other studies suggested the association between the IL1B + 3953C and the IL1RN*2 alleles to produce increased IL-1 bioactivity [34].

The main limitations of our study were the small number of analyzed polymorphisms, lack of AA genotype in IL1A polymorphism rs1800587, and lack of IL-1 measurements in blood and adipose tissue. Since the study involved healthy volunteers recruited from the students of Military University, only minimally invasive procedures were possible and so that even blood could not be taken. However, in order to further confirm the interaction between IL1B and IL-1RN which has been found in this study, more polymorphisms and larger population should be analyzed.


This study further confirmed that IL-1 family genes are involved in genetic background of obesity. Our results indicated that even if certain polymorphisms in single genes of IL-1 family cannot be defined as related to obesity in studied population, the interaction should be analyzed. We have found that IL1B rs1143634/IL1RN rs2234677 interaction for both dominant model and codominant model TC × GG was significant for increased BMI. In this interaction, rare wild allele T as a heterozygote TC of gene IL1B (rs1143634) protected the young men from high BMI. Similarly, the SNP-SNP interaction for the IL1B x IL1RN in the codominant model T/C × G/G and the recessive model T/C × G/G-A/G were related to the fat %.

Thus, the association studies may help to establish potential new genetic biomarkers for clinical molecular diagnosis of obesity before its metabolic implications

Materials and methods

Study protocol

This case-control association study included 101 volunteers—healthy male cadets 19–25 years old. They were randomly recruited via advertising on the Military University campus. The groups did not differ in sex (males only), age, and height (Table 1) as well as also nutrition and activity, as all were students of Military University, lived and ate at the university, and covered the same fitness tests. All participants were selected on the basis of the questionnaire screening for exclusion criteria such as past diseases, injuries, and related ailments and the presence of severe and chronic pain of any organ or system, both in the past and currently. To confirm the good health of the volunteers, general medical examinations were conducted, as well as electrocardiography (ECG).

The examined soldiers lived in a hall of residence on the university’s premises and ate the same meals at the students’ cafeteria. They performed the same physical activity, which is a part of academic syllabus and relates to their military service and responsibilities. The general fitness of soldiers in Polish Army is estimated based on the obligatory fitness test, performed once a year, which is also applied to the Military University students.

Before the study was conducted, participants were acquainted with the protocol and research methods. All participants received a written information sheet concerning the study, providing all pertinent information (purpose, procedures, risks, and benefits of participation). All cadets gave the written consent to participate in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of the Military Institute of Hygiene and Epidemiology—resolution number 07/2018, dated 2.02.2018.

Anthropometry and body composition

Anthropometric measurements and body composition were obtained using standard methods. Height was measured using a portable stadiometer with a precision of 0.1 cm (without shoes) (TANITA HR-001, Tanita Corporation, Japan). Body composition (including fat %) and body weight were measured using bioelectrical impedance analysis (BIA) performed by the TANITA MC-780 machine (Tanita Corporation, Japan) with accuracy to 0.1 kg according to the procedure specified in the instruction manual (lightly dressed, without shoes). The assessment of BMI values was made in accordance with the criteria set out by WHO [1, 46]

The subjects were divided into two groups depending on their BMI value. The overweight group (OVERBMI) was made up of people with BMI of ≥ 25.0, while the control group (CONBMI) consisted of people with BMI values between 20.0 and 25.0 [2]. Further analyses were made in groups divided on the basis of sharing fat mass in total subject’s weight (fat %). The control group (CONFat) comprised of subjects with fat % below 20.0% while the overweight group (OVERFat) was characterized with fat % over 20.0% [6].

Genetic analyses

The buccal cells donated by the subjects were collected using two Copan FLOQSwabs (Interpath, Australia) according to standard procedure. Genomic DNA was extracted from the buccal cells using a High Pure PCR Template Preparation Kit (Roche Diagnostics, Germany). The extraction was performed according to the manufacturer’s instructions. DNA samples of good quality and quantity were stored at − 20 °C for further analysis. The exclusion criteria were as follows: failure in DNA extraction, DNA degradation, abnormal gene detect results, and incomplete basic information.

All samples were genotyped in duplicate on a CFX Connect Real-Time PCR Detection System (BioRad, USA) using TaqMan Pre-Designed SNP Genotyping Assays for IL1B (rs1143634) C___9546517_10, IL1A (rs1800587) C___9546481_30, and IL1RN (rs2234677) C__11948096_10 single-nucleotide polymorphisms (SNPs) (Applied Biosystems, USA), which include primers and fluorescently labeled (VIC and FAM) MGBTM probes for alleles detection. Genotyping was performed according to the manufacturer’s protocol using TaqPath™ ProAmp™ Master Mix (Applied Biosystems, USA). Briefly, conditions for reaction were as following: 30 s of pre-read in 60 °C, 5 min of initial denaturation in 95 °C, cycling 15 s of denaturation in 95 °C, 1 min of primer hybridization and elongation in 60 °C, repeated in 40 cycles, 30 s of final elongation in 60°. The amplified products were visualized and primarily analyzed using CFX Maestro 4.0 Software (BioRad, USA).

Statistical analysis

IBM SPSS Statistics (version 5.1) was used to calculate differences between groups with Student t-test, and genotype frequencies were analyzed using Fisher’s exact test. Statistical analysis was done using the SNPassoc package for R (version 1.9-2, R Foundation for Statistics Computing, Single-locus analysis including SNP × SNP interaction was performed using SNPassoc (version 1.9.2) considering four genetic models (codominant, dominant, recessive, and overdominant). The models were constructed with respect to the minor allele. The p-values for SNP × SNP (epistasis) were calculated using the log-likelihood ratio test. The level of statistical significance was set at p < 0.05.

Availability of data and materials

All datasets used and/or analyzed during whole procedures are available from the corresponding author on reasonable request.



Interleukin 1


World Health Organization


Body mass index


Army Weight Control Program


Interleukin 1 receptor


Genome-wide association studies


Tumor necrosis factor alpha




Transforming growth factor beta


Leukemia inhibitory factor


Monocyte chemoattractant protein


Macrophage inflammatory protein


IL-1 receptor antagonist


Waist to hip ratio


IL-1R-associated kinases


Single-nucleotide polymorphism


IL-1R accessory protein


Single-Ig IL-1R-related molecule


Three-Ig-domain-containing IL-1R related


  1. Obesity: preventing and managing the global epidemic. Report of a WHO Consultation (WHO Technical Report Series 894).


  3. Maffetone PB, Rivera-Dominguez I, Laursen PB. Overfat and underfat: new terms and definitions long overdue. Hypothesis and Theory Front. Public Health. 2017;4:279

  4. Rajender Rao K, Nirupama Lal, Giridharan NV. Genetic & epigenetic approach to human obesity. Indian J Med Res. 2014;140:589-603

  5. Herrera BM, Lindgren CM. The genetics of obesity. Curr Diab Rep. 2010;10:498–505.

    Article  CAS  Google Scholar 

  6. Friedl KE. Body composition and military performance – many things to many people. J Strength Cond Res. 2012; 26 suppl 7, S87-S100

  7. Ballak DB, Stienstra R, Tack CJ, Dinarello CA, van Diepen JA. IL-1 family members in the pathogenesis and treatment of metabolic disease: Focus on adipose tissue inflammation and insulin resistance. Cytokine. 2015;75:280–90.

    Article  CAS  Google Scholar 

  8. Elks CE, den Hoed M, Zhao JH, Sharp SJ, Wareham NJ, Loos RJ, Ong KK. Variability in the heritability of body mass index: a systematic review and meta-regression. Front Endocrinol (Lausanne). 2012;3:29.

  9. Maes HH, Neale MC, Eaves LJ. Genetic and environmental factors in relative body weight and human adiposity. Behav Genet. 1997;27:325–51.

    Article  CAS  Google Scholar 

  10. Kershaw EE, Flier JS. Adipose tissue as an endocrine organ. J Clin Endocrinol Metabol. 2004;89:2548–56.

    Article  CAS  Google Scholar 

  11. Yu GI, Song DK, Shin DH. Associations of IL1RAP and IL1RL1 gene polymorphisms with obesity and inflammation mediators. Inflamm Res. 2020;69:191–202.

    Article  CAS  Google Scholar 

  12. Febbraio MA. Role of interleukins in obesity: implications for metabolic disease. Trends Endocrinol Metabol. 2014;25:312–9.

    Article  CAS  Google Scholar 

  13. Kern PA, Ranganathan S, Li C, Wood L, Ranganathan G. Adipose tissue tumor necrosis factor and interleukin-6 expression in human obesity and insulin resistance. Am J Physiol Endocrinol Metabol. 2001;280:E745–51.

    Article  CAS  Google Scholar 

  14. Cottam DR, Mattar SG, Barinas-Mitchell EEG, Kuller L, Kelley DE, Schauer PR. The chronic inflammatory hypothesis for the morbidity associated with morbid obesity: implications and effects of weight loss. Obesity Surg. 2004;14:589–600.

    Article  Google Scholar 

  15. Borthwick LA. The IL-1 cytokine family and its role in inflammation and fibrosis in the lung. Semin Immunopathol. 2016;38:517–34.

    Article  CAS  Google Scholar 

  16. Palomo J, Dietrich D, Praxedis M, Palmer G, Cem G. The interleukin (IL)-1 cytokine family--balance between agonists and antagonists in inflammatory diseases. Cytokine. 2015;76:25–37.

    Article  CAS  Google Scholar 

  17. Yazdi AS, Ghoreschi K. The interleukin-1 family. Adv Exp Med Biol. 2016;941:21–9.

    Article  CAS  Google Scholar 

  18. Witkin SS, Gerber S, Ledger WJ. Influence of interleukin-1receptor antagonist gene polymorphism on disease. Clin Infect Dis. 2002;34:204–9.

    Article  CAS  Google Scholar 

  19. Dursun E, Gezen-Ak D, Ertan T, Bilgiç B, Gürvit H, Emre M, et al. Interleukin-1α –889 C/T polymorphism in turkish patients with late-onset Alzheimer’s disease. Dement Geriatr Cogn Disord. 2009;27:82–7.

    Article  CAS  Google Scholar 

  20. Serretti A, Olgiati P, Politis A, Malitas P, Albani D, Dusi S, et al. Lack of Association between interleukin-1 alpha rs1800587 polymorphism and Alzheimer’s disease in two independent European samples. J Alzheimer's Disease. 2009;16:181–7.

    Article  CAS  Google Scholar 

  21. Eskola PJ, Kjaer P, Daavittila IM, Solovieva S, Okuloff A, Sorensen JS, et al. Genetic risk factors of disc degeneration among 12-14-year-old Danish children: a population study. Int J Mol Epidemiol Gen. 2010;1:158.

    CAS  Google Scholar 

  22. Carter KW, Hung J, Powell BL, Wiltshire S, Foo BTX, Leow YC, et al. Association of Interleukin-1 gene polymorphisms with central obesity and metabolic syndrome in a coronary heart disease population. Human Gen. 2008;124:199–206.

    Article  CAS  Google Scholar 

  23. Um JY, Rim HK, Kim SJ, Kim HL, Hong SH. Functional polymorphism of IL-1 alpha and its potential role in obesity in humans and mice. PloS one. 2011;6:29524.

    Article  Google Scholar 

  24. Kopelman PG. Obesity as a medical problem. Nature. 2000;404:635–43.

    Article  CAS  Google Scholar 

  25. O’Rahilly S. Human genetics illuminates the paths to metabolic disease. Nature. 2009;462:307–14.

    Article  Google Scholar 

  26. Elks CE, den Hoed M, Zhao JH, Sharp SJ, Wareham NJ, Loos RJ, et al. Variability in the heritability of body mass index: a systematic review and meta-regression. Front Endocrinol. 2012;3:29.

    Article  Google Scholar 

  27. Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Magi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518:187–96.

    Article  CAS  Google Scholar 

  28. Castillo JJ, Orlando RA, Garver WS. Gene-nutrient interactions and susceptibility to human obesity. Genes & Nutrition. 2017;12:29.

    Article  Google Scholar 

  29. Albuquerque D, Nóbrega C, Manco L, Padez C. The contribution of genetics and environment to obesity. Brit Med Bull. 2017;123:159–73.

    Article  Google Scholar 

  30. Hotamisligil GS, Shargill NS, Spiegelman BM. Adipose expression of tumor necrosis factor-alpha: direct role in obersity-linked insulin resistance. Science. 1993;259:87–91.

    Article  CAS  Google Scholar 

  31. Frühbeck G. Adipose tissue as a source of vasoactive factors. Current Medicinal Chemistry - Cardiovas Hematol Agents. 2004;2:197–208.

    Article  Google Scholar 

  32. Caër C, Rouault C, Le Roy T, Poitou C, Aron-Wisnewsky J, Torcivia A, et al. Immune cell-derived cytokines contribute to obesity-related inflammation, fibrogenesis and metabolic deregulation in human adipose tissue. Sci Rep. 2017;7:3000.

    Article  Google Scholar 

  33. Mendoza-Carrera F, Ramírez-López G, Hernández-Ramos LE, Leal-Cortés C, Portilla-de-Buen E, Castro-Martínez XH, et al. Interleukin-1 alpha polymorphisms are associated with body mass index in male but not in female adolescents. Arch Med Res. 2019;50:151–7.

    Article  CAS  Google Scholar 

  34. Strandberg L, Lorentzon M, Hellqvist A, Nilsson S, Wallenius V, Ohlsson C, et al. Interleukin-1 system gene polymorphisms are associated with fat mass in young men. J Clin Endocrinol Metab. 2006;91:2749–54.

    Article  CAS  Google Scholar 

  35. Dominici R, Cattaneo M, Malferrari G, Archi D, Mariani C, Grimaldi L, et al. Cloning and functional analysis of the allelic polymorphism in the transcription regulatory region of interleukin-1α. Immunogenetics. 2002;54:82–6.

    Article  CAS  Google Scholar 

  36. Um JY, Chung HS, Song MY, Shin HD, Kim HM. Association of interleukin-1beta gene polymorphism with body mass index in women. Clin Chem. 2004;50:647–50.

    Article  CAS  Google Scholar 

  37. Strandberg L, Mellström D, Ljunggren Ö, Grundberg E, Karlsson MK, Holmberg AH, et al. IL6 and IL1B polymorphisms are associated with fat mass in older men: the MrOS Study Sweden. Obesity. 2008;16:710–3.

    Article  CAS  Google Scholar 

  38. Moore JH. The ubiquitous nature of epistasis in determining susceptibility to common human diseases. Hum Hered. 2003;56:73–82.

    Article  Google Scholar 

  39. Moore JH, Williams SM. Traversing the conceptual divide between biological and statistical epistasis: systems biology and a more modern synthesis. Bioessays. 2005;27:637–46.

    Article  CAS  Google Scholar 

  40. Phillips PC. Epistasis – the essential role of gene interactions in the structure and evolution of genetic systems. Nat Rev Genet. 2008;9:855–67.

    Article  CAS  Google Scholar 

  41. Carlborg O, Haley CS. Epistasis: too often neglected in complex trait studies? Nat Rev Genet. 2004;5:618–25.

    Article  CAS  Google Scholar 

  42. Antonioli L, Moriconi D, Masi S, Bottazzo D, Pellegrini C, Fornai M, et al. Differential impact of weight loss and glycemic control on inflammasome signaling. Obesity (Silver Spring). 2020;28:609–15.

    Article  CAS  Google Scholar 

  43. Cordell HJ. Detecting gene-gene interactions that underlie human diseases. Nat Rev Genet. 2009;10:392–404.

    Article  CAS  Google Scholar 

  44. Tavernarakis N. Inflammation brakes mitochondrial metabolism in obesity. Nature Immunol. 2020;21:1143–5.

    Article  CAS  Google Scholar 

  45. Zhou H, Wang H, Yu M, Schugar RC, Qian W, Tang F, et al. IL-1 induces mitochondrial translocation of IRAK2 to suppress oxidative metabolism in adipocytes. Nat Immunol. 2020;21:1219–31.

    Article  CAS  Google Scholar 

  46. Defining adult overweight and obesity. Overweight & Obesity CDC. 2020,

Download references


Not applicable.


The study was financed by the Ministry of Science and Higher Education in 2020/2022 as part of the Scientific School of the Academy of Physical Education in Warsaw - SN No. 5 “Biomedical determinants of physical fitness and sports training in adult population” and by grant “KOSCIUSZKO ” II edition (for Ryszard Puchala, PhD) Decision No. 509/2017/DA from the Ministry of National Defense.

Author information

Authors and Affiliations



Conceptualization: E.M. and P.C.; methodology: E.M., P.C., B.A., A.L.D., and O.A.; formal analysis: E.M., K.H.L, M.M.S., and A.G.; investigation: E.M., A.G., B.A., O.A., E.S., and P.C.; resources: B.A, O.A., and E.M.; data curation: E.M., K.H.L, M.M.S., K.L. L.D., and B.A.; writing—original draft preparation: EM and AC; writing—review and editing, P.C., A.B., A.M., A.C., and A.L.D.; supervision: P.C. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Anna Cywińska.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Military Institute of Hygiene and Epidemiology—resolution number 07/2018, dated 2.02.2018. Informed consent statement to participate in the study was obtained from all subjects involve in the study.

Consent for publication

Not applicable

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maculewicz, E., Antkowiak, B., Antkowiak, O. et al. The interactions between interleukin-1 family genes: IL1A, IL1B, IL1RN, and obesity parameters. BMC Genomics 23, 112 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: