Skip to main content

The influence of six polymorphisms of uncoupling protein 3 (UCP3) gene and childhood obesity: a case–control study

Abstract

Background

Obesity is defined as a multifactorial disease, marked by excessive accumulation of body fat, responsible for compromising the individual’s health over the years. The energy balance is essential for the proper functioning of the body, as the individual needs to earn and spend energy in a compensatory way. Mitochondrial Uncoupling Proteins (UCP) help in energy expenditure through heat release and genetic polymorphisms could be responsible for reducing energy consumption to release heat and consequently generate an excessive accumulation of fat in the body. Thus, this study aimed to investigate the potential association between six UCP3 polymorphisms, that have not yet been represented in ClinVar®, and pediatric obesity susceptibility.

Methods

A case–control study was conducted with 225 children from Central Brazil. The groups were subdivided into obese (123) and eutrophic (102) individuals. The polymorphisms rs15763, rs1685354, rs1800849, rs11235972, rs647126, and rs3781907 were determined by real-time Polymerase Chain Reaction (qPCR).

Results

Biochemical and anthropometric evaluation of obese group showed higher levels of triglycerides, insulin resistance, and LDL-C and low level of HDL-C. Insulin resistance, age, sex, HDL-C, fasting glucose, triglyceride levels, and parents’ BMI explained up to 50% of body mass deposition in the studied population. Additionally, obese mothers contribute 2 × more to the Z-BMI of their children than the fathers. The SNP rs647126 contributed to 20% to the risk of obesity in children and the SNP rs3781907 contribute to 10%. Mutant alleles of UCP3 increase the risk for triglycerides, total cholesterol, and HDL-C levels. The polymorphism rs3781907 is the only one that could not be a biomarker for obesity as the risk allele seem to be protective gains the increase in Z-BMI in our pediatric population. Haplotype analysis demonstrated two SNP blocks (rs15763, rs647126, and rs1685534) and (rs11235972 and rs1800849) that showed linkage disequilibrium, with LOD 76.3% and D’ = 0.96 and LOD 57.4% and D’ = 0.97, respectively.

Conclusions

The causality between UCP3 polymorphism and obesity were not detected. On the other hand, the studied polymorphism contributes to Z-BMI, HOMA-IR, triglycerides, total cholesterol, and HDL-C levels. Haplotypes are concordant with the obese phenotype and contribute minimally to the risk of obesity.

Peer Review reports

Background

Obesity is a global health issue, affecting people of different age groups, sexes, and economic status with two-thirds of the obese population living in developing countries [1, 2]. Increased adiposity is responsible for changes in body metabolism, fat mass distribution, and compromising several aspects of the individual health [3]. In addition, when obesity is a childhood outcome, it may be responsible for the development of life-threatening chronic diseases throughout their lives, which in turn contributes to increase public health expenditures adding a critical component to the financing of health systems [4].

Obesity is a complex multifactorial disease. Several factors, including lifestyle such as eating habits and social behaviors, metabolic individual heterogeneity, psychological disorders, neurological changes, and individual genetic makeup are involved in the origin and maintenance of this disease throughout life [5, 6]. It has been estimated that by the year 2025, approximately 2.3 billion adults worldwide will be overweight and that 1/3 of this number will be severely obese. In 2020, 39 million children under the age of 5 were overweight or obese [7]. Regarding pediatric obesity in Brazil, around 12.9% of children between 5 and 9 years old and 7% of the adolescents aged between 12 and 17 are already obese [8].

The fundamental cause underlying an obese phenotype is an imbalance between the amount of energy ingested compared to the amount of fat consumed by the body [5, 6, 9]. All consumed energy maintained in adipocytes is metabolized in the production of ATP or, in cases of excess, is alternatively eliminated by thermoregulation, producing heat. This burning through thermoregulation is only possible due to the presence of mitochondrial uncoupling proteins (UCP), formerly known as thermogenic. There are five isoforms of UCP in mitochondrial ridges, namely UCP1 to UCP5 [10].

The genes encoding the human UCP are studied in association with obesity, as several gene variants have been related to hypertension susceptibility, diabetes development, and cardiovascular diseases [11]. The UCP3 gene, located at 11q13.4, has seven exons intervened by six intronic sequences. It is mainly expressed in both skeletal muscle and brown adipose tissue and the protein is related to cellular fatty acid metabolism. Upregulated expression of UCP3 has already been related to physical activity, fasting, and high-fat diet. On the other hand, downregulation has been reported to improve fat oxidative capacity [12, 13]. When UCP3 activity is reduced, a decrease in energy expenditure is observed, while an increase in its expression is correlated with an increase in the metabolic rate, consequently, contributing to decrease BMI. The physiological function of UCP3 is yet unclear, nevertheless, it is found abundantly in brown adipose tissue (BAT) and may interplay with UCP1 in the regulation of BAT. Although several research groups have attempted to link UCP3 genomic variations with obesity, susceptibility in human and animal models remain conflicting and current trends show an increasing association of UCP3 with obesity and diabetes [14,15,16]

The gene variants most frequently used to investigate potential association with obesity susceptibility are Single Nucleotide Polymorphisms (SNP). To date, the risk alleles of about 750 loci of well recognized SNP with minimal allele frequencies as small as 1.6% have been contributing per-allele effects as low as 0.04 kg/m2 per allele [17]. Moreover, heritable factors appear to be responsible for up to 50% of the variation in body adiposity in non-syndromic polygenic obesity [6]. Different studies have been developed to observe the influence of genetic polymorphisms in the reduction of energy expenditure and, consequently, in the increase of fat deposition in tissues. Advancement on genome sequencing and genotyping has made it possible to investigate the genetic influence on the regulation of body fat. Thus, genomic variants have been associated with as high as 6% surplus on BMI when people carry the risk alleles [18]. Among UCP3 polymorphisms, the frequently studied SNP are rs15763, rs1800849, rs647126, rs7930460, rs1685356, rs1685354, rs11235972, and rs378190719.

In the current study, a total of six UCP3’s SNP that have not yet been represented in ClinVar® – a public repository which aggregates information about genotype–phenotype correlation and its relationship to human health – were used to investigate the potential association between UCP3 genomic variation and pediatric obesity susceptibility. From the 6 SNP, 4 were categorized as mutations affecting gene transcription, namely rs1800849 (upstream variant), rs15763, rs647126, and rs1685354 (downstream variants) and 2 were in the class of intron variants, namely rs3781907 and rs11235972. The study group comprised of eutrophic and obese children from Central Brazil.

Methods

Study design and participants

The present observational case–control study included 225 children and adolescents, aged 5 to 19 years, of both sexes from Central Brazil. The participants were distributed as eutrophic or obese, corresponding to 102 and 123 people in each group, respectively. The study was performed following the ethical guidelines of the Declaration of Helsinki and was approved by the Ethics Committee on Human Research from the Pontifical Catholic University of Goiás under the registered number CAAE: 16303313.4.0000.0037. Written informed consent was voluntarily signed by parents/legal guardians of all participants.

Participants were recruited from the Endocrinology Service at the Children’s Hospital (Goiânia, Goiás, Brazil). Laboratory analyses were carried out at the genetic laboratory of Replicon Research Group/School of Medical and Life Sciences from the Pontifical Catholic University of Goiás and Mutagenesis Laboratory from Federal University of Goiás. Malnutrition, overweight, severe chronic diseases, genetic syndromes, and the use of medications that affected body weight and mass deposition were defined as exclusion criteria.

Anthropometric measurements including weight, height, and body mass index (BMI) of all participants and their biological parents were obtained by a pediatric endocrinologist with the aid of a fixed Tonelli E150-A stadiometer and a WELMY mechanical scale certified by the National Institute of Metrology, Quality and Technology (Inmetro). Inmetro is a governmental accreditation body responsible to maintaining the quality of instruments and pieces of equipment used in Brazil. The physician also performed a complete physical examination of the children and their biological parents. For all children, Z-scores of the BMI (Z-BMI) were calculated and interpreted according to the recommendation of the World Health Organization (WHO) [18]. Children were grouped according to the weight status as follows: obese Z-BMI ≥  + 2SD and eutrophic Z-BMI between < -2SD and + 1SD.

Biochemical assays

Biochemical tests were performed in an accredited laboratory in Goiânia, Goiás, Brazil. Blood samples were drawn by peripheral venipuncture following a period of overnight fasting ranging from 8 to 12 h. The assays included fasting plasma glucose, lipid panel, and insulin levels. Total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and glucose levels were run in an Abbott Architect c8000® Chemistry Analyzer. Insulin levels were determined by immunoassay using an automated Abbott Architect i2000®.

The reference values for fasting blood glucose are the same for adults and children. Thus, normoglycemic lays within the range of 3.3 and 5.6 mmol/L, pre-diabetics between 5.6 and 6.9 mmol/L and diabetics ≥ 7 mmol/L. Normal fasting insulin level is defined by the ranges 2.5 to 25 mIU/mL and < 15 mIU/mL for adults and children, respectively [19, 20]. To interpret the lipid panel, the recommendations of the Brazilian Society of Pediatrics [21] was used.

To assess insulin resistance in children, the HOMA-IR (Homeostatic Model Assessment Insulin Resistance) score was obtained using the equation: fasting serum insulin (μU/mL) x fasting plasma glucose (mmol/L)/22.5 [22]. The values obtained by calculating HOMA-IR, for the Brazilian child population, were interpreted through the parameters defined by [23].

Genotyping UCP3 polymorphisms

Genomic DNA was isolated from the peripheral blood of eutrophic and obese children using a commercial kit (AxyPrep™ BloodGenomic DNA Miniprep Kit, Axy353.gen Scientific, USA), following to the manufacture’s laboratory protocol. The concentration of genomic DNA was determined using a spectrophotometer (NanoVue Plus®, GE Healthcare, USA). Genotyping of the six SNP (rs15763, rs1685354, rs1800849, rs11235972, rs647126, and rs3781907) (Table 1) was carried out using qPCR (TaqMan® Real Time PCR in StepOnePlus™ thermocycler, Thermo Fisher Scientific, USA), following the manufacturer’s instructions.

Table 1 Probes sequences used for genotyping the UCP3 gene polymorphisms in obese and eutrophic children from Central Brazil

Statistical analyses

Analyses were carried out in the SPSS® version 25.0 (IBM Corporation, USA) and Excel® version 16.62 (Microsoft Corporation, USA) utilities. In the current study, a significance level was set at 5% (p ≤ 0.05), considering confidence interval (IC) of 95%. Shapiro–Wilk’s and Levene’s tests were applied to check for normality and homogeneity of variance between groups before proceeding with additional analyses. Student’s t-test was used to compare the means of continuous variables between the studied groups. Pearson’s chi-square test (X2) was used to assess the deviations from Hardy‐Weinberg equilibrium and potential association between categorical variables, including the distribution of genotypes as a function of the risk allele between the obese and eutrophic groups. On the other hand, the standard McNemar test was used to determine if there were differences on the proportions of patients having wildtype alleles compared to the mutated genotype (risk alleles).

Multiple Linear Regression (MLR) was applied to observe whether the qualitative dependent variables could predict the outcomes of continues variables for the entire group. For the MLR the stepwise method was applied, and the equation was generated using beta standardized coefficients. Multicollinearity statistics was considered relevant to negatively affect the MLR model if tolerance was ≤ 20%. To test the assumptions’ compliances for MLR, the normality of the residues was verified according to the P-P plot graphs and the homoscedasticity by scatterplot. Furthermore, multinomial logistic regression (NR) was used to model the contribution of UCP3 genotypes to the outcome variables. In SPSS, the dependent variable was ascribed to the aggregated genotypes containing the risk allele while the outcome variables (covariates) were Z-BMI, HOMA-IR, triglycerides, total cholesterol, and HDL-C. As the reference category, fasting glucose was used as their means were not statistically different between eutrophic and obese participants. Therefore, it would be possible to test the association for a dominant model, where the presence of a risk allele would contribute to modulate the outcome variables in the whole group of participants.

Haplotype block analysis was completed in Haploview® (Broad Institute) to test potential linkage disequilibrium (LD) for the studied SNP, using Gabriel’s confidence interval method [24]. The function multimarker haplotype tests block-specified, following the option to use logistic regression (GML), was performed for haplotype-based association analysis in Plink [25].

Results

A total of 225 children and adolescents were enrolled in the study, including 123 (54.7%) and 102 (45.3%) participants in the obese and eutrophic groups, respectively. Regarding sex distribution, among the obese children, 60 (48.8%) were female and 63 (51.2%) were male, while in the eutrophic group, 52 (51%) were female and 50 (49%) males. Table 2 summarizes the demographic, anthropometric, and biochemical variables for both eutrophic and obese participants.

Table 2 Sociodemographic, anthropometric, and biochemical parameters for the group of obese and eutrophic children included in the study investigating the potential contribution of UCP3 SNP variants on the susceptibility of pediatric obesity

It was observed that the mean values for weight, height, Z -BMI, mother’s and father’s, BMI, insulin levels, HOMA-IR, triglycerides, HDL-C, and LDL-C between obese and eutrophic groups were statistically different (p < 0.05). In the group of obese children, the mean value for HOMA-IR was 2.2 × higher than the mean of the eutrophic group, indicating insulin resistance among the obese children. Moreover, for the group of pediatric obesity, the mean of all variables revealed increased values but HDL-C, which showed an expected decrease. No statistical difference was observed for fasting glucose among the two groups. Table 2 summarizes those findings.

With respect to the normal distribution of data in the current study, few outliers were constant among all variables tested. However, Q-Q plots showed that skewness were observable only for the more extreme values related to outliers while most of the values for each variable remained adhered to the line of greatest fit as it would be expected if they truly came from a normal distribution. Consequently, because the data didn’t show heavy-tailed traits and because of the relatively large sample size of each group, non-normality presented no problems in hindering invalid further testing [26].

MLR was used to test if demographic, anthropometric, and biochemical independent variables could significantly predict the dependent Z-BMI for the whole pediatric population, consisting of obese and eutrophic participants. At first, it was identified and fixed multicollinearity among the predictors, which resulted in the removal of insulin (8,6%), height (13.4%), weight (14.7%), total cholesterol (15.6), and LDL-C (16.3%) from the model. Thus, the MLR results from the whole set of participants showed a statistically significant model (F = 347,4; p < 0.0001; R2 = 0.486), indicating that the predictors included in the model account for about 50% in the variability of Z-BMI in a pediatric population. The model only included in the equation variables that reached statistical significance:

$$Z-BMI=2.37+0.37HOMA-IR+0.3BMlm-0.25Age-0.18HDL-C+0.15Sex+0.13BMIf-0.12Glucose+0.04TG$$

The Pearson’s chi-square test showed the observed genotype distributions for the six SNP studied were not statistically significant different (p > 0.05) between obese and eutrophic groups. The chi-square test also showed the frequencies of the genotypes in the studied population complied with the Hardy–Weinberg equilibrium (p > 0.05). When evaluating the genotypic distribution as a function of the risk allele, no statistical differences were observed (Table 3). On the other hand, the McNemar test was applied to assess whether there was a difference between the proportions related to the presence of the risk allele by the obese phenotype. A statistically significant difference (p < 0.05) between the proportions of individuals who presented the risk allele and were eutrophic in relation to the obese individuals who presented wild-type alleles, with exception of the polymorphism rs647126 (Table 3).

Table 3 Allele and genotype distributions of six SNPs of UCP3 of eutrophic and obese children from a pediatric population of Central Brazil

A nominal logistic regression (NLR) was performed to verify if the genotypes aggregated for the risk alleles in the UCP3 gene were good predictors for the increase in the outcome variables. The multinomial logistic model was statistically significant (P ≤ 0,0001) for all six SNP included in the study. Table 4 includes all predictors that reached statistical significance in the model along with the betas and respective odds ratios. All SNP variants in UCP3 were implicated in predicting the outcome variables in a pediatric population. Two SNP variants were able to predict positive changes in Z-BMI (rs647126 and rs3781907), two variants associated with positive changes to HOMA-IR (rs647126, rs15763) while, on the other hand, two had a negative impact on HOMA-IR (rs3781907 and rs1685354). Additionally, positive changes in total cholesterol could be predicted by two variant genotypes (rs1800849 and rs11235972) and one SNP associated with positive increase of HDL (rs3781907). Finally, two SNP variants associated with positive changes in triglycerides plasma concentrations.

Table 4 Logistic regression analysis when genotype was used as an independent variable and Z-BMI was taken as a dependent variable of participants of the study regarding the contribution of UCP3 variants to obesity susceptibility in a pediatric cohort from Central Brazil

With respect to haplotype frequencies, Fig. 1 showed the haplotypes on UCP3 gene at 11q13.4 genotyped in association with pediatric obesity in a cohort of children and teens from Central Brazil. The first three SNPS (rs15763, rs647126, and rs1685354) were linked together in a separate haplotype block with a pairwise D’ value of 0.96, while a second block was inferred for SNP rs11235972 and rs1800849 with a pairwise D’ value of 0.97. Table 5 shows the results of plink/GML haplotype.

Fig. 1
figure 1

LD heatmap of region chr11: 74,000,400- 74,009,200, spanning over 8,800 bp, containing 6 SNPs in the UCP3 gene. Figure Legend: The map depicts association signals achieved in a cohort of eutrophic and obese pediatric population from Central Brazil. Bright red signals depict pairwise linkage disequilibrium in the UCP3 gene using six marker SNPs that associated with obesity in a pediatric cohort, indicating D’ ≥ 0.85 and LOD ≥ 3

Table 5 Haplotype-based association analysis using Plink

Discussion

Understanding pediatric non-syndromic obesity requires the comprehension of the biological and environmental factors related to its development. Only after appreciating all the factors and covariables of such important heterogeneous, multifactorial, and complex disease that risk management could be planned to mitigate the health consequences for the susceptible population. For pediatric non-syndromic obese pollutions, clinical management has largely focused on metabolic concepts to cause or prevent severe outcomes. Scientific evidence has suggested that about 80% of obese children will become obese adults [27].

Genome-wide association studies (GWAS) with adult BMI and adiposity have associated many loci and thousands of variants with this phenotype. However, the contribution and role of those variants towards pathogenicity of obesity in adult and pediatric populations remains unanswered. Moreover, much less is known about the genetic variants associated with childhood obesity in the Brazilian population, well recognized for having extensive genetic admixture [28, 29]. The present study, although based on a discreet populational sample size, helps to contribute to the bulk of knowledge regarding the genetic influence of UCP3 variants on the regulation of body fat deposition in an admixed pediatric population from Central Brazil.

In the current study, the sex of participants was not used as inclusion/exclusion criteria. Thus, random acceptance of both sexes led to the observation that pediatric obesity affected equally boy and girls. However, females seemed to be more susceptible to have larger Z-IMC than males as predicted by the Eq. 1.

Our findings showed a set of covariables that better explains the variability of the Z-BMI in the pediatric cohort. Insulin resistance, age, sex, HDL-C, fasting glucose, triglyceride levels, and parents’ BMI explained up to 50% of body mass deposition in the studied population. One relevant observation regarding the contribution of parent’s BMI to Children Z-IMC favors previous observations that the combination of genetic, environmental, and lifestyle factors, such as the obesogenic environment a child is exposed to, contributes to the obese phenotype and, consequently, undesired health outcomes [30, 31]. Additionally, obese mothers contribute 2 × more to the Z-BMI of the children than their fathers. With respect to risk mitigation, it’s important to clinicians and health support professionals to understand that complex traits are dependent on both genotypes and environment. Thus, children health outcomes are influenced by their parent’s habits during child growth and if a child shares genomic variants that even discretely increase their chances for obesity, wight managing must include a family plan. Being important to address and, whenever possible, modify parenting and lifestyles, habits, routine physical workouts, and modeling the home environment to reduce the risk of pediatric obesity [32,33,34]. It has been already postulated that when both parents are obese, the chance of a child developing obesity during childhood reaches 80%, however, when only one parent has the phenotype, the chance drops to 9% [35, 36].

The group of obese children showed higher levels of triglycerides, insulin resistance, and LDL-C and low level of HDL-C. Insulin resistance is related to the decrease of insulin capacity to uptake glucose [37] and it may contribute to becoming a risk group for the development of metabolic syndrome in children [38]. The accumulation of adipose tissue is responsible for triggering various physiological responses in the body, such as insulin resistance, increased cholesterol, and elevated blood pressure in those having the phenotype [39, 40]. Furthermore, adipose tissue is responsible for directly influencing the metabolism of glucose in the body [37,38,39,40,41] and is often related to Type 2 Diabetes Mellitus (T2DM) [42]. Therefore, special attention must be given to obese children due to the risk of developing non-communicable diseases leading to undesired and life-threatening comorbidities.

In general, genomic variants in the UCP3 gene are responsible for interfering in energy homeostasis and related to cellular fatty acid metabolism. In the current study, the studied UCP3 polymorphisms showed the variants were not causative of pediatric obesity. The genotypic distribution of all 6 SNPS agreed with the ALFA database for the Latin America available at dbSNP/NCBI. On the other hand, the proportions between the risk allele were statically different between eutrophic and obese participants, except for the SNP rs647126. Thus, although the risk alleles were not causative, they contributed to the phenotype. The SNP rs647126 and rs3781907 did contribute to the increase of Z-BMI in our pediatric populations. The first contributed to about 20% to the risk of obesity in children while the latter contribute around 10%.

More importantly, the association between UCP3 variants and HOMA-IR was found for 4 UCP3 SNP variants, namely rs647126, rs15763, rs3781907, and rs1685354, suggesting another important contribution to the UCP3 in obesity susceptibly due to insulin resistance and, consequently, higher blood sugar levels, that, in turn, it is a well stablished risk for T2DM. The findings from our study leads to future experimental designs to further understand the role of UCP3 in the development of T2DM. Moreover, increases risk for triglycerides, total cholesterol, and HDL-C levels were also associated with UCP3 mutant alleles. Although the increments were discrete, future studies would be welcome to shed some light into the role of UCP3 and cholesterol metabolism in humans (Table 4).

The association of these polymorphisms with overweight, obesity, and BMI has remained controversial. Different studies suggest the wild genotype of rs1800849 polymorphism is related to high levels of HDL-C, reduced chance to develop T2DM and decrease of BMI in adults [43,44,45,46]. A study showed no association between rs1800849 and BMI among Caucasian southern Italian children [10]. A study with Mexican adolescents did not show an association with BMI but showed a statistically significant difference with a lower waist circumference [47]. On the other hand, two other studies using Brazilian and French cohorts of adults demonstrated association of mutant alleles and T2DM [34, 48]. A Chinese study with obese adolescents did not find a significant association with rs1800849 gene polymorphisms and obesity [45].

Less is known about the rs11235972, rs647126, and rs15763 polymorphisms. A study by van Abeelen et al. reported an association between the mutant allele of the rs11235972 polymorphism and increased levels of CT and LDL-C [49]. The rs647126 polymorphism has been correlated with the increase of BMI in adult populations [49, 50]. The previous study did not observe a significant association between the rs15763 polymorphism and weight gain [51]. Additionally, there were reports of the relationship between mutant alleles of the rs1685354 polymorphism with the increase in BMI and consequently in the propensity of obesity development [49].

In the current study, rs3781907 polymorphism was the only one that could not be a biomarker for obesity as the risk allele seem to be protective against the increase in Z-BMI in a pediatric population (Table 4), despite the significant differences on the proportions of obese children having the risk alleles when compared to the wildtype genotype. Another study observed the subjects with the rs3781907-G allele had higher levels of CT and LDL-C and a higher risk of T2DM [51]. On the other hand, another study with a Chinese population did not associate this polymorphism with T2DM [50] nor with essential hypertension [52]. Further studies with obese children are still needed to elucidate the role of rs3781907 as a biomarker of obesity or related comorbidities in both pediatric and adult populations.

Finally, Fig. 1 depicts the association for all the genotyped UCP3 SNP in association with pediatric obesity in a highly admixed Brazilian pediatric population. The X axis indicates the relative physical position of each SNP within the gene. The haplotype block map for the whole span of the UCP3 gene, showing pairwise LD in D′, could be found in Table 6, which also included LOD and r2. Haplotypes and haplotype frequencies were shown in Table 7 for each defined haplotype block. The SNP blocks 1 (rs15763, rs647126, and rs1685534) and 2 (rs11235972 and rs1800849) showed linkage disequilibrium, indicating deviation of haplotype frequencies from expected values based on genotype frequencies. Thus, evidencing and observed degree of concordance of a UCP3 variants with an obesity in a cohort of children from Central Brazil.

Table 6 Linkage disequilibrium scores to estimate the degree of concordance of UCP3 SNP and pediatric obesity form a cohort of Central Brazil
Table 7 Haplotype frequencies from 6 SNP distributed within the UCP3 gene in a cohort of participants of a study on pediatric obesity from Central Brazil

The discrete population size may be considered a limitation for the current study, even though participants were recruited following precise selection criteria for both obese and eutrophic groups. At the same time, the lack of research to understand the biological role of UCP3 genetic variants on cohorts of obese Brazilian children make the current study innovative in the field. Finally, another limitation of our study relates to the lack of assessment regarding the potential obesogenic environment shared by the participant children and their parents.

In conclusion, the present study enrolling a cohort of admixed Brazilian children could not detected causality between the six UCP3 SNP and obesity. However, the studied polymorphisms contribute to Z-BMI, HOMA-IR, triglycerides, total cholesterol, and HDL-C levels. Here we reported the presence of the risk alleles increased the chances of children having obesity and related comorbidities. Due to controversial findings about the role of UCP3 polymorphisms and obesity, further studies with a larger sample size and participants of distinct ethnic groups are required to confirm the potential of UCP3 variants to obesity in humans.

We also highlighted the importance of management and follow-up of obese children to prevent metabolic disorders and comorbidity complications, such hypertension, T2DM, and cardiovascular diseases. The influence of the parents’ BMI on the Z-BMI of obese children reinforces the need for a careful look of this health problem under the assumption that childhood obesity is related to both genetic and environment factors. Thus, non-syndromic obesity should be managed and surveilled as a familial disease. The mitigation of deleterious health effects requires a holistic approach to guide the family to reduce their chance of undesired obesity outcomes. Unfortunately, however, for most non-syndromic obesity cases, a prognosis or a specific medical treatment/conduct protocols based on the contribution of multiple low risk alleles is not yet available, largely, due to the lack of knowledge regarding the multiple gene interactions and pleiotropic effect subjacent to complex traits, such as human obesity.

Availability of data and materials

The data used in this study are available from the corresponding author and could be requested by email. All SNP information are available at the links: https://www.ncbi.nlm.nih.gov/snp/rs1800849; https://www.ncbi.nlm.nih.gov/snp/rs15763; https://www.ncbi.nlm.nih.gov/snp/rs3781907; https://www.ncbi.nlm.nih.gov/snp/rs647126; https://www.ncbi.nlm.nih.gov/snp/rs1685354; https://www.ncbi.nlm.nih.gov/snp/rs11235972.

Abbreviations

UCP3:

Uncoupling Protein 3

DM2:

Type 2 Diabetes Mellitus

SNP:

Single Nucleotide Polymorphism

BMI:

Body Mass Index

OR:

Odds Ratio

HOMA-IR:

Homeostatic Model Assessment Insulin Resistance

TC:

Total cholesterol

TG:

Triglycerides

HDL-C:

High-density lipoprotein cholesterol

LDL-C:

Low-density lipoprotein cholesterol

References

  1. Sámano R, Huesca-Gómez C, López-Marure R, Hernández-Cabrera AK, Rodríguez-Ventura A, Tolentino M, et al. Association between UCP polymorphisms and adipokines with obesity in Mexican adolescents. J Pediatr Endocrinol Metab. 2018. https://doi.org/10.1515/jpem-2017-0262.

    Article  PubMed  Google Scholar 

  2. McPherson K. Reducing the global prevalence of overweight and obesity. Lancet. 2014. https://doi.org/10.1016/S0140-6736(14)60767-4.

    Article  PubMed  Google Scholar 

  3. Who.World Health Organization. Obesity: Preventing and managing the global epidemic. In: Report of a WHO Consultation on Obesity. Geneva: WHO; 1998. http://www.who.int/nutrition/publications/obesity/WHO_TRS_894/en/. Accessed 15 Jan 2021.

    Google Scholar 

  4. Wannmacher, L. Obesidade como fator de risco para morbidade e mortalidade: evidências sobre o manejo com medidas não medicamentosas. In: Uso Racional de Medicamentos: fundamentação em condutas terapêuticas e nos macroprocessos da Assistência Farmacêutica. Organização Pan-Americana da Saúde/Organização Mundial da Saúde (OPAS/OMS) no Brasil. 2016;1:7. https://www.paho.org/bra/dmdocuments/Fasciculo%207.pdf. Accessed 16 Jan 2021.

  5. Grossman DC, Bibbins-Domingo K, Curry SJ, Barry MJ, Davidson KW, Doubeni CA, et al. Screening for obesity in children and adolescents: US preventive services task force recommendation statement. JAMA. 2017. https://doi.org/10.1001/jama.2017.6803.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Kumar S, Kelly AS. Review of childhood obesity: from epidemiology, etiology, and comorbidities to clinical assessment and treatment in Mayo Clinic Proceedings. Elsevier. 2017. https://doi.org/10.1016/j.mayocp.2016.09.017.

    Article  Google Scholar 

  7. Who – World Health Organization. Obesity and overweight. 2021. https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight. Accessed 06 Aug 2022

    Google Scholar 

  8. Santos FL, Esteves SS, da Costa Pereira A, Yancy WS Jr, Nunes JP. Systematic review and meta-analysis of clinical trials of the effects of low carbohydrate diets on cardiovascular risk factors. Obes Rev. 2012. https://doi.org/10.1111/j.1467-789X.2012.01021.x.

    Article  PubMed  Google Scholar 

  9. Oussaada SM, van Galen KA, Cooiman MI, Kleinendorst L, Hazebroek EJ, van Haelst MM, et al. The pathogenesis of obesity. Metabolism. 2019. https://doi.org/10.1016/j.metabol.2018.12.012.

    Article  PubMed  Google Scholar 

  10. de Souza BM, Brondani LA, Boucas AP, Sortica DA, Kramer CK, Canani LH, ... & Crispim, D. et al. Associations between UCP1 -3826A/G, UCP2 -866G/A, Ala55Val and Ins/Del, and UCP3 -55C/T polymorphisms and susceptibility to type 2 diabetes mellitus: case-control study and meta-analysis. Plos One. 2013. https://doi.org/10.1371/journal.pone.0054259.

  11. Jin F, Qiao C. Association of maternal caffeine intake during pregnancy with low birth weight, childhood overweight, and obesity: a meta-analysis of cohort studies. Int J Obes. 2020. https://doi.org/10.1038/s41366-020-0617-4.

    Article  Google Scholar 

  12. Musa CV, Mancini A, Alfieri A, Labruna G, Valerio G, Franzese A, Buono P, et al. Four novel UCP3 gene variants associated with childhood obesity: effect on fatty acid oxidation and on prevention of triglyceride storage. Int J Obes. 2012. https://doi.org/10.1038/ijo.2011.81.

    Article  Google Scholar 

  13. Alonso A, Martí A, Corbalán MS, Martínez-González MA, Forga L, Martínez JA. Association of UCP3 gene–55C> T polymorphism and obesity in a Spanish population. Ann Nutr Metab. 2005;49(3):183–8. https://doi.org/10.1159/000086883.

    Article  CAS  PubMed  Google Scholar 

  14. Pohl EE, Rupprecht A, Macher G, et al. Important trends in UCP3 investigation. Front Physiol. 2019. https://doi.org/10.3389/fphys.2019.00470.

    Article  PubMed  PubMed Central  Google Scholar 

  15. de Luis D, Domingo JC, Izaola O, Casanueva FF, Bellido D, Sajoux I. Effect of DHA supplementation in a very low-calorie ketogenic diet in the treatment of obesity: a randomized clinical trial. Endocrine. 2016. https://doi.org/10.1007/s12020-016-0964-z.

    Article  PubMed  Google Scholar 

  16. Hesselink MK, Mensink M, Schrauwen P. Human uncoupling protein-3 and obesity: an update. Obes Res. 2003. https://doi.org/10.1038/oby.2003.192.

    Article  PubMed  Google Scholar 

  17. Loos RJ, Yeo GS. The genetics of obesity: from discovery to biology. Nature Rev Gen. 2022;23(2):120–33. https://doi.org/10.1038/s41576-021-00414-z.

    Article  CAS  Google Scholar 

  18. Who – World Health Organization. Expert Comitee on Physical Status: the use and interpretation of antropometryphisical status. Geneva: World Health Organization. 1995. https://apps.who.int/iris/bitstream/handle/10665/37003/W?sequence=1. Accessed 15 Jan 2021

    Google Scholar 

  19. Col Araz N, Nacak M, OguzkanBalci S, Benlier N, Araz M, Pehlivan S, Aynacioglu AS, et al. Childhood obesity and the role of dopamine D2 receptor and cannabinoid receptor-1 gene polymorphisms. Genet Test Mol Biomarkers. 2012. https://doi.org/10.1089/gtmb.2012.0244.

    Article  PubMed  Google Scholar 

  20. SBD - Sociedade Brasileira de Diabetes. Diretrizes da Sociedade Brasileira de Diabetes: 2019–2020. http://www.saude.ba.gov.br/wp-content/uploads/2020/02/Diretrizes-Sociedade-Brasileira-de-Diabetes-2019-2020.pdf. Accessed 23 Oct 2021.

  21. SBP – Sociedade Brasileira de Pediatria. Dislipidemia na criança e no adolescente, orientações para o pediatra. https://www.sbp.com.br/fileadmin/user_upload/22336c-GPA_-_Dislipidemia_Crianca_e_Adoles.pdf. Accessed 29 Jun 2022.

  22. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and betacell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985. https://doi.org/10.1007/BF00280883.

    Article  PubMed  Google Scholar 

  23. de Almeida CA, Pinho AP, Ricco RG, Pepato MT, Brunetti IL. Determination of glycemia and insulinemia and the homeostasis model assessment (HOMA) in schoolchildren and adolescents with normal body mass index. J PEDIATRIC. 2008. https://doi.org/10.1590/S0021-75572008000200008.

    Article  Google Scholar 

  24. Barrett JC, Fry B, Maller JDMJ, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5. https://doi.org/10.1093/bioinformatics/bth457.

    Article  CAS  PubMed  Google Scholar 

  25. Weeks JP. plink: An R package for linking mixed-format tests using IRT-based methods. J Stat Softw. 2010. https://doi.org/10.18637/jss.v035.i12.

    Article  Google Scholar 

  26. Lumley T, Diehr P, Emerson S, Chen L. The importance of the normality assumption in large public health data sets. Annu Rev Public Health. 2002. https://doi.org/10.1146/annurev.publheath.23.100901.140546.

    Article  PubMed  Google Scholar 

  27. Gepstein V, Weiss R. Obesity as the main risk factor for metabolic syndrome in children. Front Endocrinol. 2019. https://doi.org/10.3389/fendo.2019.00568.

    Article  Google Scholar 

  28. Bradfield JP, Vogelezang S, Felix JF, Chesi A, Helgeland Ø, Horikoshi M, Struan FA, et al. A trans-ancestral meta-analysis of genome-wide association studies reveals loci associated with childhood obesity. Hum Mol Genet. 2019. https://doi.org/10.1093/hmg/ddz161.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zandona MR, Sangalli CN, Campagnolo PDB, Vitolo MR, Almeida S, Mattevi VS, et al. Validation of obesity susceptibility loci identified by genome-wide association studies in early childhood in South Brazilian children. Pediatr Obes. 2017;12(1):85–92. https://doi.org/10.1111/ijpo.12113.

    Article  CAS  PubMed  Google Scholar 

  30. Hesketh K, Crawford D, Salmon JO, Jackson M, Campbell K. Associations between family circumstance and weight status of Australian children. Int J Pediatr Obes. 2007. https://doi.org/10.1080/17477160601148554.

    Article  PubMed  Google Scholar 

  31. Quattrin T, Roemmich JN, Paluch R, Yu J, Epstein LH, Ecker MA. Treatment outcomes of overweight children and parents in the medical home. Pediatrics. 2014. https://doi.org/10.1542/peds.2013-4084.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Ochoa A, Berge JM. Home environmental influences on childhood obesity in the Latino population: a decade review of literature. J Immigr Minor Health. 2017. https://doi.org/10.1007/s10903-016-0539-3

  33. Yang M, Xu Y, Liang L, Fu J, Xiong F, Liu G, Zhu Y, et al. The effects of genetic variation in FTO rs9939609 on obesity and dietary preferences in Chinese Han children and adolescents. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0104574.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, Powell C, Vedantam S, Buchkovich ML, Yang J, Croteau-Chonka DC. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206. https://doi.org/10.1038/nature14177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Börjeson M. The aetiology of obesity in children A Study of 101 Twin Pairs. Acta paediatrica. 1976;65(3):279–87. https://doi.org/10.1111/j.1651-2227.1976.tb04887.x.

    Article  Google Scholar 

  36. Wisniewski AB, Chernausek SD. Gender in childhood obesity: family environment, hormones, and genes. Gend Med. 2009. https://doi.org/10.1016/j.genm.2008.12.001.

    Article  PubMed  Google Scholar 

  37. DeFronzo RA, Tripathy D. Skeletal muscle insulin resistance is the primary defect in type 2 diabetes. Diabetes Care. 2009. https://doi.org/10.2337/dc09-S302.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Zafar U, Khaliq S, Ahmad HU, Manzoor S, Lone KP. Metabolic syndrome: an update ondiagnosticcriteria, pathogenesis, andgenetic links. Hormones. 2018. https://doi.org/10.1007/s42000-018-0051-3.

    Article  PubMed  Google Scholar 

  39. Oliveira LCD, West LEM, Araújo EA, Brito JS, Nascimento CL. Prevalência de adiposidade abdominal em adultos de São Francisco do Conde, Bahia, Brasil, 2010. Epidemiologia e Serviços de Saúde. 2015. https://doi.org/10.5123/S1679-49742015000100015.

    Article  Google Scholar 

  40. Morais JB, Severo JS, Beserra JB, Oiveira AR, Cruz KJ, de Sousa Melo SR, do Nascimento GV, de Macedo GF, do Nascimento Marreiro D. Association between cortisol, insulin resistance and zinc in obesity: a mini-review. Biol Trace Element Res. 2019;191:323–30. https://doi.org/10.1007/s12011-018-1629-y.

    Article  CAS  Google Scholar 

  41. Petersen KF, Dufour S, Savage DB, Bilz S, Solomon G, Yonemitsu S, Shulman GI, et al. The role of skeletal muscle insulin resistance in the pathogenesis of the metabolic syndrome. Proc Natl Acad Sci. 2007. https://doi.org/10.1073/pnas.0705408104.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Domínguez-Cruz MG, de Lourdes Muñoz M, Totomoch-Serra A, García-Escalante MG, Burgueño J, Valadez-González N, Díaz-Badillo A, et al. Maya gene variants related to the risk of type 2 diabetes in a family-based association study. G E N. 2020. https://doi.org/10.1016/j.gene.2019.144259.

    Article  Google Scholar 

  43. Xu X, Su S, Barnes VA, De Miguel C, Pollock J, Ownby D, Wang X, et al. A genome-wide methylation study on obesity: differential variability and differential methylation. Epigenetics. 2013. https://doi.org/10.4161/epi.24506.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Meirhaeghe A, Helbecque N, Cottel D, Amouyel P. Impact of polymorphisms of the human β 2-adrenoceptor gene on obesity in a French population. Int J Obes. 2000. https://doi.org/10.1038/sj.ijo.0801168.

    Article  Google Scholar 

  45. Hamada T, Kotani K, Fujiwara S, Sano Y, Domichi M, Tsuzaki K, Sakane N. The common− 55 C/T polymorphism in the promoter region of the uncoupling protein 3 gene reduces prevalence of obesity and elevates serum high-density lipoprotein cholesterol levels in the general Japanese population. Metabolism-Clinical and Experimental. 2008. https://doi.org/10.1016/j.metabol.2007.10.019.

    Article  PubMed  Google Scholar 

  46. Su M, Chen X, Chen Y, Wang C, Li S, Ying X, Fu C, et al. UCP2 and UCP3 variants and gene-environment interaction associated with prediabetes and T2DM in a rural population: a case control study in China. BMC Med Genet. 2018. https://doi.org/10.1186/s12881-018-0554-4.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Muñoz AM, Velásquez CM, Agudelo GM, Uscátegui RM, Estrada A, Patiño FA, Bedoya G, et al. Examining for an association between candidate gene polymorphisms in the metabolic syndrome components on excess weight and adiposity measures in youth: a cross-sectional study. Genes Nutr. 2017. https://doi.org/10.1186/s12263-017-0567-1.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Schnor NPP, Verlengia R, Novais PFS, Crisp AH, Leite CVDS, Rasera-Junior I, Oliveira MRMD. Association of 5-HT2C (rs3813929) and UCP3 (rs1800849) gene polymorphisms with type 2 diabetes in obese women candidates for bariatric surgery. Archives of endocrinology and metabolismo. 2017. https://doi.org/10.1590/2359-3997000000260.

    Article  Google Scholar 

  49. van Abeelen AF, de Krom M, Hendriks J, Grobbee DE, Adan RA, van der Schouw YT. Variations in the uncoupling protein-3 gene are associated with specific obesity phenotypes. Eur J Endocrinol. 2008. https://doi.org/10.1530/EJE-07-0834.

    Article  PubMed  Google Scholar 

  50. Su, M., Chen, X., Chen, Y., Wang, C., Li, S., Ying, X., Fu, C. et al. UCP2 and UCP3 variants and gene-environment interaction associated with prediabetes and T2DM in a rural population: a case control study in China. BMC Medical Gen. https://doi.org/10.1186/s12881-018-0554-4.

  51. Salopuro T, Pulkkinen L, Lindström J, Kolehmainen M, Tolppanen AM, Eriksson JG, Uusitupa M, et al. Variation in the UCP2 and UCP3 genes associates with abdominal obesity and serum lipids: the Finnish Diabetes Prevention Study. BMC Med Genet. 2009. https://doi.org/10.1186/s12881-018-0554-4.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Sun M, Feng W, Wang F, Li P, Li Z, Li M, ... Tse LA, et al. Meta‐analysis on shift work and risks of specific obesity types. Obes Rev. 2018. https://doi.org/10.1111/obr.12621.

Download references

Acknowledgements

The authors wish to thank all participants and the parents/guardians for participating in the study. We also express our gratitude to Mr. Sean M. Quail, a specialist on teaching English as foreign language, who proofread the draft. We acknowledged JSF’s scholarship from the Coordination for the Improvement of Higher Education Personnel (CAPES). ADC and DMS are CNPq fellows. The project was supported by Fundação de Amparo à Pesquisa do Estado de Goiás (FAPEG) under the call PPSUS No. 04/2017.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

JSF, RMP, ADC, RSC, and LBM conceived and designed the study. JSF, RFS, RMP, FRG, and DMS collected and processed the biological samples and performed the Real Time PCR experiments. HPLF, and ADC conducted the statistical analysis. RMP was the patient’s attending physician. All authors drafted, read, and approved the manuscript. ADC and LBM reviewed and edited the manuscript.

Corresponding author

Correspondence to Lysa B. Minasi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee on Human Research from the Pontifical Catholic University of Goiás under the number CAAE: 16303313.4.0000.0037. Written informed consent was signed by parents of all children. Involvement in the study was voluntary and there were no repercussions for non-participation. This study was performed in accordance with the declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Fortes, J.S., Pinto, R.M., de Souza, R.F. et al. The influence of six polymorphisms of uncoupling protein 3 (UCP3) gene and childhood obesity: a case–control study. BMC Pediatr 23, 87 (2023). https://doi.org/10.1186/s12887-023-03905-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12887-023-03905-6

Keywords