Abstract

Background. The PPARG2 Pro12Ala (rs1801282) and IL6 -174G >C (rs1800795) have important function in body weight regulation and a potential role in obesity risk. We aimed to investigate the association between PPARG2 Pro12Ala and IL6 -174G >C variants and the genotypes interaction with body composition, metabolic markers, food consumption, and physical activity in severely obese patients. Methods. 150 severely obese patients (body mass index (BMI) ≥ 35 kg/m2) from Central Brazil were recruited. Body composition, metabolic parameters, physical activity, and dietary intake were measured. The genotype was determined by the qPCR TaqMan Assays System. Multiple linear regression and multiple logistic regression models were fitted adjusting for confounders. Results. Ala carriers of the Pro12Ala polymorphism had higher adiposity measures (BMI: , and fat mass: ) and systolic blood pressure () compared to Pro homozygotes. We found no important associations between the -174G >C polymorphism and obesity phenotypes. When genotypes were combined, individuals with genotypes ProAla + AlaAla and GC + CC presented higher BMI () and higher polyunsaturated fatty acids (PUFAs) consumption () compared to the ones with genotypes ProPro and GG, and individuals carriers of the PPARG2 Ala allele only (genotype ProAla + AlaAla and GG) had higher fat mass and systolic and diastolic blood pressure compared to the ones with genotypes ProPro and GG. Conclusions. Severely obese individuals carrying the Ala allele of the PPARG2 Pro12Ala polymorphism had higher measures of adiposity and blood pressure, while no important associations were found for the IL6 -174G >C polymorphism.

1. Introduction

Studies about the genetic obesity susceptibility have investigated polymorphisms related to genes encoding factors involved in food and energy intake regulation, energy expenditure, and adipogenesis control [1]. In this context, the polymorphisms PPARG2 Pro12Ala (rs1801282) and IL6 -174G >C (rs1800795) have received attention for their possible influence in body weight regulation [13].

The PPARG2 has an important role in modulating the expression of genes involved in the regulation of adipose tissue differentiation and lipid metabolism [4, 5]. The mostly studied polymorphism of PPARG2, the Pro12Ala, may promote lower affinity of the PPAR-γ for the response element and about 50% lower transcriptional capacity [6]. The Ala allele has been associated with lower body mass index (BMI) [6, 7]; nevertheless, some meta-analyses have found contradictory results showing higher BMI in Ala carriers, especially in severely obese individuals [811]. Considering this potential association of the Ala allele and increased BMI, it is important to investigate whether this variant influences other parameters in severely obese individuals.

Regarding interleukin-6 (IL-6), one of the major proinflammatory cytokines, it is positively related to increased BMI [3, 12]. The association studies between obesity and the IL6 -174G >C polymorphism have shown higher BMI in the presence of C allele in cross-sectional and cohort studies [1317], but it was not confirmed by meta-analysis [18, 19]. Thus, the role of Pro12Ala and -174G >C in modulating BMI is still to be confirmed and also the information regarding other body composition parameters such as fat mass, fat-free mass, percentage body fat, and lean mass.

Interaction of genes with environmental factors, such as diet and physical activity, may be involved in the discrepancies of associations between studies [1, 20, 21]. Gene-diet interaction studies with Pro12Ala and -174G >C polymorphisms have found that the energy content and composition of the diet may affect obesity phenotypes, showing the importance to assess the interactions among genotypes and potentially modifiable lifestyle factors [20, 22, 23]. Considering that the PPARG2 Pro12Ala and IL6 -174G >C polymorphisms have important functions in body weight regulation with a potential role in obesity risk and the alarming increase in severe obesity worldwide [24], they are promising single nucleotide polymorphisms (SNPs) for association studies of obesity phenotypes. Thus, we aimed to investigate the association between PPARG2 Pro12Ala and IL6 -174G >C variants and the genotypes interaction with body composition, metabolic markers, food consumption, and physical activity in severely obese patients.

2. Methods

2.1. Subjects

This study is an analysis of baseline data from participants of the clinical trial “Effect of Nutritional Intervention and Olive Oil in Severe Obesity: Randomized Controlled Trial” (DieTBra Trial) (registered at ClinicalTrials.gov: NCT02463435). A total of 150 severely obese patients (BMI ≥ 35 kg/m2) aged 18 to 65 years were recruited from primary care of the Brazilian Unified Health System at Goiânia, Goiás State, in Central Brazil. Patients were referred to the Nutrition in Severe Obesity Outpatient Clinic, and the study took place at the Clinical Research Unit of the Clinical Hospital/Federal University of Goiás. The study excluded individuals that had already underwent bariatric surgery, under actual nutritional treatment for weight loss or in the previous 2 years, using antiobesity or anti-inflammatory drugs, having HIV/AIDS, as well as heart/kidney/hepatic insufficiency, chronic obstructive pulmonary disease, cancer, and pregnancy. Patients were recruited from June to November 2015. All patients who agreed to participate in the study gave written consent. The study was approved by the Research Ethics Committee of Clinical Hospital of Federal University of Goiás (protocol number 747.792).

2.2. Anthropometric and Body Composition Measurements

Body weight and height were measured using standardized procedures [25]. BMI (kg/m2) was calculated dividing the body mass (kg) by the squared height (m2). Severe obesity was defined as BMI ≥ 35 kg/m2 [24].

Fat mass (kg), fat-free mass (kg), percentage body fat, and lean mass were measured using multifrequency bioelectrical impedance analysis (BIA). The measurement was performed with the InBody S10 device (Biospace Co., Ltd., Seoul, Korea) by using different frequencies (1, 5, 50, 250, 500, and 1000 kHz) at each segment (right arm, left arm, trunk, right leg, and left leg).

For BIA assessment, patients were instructed to fast for 12 h and avoid strenuous physical activity and alcohol, as well as food and drinks containing caffeine on the previous day [26]. BIA assessment was conducted according to the manufacturer’s guidelines.

2.3. Dietary Intake

Food consumption was assessed using three 24 h records collected within seven days, being two face to face and one by phone. Data were assessed by trained registered dieticians. We used the multiple pass method (MPM) to collect the 24 h records [27], and the nutritional analysis was performed using Avanutri Online® (Avanutri Equipamentos de Avaliação Ltda, Rio de Janeiro, BR). Energy (kcal), proteins (%), carbohydrates (%), lipids (%), saturated fatty acids (SAFs) (%), monounsaturated fatty acids (MUFAs) (%), polyunsaturated fatty acids (PUFAs) (%), polyunsaturated : saturated fatty acids ratio (P : S ratio), cholesterol (g), and fiber (g) were obtained calculating the mean of the three 24 h records.

2.4. Physical Activity Assessment

Physical activity was assessed using a triaxial accelerometer ActiGraph wGT3X (ActiGraph, Pensacola, FL, USA) for movement registration. Patients were instructed to wear the accelerometer 24 h a day for six consecutive days over the nondominant wrist, even during shower and water activities, as the device was waterproof. The sampling frequency of the accelerometer was set at 30 Hz, and the data collection interval was set at one min. Accelerometers were set up and downloaded at ActiLife 6 software. Output data were processed using the R-package GGIR (http://cran.r-project.org). The outcome measures used in the present study were moderate-to-vigorous physical activity (MVPA) (>100 mg) defined as estimated time spent in ≥ 10 min per bout during a week and the sedentary time (<50 mg, without bouts) in min per day.

2.5. Blood Pressure and Comorbidities

Systolic and diastolic blood pressures were measured with the patient in the sitting position after resting for at least 5 min. Two measures were taken within the 3 min interval using the Omron HEM 742INT (Omron Healthcare Inc., Kyoto, Japan) automatic blood pressure monitor with an appropriately sized cuff, and the mean was calculated.

The presence of comorbidities was analyzed as a dichotomous variable (presence/absence). Subjects with systolic/diastolic blood pressure higher than 140/90 mmHg or under antihypertensive therapy were considered hypertensive [28]. Subjects with fasting glucose ≥ 126 mg/dL or under glucose-lowering therapy were considered diabetic [29]. Individuals with low-density lipoprotein (LDL) cholesterol ≥ 160 mg/dL and/or triglycerides ≥ 150 mg/dL or high-density lipoprotein (HDL) cholesterol < 40 mg/dL for men and < 50 mg/dL for women were classified as dyslipidaemic [30].

2.6. Laboratory Tests

Blood samples were collected for metabolic markers and genomic DNA extraction after 12 h overnight fasting. Serum glucose, total cholesterol, HDL cholesterol, LDL cholesterol, and triglycerides were measured by enzyme-colorimetric methods. Serum insulin was measured by chemiluminescence, and hemoglobin A1c (HbA1c) was measured by liquid chromatography. The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated following the formula derived by Matthews et al. [31].

2.7. DNA Extraction and Genotyping

Genomic DNA was extracted from whole blood with the PureLink™ Genomic DNA Mini Kit (Invitrogen, Carlsbad, CA, USA). DNA concentration and purity were evaluated by spectrophotometric determination of the A260/280 ratio with NanoDrop® 2000c (Thermo Fisher Scientific, Waltham, MA, USA), and DNA quality was checked using agarose gel electrophoresis. Genotyping was performed using custom TaqMan SNP genotyping assays—ID c__1129864_10 for PPARG2 Pro12Ala (rs1801282) and ID c__1839697_20 for IL6 -174G >C (rs1800795) (Applied Biosystems, Foster City, CA, USA)—on a StepOnePlus™ Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, USA). The standard real-time polymerase chain reaction (qPCR) was carried out using the TaqMan GTXpress™ Master Mix (Thermo Fisher Scientific, Waltham, MA, USA) reagent kit in a 21 µL volume according to the manufacturer’s instructions. Although DNA samples were extracted for all study participants, the qPCR amplification was only conducted for the PPARG2 Pro12Ala polymorphism on samples from 146 individuals and for the IL6 -174G >C polymorphism on samples from 148 individuals.

2.8. Statistical Analysis

The dataset was structured using EpiData 3.1, and double entry typing with validation was performed. The data were presented as mean ± SD for continuous variables and frequencies for categorical variables. The chi-squared test was used to analyze the agreement of genotype frequencies with the Hardy–Weinberg equilibrium expectation. Allele frequency was determined by manual counting. Normal distribution was tested for all measured variables using the Kolmogorov–Smirnov test and histograms; skewed variables were normalized by log transformation and then backtransformed for results presentation.

Individual genotype analysis and the combination of the two genotypes analysis (PPARG2 + IL6) were performed. For the combination of genotypes, patients were grouped as follows: no variants (genotypes ProPro and GG, ), IL6 only (genotypes ProPro and GC + CC, ), PPARG2 only (genotypes ProAla + AlaAla and GG, ), and both variants (genotypes ProAla + AlaAla and GC + CC, ). Comparisons were performed using Student’s t-test or ANOVA and chi-squared test or Fisher’s exact test. Due to the low frequency of the variant allele, we compared carriers versus noncarriers of the minor allele. We fitted multiple linear regression models to adjust the analysis for potential confounders (age, sex, BMI, sedentary time, and diabetes). For binary variables, odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated and multiple logistic regression models were fitted adjusting for the same confounders. Statistical analyses were performed in Stata 12.

3. Results

All the study participants () had DNA extracted; however, genetic samples were viable in 146 individuals for the PPARG2 Pro12Ala polymorphism and in 148 individuals for the IL6 -174G >C polymorphism. The genotype distribution for PPARG2 Pro12Ala was 86.9%, 12.4%, and 0.7% for ProPro, ProAla, and AlaAla, respectively. The minor allele frequency of the Ala allele was 0.065. For IL6 -174G >C, the frequencies of GG, GC, and CC genotypes were 65.3%, 31.3%, and 3.4%, respectively. The minor allele frequency of the C allele was 0.193. Observed genotype frequencies were in agreement with the Hardy–Weinberg equilibrium ( for PPARG2 Pro12Ala and for IL6 -174G >C).

The characteristics of the study participants according to the PPARG2 Pro12Ala polymorphism are displayed in Table 1. Analysis of the Pro12Ala polymorphism showed higher BMI () and fat mass () for Ala carriers, even after adjustment for age, sex, sedentary time, and diabetes. Ala carriers presented significantly higher SBP and DBP, but after adjustments, only SBP () remained associated (Table 1).

The characteristics of the study participants according to the IL6 -174G >C polymorphism are displayed in Table 2. Sex and MVPA were associated with the -174G >C polymorphism after adjustments ( and , respectively). Males had triple probability (OR: 3.60; 95% CI: 1.04–12.48) to be C carriers, and the C carriers spent lower amount of time in MVPA (Table 2).

For combined genotypes, participants were grouped as follows: no variants (), IL6 variant only (), PPARG2 variant only (), and both variants (). The same variables of the individual genotype analysis were tested for the combined genotypes, but the results were presented in figures only for the variables with a significant statistical difference. Analysis of the combined effects of the two genotypes showed association with BMI, fat mass, SBP, DBP, and polyunsaturated fat consumption after adjustments. Individuals with both variants had higher BMI () compared to the ones with no variants. Fat mass, SBP, and DBP were higher for participants with the PPARG2 variant only compared to those with no variants (, , and , respectively). Individuals with both variants presented higher consumption of PUFA compared to the ones with no variants () (Figure 1). The analysis of BMI in categories did not show association with the polymorphisms or the genotype combination (Table 3).

4. Discussion

The Ala allele of PPARG2 Pro12Ala (rs1801282) is potentially associated with higher BMI; however, few studies have investigated other key factors that may interact with the variant leading to severe obesity susceptibility such as metabolic, food consumption, and clinical parameters. To our knowledge, this is the first study to find association of the PPARG2 Pro12Ala (rs1801282) polymorphism with an obesity phenotype other than BMI and blood pressure in severely obese patients. We found higher measures of adiposity (BMI and fat mass) and higher SBP in Ala carriers of the Pro12Ala polymorphism compared to the ProPro individuals. These results may contribute to a better understanding of the pathophysiology of severe obesity and translate into more effective preventive and treatment measures to halt its increasing prevalence.

For the PPARG2 Pro12Ala polymorphism, we observed a genotype distribution (87% for ProPro and 13% for Ala carriers) similar to that in the few studies with Brazilians [32, 33]. Among studies with morbidly obese individuals (BMI ≥ 40 kg/m2), frequencies range from 72.5% to 86.0% for the ProPro genotype and from 14.0% to 27.5% for Ala carriers [3438]. For the -174G >C polymorphism, frequencies vary widely. In the current study, we observed 64.9% of GG and 35.1% of GC + CC, while French morbidly obese individuals presented frequencies of 43.1% for GG and 56.9% for GC + CC [39].

Despite controversial results regarding the PPARG2 Pro12Ala polymorphism, meta-analyses have shown association of the Ala allele with higher BMI and fat mass, especially in severely obese individuals, corroborating our results [68, 10, 11, 40]. The controversies observed in other studies may be explained by a suggested interaction between dietary fat intake and the Pro12Ala polymorphism, showing that when the dietary P : S ratio is similar or lower than 0.66, Ala carriers present higher BMI than Pro homozygotes, while when the P : S ratio is higher than 0.66, the opposite occurs [41]. Our patients had the mean P : S ratio ≤ 0.66, possibly explaining the higher BMI among Ala carriers and demonstrating the effect of dietary fat intake on the phenotype determination.

PPAR-γ2 has a crucial role in modulating lipid metabolism and adipose tissue accumulation [4, 5]. The PPARG2 Pro12Ala polymorphism has been suggested to have a role in the variance of fat mass among analyses of candidate genes for adiposity changes [42]. Few studies have investigated the influence of the Pro12Ala polymorphism on body composition parameters other than body mass, BMI, and waist circumference. Association between the Ala allele and higher fat mass was found in Caucasians from the Québec Family Study [43] and in Italian women [44], corroborating the results of the current study. Thus, more detailed measurements of body composition are important to be addressed in studies of the Pro12Ala polymorphism.

We found an association between the Ala allele of the Pro12Ala polymorphism and higher blood pressure in severely obese individuals. Although the association between the Pro12Ala variants and blood pressure is controversial, a meta-analysis of eight studies with 3281 individuals (1865 cases and 1416 controls) suggested that the Ala allele might be protective for hypertension among East Asians, but not among Caucasians [45]. The mechanisms behind this association are unclear. Studies have suggested that the activation of the PPAR-γ2 inhibits processes related to inflammation and hypertension. The Ala allele has been associated with impaired function of the PPAR-γ2, leading to increased blood pressure [4648].

We found higher probability of severely obese males to be C carriers of the IL6 -174G >C polymorphism after adjustments. Studies comparing elderly (octa/nonagenarians and centenarians) with younger individuals have found a tendency in the reduction of GG genotypes in males [4952]. Regarding MVPA, we did not find any evidence that justifies this association. More investigation is needed to clarify the association of the higher frequency of C alleles in severely obese males and the -174G >C polymorphism and also the association of the C allele with lower time spent in MVPA.

When the Pro12Ala and the -174G >C genotypes were combined, association was found between the PPARG2 variant only and fat mass, SBP, and DBP compared to the ones with no variants, as demonstrated in the independent analysis of genotypes. We also observed higher BMI and higher PUFA consumption in the presence of both variants compared to the ones with no variants. Nevertheless, this association of the combined genotypes seems to be determined by the PPARG2 Pro12Ala polymorphism as the IL6 -174G >C polymorphism was not associated with adiposity indexes in the independent analysis of genotypes and also due to the low number of individuals with both variants.

Our study has limitations such as the small sample size, especially in the analysis of combined genotypes, and the impossibility to demonstrate causality due to the study design. However, we may address some strengths, such as the analysis of the association between SNPs and also the combination of genotypes, with a variety of parameters, such as anthropometrical, clinical, biochemical, physical activity, and dietary intake, in severely obese patients; the adjusted analysis for potential confounders, such as sociodemographic, anthropometrical, and physical activity measurements; and the high-quality methods used to assess physical activity (triaxial accelerometer) and body composition (BIA).

5. Conclusions

In summary, we found that severely obese individuals carrying the Ala allele of the PPARG2 Pro12Ala polymorphism have higher measures of adiposity (BMI and fat mass) and blood pressure, while no important associations were found for the IL6 -174G >C polymorphism. Further studies on gene-diet/gene-environment interactions are necessary to clarify the associations and underlying mechanisms between SNPs and severe obesity. This information may eventually be used to develop tailored interventions specific to the individual’s genotype and enable more effective prevention and treatment measures to address this increasing public health problem.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

The funder had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

We thank the patients who participated in the study. We also would like to thank the Clinical Research Unit of Clinical Hospital from Universidade Federal de Goiás and its technical staff for lending us the facilities to conduct this research and giving technical support. This study was partially funded by the Goiás State Research Support Foundation (FAPEG) (grant number 201310267000003).