Abstract

The high prevalence of childhood obesity in Mexico is alarming in the health-science field. We propose to investigate the contribution of adipokines and cytokines polymorphisms and common BMI/obesity-associated loci, revealed in genome-wide association studies in Caucasian adult cohorts, with childhood obesity. This study included 773 Mexican-Mestizo children (5-15 years old) in a case-control study. The polymorphisms included were ADIPOQ (rs6444174), TNF-α (rs1800750), IL-1β (rs1143643), IL-6 (rs1524107; rs2069845), NEGR1 (rs34305371), SEC16B-RASAL2 (rs10913469), TMEM18 (rs6548238; rs7561317), GNPDA2 (rs16857402), LEP (rs2167270), MTCH2 (rs10838738), LGR4-LIN7C-BDNF (rs925946), BCDIN3D-FAIM2 (rs7138803), FTO (rs62033400), MC4R (rs11872992), MC4R (rs17782313), and KCTD15 (rs29942). No significant contribution was found with adipokines and cytokines polymorphisms in this study. Only both TMEM18 (rs6548238; rs7561317) polymorphisms were found associated with obesity (OR=0.5, P=0.008) and were in linkage disequilibrium (r2=0.87). The linear regression showed that the rs7561317 polymorphism of TMEM18 is negatively associated with obesity. This report highlights the influence of TMEM18 in Mexican-Mestizo children obesity, while adipokine and cytokine polymorphisms were not associated with it.

1. Introduction

Mexico has the highest prevalence of childhood obesity worldwide [1]. Children obesity is a predictive factor for adult obesity [2] and leads to the development of insulin resistance, diabetes, hypertension, cardiovascular disease, or cancer [3].

Children obesity develops in response to several factors, where nutrition and energy expenditure play a determinant role [4]. Adipose tissue is an endocrine organ that secretes hormones (adipokines) and proinflammatory cytokines [5]. Adiponectin is secreted by adipocytes with antidiabetic, antiatherogenic, and anti-inflammatory functions [6]. During obesity adiponectin decreases and serum levels are negatively correlated with body mass index (BMI) in prepubertal children [7]. An increase in proinflammatory cytokines like Tumor Necrosis Factor-alpha (TNF-α), interleukin-1 beta (IL-1β), and interleukin-6 (IL-6) suppresses adiponectin expression [8], which is regulated by the ADIPOQ gene, and single nucleotide polymorphisms (SNPs) in this gene have been associated with obesity-related traits [9].

Also, genetic factors also have a significant contribution to obesity as evidenced by the large disparities observed between ethnic groups. Results of large cohort studies showed that American Indian, Mexican, and Hispanic children have a higher risk to develop obesity than Africans children [10]. It has been estimated that obesity heritability goes from 6% to 85% among several populations. Genome Wide Association Studies (GWAS) have revealed several loci highly associated with body mass index (BMI), waist circumference, and body fat percentage (BF%) among populations [1113]. Fat mass and obesity-associated gene (FTO) has been widely linked to obesity among populations, and posterior replications have shown an additional signal in the melanocortin 4 receptor (MC4R) locus. From the Genetic Investigation of Anthropometric Traits (GIANT) consortium, additional loci were associated with obesity, like transmembrane protein 18 (TMEM18), potassium channel tetramerization domain containing 15 (KCTD15), glucosamine-6-phosphate deaminase 2 (GNPDA2), Src homology 2B (SH2B) family member 1 (SH2B1), mitochondrial carrier 2 (MTCH2), the melanocortin 4 receptor (MC4R), and the neuronal growth regulator 1 (NEGR1) genes [14]. Thus, the aim of this study was to evaluate the association between polymorphisms previously associated in GWAS in Caucasian adult population, adipokines, and proinflammatory cytokines with obesity in Mexican-Mestizo children.

2. Methodology

2.1. Study Design

A case of control study was conducted in 773 children ranging from 5 to 15 years old, from the children obesity cohort COIPIS (Cohorte de Obesidad Infantil-Proyecto Infancia Saludable) from Hospital Regional Lic, Adolfo López Mateos, ISSSTE (Instituto de Seguridad y Servicios Sociales de los Trabajadores del Estado) [15].

After parents signed an informed consent to include their children in this study, anthropometric measurements (weight, height, BMI, and BF%) were determined. BMI was calculated as weight (kg) divided by height (m) squared. Obesity was defined according to the World Health Organization (WHO) Child Growth Standard for children as BMI-for-age value over +2 SD [16] using the 95 percentile. The BF% was evaluated using a Body composition analyzer with an accurate stadiometer (InBody J10). A whole blood sample was taken to evaluate clinical parameters at fasting: glucose, glycated hemoglobin (HbA1c), creatinine, triglycerides (TAG), total cholesterol (CHOL), and high-density lipoprotein cholesterol (HDL).

An exclusion criterion was related to children with foreign parents and grandparents. This study was conducted in accordance with the Declaration of Helsinki and approved by the Hospital Regional Lic, Adolfo López Mateos Research, Ethics, and Biosafety Committees (registration number 245.2012).

2.2. DNA Extraction and Genotyping

Genomic DNA was obtained from 200 μl of the whole blood-EDTA with the InviMag Blood DNA Mini Kit, using an automated system (InviGenius, Stratec). For genotyping, we selected 18 SNPs (Table 1) which were amplified using predesigned 5′ exonuclease TaqMan genotyping assays on a 7500 series Real-Time PCR system (Applied Biosystems, Foster City, CA, USA).

3. Statistical Analysis

The results of continue variables are presented as mean ± standard deviation for Gaussian distribution and analyzed by the Student t test. For non-Gaussian distribution, the continuous variables are presented as median and interquartile range and analyzed by Mann-Whitney test. A chi-squared test was performed to estimate the Hardy-Weinberg equilibrium in controls. Pearson correlations (R) were used to assess linkage disequilibrium (LD) between SNPs in TMEM18, IL-6, and MC4R genes using Haploview® Software (Broad Institute, Cambridge, MA, USA). Statistical power of a case-control outcome design was carried out including 17% of childhood obesity prevalence [17], gene only and allele frequencies from 0.02 to 0.5. To evaluate the effect of the SNPs on children obesity, a logistic regression was used, including dominant and additive inheritance models. In addition, for further assessment a linear model was used to test the significant loci (resulted from logistic regression) on BMI, BM%, and clinical parameters. All statistical analyses were adjusted by gender, age, and BMI (in the case of logistic regression) and were performed using STATA12 (StataCorp, Texas). The statistical power calculation was performed using Quanto® Software (USC Biostats, California, USA). Data (.dta) used to support the findings of this study are available from the corresponding author upon request.

4. Results

Clinical and anthropometrical characteristics are summarized in Table 2. This study included 204 children with obesity and 569 children without obesity. BF%, glucose, creatinine, TAG, and CHOL levels were found higher in children with obesity. However, the levels of these clinical outcomes in the obesity group fall in normal ranges. Only HDL levels decreased in this group. There was no difference in HbA1c concentration between groups.

All the SNPs analyzed were in Hardy-Weinberg equilibrium. The allelic distribution of the SNPs is shown in Table 3. Only rs6548238 and rs7561317 SNPs in TMEM18 were associated in the dominant and additive inheritance models. Both polymorphisms confer significant protection against obesity. These polymorphisms were in linkage disequilibrium in our studied population (r2=0.87).

A marginal association was found between rs925946 in LGR4-LIN7C-BDNF and rs1800750 in TNF-α with obesity but does not reach statistical significance. In this study we did not identify any association between ADIPOQ or any cytokine polymorphism with obesity (Table 3).

To evaluate the effects of the associated polymorphisms on BMI and BF%, a linear regression was performed polymorphisms (Table 4). Both TMEM18 rs6548238 and rs7561317 loci were only associated BF% after adjustment by age and gender. However, the rs7561317 polymorphism was significantly associated with BMI. There was no association between these polymorphisms and clinical parameters.

5. Discussion

It is well recognized that adipokines and cytokines homeostasis plays an important role in obesity condition. Additionally, GWAS in Caucasian populations revealed many common loci associated with adult obesity [18]. Here, we analyzed the effect of several of these polymorphisms, as well as adipokines (ADIPOQ) and cytokines (TNF-α, IL-1β, and IL-6).

We observed a significant association between TMEM18 (rs7561317) and children obesity; however, we did not observe significant association between adipokines and obesity in our study. A meta-analysis of the influence of adipokine polymorphisms in adipokine genes (LEP, ADIPOQ, IL-1β, IL-6, and TNF-α) in obesity susceptibility showed that there was not association with the risk of obesity and LEP variants in adults. Interestingly, the results of this meta-analysis suggested that polymorphisms in the adipokine genes, ADIPOQ, IL-1β, IL-6, and TNF-α, increase the risk of obesity [19]. Our results are consistent with three studies in Mexican children that showed no contribution of several polymorphisms in ADIPOQ gene (rs2241766, rs182052, rs266729, and rs822393) in children obesity [9, 20, 21]. A similar effect was observed in TNF-α and interleukin-10 [22], suggesting that genetic variants in adipokines do not have a significant effect in Mexican-Mestizo children obesity.

Polymorphisms in FTO and MC4R have been associated with obesity in children. A recent meta-analysis which included 13 articles published between 2011 and 2015, comprising 15,613 participants, shows that the polymorphism rs9939609 in FTO was significantly associated with an increased risk of obesity in children and adolescents in several populations, except Mexicans [23]. Another recent report exposed that the rs9939609 polymorphism interacts with the Native American-origin polymorphism rs9282541, in the ATP-binding cassette transporter (ABCA1). Also, in the presence of the rs9282541 risk allele, the rs9939609 variant was not associated with BMI in Mexican adults, but this interaction was not observed in children [24]. This lack of significance was attributed to the to sample size; however the study of Mejía-Benitez (REF) showed no association between several polymorphisms in FTO with obesity in Mexican children, using a larger sample size. We propose that large cohorts during extensive periods could help to understand the real contribution of FTO variants in Mexican children.

The TMEM18 rs7561317 variant has been widely studied. This polymorphism, together with rs10913469 in SEC16B gene, was associated with adult obesity in a Japanese population [25] and Caucasians children [26]. On the contrary, the minor allele of TMEM18 rs7561317 was related to underweight in Czech adolescents [27]. This polymorphism has also been associated with BMI in Mexican children [28]. In our study, after adjustment by age, sex, and BMI, we observed the same association. Thus, this report indicates that TMEM18 rs7561317 has an effect on obesity development at all stages of life development.

6. Conclusions

Here we report no association of adipokines and cytokines polymorphisms with children obesity. However, the rs7561317 polymorphism in TMEM18 gene is associated with BMI and BF% in our sample of Mexican-Mestizo children.

Data Availability

Data (.dta) used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This research was supported by ISSSTE-E015, project number 245.2012.