Adenosine deaminase and interleukin-1 receptor antagonist genetic polymorphisms among obese children with versus without metabolic dysfunction-associated fatty liver disease

Article information

Clin Exp Pediatr. 2025;68(10):808-818

Publication date (electronic) : 2025 May 29

doi : https://doi.org/10.3345/cep.2025.00731

Hala M. Sakhr, MD^,1, Mohammed H. Hassan, MD^,2, Azza Mohamed Taha, MD¹, Ali Helmi Bakri, MD¹

¹Department of Pediatrics, Faculty of Medicine, South Valley University, Qena, Egypt

²Department of Medical Biochemistry, Faculty of Medicine, South Valley University, Qena, Egypt

Corresponding author: Mohammed H. Hassan, MD. Department of Medical Biochemistry, Faculty of Medicine, South Valley University, Qena, Egypt Email: mohammedhosnyhassaan@yahoo.com, mohammedhosnyhassaan@med.svu.edu.eg

Co-corresponding author: Hala M. Sakhr. Department of Pediatrics, Faculty of Medicine, South Valley University, Qena , Egypt Email: hala.sakhr@med.svu.edu.eg, halasakhr@yahoo.com

Received 2025 March 28; Revised 2025 April 22; Accepted 2025 May 8.

Abstract

Introduction

Metabolic-associated fatty liver disease (MAFLD), previously known as nonalcoholic fatty liver disease (NAFLD), has become the most prevalent cause of chronic liver disease among children and adolescents [1]. The global prevalence of MAFLD in this demographic has risen alarmingly, paralleling the increase in pediatric obesity rates. Notably, obesity is a major contributor to MAFLD, with prevalence rates of 34% among overweight and obese children and adolescents aged 1-19 years [2].

Beyond environmental factors, genetic predispositions play a significant role in MAFLD susceptibility. Interleukin- 1 receptors (IL-1Rs) are integral to the inflammatory response, mediating the effects of interleukin-1 cytokines. Research indicates that IL-1β, a key cytokine interacting with IL-1Rs, plays a significant role in obesity-related liver conditions. In studies involving IL-1β knockout mice subjected to high-fat diets, these mice exhibited increased adipose tissue expansion and reduced liver weight compared to wild-type controls. This suggests that IL-1β contributes to adipose tissue inflammation and limits its expandability, leading to ectopic fat deposition in the liver, thereby promoting hepatic steatosis [3].

Interleukin-1 receptor antagonist (IL-1Ra), encoded by the IL1RN gene, plays a pivotal role in modulating inflammatory responses by inhibiting the activities of interleukin-1 cytokines. IL-1 family cytokines, including IL-1Ra, are integral to liver inflammation. Dysregulation of IL-1 signaling contributes to the progression of NAFLD, now commonly referred to as MAFLD. Elevated IL-1Ra levels have been observed in individuals with MAFLD, implicating its role in disease pathogenesis [4]. Therapeutic Potential targeting IL-1 pathways, particularly using IL-1 inhibitors, has been proposed as a therapeutic strategy for MAFLD. Modulating IL-1Ra levels may ameliorate liver inflammation and improve metabolic outcomes in affected individuals [5].

Adenosine deaminase (ADA) is crucial for purine metabolism, converting adenosine into inosine. Its activity influences extracellular adenosine levels, which interact with adenosine receptors affecting various metabolic pathways. Alterations in adenosine signaling have been implicated in obesity and liver diseases. For instance, activation of A1 adenosine receptors stimulates adipogenesis and lipid accumulation, promoting weight gain, whereas A2B receptor activation reduces adipogenesis and insulin resistance [6]. Furthermore, studies have shown that adenosine generated by ethanol metabolism contributes to hepatic steatosis via both A1 and A2B receptors, suggesting that targeting adenosine receptors may be effective in preventing alcohol-induced fatty liver [7].

The aim of the current study was to assess the relationship between the development of metabolically related fatty liver and the ADA G22A and IL-1NR gene polymorphisms in obese Egyptian children and adolescents.

Methods

3. Laboratory investigations

1) Routine laboratory investigations

Five milliliters of venous blood were obtained from all participants and separated into 2 samples; the first (3 mL) was collected in a plain vacotainer tube, centrifuged, and the resultant sera were utilized for routine biochemical analyses in the form of liver enzymes (alanine transaminase [ALT] and aspartate transaminase [AST], lipid profile (total cholesterol, TG, very low-density lipoprotein [VLDL], HDL, and low-density lipoprotein [LDL]; VLDL/ HDL-C and TG/HDL-C ratios were calculated), glycemic profile (random blood glucose, fasting blood glucose, 2-hour postprandial glucose levels, glycated hemoglobin [HbA1c], homeostatic model assessment for insulin resistance [HOMA-IR]), and serum cortisol levels. While the second (2 mL) was collected in a tube containing the anticoagulant ethylenediaminetetraacetic acid (EDTA) and employed for complete blood count.

2) Genetic analysis

(1) DNA Extraction and Genotyping of ADA G22A Polymorphism

Venous blood samples (2 mL) were collected in EDTA tubes and stored at -80°C until genetic analysis. DNA was extracted using the G-spin Total DNA Extraction Kit (iNtRON Biotechnology, Korea) and stored at -80°C. Genotyping of the ADA G22A polymorphism was performed using restriction fragment length polymorphism-polymerase chain reaction (PCR). The PCR reaction (25 μL) contained 12.5-μL PCR master mix (catalog no. 25028, iNtRON Biotechnology), 1 μL of forward primer (5′-GCCCGGCCCGTTAAGAAGAGC-3′), 1 μL of reverse primer (5′-GGTCAAGTCAGGGGCAGAAGCAGA-3′) [13,14], 8.5 μL nuclease-free water, and 2 μL of extracted DNA. Thermal cycling conditions included an initial denaturation at 94°C for 15 minutes, followed by 36 cycles of denaturation (94°C, 40 seconds), annealing (66°C, 80 seconds), and extension (72°C, 80 seconds), with a final extension at 72°C for 8 minutes using an Applied Biosystems thermal cycler (Perkin-Elmer 9600, USA).

PCR products were digested with TaqI enzyme (Cat. No. ER0671, Thermo Fisher Scientific, USA) in a reaction containing 1-μL enzyme, 2-μL buffer, 7-μL nuclease-free water, and 10-μL PCR product, incubated at 65℃ for 1 hour. The digested fragments were resolved on a 2% agarose gel stained with ethidium bromide and visualized under a UV transilluminator. A 50-bp DNA ladder (Catalog No. 10416014, Thermo Fisher Scientific) was used for band size determination. The ADA G22A genotype was identified based on banding patterns: ADA1/ADA1 (GG) showed 2 bands (245 and 152 bp), while ADA1/ADA2 (GA) displayed 3 bands (397, 245, and 152 bp) (Fig. 1).

Fig. 1.

ADA G22A and IL-1RN single nucleotide polymorphisms using PCR techniques. (A) Gel electrophoresis of ADA G22A polymorphism using the RFLP-PCR method; the numbers refer to lanes. Lane 1 shows the 50-bp DNA ladder; lanes 2–6 and 8 represent heterozygotes GA genotypes (397 bp, 245 bp, and 152 bp). Lane 7 represents homozygotes GG (245 bp and 152 bp). (B) Gel electrophoresis of IL-1RN polymorphism using conventional PCR method; the numbers refer to lanes. Lane 1 shows 50-bp DNA ladder, lanes 2–7 represents homozygotes *1/*1, and lane 8 represents homozygotes *2/*2. ADA, adenosine deaminase; IL-1RN, interleukin-1 receptor antagonist; RELP-PCR, restriction fragment length polymorphism-polymerase chain reaction.

(2) Genotyping of interleukin-1 receptor antagonist (IL-1RN) gene polymorphism

Genotyping of IL-1RN gene polymorphism was performed using conventional PCR without restriction enzyme digestion. The primers used (GeneBank accession no. X64532) were: sense (5′-CTCAGCCAACACTCCTAT-3′) and antisense (5′-TCCTGGTCTGCAGGTAA-3′), as described by Swellam et al. [15]. The PCR reaction (25 μL) included 12.5-μL PCR master mix (catalog no. 25028, iNtRON Biotechnology), 1 μL each of forward and reverse primers, 8.5-μL nuclease-free water, and 2 μL of extracted DNA. Thermal cycling conditions consisted of an initial denaturation at 94°C for 3 minutes, followed by 35 cycles of denaturation (94°C, 1 minute), annealing (°C, 1 minute), and extension (72°C 1 minute), with a final extension at 72°C for 10 minutes. Reactions were conducted using a Biometra thermal cycler (Germany, serial number 2603204). PCR products were resolved on a 2% agarose gel stained with ethidium bromide, along with a 50 bp DNA ladder (Catalog No. 10416014, Thermo Fisher Scientific). Electrophoresis was performed at 120 V for 15–20 minutes, and bands were visualized under UV light using a gel documentation system. The IL-1RN alleles were identified based on product sizes: IL-1RN0 (150 bp), IL-1RN1 (410 bp), IL-1RN2 (240 bp), IL-1RN3 (325 bp), IL-1RN4 (500 bp), and IL-1RN5 (595 bp).

Results

This study compared 100 obese children with 50 ageand sex-matched healthy controls (median age: 12.0 years vs. 12.1 years; sex distribution: 38% vs. 40% males, 62% vs. 60% females). Obese children had significantly higher weight (73 kg vs. 54.5 kg, P<0.001), weight z score (2.52 vs. 0.41, P<0.001), BMI (33.69 kg/m² vs. 20.82 kg/m², P<0.001), BMI z score (2.5 vs. 0.5, P<0.001), and waist circumference (89 cm vs. 67.5 cm, P<0.001). These findings confirm a strong association between obesity and central adiposity. Lower hemoglobin levels in obese children (11.06±1.32 g/dL vs. 11.78±0.97 g/dL, P<0.001). Higher insulin resistance (HOMA-IR: 2.6 vs. 1.01, P<0.001), strongly indicating metabolic dysfunction. Elevated post meridiem (PM) cortisol levels (P<0.001), suggesting altered diurnal cortisol regulation. AST and ALT are significantly higher in patients (P=0.032 and P<0.001, respectively) and AST/ALT ratio is significantly lower in obese children (0.91 vs. 1.28, P<0.001), indicating altered liver enzyme balance. Obese children had significantly higher total cholesterol, LDL, VLDL/ HDL ratio, and TGs/ HDL ratio (P<0.001). Obese children had significantly lower VLDL and HDL (Table 1).

Table 1.

Clinical and laboratory data of study participants

Table 2 compares obese children with MAFLD to those without, highlighting greater obesity severity, metabolic dysfunction, inflammation, and liver injury in the MAFLD group. Higher weight (80 kg vs. 62 kg, P<0.001) and BMI (36.15 kg/m² vs. 31.12 kg/m², P<0.001) in the MAFLD group. Higher weight z score & BMI z score (P<0.001), confirming greater obesity severity. Larger waist circumference (89.5 cm vs. 85 cm, P<0.001), indicating stronger central obesity association. laboratory data showed a higher white blood cell (WBC) count (P=0.001), suggesting increased inflammation, and higher platelet count (P<0.001), possibly linked to liver dysfunction. Greater insulin resistance (HOMA-IR, P=0.016) and higher HbA1C levels (P=0.025), indicating metabolic disturbances. Significant higher AST and ALT levels in the MAFLD group. Lipid profile showed higher total cholesterol and LDL (P<0.001), reflecting a more atherogenic profile. Higher VLDL/HDL-C ratio (P=0.002) and TG/HDL-C ratio (P<0.001), suggesting increased cardiovascular risk. Lower HDL (P<0.001), increasing cardiovascular risk.

Table 2.

Clinical and laboratory data by study subgroup

Table 3 examines the genetic distribution of ADA G22A polymorphisms among patients and controls, as well as between obese children with and without MAFLD. GG genotype is significantly more frequent in patients (38%) than in controls (14%) (P=0.002), suggesting a strong association with disease susceptibility. GA genotype is more common in controls (86%) than in patients (62%), implying a potential protective effect of the GA variant. G allele is significantly more frequent in patients (69%) versus controls (57%) (P=0.040), indicating an increased disease risk associated with the G allele. GG genotype (odds ratio [OR], 3.765; 95% confidence interval [CI], 1.538–9.215), meaning patients are ~3.76 times more likely to have this genotype and G allele (OR, 1.679; 95% CI, 1.022–2.759) suggesting a moderate association with disease risk. No significant differences in ADA G22A genotype or allele distributions between the 2 groups with and without MAFLD.

Table 3.

Genetic profile of ADA G22A polymorphisms among the study participants

Table 4 presents the genetic profile of IL-1RN polymorphisms among different studied groups, focusing on genotype and allele distributions. Comparison between patients and controls showed that the *1/*1 genotype is more frequent in controls (66%) than in patients (61%). The *1/*2 genotype appears in 25% of patients versus (4%) in controls. The *2/*2 genotype is found in 14% of patients versus 30% in controls. The genotype distribution shows a statistically significant difference (P<0.001), suggesting an association between IL-1RN polymorphism and disease presence. Comparison between obese children with and without MAFLD revealed that the *1/*1 genotype is less frequent in obese children with MAFLD (58%) compared to those without MAFLD (64%). The *1/*2 genotype is more common in obese children with MAFLD (38%) than those without (12%). The *2/*2 genotype is more frequent in obese children without MAFLD (24%) than in those with MAFLD (4%). The genotype distribution is significantly different (P<0.001), indicating a potential role of IL-1RN polymorphism in MAFLD development. The allele frequencies (ILRN1 vs. ILRN2) show no significant difference (P=0.262), suggesting that specific genotypes, rather than alleles alone, might be influencing MAFLD risk.

Table 4.

Genetic profile of IL-1RN polymorphisms among study participants

Table 5 presents the mean lipid profile parameters for different genotypes in relation to the 2 gene polymorphism. VLDL is significantly higher in GA (26.5 mg/dL) than in GG (21 mg/dL) (P<0.001) and HDL is also higher in GA (35.05 mg/dL) than in GG (31.95 mg/dL) (P=0.049). No significant differences in other lipid profile. Lipid profile analysis in relation to IL-1RN Genotype showed that HDL is highest in *2/*2 (39.1 mg/dL), followed by *1/*1 (33 mg/dL) and *1/*2 (32.8 mg/dL) (P=0.004), suggesting *2/*2 is linked to better HDL levels. LDL is highest in *1/*2 (155 mg/dL), lowest in *2/*2 (31.47 mg/dL) (P=0.045), indicating a possible protective effect of *2/*2. No significant differences in total cholesterol, TGs, VLDL, VLDL/HDL ratio, or TGs/HDL ratio.

Table 5.

Mean lipid profile parameters of different genotypes and alleles of ADA G22A polymorphisms among included patients (N=100)

Fig. 2 showed that VLDL/HDL-C ratio is fair (area under the curve [AUC]=0.68) to discriminate Obese children with MAFLD from those without MAFLD among the included patients at cutoff value of >0.6308 with sensitivity (80%), specificity (58%), positive predictive value (PPV; 65.6%), negative predictive value (NPV; 74.4%), and accuracy (69%). TG/HDL-C ratio at cutoff value >3.0685 is a moderately effective tool (AUC=0.752) for identifying MAFLD in obese children, with good specificity (88%) and PPV (84.2%) but moderate sensitivity (64%), NPV (71%), and overall accuracy=76%.

Fig. 2.

Receiver operating characteristic curve showing validity of VLD-C-to-HDL-C and TG-to-HDL-C ratios at discriminating between obese children with versus without MAFLD. HDL-C, high-density lipoprotein cholesterol; MAFLD, metabolic dysfunction-associated fatty liver disease; TG, triglyceride; VLD-C, very-low-density cholesterol.

Multivariate logistic regression analyses identified BMI, waist circumference, TG/HDL ratio, and LDL as significant predictors of MAFLD in obese children.

Discussion

Pediatric obesity is a rising global health concern, with a significant increase in prevalence among children and adolescents. The long-term consequences include type 2 diabetes, cardiovascular and respiratory diseases, and autoimmune disorders [16,17]. While pediatric MAFLD pathogenesis remains unclear, studies suggest that obesity [18,19] and insulin resistance [20,21] are key risk factors. However, susceptibility varies among children [22,23].

In tandem with the rising incidence of childhood obesity, fatty liver is anticipated to grow among the leading causes of end-stage liver disease in young adults and children [24]. The magnitude of obesity among Egyptian children in various governorates of Egypt revealed that one out of seven of 6–12 years old children in Qena governorate were obese with female predominance due to bad dietary habits with lack of physical activity [25]. The overall prevalence of obesity among primary school children in Port Said governorate was 20.6% [26], while in Assiut governorate it was 12.28% [27], in Alexandria governorate it was 9% [28]. A country-based study from the Eastern Mediterranean region regarding the prevalence of obesity among school aged children in Egypt revealed that the overall prevalence was 10.6%, highest in Lower Egypt (10.0%) and urban governorates (9.2%) compared to Upper Egypt governorates (6.3%) [29]. The findings of the current work revealed female sex predominance among obese children with significantly higher BMI and waist circumference among obese children particularly those with MAFLD. These findings were in agreement with previous researches [25-28].

In the current study, patients had lower hemoglobin levels, higher insulin resistance, altered cortisol levels, and significant lipid profile changes compared to controls. Anemia in children may result from nutritional deficiencies or chronic diseases. Studies suggest a link between iron metabolism and lipid profiles in children. Zhu et al. [20] observed that children with dyslipidemia had lower serum ferritin and transferrin levels, indicating a connection between iron deficiency and lipid abnormalities, especially in obese children. This suggests that iron deficiency is more prevalent in overweight or obese individuals.

Cortisol, produced by the adrenal glands, plays a crucial role in glucose metabolism and stress response. Dysregulated cortisol levels are linked to metabolic syndrome, as chronic stress elevates cortisol, leading to visceral adiposity, insulin resistance, dyslipidemia, and hypertension, highlighting the The hypothalamic-pituitary-adrenal axis's role in metabolic health [30,31].

A study by Chiarelli and Marcovecchio [31] found that insulin resistance in obese children is associated with hypertriglyceridemia, hypercholesterolemia, and low HDL-C. Similarly, Chang et al. [32] examined inflammatory markers, lipid profiles, and insulin sensitivity, reinforcing their interconnection in pediatric obesity.

Comparing obese children with and without MAFLD reveals significant clinical and laboratory differences, emphasizing the strong link between MAFLD and metabolic syndrome in children. Higher weight z scores, BMI z scores, and waist circumference in the MAFLD group align with previous studies, supporting central obesity’s role in liver fat accumulation [33,34].

Consistent with this study, Xing et al. [35] reported elevated WBC counts in MAFLD patients, reinforcing the association between inflammation and metabolic disturbances. Reported also that children with MAFLD had higher fasting plasma glucose, HOMA-IR levels, ALT levels, and dyslipidemia, reinforcing the association between MAFLD and metabolic dysfunction. Insulin resistance both contributes to and results from hepatic steatosis, creating a vicious cycle that worsens metabolic dysfunction [36].

A study on Danish children and adolescents (6–18 years old) found that overweight and obese individuals had higher ALT levels across all age groups. The study established ALT cutoff values for hepatic steatosis diagnosis at 24.5 U/L for girls and 34.5 U/L for boys [37]. Another study using ultrasound and ALT levels to assess NAFLD found that 23.61% of overweight and obese children showed fatty liver features on ultrasound, yet 72.73% of them had normal ALT levels (<42 U/L), indicating that ALT alone may underestimate hepatic steatosis [38]. Additionally, a study on prepubertal children with obesity reported elevated liver enzymes, leptin, insulin resistance markers, inflammation, and endothelial dysfunction, further linking obesity to metabolic disturbances [39].

The current study observed dyslipidemia in the MAFLD group, indicating an atherogenic lipid imbalance that increases cardiovascular risk. This aligns with studies identifying dyslipidemia as a common metabolic abnormality in adolescents with metabolic syndrome [40-42].

ADA plays a key role in purine metabolism, immune function, and oxidative stress regulation. The current study found that the GG genotype is significantly associated with obesity susceptibility, aligning with prior research indicating that the ADA G22A polymorphism affects immune and metabolic pathways [43]. A study by Jadhav and Jain [44] reported elevated ADA activity in overweight and obese Indian individuals, suggesting a link between increased ADA activity and obesity. Additionally, Park et al. [45] investigated the A2a adenosine receptor (A2aAR) signaling in hepatic fibrosis and found that inhibiting A2aAR signaling reduced hepatic fibrosis and steatohepatitis in a mouse model, implicating adenosine pathways in liver disease progression. Given that obesity and insulin resistance are major risk factors for MAFLD, the association of the ADA G22A polymorphism with metabolic disorders such as gestational diabetes mellitus (GDM) and type 2 diabetes mellitus (T2DM) may suggest an indirect link to MAFLD [46,47].

Our findings reveal a significant association between IL-1RN polymorphisms, obesity, and MAFLD in children. These results are consistent with previous research on IL-1RN polymorphisms in metabolic disorders. Strandberg et al. [48] demonstrated that carriers of the IL-1RN*2 allele exhibit increased total fat, trunk fat, and serum leptin levels. Additionally, Mikhailova and Ivanoshchuk [49] reported that IL-1RN, an anti-inflammatory cytokine, is secreted by white adipose tissue, with serum levels correlating positively with BMI in obese individuals, indicating a link between IL-1RN expression and obesity severity.

Knockout studies in mice suggest that the absence of the IL1RN gene confers resistance to obesity, implicating IL1RN in fat accumulation and energy balance. Moreover, Bojarczuk et al. [50] highlighted that interleukin gene polymorphisms, including IL1RN, influence cytokine production and may impact energy balance, fat storage, and metabolic regulation, contributing to obesity susceptibility.

IL-1RN polymorphisms play a role in liver diseases such as MAFLD, as inflammatory cytokines, including those in the interleukin-1 family, contribute to the onset and progression of NAFLD, which shares features with MAFLD. Variations in IL-1RN may influence liver inflammation and disease progression [51]. A study by Maculewicz et al. [1] reported that interactions between IL-1 family genes (IL1B and IL-1RN) are linked to BMI and fat percentage, with specific alleles potentially protecting against higher BMI and fat accumulation, influencing MAFLD risk.

Regarding lipid ratios and liver disease severity, Nobili et al. [52] found that total cholesterol/HDL-C and LDL/HDL-C ratios correlate with liver injury severity in children with NAFLD, making them potential markers for disease progression.

The TG/HDL-C ratio has been explored as a marker for MAFLD in obese children. Iwani et al. [53] associated a higher TG/HDL-C ratio with increased insulin resistance and hepatic steatosis. Similarly, Demiral et al. [54] reported a significant link between TG/HDL-C ratios, insulin resistance, and hepatosteatosis in overweight and obese children. Additionally, Ting et al. [55] identified TG/HDL-C ratio as an independent predictor of liver fibrosis in pediatric NAFLD, emphasizing its utility in identifying children at risk for advanced liver disease.

The relatively small sample size, evaluation of the effect of possible early lifestyle interventions with lack of long-term follow-up of the included pediatric patients with MAFLD to assess the possible outcomes were the main limitations of the current work.

Conclusion

Children with MAFLD tend to have higher weight, BMI, central obesity, and significant metabolic disturbances, including higher insulin resistance, dyslipidemia, and liver enzyme elevations. The current study highlights a potential genetic association between the ADA G22A variant and childhood obesity, while also emphasizing the significant role of the IL-1RN SNP in the development of MAFLD among obese children. Additionally, the TG-to-HDL ratio appears to be a more effective predictor of pediatric MAFLD than the VLDL/HDL ratio. Understanding these genetic influences can help in early identification of high-risk individuals and the development of targeted therapeutic approaches to prevent and manage pediatric MAFLD more effectively.

Article information Continued