Obese children and adolescents have an impaired metabolic and anthropometric profile
The anthropometric, body composition, and biochemical variables of the studied individuals stratified by BMI percentile status are shown in Table 1. Age and gender of the participants were not significantly different among the control and obese group, the median age was 11 years old, and 48% were girls. All mean values of skinfold thickness, waist circumference (WC), arm circumference (AC), waist-to-height ratio (WHR), arm muscle circumference, and neck circumference were significantly increased in participants with obesity in comparison with healthy weight children and adolescents (Table 1). Regarding body composition, the percentage of arm, leg, trunk, android, gynoid, and body fat were significantly elevated in obese patients compared to control subjects (Table 1). In obese children, fat distribution was similar between different body areas, and there were no significant differences in mean values of lean mass between groups.
Table 1
Anthropometric, body composition and clinical cardiometabolic risk factors in children and adolescents†
| Control | Obesity | p value |
| n = 19 | n = 18 | |
Gender | | | |
Female | 8 | 10 | |
Male | 11 | 8 | |
Age (years) | 10,47 ± 2,87 | 11,33 ± 2,93 | 0,3515 |
Tanner | 2,88 ± 1,72 | 2,72 ± 1,32 | 0,8826 |
BMI | 49,92 ± 28,03 | 96,37 ± 8,50 | < 0,0001**** |
Arm circumference (cm) | 20,41 ± 3,23 | 27,45 ± 3,50 | < 0,0001**** |
Arm muscle circumference (cm) | 16,77 ± 2,83 | 20,64 ± 3,19 | 0,0008*** |
Neck circumference (cm) | 28,29 ± 3,03 | 32,47 ± 2,82 | 0,0003*** |
Triceps skinfold (cm) | 11,68 ± 4,61 | 21,64 ± 4,94 | < 0,0001**** |
Subscapular skinfold (cm) | 6,85 ± 2,27 | 21,17 ± 5,40 | < 0,0001**** |
Suprailiac skinfold (cm) | 10,91 ± 4,78 | 26,61 ± 5,79 | < 0,0001**** |
Abdominal skinfold (cm) | 10,53 ± 5,01 | 27,64 ± 6,17 | < 0,0001**** |
Waist circumference (WC) (cm) | 63,75 ± 7,53 | 83,75 ± 6,94 | < 0,0001**** |
Waist/Height ratio (WHR) | 0,45 ± 0,03 | 0,58 ± 0,04 | < 0,0001**** |
Arm fat (%) | 33,07 ± 9,44 | 45,91 ± 7,73 | < 0,0001**** |
Leg fat (%) | 32,82 ± 7,74 | 44,41 ± 5,82 | < 0,0001**** |
Trunk fat (%) | 21,06 ± 8,95 | 42,9 ± 6,79 | < 0,0001**** |
Android fat (%) | 18,58 ± 9,93 | 44,33 ± 7,53 | < 0,0001**** |
Gynoid fat (%) | 28,88 ± 8,21 | 44,68 ± 5,96 | < 0,0001**** |
Body fat (%) | 26,28 ± 7,35 | 39,9 ± 5,89 | < 0,0001**** |
Lean mass (kg) | 32,83 ± 10,31 | 22,84 ± 11,25 | 0,0073** |
Fasting blood glucose (mg/dL) | 87,61 ± 5,55 | 87,59 ± 6,25 | 0,936 |
2 hours glucose (mg/dL) | 91,74 ± 13 | 94,23 ± 13,58 | 0,7617 |
Insulin (mcU/mL) | 7,77 ± 4,78 | 19,77 ± 8,02 | < 0,0001**** |
Total cholesterol (mg/dL) | 150,15 ± 21,83 | 158,29 ± 24,91 | 0,3758 |
Triglycerides (mg/dL) | 76,48 ± 18,47 | 117,27 ± 31,64 | 0,0001*** |
HDL (mg/dL) | 50,34 ± 8,22 | 41,39 ± 9,51 | 0,0077** |
LDL (mg/dL) | 84,55 ± 18,04 | 94,70 ± 20,67 | 0,2547 |
VLDL (mg/dL) | 15,29 ± 3,69 | 23,45 ± 6.33 | 0,0001*** |
HOMA-IR (score) | 1,50 ± 1,19 | 4,28 ± 1,72 | < 0,0001**** |
TG/HDL ratio | 1,60 ± 0,61 | 3,02 ± 1,24 | < 0,0001**** |
Android-gynoid % fat ratio | 0,61 ± 0,17 | 0,99 ± 0,08 | < 0,0001**** |
GLP-1 (ng/mL) | 15,98 ± 0,69 | 17,70 ± 0,66 | 0,0725 |
GIP (pg/mL) | 2637 ± 486,89 | 3137 ± 512,44 | 0,6299 |
†BMI: body mass index. HDL: high-density lipoprotein. VLDL: very low-density lipoprotein. LDL: low-density lipoprotein. HOMA-IR: homeostatic model of assessment for insulin resistance. GLP-1: glucagon-like peptide 1. GIP: gastric inhibitory polypeptide. Values are presented as percentage or mean ± standard deviation (SD). Differences between groups were evaluated using Wilcoxon signed-rank tests. * p < 0.05, ** p < 0.005, *** p < 0.0001. |
The mean of fasting glucose and 2-hours glucose fell within the normal ranges, as recommended by the American Academy of Pediatrics (Table 1) [20]. However, insulin was significantly higher in children with obesity than in the standard weight group. Accordingly, the HOMA score was increased in obese children; this score overcomes the 3,4 cut-off indicating IR (Table 1) [21]. TG and VLDL levels were higher in obese children than those of lean subjects, whereas HDL-c levels were significantly reduced in obese subjects in contrast with the control group (Table 1). The levels of total cholesterol were similar among groups (Table 1).
TG/HDL and A/G ratios are frequently used as indicators of cardiometabolic risk, and here, those ratios were significantly elevated in the obese group in comparison with the lean patients. Based on this data, obese children and adolescents have an increased risk of developing cardiovascular disease and T2D. Incretins gastric inhibitory polypeptide (GIP) and GLP-1 did not show significant differences between the study groups (Table 1).
Triceps skinfold displayed a positive correlation with body composition parameters
Relationships between fat body distribution and anthropometric data are shown in Table 2. Triceps skinfold displayed a positive correlation with all body composition parameters, the stronger correlation was observed with total body fat percentage (r = 0,6554, p = 0,0071), and a weaker association with Android fat percentage (r = 0,4974, p = 0,0516) (Table 2). Furthermore, WHT ratio showed a positive association with arm (r = 0,691, p = 0,0038), leg (r = 0,5245, p = 0,0388), android (r = 0,5528, p = 0,0282), and total body fat percentages (r = 0,5379, p = 0,0334) (Table 2).
Table 2
Correlation of DEXA fat measurements with anthropometric data†
| Arm fat (%) | Leg fat (%) | Trunk fat (%) | Android fat (%) | Gynoid fat (%) | Body fat (%) |
BMI | 0,060 | 0,131 | -0,013 | -0,063 | -0,028 | 0,081 |
Arm circumference | -0,118 | 0,014 | -0,167 | -0,217 | -0,179 | -0,046 |
Arm muscle circumference | -0,340 | -0,202 | -0,437 | -0,419 | -0,425 | -0,299 |
Neck circumference | -0,449 | -0,278 | -0,411 | -0,496 | -0,391 | -0,355 |
Triceps skinfold | 0,590a | 0,627a | 0,604a | 0,497a | 0,579a | 0,655a |
Subscapular skinfold | 0,148 | 0,022 | 0,294 | 0,186 | 0,179 | 0,171 |
Suprailiac skinfold | -0,041 | -0,186 | 0,099 | 0,140 | 0,031 | -0,027 |
Abdominal skinfold | 0,189 | -0,029 | 0,117 | 0,289 | -0,004 | 0,119 |
Waist circumference (WC) | -0,105 | -0,110 | -0,255 | -0,217 | -0,234 | -0,149 |
Waist/Height ratio (WHR) | 0,691b | 0,525a | 0,442 | 0,553a | 0,482 | 0,538a |
† BMI: body mass index. Spearman's rank correlation test, a p < 0.05 and b p < 0.005. |
Trunk fat percentage presented a significant association with A/G ratio (r = 0,558, p = 0,027). Likewise, Android fat percentage also displayed a positive correlation with A/G ratio (r = 0,558, p = 0.006). Arm circumference parameter correlated negatively with LDL (r=-0,572, p = 0,013) and no-HDL (r=-0,502, p = 0,034) whereas arm muscle circumference exhibited a weak negative correlation with LDL (r=-0,486, p = 0,041) (Supplementary table 2).
The correlation of different variables with the suprailiac skinfold indicates a positive association with Total Cholesterol (r = 0,658, p = 0,003), and no-HDL (r = 0,724, p = 0,001) while showed a weaker association with LDL (r = 0,550, p = 0,018) (Supplementary Table 2). There was a positive correlation between Abdominal skinfold with HOMA-IR (r = 0,459, p = 0,055).
PPARα and GLP-1R gene expression is reduced in leukocytes from obese subjects
To determine the expression of PPAR isotypes in the leukocytes, the expression of each isotype was analyzed by qPCR. PPAR-α expression showed a significant reduction in the leukocytes from obese patients in comparison with the control group, with a 50% reduction in the expression of this gene (p = 0,0484) (Fig. 1a). PPAR-β did not show differential expression between the study groups (Fig. 1b). In contrast, PPAR-γ expression was undetectable in leukocytes from all the samples. To corroborate this finding, we used adipose tissue cDNA as a positive control, since this isotype is mainly present in adipose cells. This sample presented a positive amplification and a single peak in the dissociation curve, indicating that the reaction was specific (data not shown). Likely, PPAR-γ transcripts were not expressed on leukocytes or are expressed at low levels but just in some cell lineages of these samples.
Besides, analysis from mRNA expression levels of incretin receptors was also evaluated in the studied individuals. GLP-1R expression was significantly reduced in the obese group as compared to the healthy weight participants (p = 0,1358) (Fig. 2a), whereas no difference was observed in GIPR expression between the two groups (Fig. 2b).
Obese children and adolescents showed a proinflammatory profile of adipokines and cytokines
To identify the inflammatory serum profile in obese children and adolescents, the levels of several adipokines and cytokines were measured by flow cytometry using a multiplex assay. Obese subjects showed a significant increase in the levels of IL-8 (p = 0,0081), IL-6 (p = 0,0005), TNF-α (p = 0,0004) IFN-γ (p = 0,0110), and MCP-1 levels (p = 0,0452) in comparison with the normal weight group (Fig. 3a and b). Serum levels of IL-10 and IP-10 were detected but did not differ significantly between the groups.
The serum adiponectin concentration presented a significant reduction in the obese group (p = 0,0452) compared to the control group (Fig. 4). Conversely, adipsin displayed a concentration significantly higher (p = 0,0397) in the obese group in comparison with the control group, which had a concentration of 7,688 ± 6,116 pg/mL. (Fig. 4). Resistin showed similar levels between the study groups, with concentrations of 13,102 ± 8,342 pg/mL and 16,792 ± 8,302 pg/mL for the obese and control participants, respectively (Fig. 4). Leptin was also analyzed; however, it presented considerably high levels in the group with obesity, which exceeded the concentrations determined in the calibration curve, which did not allow to quantify it using the test used. Finally, RBP4 did not present a representative concentration in any of the study groups, indicating that it is not secreted in serum at detectable levels under the evaluated conditions.
PPARα and PPARβ/δ correlate negatively with proinflammatory markers
Analysis of the relationships between PPAR-α, PPAR-β, and GLP-1R expression and the levels of cytokines, chemokines, adipokines and anthropometric parameters of the obese subjects showed that PPAR-α transcript levels had a significant negative correlation with TNF-α levels (r = -0,583, p = 0,03883) (Table 3), as well as with abdominal skinfold (r = -0,712, p = 0,0016). Besides PPAR-β showed a significant negative correlation with IL-8 levels (r = -0,667, p = 0,0085) and arm fat percentage (r = -0,651, p = 0,0132) (Tables 3 and 4).
Table 3
Correlation between hormones and cytokines with PPAR-α, PPAR-β and GLP-1R expression in children and adolescents†
| PPAR-α | PPAR-β | GLP-1R |
Adiponectin | 0,496 | -0,330 | -0,084 |
Adipsin | 0,067 | -0,667 | 0,357 |
MCP-1 | 0,324 | 0,011 | 0,235 |
IP-10 | 0,370 | -0,300 | 0,033 |
IL-10 | -0,125 | 0,000 | -0,136 |
IL-8 | 0,233 | -0,664b | -0,007 |
IL-6 | -0,343 | -0,200 | 0,317 |
IFN-γ | 0,130 | -0,411 | -0,130 |
Resistin | 0,172 | -0,311 | 0,323 |
TNF-α | -0,583a | -0,059 | 0,304 |
† Spearman's rank correlation test. a p = 0,0383, b p = 0,0085 |
Table 4
Correlation between anthropometric and biochemical parameters with gene expression of metabolic markers in children and adolescents†
| PPAR-α | p value | PPAR-β | p value | GLP-1R | p value |
BMI | 0,049 | 0,854 | -0,257 | 0,354 | -0,086 | 0,773 |
Arm fat (%) | 0,052 | 0,850 | -0,651 | 0,013 | -0,364 | 0,246 |
Leg fat (%) | 0,208 | 0,438 | -0,515 | 0,061 | -0,266 | 0,404 |
Trunk fat (%) | 0,096 | 0,723 | -0,436 | 0,118 | -0,322 | 0,309 |
Android fat (%) | -0,028 | 0,915 | -0,339 | 0,231 | -0,420 | 0,177 |
Gynoid fat (%) | 0,108 | 0,690 | -0,513 | 0,062 | -0,189 | 0,546 |
Body fat (%) | 0,096 | 0,723 | -0,511 | 0,063 | -0,343 | 0,276 |
Arm circumference | -0,094 | 0,713 | 0,306 | 0,266 | -0,147 | 0,608 |
Arm muscle circumference | -0,173 | 0,501 | 0,411 | 0,128 | 0,042 | 0,892 |
Waist circumference (WC) | -0,308 | 0,226 | 0,050 | 0,860 | 0,059 | 0,844 |
Neck circumference | -0,048 | 0,852 | 0,190 | 0,496 | 0,090 | 0,758 |
Waist/Height ratio (WHR) | -0,262 | 0,291 | -0,448 | 0,090 | -0,326 | 0,241 |
Triceps skinfold | 0,210 | 0,415 | -0,390 | 0,146 | -0,277 | 0,330 |
Subscapular skinfold | 0,142 | 0,586 | -0,225 | 0,413 | 0,183 | 0,527 |
Suprailiac skinfold | -0,429 | 0,084 | 0,285 | 0,301 | -0,236 | 0,405 |
Abdominal skinfold | -0,712 | 0,002 | 0,379 | 0,163 | -0,678 | 0,009 |
Fasting blood glucose | -0,364 | 0,149 | 0,340 | 0,214 | -0,279 | 0,333 |
2 hours glucose | 0,042 | 0,876 | -0,412 | 0,125 | -0,051 | 0,868 |
Insulin | -0,185 | 0,471 | 0,293 | 0,287 | -0,288 | 0,318 |
Total cholesterol | -0,021 | 0,934 | 0,075 | 0,790 | 0,209 | 0,473 |
Triglycerides | -0,247 | 0,335 | -0,064 | 0,815 | -0,125 | 0,671 |
HDL | 0,416 | 0,098 | -0,054 | 0,845 | 0,354 | 0,215 |
VLDL | -0,248 | 0,333 | -0,064 | 0,815 | -0,117 | 0,685 |
LDL | 0,033 | 0,900 | -0,086 | 0,756 | 0,077 | 0,797 |
HOMA-IR | -0,367 | 0,146 | 0,343 | 0,209 | -0,341 | 0,234 |
TG/HDL ratio | -0,396 | 0,115 | 0,097 | 0,731 | -0,182 | 0,532 |
† BMI: body mass index. HDL: high-density lipoprotein. VLDL: very low-density lipoprotein. LDL: low-density lipoprotein. HOMA-IR: homeostatic model of assessment for insulin resistance. Spearman's rank correlation test. |
In contrast, GLP-1R expression did not correlate with any inflammatory parameters (Table 3). However, GLP-1R showed a negative correlation with Abdominal skinfold (r= -0,678, p = 0,009) (Table 4).