Next Article in Journal
Evolution of NIN and NIN-like Genes in Relation to Nodule Symbiosis
Previous Article in Journal
Evolutionary History of the Risk of SNPs for Diffuse-Type Gastric Cancer in the Japanese Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pro12Ala PPAR-γ2 and +294T/C PPAR-δ Polymorphisms and Association with Metabolic Traits in Teenagers from Northern Mexico

by
Martín A. Carrillo-Venzor
1,
Nancy R. Erives-Anchondo
1,
Janette G. Moreno-González
1,
Verónica Moreno-Brito
1,
Angel Licón-Trillo
1,
Everardo González-Rodríguez
1,
Pilar del Carmen Hernández-Rodríguez
2,
Sandra A. Reza-López
1,
Verónica Loera-Castañeda
3 and
Irene Leal-Berumen
1,*
1
Faculty of Medicine and Biomedical Sciences, Autonomous University of Chihuahua, Circuito Universitario, Campus II, Chihuahua 31109, Mexico
2
Faculty of Chemical Sciences, Autonomous University of Chihuahua, Circuito Universitario, Campus II, Chihuahua 31109, Mexico
3
CIIDIR-IPN, Durango 34220, Mexico
*
Author to whom correspondence should be addressed.
Genes 2020, 11(7), 776; https://doi.org/10.3390/genes11070776
Submission received: 12 May 2020 / Revised: 1 July 2020 / Accepted: 7 July 2020 / Published: 10 July 2020
(This article belongs to the Section Human Genomics and Genetic Diseases)

Abstract

:
Peroxisome proliferator-activated receptors (PPARs) play roles in glucose and lipid metabolism regulation. Pro12Ala PPAR-γ2 and +294T/C PPAR-δ have been associated with dyslipidemia, hyperglycemia and high body mass index (BMI). We compared metabolic traits and determined associations with Pro12Ala PPAR-γ2 or +294T/C PPAR-δ polymorphism among teenagers from different ethnicity. Four hundred and twelve samples with previous biochemical and biometric measurements were used. Genomic DNA from peripheral blood was extracted and analyzed by end-point PCR for Pro12Ala PPAR-γ2. The +294T/C PPAR-δ PCR product was also digested with Bsl I. Two genotype groups were formed: major allele homozygous and minor allele carriers. Pro12Ala PPAR-γ2 G minor allele frequencies were: 10% in Mestizo-1, 19% in Mestizo-2, 23% in Tarahumara, 12% in Mennonite, and 17% in the total studied population. The +294T/C PPAR-δ C minor allele frequencies were: 18% in Mestizo-1, 20% in Mestizo-2, 6% in Tarahumara, 13% in Mennonite, and 12% in the total studied population. Teenagers with PPAR-γ2 G allele showed a greater risk for either high waist/height ratio or low high-density lipoprotein; and, also had lower total cholesterol. Whereas, PPAR-γ2 G allele showed lower overweight/obesity phenotype (BMI Z-score) frequency, PPAR-δ C allele was a risk factor for it. Metabolic traits were associated with both PPAR polymorphisms.

1. Introduction

The worldwide prevalence of obesity in 2015 was 5% in children and 12% in adults [1]. According to the Organization for Economic Co-operation and Development (OECD) Heavy Burden of Obesity 2018, close to 60% of the population from the 36 countries members of the OECD have overweight and nearly 25% have obesity. Obesity and related diseases will reduce life expectancy 0.9–4.2 years in the future and 92 million premature deaths are expected due to obesity-related diseases by 2050 [2]. In the US, the obesity rate has increased from 10.5 to 20.6% in 12–19 year-old teenagers between 1988 and 2014 [3]. Obesity has also increased drastically in Mexico. According to the Encuesta Nacional de Salud y Nutrición 2018, the overweight/obesity prevalence in Mexico was 38.4% for 12–19 year-olds [4].
PPARs play roles in physiology, pathology, and whole-body energy regulation, lipid and glucose metabolism. Pro12Ala PPAR-γ2 and +294T/C PPAR-δ have been associated with dyslipidemia, hyperglycemia, and overweight/obesity. PPARs are peroxisome proliferator-activated receptors and members of the nuclear hormone receptor superfamily, located in a wide variety of tissues and cells. They are ligand-activated nuclear transcription factors that interact with complex metabolic networks [5,6]. The PPAR-γ gene is located on chromosome 3p25 (OMIM number 601487) and it encodes three proteins isoforms [7,8]. The PPAR-γ1 isoform is expressed in most tissues; the γ3 is seen in macrophages, colon epithelium, and adipose tissue [9]. However, the γ2 is specific for adipose tissue and it can be activated by fatty acids including prostanoids, thiazolidinediones, and other insulin-sensitizing antidiabetic agents [10,11]. The Pro12Ala polymorphism is a missense coding variant (CCA to GCA) in γ2, which results in a proline (major allele) to alanine (minor allele) substitution (rs1801282) [12]. The presence of this single nucleotide polymorphism (SNP) has been associated with type 2 diabetes (T2D) related traits, hypertension, abnormal lipid profile, insulin resistance, and greater body mass index (BMI) [9,13]. Several meta-analysis studies have reported between 13–16% increased risk of T2D in adult population from Europe, North America, and East Asia with a Pro allele [8,14]. However, different results have been reported regarding BMI and Pro12Ala polymorphism in other ethnic groups [15,16,17,18,19]. The prevalence of the modified Ala12 allele varies from about 4% in Asian populations to 28% in Caucasians [20].
The PPAR-δ gene is located in chromosome 6p21.3 (OMIM number 600409) with 13 exons [21]. It has ubiquitous expression; however, its greatest expression occurs in tissues with high lipid metabolism, including the small intestine, heart, adipose tissue, and skeletal muscle [22]. PPAR-δ can be activated by polyunsaturated fatty acids and synthetic compound GW501516 [23]. The +294T/C polymorphism (rs2016520) is also called −87 T/C or +15C/T and is a T (major allele)–C (minor allele) base exchange in the 5′-UTR region of the gene [24,25]. The presence of this SNP has been related to several metabolic traits. It is mostly associated with alterations in the blood levels of HDL, LDL, VLDL, and triglycerides [24,25,26,27,28]. There is also an association between the C allele and BMI, glucose levels, height, and coronary heart disease (CHD) [24,26,27,29,30,31,32,33,34]. The prevalence of the C allele varies among different populations: 22.5% to 31.04% in Chinese people, 19.2% in Germans, 20.1% in French Canadians and 26% in Mexicans [22,25,28,35,36].
A few studies in Mexican teenagers have shown association with Pro12Ala PPAR-γ2 polymorphism with obesity-related traits. Stryjecki et al. reported an association between this polymorphism with insulin resistance in Mexican children [11]. No studies about +294T/C PPAR δ polymorphism and metabolic traits in Mexican children were found. The aim of this study is to determine the association between Pro12Ala PPAR-γ2 or +294T/C PPAR-δ polymorphisms with metabolic traits in Mexican teenagers from different ethnic populations in Chihuahua, Mexico. We focused on teenagers to investigate whether polymorphisms and clinical metabolic traits relate, since the early detection of inherited polymorphisms, which cannot be modified, will allow us to promote preventive behaviors as eating healthy food, eating schedules, and physical activity.

2. Materials and Methods

2.1. Study Population

The Facultad de Medicina y Ciencias Biomédicas, Ethical Research Committee Board (Protocol registration numbers FM-FM-A269/12 and CI-040-17) from the Universidad Autónoma de Chihuahua approved the study protocol; written authorization was obtained from one of the parents or tutors. There were 412 blood samples collected from unrelated adolescents aged 12–18 years. Samples were collected from a previous scholar health care study in three different populations (Mestizo, Tarahumara, and Mennonite) from Chihuahua, Mexico. However, two Mestizo groups were considered (Mestizo-1 and Mestizo-2), since they were from different geographic locations. Moreover, we found significant differences when clinical metabolic traits were compared between them. On the other hand, a parallel study in our lab (same mestizos) showed significant differences in mitochondrial haplotype frequencies between these two groups (results not shown). Recruitment was done in collaboration with public high schools not randomly selected. The data used to support the findings of this study are restricted by the Ethics Board “Comité de Ética en Investigación de la Facultad de Medicina y Ciencias Biomédicas de la Universidad Autónoma de Chihuahua”, in order to protect participants’ confidentiality. Data are available from Irene Leal-Berumen, [email protected], for the researchers who meet the criteria for access to confidential data.

2.2. Phenotyping

All of the participants were weighed while using a digital scale (Tanita BC-418, Ilinois, USA), and height was measured with a portable stadiometer (HM200P, Taichung, Taiwan). The waist circumference (WC) was measured at the midpoint between the lowest rib and the iliac crest after a normal exhalation with students in the standing position, and percentile waist circumference was classified as <p90 normal or ≥p90 increased [37,38]. The waist-to-height ratio (WHR) is a novel abdominal overweight indicator < 0.45 is normal, ≥0.45 is overweight, and ≥0.5 is obese [39]. Body mass index was calculated as weight (kg)/height(m2) and individual values were then compared to BMI-for age and sex-specific reference charts in order to obtain their percentile and the Z-score value [40,41,42]. The percentile BMI (pBMI) was classified as p5 to <p85 as normal, p85 to <p95 overweight, and ≥p95 obese according to CDC growth reference charts, while for Z-score > 1.04 [40]. Blood samples were obtained after 8–12 h fast to measure the following via an automated clinical analyzer (Prestige 24i; Tokio, Boeki Medical System LTD, Japan): fasting glucose, total cholesterol (TC), high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), and triglycerides (TG). The metabolic traits cut-off points included high blood pressure, defined as systolic and/or diastolic blood pressure ≥ p90 (SBP/DBP); waist circumference ≥ p90; TG ≥ 110 mg/dL; HDL ≤40 mg/dL; and, high blood glucose ≥ 110 mg/dL according to Cook et al. [37]. Other metabolic traits risk values were included TC ≥ 150 mg/dL, LDL ≥ 110 mg/dL, and the atherogenic index (AI = TC/HDL; male > 4.0 and female > 3.5) [43,44].

2.3. Genotyping

Genomic DNA was isolated from a peripheral blood buffy coat while using a Master Pure Epicentre kit (Thermo Scientific, Madison, WI, USA). The Pro12Ala PPAR-γ2 polymorphism was genotyped by mutagenically separated PCR with two different length allele-specific downstream primers (P1 and P2) and a common upstream primer (P3). The P1 (5′-GTGTATCAGTGAAGGAATCGCTTTCTTG-3′) was specific for the C allele (Pro); and P2 (5′-TTGTGATATGTTTGCAGACAAGGTATCAGTGAAGGAATCGCTTTGTGC-3′) bound to the G allele (Ala). The P3 (upstream primer) was 5′-TTTCTGTGTTTATTCCCATCTCTCCC-3′. The bases underlined and in bold type indicate the location of mismatches to maintain the specificity of the two separate amplification reactions. Here, DNA (30 ng) was added to a 25 mL reaction mix containing 3 mM MgCl2, 5 pmol P1, 5 pmol P2, 5 pmol P3, 200mM dNTP´s, and 1 U Taq polymerase (Invitrogen, Carlsbad, CA, USA). The PCR conditions used an initial denaturation of 3 min. at 94 °C followed by 35 cycles of denaturation at 94 °C for 45 s, annealing at 62 °C for 45 s, and extension at 72 °C for 45 s. The final extension step was 5 min. at 72 °C (Agilent SureCycler 8800, Santa Clara, CA, USA). A 230 bp product identified the Pro allele, and the Ala-specific product was 250 bp. Electrophoresis used 3.5% agarose (Figure S1) [9]. The +294T/C PPAR-δ polymorphism genotyping used endpoint PCR with 5′-CATGGTATAGCACTGCAGGAA-3′ (forward) and 5′-CTTCCTCCTGTGGCTGCTC-3′ (reverse) primers. The 25-mL mix reaction contained 1.5 mM MgCl2, 0.2 μM forward primer, 0.2 μM reverse primer, 200 μM dNTPs, and 1 U Taq polymerase (Invitrogen, Carlsbad, CA, USA). The PCR conditions were the same as in the Pro12Ala polymorphism. The 269 bp resulting PCR product was digested with 5 U fast Bsl I (Thermo Scientific, Madison, WI, USA) restriction enzyme for 1 h at 37 °C. The 12% polyacrylamide electrophoresis was used to identify the three different genotypes: single 269 bp fragment for TT, three fragments for TC (269, 167, and 102 bp), and two fragments for CC (167 and 102 bp) (Figure S1) [25].

2.4. Statistical Analysis

An exploratory analysis was performed to verify the data quality and observe variable distribution and frequencies. Normality tests were conducted for variables measured in ratio scale. The mean and standard deviation or median and interquartile range were used to describe normally (SBP, glucose, total cholesterol, and LDL) or non-normally distributed variables (the other metabolic traits), respectively. Metabolic traits measured in ratio scale were analyzed as continuous or categorical variables with the cutoffs defined in Section 2.2. For comparisons among ethnic groups or polymorphism groups, we used either ANOVA or t-test for normally distributed variables, and Kruskal–Wallis or Wilcoxon rank sum test (non-parametric tests) for non-normally distributed variables. For linear regression analyses we log-transformed non-normally distributed variables, as described in the following paragraphs.
To compare metabolic traits measured in continuous scale among ethnic groups we used ANOVA—for normally distributed—or Kruskal–Wallis test—for non-normally distributed—variables, followed by post hoc Bonferroni´s test or Dunn´s test with Bonferroni adjustment, respectively, for pairwise comparisons. Significant differences between ethnic groups were indicated by different superscripts (a–e). When groups were not statistically different, they share the same letter; therefore, some groups may have more than one superscript letter. The frequency and percentage were used to describe categorical variables. The Chi2 test was used to compare groups and to verify Hardy-Weinberg equilibrium (HWE).
We formed two groups according to PPAR-γ2 and PPAR-δ genotypes (major vs. minor allele carriers). Comparison of metabolic traits between the two allelic groups was performed with t-test or Wilcoxon´s rank sum test for normal or non-normally distributed variables, respectively. The α value was adjusted for multiple comparisons by the Holm–Bonferroni method, to correct for family-wise error rate.
Multiple linear regression models were used to adjust the relation between metabolic traits (as continuous variables) and genotype groups for potential confounders and effect modifiers using a forward inclusion of variables. Metabolic trait variables with a non-normal distribution were log-transformed (TG, HDL, VLDL, and atherogenic index) to reach normality, only in the case of the diastolic blood pressure, we left the original values, because the distribution did not reach normality with any of the tested transformations. The model residuals were analyzed and a heteroskedasticity test was performed to verify linear regression assumptions. In further analyses, we included ethnic groups as dummy variables and as fixed effects, for the analyses performed in the entire sample. Because we found interaction effects between the polymorphism and ethnic groups, we conducted a stratified analysis, by ethnic group.
Dichotomous variables of metabolic traits were compared by Pearson’s Chi2 or Fisher exact test. Logistic regression models were used to calculate the odds ratio (OR) of metabolic traits (dichotomous) according to the genotype groups. The reference group was homozygous for the major allele for both PPAR polymorphisms. Logistic regression analysis was performed on the total studied population only due to the limited sample size in some ethnic groups, which led to very low frequencies in some categories. A p-value < 0.05 was considered to be statistically significant, unless otherwise specified. All of the statistical analyses were performed with STATA software (v. 11.0, StataCorp, College Station, TX, USA).

3. Results

3.1. Population Characteristics and Metabolic Traits among Them

Three main ethnic populations were studied, Mestizo, Tarahumara, and Mennonite. However, the Mestizo group was separated in two: Mestizo-1 and Mestizo-2, since they were from a different location and showed significant differences between them. The median age in the total studied population was 14 years and 60.92% were females. The Mestizo-2 group had the highest overweight/obese phenotype proportion according to their pBMI (Table 1).
Metabolic traits were compared among the ethnic populations (Table 2). Relevant differences were observed among them. The Mestizo-2 group showed greater values for BMI Z-score, waist circumference, and waist/height ratio, but lower glucose levels when compared to the other populations. The Tarahumara and Mestizo-2 populations had higher triglycerides values compared to Mestizo-1 and Mennonite, whereas the Tarahumara group had the lowest HDL levels with greater atherogenic index as compared to the others. The prevalence of altered metabolic traits in total studied population included overweight/obesity 20.15% (by BMI Z-score), high blood pressure 32.28%, high waist circumference 3.4%, high waist/height ratio 39.56%, hyperglycemia 1.21%, hypertriglyceridemia 25.73%, high cholesterol 12.14%, low HDL 32.11%, high LDL 3.88%, high VLDL 7.04%, and high atherogenic index 31.34%.

3.2. Pro12Ala PPAR-γ2 and +294T/C PPAR-δ

The genotype and allelic frequency distribution were compared among populations. For PPAR-γ2, no significant difference was observed among Mestizo-1, Mestizo-2, and Mennonite populations; however, the genotype distribution in Tarahumara was different between Mestizo-1 and Mennonite, but not with Mestizo-2 (Table 3). Similar results were observed for PPAR-δ genotype distribution—only the Tarahumara population showed a different genotype distribution (Table 4). The polymorphisms Pro12Ala PPAR-γ2 and +294T/C PPAR-δ were in Hardy–Weinberg equilibrium in all of the ethnic groups and total studied population.

3.3. Metabolic Traits by Pro12Ala PPAR-γ2 and +294T/C PPAR-δ Polymorphism Genotype: Major Allele Homozygous vs. Minor Allele Carriers in the Total Studied Population

Statistical analysis showed significant differences in waist/height ratio and HDL median values between PPAR-γ2 alleles. The α value was adjusted for multiple comparisons by the Holm–Bonferroni method, to correct for family-wise error rate. This association remained significant after adjusting for multiple testing (cut off p value = 0.0125), whereas the difference in waist/height ratio did not (cut off p value = 0.006). Individuals showed significant differences in TG, HDL, and VLDL median values, according with their PPAR-δ genotype. However, the association between PPAR-δ and HDL did not remain significant after correcting for multiple testing, while significance remained for TG and VLDL (Table 5).

3.4. Association between Pro12Ala PPAR-γ2 and Metabolic Traits

Among the clinical variables, sex was significantly associated (p < 0.05) with all metabolic traits, except for diastolic blood pressure and age was also associated with all, except for LDL and atherogenic index. Therefore, we included them as covariates in multivariate models. Associations between metabolic traits and the polymorphism Pro12Ala PPAR-γ2 were found with glucose, HDL, and atherogenic index in the total studied population, adjusting for clinical variables (age, sex, and or waist circumference). Minor allele carries had lower HDL and glucose levels and higher atherogenic index, than major allele homozygous. However, analysis according to ethnic populations only showed a significant association in the Mestizo-1 group with SBP and HDL. The minor allele carriers had ~5.59 mmHg lower SBP (adjusted by waist circumference, sex, and age) and ~0.90 mg/dL lower HDL than major allele homozygotes (log(−0.11), adjusted by waist circumference and sex) (Table 6). After also adjusting for ethnic group, only the association of atherogenic index remained significant. The interaction effect of sex and age for SBP showed that men had a slight elevation in SPB as age increased. The PPAR-γ2 and atherogenic index relation was modified by the waist circumference. Atherogenic index increased as waist circumference increased, but this relation differed by the studied polymorphism (Table 6 footnote). After also adjusting for ethnic group (last column), only the association of atherogenic index remained significant. An interaction effect between the polymorphism and the ethnic group was also observed for LDL.

3.5. Association between +294T/C PPAR-δ and Metabolic Traits

Associations between +294T/C PPAR-δ and metabolic traits in the total studied population were observed with glucose, TG, TC, and VLDL. C allele carries had 2.57 mg/dL less glucose (adjusted by waist circumference, sex and age), 0.89 mg/dL less TG (log(−0.12), adjusted by sex and age), 0.89 mg/dL less VLDL (log(−0.12), adjusted by waist circumference, sex and age), and 79.02 mg/dL more TC (adjusted by sex and age) than major allele homozygotes. Among ethnic groups, the only association was found in the Tarahumara population with HDL where C allele carriers had 0.90 mg/dL less HDL (log(−0.10), adjusted by waist circumference, sex, and age) than major allele homozygotes (Table 7). The relation between PPAR-δ allele and HDL was modified by sex. In male teenagers, the HDL levels were higher in C allele carriers versus major allele homozygotes; no differences were noted in females (Figure 1A). In the total studied population, PPAR-δ and LDL relation was also modified by sex. Males with the major allele had the lowest values of LDL versus females and males with minor allele carriers (Figure 1B). The same interaction was found in Mestizo-1 population (Table 7 footnote). After adjusting for ethnic group (Table 7, last column), the association of the polymorphism and glucose levels remained significant, as well as the interaction effect between the polymorphism and sex for LDL. We also observed a significant interaction effect between ethnic group and this polymorphism in its relation to HDL and atherogenic index.

3.6. Logistic Regression

In the total studied population, Pro12Ala PPAR-γ2 minor allele carriers had greater likelihood of having an elevated waist-height ratio (OR = 1.79, adjusted by sex) and 2.5-fold greater chance of having low HDL values (adjusted by waist circumference, sex, and age) than major allele homozygous. In contrast, they had a nearly two-fold lower likelihood of having a BMI Z-score suggestive of overweight/obesity (OR = 0.46) and increased TC (OR = 0.47, adjusted by sex and age) than major allele homozygous. The PPAR-δ minor allele carriers had greater likelihood (OR = 1.78) of having a BMI Z-score that was suggestive of overweight/obesity (Table 8). No significant interactions were seen between both minor alleles.

4. Discussion

In this study, we analyzed the association between the Pro12Ala PPAR-γ2 or +294T/C PPAR-δ polymorphisms and metabolic traits in teenagers from the North of Mexico, including Mestizo, Tarahumara, and Mennonite populations. Hispanic populations have a higher tendency for obesity than other Caucasians due to their genetic origin and cultural characteristics [45,46,47]. Mexican-origin population have Native American and Spanish admixture. However, Salzano et al. reported that Amerindian ancestry is most prevalent in the general population of Mexico followed by European ancestry [48]. The Mexican territory is extensive, it borders with countries of North and Central America. In addition, Mexico has a great history of Spanish colonization and the arrival of a large number of immigrants of African and European origin. Therefore, it is not surprising to find diversity in the population in this country. The interest of this study was to include the three most representative populations of the Chihuahua state, Northern of Mexico.
The overweight/obesity prevalence in the total studied population was 20.15%, which is lower than the national prevalence reported among Mexican teenagers in 2018 (38.4%) [4]. However, it was similar to the US national report in 2015 (20.6%) [49]. Mestizo-2 group had the highest overweight/obesity prevalence (70%) when compared to the other ethnic populations—this might be because they are from an urbanized location. Benitez et al. found that urban Tarahumara children were 10% more overweight than rural Tarahumara children in Chihuahua, Mexico [50]. The high blood pressure prevalence in the total studied population was 32.28% similar to the frequency that was reported by Salcedo-Rocha et al. (31.1%) and higher than that found by Cardoso-Saldaña et al. (19.2%) in Mexican adolescents [51,52]. The results in Mexicans are significantly higher than values reported by the NHANES 2015 in American teenagers (4.2%) [53]. In contrast, the hypertriglyceridemia and HDL prevalence in this study was lower than most reported in Mexican teenagers [54,55]. The prevalence of hyperglycemia was low (1.12%) similar to the reported by Cardoso-Saldaña et al. (1.7%), but different the reported by Camarillo-Romero et al. (4.9%) [52,56]. Among the four ethnic groups analyzed in this study, the Mennonites were the healthiest, whereas the Mestizo-2 population had the greater number of metabolic traits above normal values. The metabolic indicators we chose are considered metabolic syndrome risk factors for adults. Even though the study population appeared to be in healthy individuals—no clinical setting—there was a percentage of teenagers with values outside the expected for healthy populations. This suggests that at these early ages, some metabolic risks could be detected, thus providing an opportunity window for intervention targeted to modifiable factors of obesity and chronic non-transmissible diseases.
PPARs play complex, overlapping, and specific roles in metabolic networks by regulating cellular energy homeostasis during lipid and carbohydrate metabolism. Therefore, PPARs have been considered as potential therapeutic targets [6]. However, the allele frequencies vary among different populations. The Pro12Ala PPAR-γ2 G minor allele frequency in the literature ranges from 9% to 25% [9,15,16,17,57,58]. In this study, the frequency in the total studied population was 17%, while other studies in Mexican population report 10–21% [11,17]. Similarly, the +294T/C PPAR-δ C minor allele frequencies vary from 12% to 31% among different populations [22,24,25,27,28,29,30,35,59,60,61]. We found that the +294T/C PPAR-δ C allele frequency in the total studied population was 12%, which is lower than what was found in other Mexican population (26%) by Rosales-Reynoso et al. [36]. Differences in both PPAR-γ2 and PPAR-δ allelic and genotype frequencies were only observed in the Tarahumara group. This might be explained by differences in ethnic admixture; however, further studies are needed.
The Pro12Ala PPAR-γ2 G minor allele has been associated with hypertension, insulin resistance, greater BMI, and lipemia profiles, mainly in adult population. With the multiple linear regression analysis, we found significant associations with several metabolic traits: systolic blood pressure, glucose, HDL, and atherogenic index. However, except for the atherogenic index, these associations were no longer significant after adjusting for ethnic group. Hasan et al. in Egyptian populations reported that the polymorphism was related to lower SBP and DBP values in patients with T2D and CHD [62]. However, most literature reports an opposite association with higher blood pressure values (SBP and DBP) in individuals from the USA, Finland, and Spain [63,64,65]. We also found a relation between the Pro12Ala PPAR-γ2 G allele and lower glucose values (adjusted by waist circumference, sex, and age). Most studies showed no association between the polymorphism and glucose levels in other populations, such as Tunisian (adjusted by age and BMI), Cypriot, Egyptian, and Emirati populations [18,58,62,66]. Other studies have reported an association of minor allele with increased insulin [15,16]. It would have been interesting to have insulin data in this study.
In terms of lipemia, we only found an association between the Pro12Ala PPAR-γ2-minor allele and lower HDL values in the total studied population and the Mestizo-1 group (adjusted by waist circumference and sex). After the adjustment for ethnic group, this relation was no significant. In contrast, to the higher HDL levels reported in Finnish population, Becer et al. found no association between the SNP and lipemia traits in non-obese Cypriot subjects; however, obese subjects with homozygous minor genotype had increased TG values [15,58]. This suggests that obesity may interact with Pro12Ala PPAR-γ2 G minor allele and lipemia traits. As expected, we found an association between the Pro12Ala PPAR-γ2 G allele and higher atherogenic index.
The +294T/C PPAR-δ polymorphism has been related to HDL, LDL, VLDL, and TG. With a multiple linear regression analysis, we found significant associations between the PPAR-δ polymorphism and glucose, TG, total cholesterol, HDL, and VLDL. The presence of the minor allele was associated with lower glucose values in this study (adjusted by waist circumference, sex and age). In contrast, Hu et al. found that the minor allele associated with higher glucose concentrations in Shanghai, China population (adjusted by age and sex) [67]. Whereas, no association between the C allele and glucose was observed in other Shanghai study (adjusted by age) and in Denmark subjects (adjusted by age, sex, and BMI) [30,31].
The +294T/C PPAR-δ C minor allele was also associated with lower TG (adjusted by waist circumference, sex and age), higher total cholesterol (adjusted by sex and age), and lower VLDL (adjusted by waist circumference, sex, and age) values in the total studied population. However, an association with lower HDL values was also found in the Tarahumara group (adjusted by waist circumference, sex, and age). This polymorphism has been associated with HDL since its discovery by Skogsberg et al. who observed its relation with lower HDL in Scottish population [24]. Similar results were reported by Aberle et al. in German woman, where +294T/C PPAR-δ-minor allele carriers had lower HDL and VLDL values (adjusted by age, smoking, and BMI) [26]. In contrast, association between the SNP and higher HDL values has been reported in Scottish and Canadian studies [27,28]. No relation was found with SBP in this study, while others have related the minor allele with lower blood pressure in Chinese population [32,68].
The logistic regression analysis with dichotomous metabolic traits values indicated that individuals with the Pro12Ala PPAR-γ2 G allele had greater risk of having an elevated waist/height ratio (adjusted by sex). Similar results were reported by Sözen et al. in Turkish populations with obesity, where the minor allele was associated with higher values of waist/height ratio; however, the same SNP was also associated with lower values in the non-obese group [69]. Whereas, Meirhaeghe et al. found no association [70]. In this study the Pro12Ala G allele carriers had two-fold less possibilities of being overweight/obese (adjusted by sex and age). In contrast to most studies in the literature, such as different meta-analysis that have shown a BMI increment in minor allele carriers [71,72]. In this study, we also found that Pro12Ala-minor allele carriers had 2.5-fold more chances of having low HDL values when adjusted by waist circumference, sex, and age and two-fold more total cholesterol (adjusted by sex and age). Regarding +294T/C PPAR-δ, we observed that C minor allele carriers had greater probability of having an overweight/obese phenotype. In contrast, this polymorphism has been related with lower BMI in Scottish males, as well as German and Chinese adult populations [27,59,73]. Nevertheless, an association was found between +294T/C minor allele and BMI in Greek toddlers [59]. These analyses were not adjusted for ethnic group, due to the low frequency in several categories. Taking into account the observed effects in the metabolic traits analyzed as continuous variables, these results may vary across the groups.
Metabolic diseases are of public health concern. For years, health programs have been focused on how to maintain normal values of clinical indicators. Unfortunately, most actions are taken in adulthood when it is already difficult to change a person’s habits. However, it should be considered that metabolic diseases depend on modifiable and non-modifiable factors. The latter are those that are related to genetic load. Ideally, the era of “Omics” aims to create networks between the different levels of information that an individual has (genome, transcriptome, metabolome, proteome, etc.). We consider that the detection of risk components for metabolic diseases at an early age would allow us to improve preventive programs for health care. To the best of our knowledge, this is the first study to report differences among ethnic populations that live in Chihuahua State from Mexico regarding PPAR-γ2 and PPAR-δ and their association/interaction with metabolic traits. Population admixtures, habits, and cultural aspects need to be considered in public health strategies. Networking disciplines, including genetics, may help to identify new therapeutic answers, and PPARs are potential therapeutic targets to control metabolic diseases. Further studies in the same ethnic populations with larger sample size and including adults will be important to confirm our findings.

5. Study Limitations

The cross-sectional nature of this study does not allow to follow up if participants would develop a metabolic disease in the future associated with the studied polymorphisms. The sample size by ethnic group was small for logistic regression analysis and the variable diastolic blood pressure did not reach normality with any transformation, then we used the original scale. On the other hand, the initial study was focused on early detection of metabolic syndrome components in teenagers, thus most of the studied population was in good health, contrary to most studies that include adults with specific inclusion criteria. Moreover, metabolic diseases, such as diabetes type 2, hypertension, and obesity, are multifactorial, meaning that there are external factors, such as environmental conditions, diet, and exercise, which can be modified. However, we did not explore those issues. Multifactorial diseases are also multigenic, whereas, in this study, we just included two PPAR polymorphisms and multiple SNPs analysis could still give more information. Finally, the schools were not randomly selected to assure the sample representativity.

6. Conclusions

The results of this study showed that Pro12Ala PPAR-γ2 and +294T/C PPAR-δ polymorphisms are associated with metabolic indicators among ethnic populations. The total studied population showed a significant relationship between the PPAR-γ2 and PPAR-δ polymorphisms with some metabolic traits. This PPAR polymorphism study seems to be the first that includes different ethnic teenagers from the North of Mexico and, to our knowledge, it is the first report to estimate +294T/C PPAR-δ minor allele frequencies among them. An early metabolic syndrome components detection will provide the opportunity to implement prophylactic measures on modifiable factors. Further research is warranted in order to confirm the role of other PPAR polymorphisms on metabolic diseases.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/7/776/s1, Figure S1: Electrophoresis example for both polymorphisms.

Author Contributions

Conceptualization, V.L.-C. and I.L.-B.; Data curation, J.G.M.-G.; Formal analysis, S.A.R.-L., M.A.C.-V., I.L.-B.; Funding acquisition, I.L.-B.; Investigation, I.L.-B.; Methodology, J.G.M.-G., M.A.C.-V., and N.R.E.-A.; Project administration, I.L.-B.; Resources, P.d.C.H.-R.; Supervision, V.M.-B. and A.L.-T.; Visualization, E.G.-R.; Writing–original draft, M.A.C.-V.; Writing–review & editing, I.L.-B., S.A.R.-L., M.A.C.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This research study was supported by the Universidad Autónoma de Chihuahua (Facultad de Medicina y Ciencias Biomédicas) with PFCE-2018-08MSU0017H-05 and publication expenses will be covered by Universidad Autónoma de Chihuahua.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

References

  1. Afshin, A.; Forouzanfar, M.H.; Reitsma, M.B.; Sur, P.; Estep, K.; Lee, A.; Marczak, L.; Mokdad, A.H.; Moradi-Lakeh, M.; Naghavi, M.; et al. Health effects of overweight and obesity in 195 countries over 25 years. New Engl. J. Med. 2017, 377, 13–27. [Google Scholar] [PubMed]
  2. OECD. The Heavy Burden of Obesity; OECD Health Policy Studies, Ed.; OECD Publishing: Paris, France, 2019. [Google Scholar]
  3. Ogden, C.L.; Carroll, M.D.; Lawman, H.G.; Fryar, C.D.; Kruszon-Moran, D.; Kit, B.K.; Flegal, K.M. Trends in obesity prevalence among children and adolescents in the United States, 1988-1994 through 2013-2014. JAMA 2016, 315, 2292–2299. [Google Scholar] [CrossRef] [PubMed]
  4. Instituto Nacional de Salud Pública. Encuesta Nacional de Salud y Nutrición 2018. México: Instituto Nacional de Salud Pública. 2018. Available online: https://www.inegi.org.mx/contenidos/programas/ensanut/2018/doc/ensanut_2018_diseno_conceptual.pdf (accessed on 19 June 2020).
  5. Willson, T.M.; Brown, P.J.; Sternbach, D.D.; Henke, B.R. The PPARs: From orphan receptors to drug discovery. J. Med. Chem. 2000, 43, 527–550. [Google Scholar] [CrossRef]
  6. Hong, F.; Pan, S.; Guo, Y.; Xu, P.; Zhai, Y. PPARs as Nuclear Receptors for Nutrients and Energy Metabolism. Molecules 2019, 24, 2545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Chen, Y.; Jiménez, A.R.; Medh, J.D. Identification and regulation of novel PPAR-γ splice variants in human THP-1 macrophages. Biochim. Biophys. Acta 2006, 1759, 32–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gouda, H.N.; Sagoo, G.S.; Harding, A.H.; Yates, J.; Sandhu, M.S.; Higgins, J.P.T. The association between the peroxisome proliferator-activated receptor-γ2 (PPARG2) Pro12Ala gene variant and type 2 diabetes mellitus: A HuGE review and meta-analysis. Am. J. Epidemiol. 2010, 171, 645–655. [Google Scholar] [CrossRef] [PubMed]
  9. Swarbrick, M.M.; Chapman, C.M.L.; McQuillan, B.M.; Hung, J.; Thompson, P.L. A Pro12Ala polymorphism in the human peroxisome proliferator-activated receptor-γ2 is associated with combined hyperlipidaemia in obesity. Eur. J. Endocrinol. 2001, 144, 277–282. [Google Scholar] [CrossRef] [Green Version]
  10. Stumvoll, M.; Häring, H. The Peroxisome Proliferator-Activated Receptor-γ2 Pro12Ala Polymorphism. Perspect. Diabetes 2002, 51, 2341–2347. [Google Scholar]
  11. Stryjecki, C.; Peralta-Romero, J.; Alyass, A.; Karam-Araujo, R.; Suarez, F.; Gomez-Zamudio, J.; Burguete-Garcia, A.; Cruz, M.; Meyre, D. Pro12Ala genotype and insulin resistance is modified by circulating lipids in Mexican children. Sci. Rep. 2016, 6, 1–7. [Google Scholar] [CrossRef] [Green Version]
  12. Yen, C.; Beamer, B.A.; Negri, C.; Silver, K.; Yarnall, D.P. Molecular Scanning of the Human Peroxisome Proliferator Activated Receptor γ (hPPARγ) Gene in Diabetic Caucasians: Identification of a Pro12Ala γ2 Missense Mutation. Biochem. Biophys. Res. Commun. 1997, 274, 270–274. [Google Scholar] [CrossRef]
  13. Galbete, C.; Toledo, E.; Martínez-González, M.A.; Martínez, J.A.; Guillén-Grima, F.; Marti, A. Pro12Ala variant of the PPARG2 gene increases body mass index: An updated meta-analysis encompassing 49,092 subjects. Obesity 2013, 21, 1486–1495. [Google Scholar] [CrossRef]
  14. Morris, A.; Voight, B.; Teslobich, T.; Ferreira, T.; Segre, A.; Steinthorsdottir, B.; Strawbridge, R.J.; Khan, H.; Grallert, H.; Mahajan, A.; et al. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat. Genet. 2012, 44, 981–990. [Google Scholar] [PubMed]
  15. Deeb, S.S.; Fajas, L.; Nemoto, M.; Pihlajamäki, J.; Mykkänen, L.; Kuusisto, J.; Laakso, M.; Fujimoto, W.; Auwerx, J. A Pro12Ala substitution in PPARγ2 associated with decreased receptor activity, lower body mass index and improved insulin sensitivity. Nat. Genet. 1998, 20, 284–287. [Google Scholar] [CrossRef]
  16. Muller, Y.L.; Bogardus, C.; Beamer, B.A.; Shuldiner, A.R.; Baier, L.J. A Functional Variant in the Peroxisome Proliferator—Activated Receptor γ2 Promoter Is Associated with Predictors of Obesity and Type 2 Diabetes in Pima Indians. Diabetes 2003, 52, 1864–1871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Canizales-Quinteros, S.; Aguilar-Salinas, C.A.; Ortiz-López, M.G.; Rodríguez-Cruz, M.; Villarreal-Molina, M.T.; Coral-Vázquez, R.; Huertas-Vázquez, A.; Hernández-Caballero, A.; López-Alarcón, M.; Brito-Zurita, O.R.; et al. Association of PPARG2 Pro12Ala Variant with Larger Body Mass Index in Mestizo and Amerindian Populations of Mexico. Hum. Biol. 2007, 79, 111–119. [Google Scholar] [CrossRef] [PubMed]
  18. Ben Ali, S.; Ben Yahia, F.; Sediri, Y.; Kallel, A.; Ftouhi, B. Gender-specific effect of Pro12Ala polymorphism in peroxisome proliferator-activated receptor γ-2 gene on obesity risk and leptin levels in a Tunisian population. Clin. Biochem. 2009, 42, 1642–1647. [Google Scholar] [CrossRef] [PubMed]
  19. Baldani, D.P.; Skrgatic, L.; Cerne, J.Z.; Ferk, P.; Simunic, V.; Gersak, K. Association of PPARG Pro12Ala polymorphism with insulin sensitivity and body mass index in patients with polycystic ovary syndrome. Biomed. Rep. 2014, 2, 199–206. [Google Scholar] [CrossRef] [Green Version]
  20. Tönjes, A.; Stumvoll, M. The role of the Pro12Ala polymorphism in peroxisome proliferator-activated receptor γ in diabetes risk. Lippincott Williams & Wilkins. Curr. Opin. Clin. Nutr. Metab. Care 2007, 10, 410–414. [Google Scholar]
  21. Information NC for B. PPARD peroxisome proliferator activated receptor delta [Homo sapiens (human)]. 2018. Available online: https://www.ncbi.nlm.nih.gov/gene/5467#gene-expression (accessed on 2 November 2018).
  22. Gouni-Berthold, I.; Giannakidou, E.; Faust, M.; Berthold, H.K.; Krone, W. The peroxisome proliferator-activated receptor delta +294T/C polymorphism in relation to lipoprotein metabolism in patients with diabetes mellitus type 2 and in non-diabetic controls. Atherosclerosis 2005, 183, 336–341. [Google Scholar] [CrossRef]
  23. Sznaidman, M.L.; Haffner, C.D.; Maloney, P.R.; Fivush, A.; Chao, E.; Goreham, D.; Sierra, M.L.; LeGrumelec, C.; Xu, H.E.; Montana, V.G.; et al. Novel selective small molecule agonists for peroxisome proliferator-activated receptor δ (PPARδ) - Synthesis and biological activity. Bioorg Med. Chem. Lett. 2003, 13, 1517–1521. [Google Scholar] [CrossRef]
  24. Skogsberg, J.; McMahon, A.D.; Karpe, F.; Hamsten, A.; Packard, C.J.; Ehrenborg, E. Peroxisome proliferator activated receptor delta genotype in relation to cardiovascular risk factors and risk of coronary heart disease in hypercholesterolaemic men. J. Intern. Med. 2003, 254, 597–604. [Google Scholar] [CrossRef] [PubMed]
  25. Luo, C.Y.; Liu, C.W.; Ge, L.; Pang, G.F.; Yang, M.; Hu, C.Y.; Ze-Ping, L.; Ning-Yuan, C.; Hai-Yan, L.; Hua-Yu, W.; et al. PPARD +294C overrepresentation in general and long-lived population in China Bama longevity area and unique relationships between PPARD +294T/C polymorphism and serum lipid profiles. Lipids Health Dis. 2015, 14, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Aberle, J.; Hopfer, I.; Beil, F.U.; Seedorf, U. Association of the T+294C polymorphism in PPAR δ with low HDL cholesterol and coronary heart disease risk in women. Int. J. Med Sci. 2006, 3, 108–111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Burch, L.R.; Donnelly, L.A.; Doney, A.S.F.; Brady, J.; Tommasi, A.M.; Whitley, A.L.; Goddard, C.; Morris, A.D.; Hansen, M.K.; Palmer, C.N.A. Peroxisome proliferator-activated receptor-δ genotype influences metabolic phenotype and may influence lipid response to statin therapy in humans: A genetics of diabetes audit and research tayside study. J. Clin. Endocrinol. Metab. 2010, 95, 1830–1837. [Google Scholar] [CrossRef] [Green Version]
  28. Robitaille, J.; Gaudet, D.; Pérusse, L.; Vohl, M.-C. Features of the metabolic syndrome are modulated by an interaction between the peroxisome proliferator-activated receptor-delta -87T>C polymorphism and dietary fat in French-Canadians. Int. J. Obes. 2007, 31, 411–417. [Google Scholar]
  29. Aberle, J.; Hopfer, I.; Beil, F.U.; Seedorf, U. Association of peroxisome proliferator-activated receptor delta +294T/C with body mass index and interaction with peroxisome proliferator-activated receptor alpha L162V. Int. J. Obes. 2006, 30, 1709–1713. [Google Scholar] [CrossRef] [Green Version]
  30. Grarup, N.; Albrechtsen, A.; Ek, J.; Borch-Johnsen, K.; Jørgensen, T.; Schmitz, O.; Pedersen, O. Variation in the peroxisome proliferator-activated receptor δ gene in relation to common metabolic traits in 7,495 middle-aged white people. Diabetologia 2007, 50, 1201–1208. [Google Scholar] [CrossRef] [Green Version]
  31. Villegas, R.; Williams, S.; Gao, Y.; Cai, Q.; Li, H.; Elasy, T.; Cai, H.; Edwards, T.; Xiang, Y.B.; Zheng, W.; et al. Peroxisome proliferator-activated receptor (PPAR) delta genetic variation and type 2 diabetes in middle age Chinese women. Ann. Hum. Genet. 2011, 75, 621–629. [Google Scholar] [CrossRef] [Green Version]
  32. Tang, L.; Lü, Q.; Cao, H.; Yang, Q.; Tong, N. PPARD rs2016520 polymorphism is associated with metabolic traits in a large population of Chinese adults. Gene 2016, 585, 191–195. [Google Scholar] [CrossRef]
  33. Qian, Y.; Li, P.; Zhang, J.; Shi, Y.; Chen, K.; Yang, J.; Wu, Y.; Ye, X. Association between peroxisome proliferator-activated receptor-alpha, delta, and gamma polymorphisms and risk of coronary heart disease: A case-control study and meta-analysis. Medicine 2016, 95, 1–9. [Google Scholar] [CrossRef]
  34. Jguirim-Souissi, I.; Jelassi, A.; Hrira, Y.; Najah, M.; Slimani, A.; Addad, F.; Hassine, M.; Hamda, K.B.; Maatouk, F.; Rouis, M.; et al. +294T/C polymorphism in the PPAR-δ gene is associated with risk of coronary artery disease in normolipidemic Tunisians. Genet. Mol. Res. 2010, 9, 1326–1333. [Google Scholar] [CrossRef]
  35. Miao, L.; Yin, R.-X.; Wu, D.-F.; Cao, X.-L.; Li, Q.; Hu, X.-J.; Yan, T.-T.; Aung, L.H.H.; Yang, D.-Z.; Lin, W.-X.; et al. Peroxisome proliferator-activated receptor delta +294 > C polymorphism and serum lipid levels in the Guangxi Bai Ku Yao and Han populations. Lipids Health Dis. 2010, 9, 145–157. [Google Scholar]
  36. Rosales-Reynoso, M.A.; Wence-Chavez, L.I.; Arredondo-Valdez, A.R.; Dumois-Petersen, S.; Barros-Núñez, P.; Gallegos-Arreola, M.P.; Flores-Martínez, S.E.; Sánchez-Corona, J. Protective role of +294 T/C (rs2016520) polymorphism of PPARD in Mexican patients with colorectal cancer. Genet. Mol. Res. 2017, 16, 1–10. [Google Scholar] [CrossRef] [PubMed]
  37. Cook, S.; Weitzman, M.; Auinger, P.; Michael, N.; Dietz, W.H. Prevalence of a Metabolic Syndrome Phenotype in Adolescents. Arch. Pediatr. Adolescent Med. 2003, 157, 821–827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Mederico, M.; Paoli, M.; Zerpa, Y.; Camacho, N.; Cichetti, R.; Molina, Z.; Camachob, N.; Cichettib, R.; Molinab, Z.; Mora, Y.; et al. Valores de referencia de la circunferencia de la cintura e índice de la cintura/cadera en escolares y adolescentes de Mérida, Venezuela: Comparación con referencias internacionales. Endocrinol. Y Nutr. 2013, 60, 235–242. [Google Scholar] [CrossRef] [PubMed]
  39. Mokha, J.S.; Srinivasan, S.R.; DasMahapatra, P.; Fernandez, C.; Chen, W.; Xu, J.; Berenson, G.S. Utility of waist-to-height ratio in assessing the status of central obesity and related cardiometabolic risk profile among normal weight and overweight/obese children: The Bogalusa Heart Study. BMC Pediatr. 2013, 10, 73. [Google Scholar]
  40. CDC. Defining Childhood Obesity. 2018. Available online: https://www.cdc.gov/obesity/childhood/defining.html (accessed on 19 June 2020).
  41. Must, A.; Anderson, S.E. Body mass index in children and adolescents: Considerations for population-based applications. Int. J. Obes. 2006, 30, 590–594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. World Health Organization. Growth Reference 5-19 years. 2007. Available online: https://www.who.int/growthref/who2007_bmi_for_age/en/ (accessed on 19 June 2020).
  43. National Heart Lung and Blood Institute. Expert Panel on Integrated Guidelines for Cardiovascular Health and Risk Reduction in Children and Adolescents. 2012. Available online: http://www.nhlbi.nih.gov/files/docs/guidelines/peds_guidelines_full.pdf (accessed on 19 June 2020).
  44. Millán, J.; Pintó, X.; Muñoz, A.; Zúñiga, M.; Rubiés-Prat, J.; Pallardo, L.F.; Masana, L.; Mangas, A.; Hernández-Mijares, A.; González-Santos, P.; et al. Cocientes lipoproteicos: Significado fisiológico y utilidad clínica de los índices aterogénicos en prevención cardiovascular. Clin. E Investig. En Arterioscler. 2010, 22, 25–32. [Google Scholar] [CrossRef]
  45. Cossrow, N.; Falkner, B. Race/ethnic issues in obesity and obesity-related comorbidities. J. Clin. Endocrinol. Metab. 2004, 89, 2590–2594. [Google Scholar] [CrossRef] [Green Version]
  46. Anderson, S.E.; Whitaker, R.C. Prevalence of obesity among US preschool children in different racial and ethnic groups. Arch. Pediatrics Adolesc. Med. 2009, 163, 344–348. [Google Scholar] [CrossRef] [Green Version]
  47. Taveras, E.M.; Gillman, M.W.; Kleinman, K.; Rich-Edwards, J.W.; Rifas-Shiman, S.L. Racial/Ethnic Differences in Early-Life Risk Factors for Childhood Obesity. Pediatrics 2010, 125, 686–695. [Google Scholar] [CrossRef] [Green Version]
  48. Salzano, F.M.; Sans, M. Interethnic admixture and the evolution of Latin American populations. Genetics Mol. Biol. 2014, 37 (Suppl. 1), 151–170. [Google Scholar] [CrossRef] [Green Version]
  49. Hales, C.M.; Carroll, M.D.; Fryar, C.D.; Ogden, C.L. Prevalence of Obesity Among Adults and Youth: United States, 2015–2016; NCHS data brief; 2017; 288, pp. 1–8. Available online: https://www.cdc.gov/nchs/products/databriefs/db288.htm (accessed on 19 June 2020).
  50. Benítez-Hernández, Z.P.; Hernández-Torres, P.; Cabañas, M.D.; De La Torre-Díaz, M.D.L.; López-Ejeda, N.; Marrodán, M.D.; Cervantes-Borunda, M. Composición corporal, estado nutricional y alimentación en escolares Tarahumaras urbanos y rurales de Chihuahua, México. Nutr. Clínica Y Dietética Hosp. 2014, 34, 71–79. [Google Scholar]
  51. Salcedo-Rocha, A.L.; García de Alba, J.E.; Contreras-Marmolejo, M. Presión arterial en adolescentes mexicanos: Clasificación, factores de riesgo e importancia. Rev. De Salud Pública 2010, 12, 612–622. [Google Scholar]
  52. Cardoso-Saldaña, G.C.; Yamamoto-Kimura, L.; Medina-Urrutia, A.; Posadas-Sánchez, R.; Caracas-Portilla, N.A.; Posadas-Romero, C. Exceso de peso y síndrome metabólico en adolescentes de la Ciudad de México. Archivos de Cardiología de México 2010, 80, 12–18. [Google Scholar]
  53. Jackson, S.L.; Zhang, Z.; Wiltz, J.L.; Loustalot, F.; Ritchey, M.D.; Goodman, A.B.; Yang, Q. Hypertension Among Youths—United States, 2001–2016. Morb. Mortal. Wkly. Rep. 2018, 67, 758–762. [Google Scholar] [CrossRef]
  54. Romero-Velarde, E.; Campollo-Rivas, O.; De La Rosa, A.C.; Vásquez-Garibay, E.M.; Castro-Hernández, J.F.; Cruz-Osorio, R.M. Factores de riesgo de dislipidemia en niños y adolescentes con obesidad. Salud Pública De México 2007, 49, 103–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Perea-Martínez, A.; Bárcena-Sobrino, E.; Rodríguez-Herrera, R.; Greenawalt-Rodríguez, S.; Carbajal-Rodríguez, L.; Zarco-Román, J. Obesidad y comorbilidades en niños y adolescentes asistidos en el Instituto Nacional de Pediatría. Acta Pediátrica de México 2009, 30, 167–174. [Google Scholar]
  56. Camarillo-Romero, E.; García, M.V.D.; Amaya-Chávez, A.; Huitrón-Bravo, G.; Majluf-Cruz, A. Dificultades en la clasificación del síndrome metabólico. El ejemplo de los adolescentes en México. Salud Pública de México 2010, 52, 524–527. [Google Scholar] [PubMed]
  57. Almeida, S.M.; Furtado, J.M.; Mascarenhas, P.; Ferraz, M.E.; Ferreira, J.C.; Monteiro, M.P.; Vilanova, M.; Ferraz, F.P. Association between LEPR, FTO, MC4R, and PPARG-2 polymorphisms with obesity traits and metabolic phenotypes in school-aged children. Endocrine 2018, 60, 466–478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Becer, E.; Çlrakoǧlu, A. Effect of the Pro12Ala polymorphism of the peroxisome proliferator-activated receptor γ2 gene on lipid profile and adipokines levels in obese subjects. Balk. J. Med Genet. 2017, 20, 71–79. [Google Scholar] [CrossRef] [Green Version]
  59. Lagou, V.; Scott, R.A.; Manios, Y.; Chen, T.L.J.; Wang, G.; Grammatikaki, E.; Kortsalioudaki, C.; Liarigkovinos, T.; Moschonis, G.; Roma-Giannikou, E.; et al. Impact of peroxisome proliferator-activated receptors γ and δ on adiposity in toddlers and preschoolers in the GENESIS study. Obesity 2008, 16, 913–918. [Google Scholar] [CrossRef]
  60. Akhmetov, I.; Astranenkove, I.; Togozkin, V. Association of PPARD gene polymorphism with human physical performance. Mol. Biol. 2007, 41, 852–857. [Google Scholar]
  61. Shin, H.D.; Park, B.L.; Kim, L.H.; Jung, H.S.; Cho, Y.M.; Moon, M.K.; Park, Y.J.; Lee, H.K.; Park, K.S. Genetic Polymorphisms in Peroxisome Proliferator-Activated Receptor δ Associated with Obesity. Diabetes 2004, 53, 847–851. [Google Scholar] [CrossRef] [Green Version]
  62. Hasan, N.S.; Kamel, S.A.; Hamed, M.; Awadallah, E.; Rahman, A.H.A.; Musa, N.I.; Hussein, G.H.S. Peroxisome proliferator-activated receptor-γ polymorphism (rs1801282) is associated with obesity in Egyptian patients with coronary artery disease and type 2 diabetes mellitus. J. Genet. Eng. Biotechnol. 2017, 15, 409–414. [Google Scholar] [CrossRef]
  63. Hasstedt, S.J.; Ren, Q.F.; Teng, K.; Elbein, S.C. Effect of the peroxisome proliferator-activated receptor-γ2 Pro12Ala variant on obesity, glucose homeostasis, and blood pressure in members of familial type 2 diabetic kindreds. J. Clin. Endocrinol. Metab. 2001, 86, 536–541. [Google Scholar] [CrossRef] [Green Version]
  64. Douglas, J.A.; Erdos, M.R.; Watanabe, R.M.; Braun, A.; Johnston, C.L.; Oeth, P.; Mohlke, K.L.; Valle, T.T.; Ehnholm, C.; Buchanan, T.A.; et al. The Peroxisome Proliferator-Activated Receptor-γ2 Pro12Ala Variant Association with Type 2 Diabetes and Trait Differences. Diabetes 2001, 50, 886–890. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Rodríguez-Esparragón, F.J.; Rodríguez-Pérez, J.C.; Macías-Reyes, A.; Alamo-Santana, F. Peroxisome proliferator-activated receptor-gamma2-Pro12Ala and endothelial nitric oxide synthase-4a/b gene polymorphisms are associated with essential hypertension. J. Hypertens. 2003, 21, 1649–1655. [Google Scholar] [CrossRef] [PubMed]
  66. Al-Safar, H.; Hassoun, A.; Almazrouei, S.; Kamal, W.; Afandi, B.; Rais, N. Association of the genetic polymorphisms in transcription factor 7-like 2 and peroxisome proliferator-activated receptors-γ2 with type 2 diabetes mellitus and its interaction with obesity status in Emirati population. J. Diabetes Res. 2015, 2015, 129695. [Google Scholar] [CrossRef] [Green Version]
  67. Hu, C.; Jia, W.; Fang, Q.; Zhang, R.; Wang, C.; Lu, J.; Xiang, K. Peroxisome proliferator-activated receptor (PPAR) delta genetic polymorphism and its association with insulin resistance index and fasting plasma glucose concentrations in Chinese subjects. Diabet. Med. 2006, 23, 1307–1312. [Google Scholar] [CrossRef]
  68. Ye, H.; Li, Y.; Hong, Q.; Zhou, A.; Zhao, Q.; Xu, L.; Xu, M.Q.; Xu, X.T.; Tang, L.L.; Dai, J.; et al. Positive association between PPARD rs2016520 polymorphism and coronary heart disease in a Han Chinese population. Genet. Mol. Res. 2015, 14, 12254–12261. [Google Scholar] [CrossRef] [PubMed]
  69. Sözen, M.A.; Özcan, M.U.; Çıldır, M.; Doğru, I.H.; Aygök, A.G.; Balkan, K. Association of the human PPARγ2 Pro12Ala polymorphism with obesity in a population from Turkey. Acta Endocrinol. 2018, 14, 459–465. [Google Scholar]
  70. Meirhaeghe, A.; Fajas, L.; Helbecque, N.; Cottel, D.; Auwerx, J.; Deeb, S.S.; Amouyel, P. Impact of the Peroxisome Proliferator Activated Receptor γ2 Pro12Ala polymorphism on adiposity, lipids and non-insulin-dependent diabetes mellitus. Int. J. Obes. 2000, 24, 195–199. [Google Scholar] [CrossRef] [Green Version]
  71. Masud, S.; Ye, S. Effect of the peroxisome proliferator activated receptor-y gene Pro12Ala variant on body mass index: A meta-analysis. J. Med. Genetics 2003, 40, 773–780. [Google Scholar] [CrossRef] [Green Version]
  72. Mansoori, A.; Amini, M.; Kolahdooz, F.; Seyedrezazadeh, E. Obesity and Pro12Ala polymorphism of peroxisome proliferator-activated receptor-gamma gene in healthy adults: A systematic review and meta-analysis. Ann. Nutr. Metab. 2015, 67, 104–118. [Google Scholar] [CrossRef]
  73. Luo, W.; Guo, Z.; Wu, M.; Hao, C.; Hu, X.; Zhou, Z.; Zhou, Z.; Yao, X.; Zhang, L.; Liu, J. Association of Peroxisome Proliferator-Activated Receptor α/δ/γ with Obesity, and Gene-Gene Interaction, in the Chinese Han Population. J. Epidemiol. 2013, 23, 187–194. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Interaction between PPAR-δ mutated allele and sex: effect in HDL levels. *β = -0.13±0.06, p = 0.029. (B) Interaction between PPAR-δ mutated allele and sex: effect in LDL levels. *β = 9.71±4.28, p = 0.024. Genes 11 00776 i001 female, Genes 11 00776 i002 male.
Figure 1. (A) Interaction between PPAR-δ mutated allele and sex: effect in HDL levels. *β = -0.13±0.06, p = 0.029. (B) Interaction between PPAR-δ mutated allele and sex: effect in LDL levels. *β = 9.71±4.28, p = 0.024. Genes 11 00776 i001 female, Genes 11 00776 i002 male.
Genes 11 00776 g001
Table 1. General characteristics of studied populations.
Table 1. General characteristics of studied populations.
CharacteristicMestizo-1 n = 96 (23%)Mestizo-2 n = 37 (9%)Tarahumara n = 173 (42%)Mennonite n = 106 (26%)Total (n = 412)
Md (IQR) Md (IQR)Md (IQR) Md (IQR)Md (IQR)
Age (years)13 (13–14) a13 (12–13) a14 (13–15) b15 (14–15) c14 (13–15)
Weight (kg)47 (43–56) a59 (46–69) b,d48 (42–55) a57 (51–66) c,d51 (44–59)
Height (cm)158 (152–164) a157 (152–162) a152 (148–158) b171 (164–175) c158 (151–166)
BMI (kg/m2)19 (17–21) a24 (20–25) b20 (18–23) c,d20 (18–22) a,d20 (18–23)
n (%)n (%)n (%)n (%)n (%)
Sex
Male41 (43)a21 (57)a52 (30)b47 (44)a161 (39)
Female55 (57)16 (43)121 (70)59 (56)251 (61)
pBMI ≥ 85th (kg/m2)17 (18)a26 (70)b30 (17)a12 (11)a85 (21)
Ethnic group characteristics were compared by Kruskal-Wallis test, followed by post hoc Dunn’s test with Bonferroni adjustment a–d Different letters indicate significant differences between ethnic groups by pairwise comparisons. Abbreviations: Md = Median, IQR = Interquartile range, BMI (Body mass index, pBMI (Body mass index percentile).
Table 2. Metabolic traits comparison among the studied populations.
Table 2. Metabolic traits comparison among the studied populations.
CharacteristicMestizo-1Mestizo-2TarahumaraMennoniteTotal
Mean ± SD/Md (IQR)Mean ± SD/Md (IQR)Mean ± SD/Md (IQR)Mean ± SD/Md (IQR)Mean ± SD/Md (IQR)
BMI (Z-score)0.05 (−0.66–0.70) a1.45 (.75–1.64) b0.41 (-0.19–0.88) a,c-0.05 (-0.63–0.77) a,d0.25 (-0.50–0.92)
SBP (mmHg)110 ± 11 a115 ± 11 a108 ± 10 a108 ± 10 a112 ± 10
DBP (mmHg)70 (63–80) a70 (65–80) a,c70 (60–75) a75 (68–80) b,c70 (64–80)
WC (cm)64 (61–70) a80 (72–85) b,e70 (65–75) c68 (64–74) d,e68 (63–74)
WHR 0.41 (0.38–0.44) a0.52 (0.47–0.54) b0.46 (0.43–0.49) c0.40 (0.38–0.42) a0.43 (0.40–0.48)
Glucose (mg/dL)87 ± 10 a75 ± 6 b83 ± 11 a78 ± 9 c82 ± 10
TG (mg/dL)77 (57–109) a93 (72–129) a,c92 (75–115) b,c75 (60–92) a,d83 (64–111)
TC (mg/dL)145 ± 22 a156 ± 26 a133 ± 25 a144 ± 21 a140 ± 24
HDL (mg/dL)48 (40–55) a51 (44–56) a37 (32–42) b50 (43–56) a43 (36–52)
LDL (mg/dL)79 ± 17 a84 ± 19 a76 ± 20 a77 ± 16 a78 ± 18
VLDL (mg/dL)15 (11–22) a19 (14–26) a,c18 (15–23) b,c15 (12–18) a,d17 (13–22)
AI (Index)3.01 (2.59–3.48) a3.02 (2.55–3.74) a3.59 (3.12–4.20) b2.90 (2.56–3.19) a3.17 (2.76–3.69)
Ethnic group metabolic traits were compared by ANOVA or Kruskal-Wallis test for normally or non-normally distributed data, followed by post hoc Bonferroni or Dunn’s tests, respectively. a–e Different letters in the superscript indicate significant difference between ethnic groups by pairwise comparisons. Abbreviations: Md = Median, IQR = Interquartile range, SD = Standard deviation, BMI (body mass index), SBP (systolic blood pressure), DBP (diastolic blood pressure), WC (waist circumference), WHR (waist/height ratio), TC (total cholesterol), TG (triglycerides), VLDL (very density lipoproteins cholesterol), LDL (low density lipoproteins cholesterol), HDL (high low density lipoproteins cholesterol), and AI (atherogenic index).
Table 3. Pro12Ala PPAR-γ2 genotype and allelic frequencies comparison among populations.
Table 3. Pro12Ala PPAR-γ2 genotype and allelic frequencies comparison among populations.
PPAR-γ2Mestizo-1 n (%)Mestizo-2 n (%)Tarahumara n (%)Mennonite n (%)Total n (%)
CC79 (82) a25 (68) a,c103 (60) b,c83 (78) a290 (70)
CG14 (15)10 (27)61 (35)20 (19)105 (26)
GG3 (3)2 (5)9 (5)3 (3)17 (4)
HWE p0.090.7510.450.18
C172 (90) a60 (81) a,c267 (77) b,c186 (88) a685 (83)
G20 (10)14 (19)79 (23)26 (12)139 (17)
Chi2 test. a–c Different letters indicate significant difference in genotype and allelic frequencies among groups p < 0.05. HWE p (Hardy-Weinberg Equilibrium p). G = minor allele. The Tarahumara (b,c) group differed from Mestizo-1 (a) and Mennonite (a), but not from Mestizo-2 (a,c). Mestizo-1 (a), Mestizo-2 (a) and Mennonite (a) were not different. The difference between Mestizo-1 and Mestizo-2 showed a p value = 0.06.
Table 4. +294T/C PPAR-δ genotype and allelic frequencies comparison among populations.
Table 4. +294T/C PPAR-δ genotype and allelic frequencies comparison among populations.
PPAR-δMestizo-1 n (%)Mestizo-2 n (%)Tarahumara n (%)Mennonite n (%)Total n (%)
TT62 (65) a23 (62) a153 (88) b81 (76) a319 (77)
TC33 (34)13 (35)20 (12)23 (22)89 (22)
CC1 (1)1 (3)0 (0)2 (2)4 (1)
HWE p0.320.880.710.970.72
T157 (82) a59 (80) a326 (94) b185 (87) a727 (88)
C35 (18)15 (20)20 (6)27 (13)97 (12)
Chi2 test. a,b Different letters indicate significant differences in genotype and allelic frequencies among groups p < 0.05. HWE p (Hardy-Weinberg Equilibrium p). C = minor allele. The Tarahumara group differed from the other ethnic groups.
Table 5. Metabolic traits comparison between Pro12Ala PPAR-γ2 and +294T/C PPAR-δ genotype (major allele homozygous vs. minor allele carriers) in the total studied population.
Table 5. Metabolic traits comparison between Pro12Ala PPAR-γ2 and +294T/C PPAR-δ genotype (major allele homozygous vs. minor allele carriers) in the total studied population.
Metabolic TraitsPro12Ala PPAR-γ2+294T/C PPAR-δ
CC (n = 290)CG/GG (n = 122) TT (n = 319)TC/CC (n = 93)
Mean ± SD/Md (IQR)Mean ± SD/Md (IQR)pMean ± SD/Md (IQR)Mean ± SD/Md (IQR)p
BMI (Z-score)0.22 (−0.54–0.99)0.28 (-0.37–0.77)0.46 #0.21 (-0.52–0.87)0.39 (−0.41–1.18)0.06 #
SBP (mmHg)111 ± 11109 ± 100.08 *110 ± 11112 ± 110.15 *
DBP (mmHg)70 (65–80)70 (60–78)0.10 #70 (60–80)70 (65–78)0.78 #
WC (cm)68 (63–75)69 (65–74)0.37 #68 (63–74)69 (65–77)0.30 #
WHR0.42 (0.39–0.48)0.44 (0.41–0.48)0.04 #0.43 (0.40–0.48)0.42 (0.40–0.47)0.82 #
Glucose (mg/dL)82 ± 1082 ± 100.90 *82 ± 1180 ± 80.07 *
TG (mg/dL)83 (64–112)87 (66–110)0.51 #87 (68–112)75 (59–109)0.01 #
TC (mg/dL)142 ± 24138 ± 240.14 *140 ± 25143 ± 210.25 *
HDL (mg/dL)45 (37–52)39 (35–48)<0.01 #43 (36–51)48 (36–54)0.03 #
LDL (mg/dL)78 ± 1976 ± 170.37 *77 ± 1980 ± 160.24 *
VLDL (mg/dL)17 (13–22)17 (13–22)0.51 #17 (14–22)15 (12–22)0.01 #
AI (Index)3.14 (2.70–3.66)3.23 (2.82–3.77)0.15 #3.18 (2.81–3.68)3.12 (2.66–3.70)0.37 #
Groups were compared by * T-Student test or # Wilcoxon’s rank sum test for normally or non-normally distributed data, respectively. The association between PPAR-γ2 polymorphism and HDL remained significant after adjusting for multiple testing by the Holm-Bonferroni method (cut off p value = 0.0125) but the one of WHR did not (cut off p value = 0.006). The association between PPAR-δ and HDL did not remain significant after correcting for multiple testing; while significance remained for TG and VLDL. Abbreviations: Md = Median, IQR = Interquartile range, SD = Standard deviation, BMI (body mass index), SBP (systolic blood pressure), DBP (diastolic blood pressure), WC (waist circumference), WHR (waist/height ratio), TC (total cholesterol), TG (triglycerides), VLDL (very density lipoproteins cholesterol), LDL (low density lipoproteins cholesterol), HDL (high low density lipoproteins cholesterol), AI (atherogenic index), pBP (blood pressure percentile), and pWC (waist circumference percentile).
Table 6. Association/Interaction between Pro12Ala PPAR-γ2 polymorphism and metabolic traits by population groups.
Table 6. Association/Interaction between Pro12Ala PPAR-γ2 polymorphism and metabolic traits by population groups.
Metabolic TraitMestizo-1Mestizo-2TarahumaraMennoniteTotalTotal (Adjusted)
n = (96)n = (37)n = (173)n = (106)n = (412)n = (412)
β ± SEβ ± SEβ ± SEβ ± SEβ ± SEβ ± SE
(p)(p)(p)(p)(p)(p)
SBP (mmHg) *−5.59 ± 2.60−0.62 ± 3.03−1.37 ± 1.440.61 ± 2.64−2.05 ± 1.07−1.41 ± 1.08 (0.191) *, ++
(0.035)(0.840)(0.341)(0.817)(0.057) +
DBP (mmHg) *−2.14 ± 2.34−0.70 ± 3.36−0.34 ± 1.371.25 ± 2.10−1.68 ± 0.98−0.73 ± 0.98 (0.455)
(0.363)(0.835)(0.802)(0.554)(0.086)
Glucose (mg/dL) *−4.77 ± 2.711.74 ± 2.200.73 ± 1.531.46 ± 1.81−18.77 ± 9.20−13.82 ± 8.48 (0.104) *, --
(0.082)(0.436)(0.632)(0.423)(0.042)-
TG (log mg/dL) *−0.02 ± 0.11−0.07 ± 0.14−0.03 ± 0.050.10 ± 0.090.03 ± 0.04-0.008 ± 0.04 (0.84) §
(0.863)(0.654)(0.586)(0.276)(0.475)
TC (mg/dL) #−4.07 ± 5.85−10.67 ± 9.27−1.66 ± 3.748.01 ± 4.57−3.91 ± 2.59−1.48 ± 2.53 (0.558)
(0.489)(0.258)(0.658)(0.083)(0.132)
HDL (log mg/dL) ^−0.11 ± 0.05<−0.01 ± 0.080.03 ± 0.030.01 ± 0.04−0.06 ± 0.03−0.005 ± 0.02 (0.825)
(0.033)(0.961)(0.250)(0.815)(0.031) ~
LDL (mg/dL) *2.07 ± 4.71−9.94 ± 6.74−2.77 ± 2.874.98 ± 3.66−1.94 ± 1.92−9.63 ± 6.19 (0.119) ^
(0.660)(0.150)(0.336)(0.176)(0.312)
VLDL (log mg/dL) *−0.02 ± 0.11−0.07 ± 0.14−0.03 ± 0.050.10 ± 0.090.03 ± 0.04−0.28 ± 0.78 (0.720) §
(0.863)(0.654)(0.586)(0.276)(0.475)
AI (log) ^0.08 ± 0.05−0.07 ± 0.07−0.05 ± 0.030.05 ± 0.040.40 ± 0.190.39 ± 0.17 (0.021) ^, //
(0.114)(0.367)(0.084)(0.244)(0.033) /
Analyses on total studied population were adjusted for waist circumference and/or sex, and/or age, as indicated in variable names and ethnic groups. Total studied population adjusted models included ethnic groups as fixed effects, in addition to those indicated by the superscript. Variables with non-normal distribution were log transformed (HDL, TG, VLDL and AI). Coefficients for interaction terms are described below and indicated by symbols in the table. * Adjusted for waist circumference, sex and age. # adjusted for sex and age. ^ adjusted for waist circumference and sex. adjusted for waist circumference. § adjusted for waist circumference and age. adjusted for sex. + Interaction sex X age (β = −2.10 ± 0.72, p = 0.004), ++ Interaction sex X age (β = 1.84 ± 0.72, p = 0.011). - interaction polymorphism X waist circumference (β = 0.27 ± 0.13, p = 0.040), -- interaction polymorphism X waist circumference (β = 0.20 ±0.12, p = 0.10). ~ interaction sex X waist circumference (β = 0.007 ± 0.003, p = 0.008). / interaction polymorphism X waist circumference (β = −0.005 ± 0.003, p = 0.047), // interaction polymorphism X waist circumference (β = −0.006 ± 0.002, p = 0.019). Abbreviations: SBP (systolic blood pressure), DBP (diastolic blood pressure), TC (total cholesterol), TG (triglycerides), HDL (high density lipoproteins cholesterol), LDL (low density lipoproteins cholesterol), VLDL (very low-density lipoproteins cholesterol), and AI (atherogenic index).
Table 7. Association between +294T/C PPAR-δ polymorphism and metabolic traits in studied populations.
Table 7. Association between +294T/C PPAR-δ polymorphism and metabolic traits in studied populations.
Metabolic TraitMestizo-1Mestizo-2TarahumaraMennoniteTotalTotal (Adjusted)
n = (96)n = (37)n = (173)n = (106)n = (412)n = (412)
β ± SEβ ± SEβ ± SEβ ± SEβ ± SEβ ± SE
(p)(p)(p)(p)(p)(p)
SBP (mmHg) *2.07 ± 2.090.74 ± 2.971.71 ± 2.26−1.54 ± 2.551.94 ± 1.191.34 ± 1.20
(0.326)(0.806)(0.451)(0.546)(0.103)(0.263) *
DBP (mmHg) #−0.25 ± 1.84−0.38 ± 3.250.29 ± 2.10−0.67 ± 2.030.85 ± 1.080.01 ± 1.08
(0.894)(0.907)(0.889)(0.742)(0.432)(0.994) #
Glucose (mg/dL) *−2.62 ± 2.150.04 ± 2.17−3.91 ± 2.39−0.75 ± 1.75−2.57 ± 1.18−2.22 ± 1.12
(0.048) *
(0.227)(0.984)(0.103)(0.669)(0.031)
TG (log mg/dL) *−0.10 ± 0.09−0.12 ± 0.14−0.07 ± 0.08−0.15 ± 0.09−0.12 ± 0.05 −0.08 ± 0.05
(0.08) #
(0.273)(0.415)(0.374)(0.086)(0.008)
TC (mg/dL) ^−0.002 ± 4.602.12 ± 9.35−4.89 ± 5.87−5.76 ± 4.4879.02 ± 29.6149.8 ± 28.9
(0.085) &, ++
(1)(0.822)(0.406)(0.202)(0.008) +
HDL (log mg/dL) *0.02 ± 0.040.11 ± 0.07−0.10 ± 0.050.01 ± 0.04<0.01 ± 0.040.08 ± 0.07
(0.256)
(0.594)(0.144)(0.046)(0.685)(0.904) -
LDL (mg/dL) *−6.75 ± 4.642.66 ± 6.60.47 ± 4.52−4.15 ± 3.49−1.67 ± 2.66−2.89 ± 2.68
(0.281) ‡,~~
(0.149) /(0.688)(0.917) ^(0.238) &(0.531) ‡,~
VLDL (log mg/dL) *0.10 ± 0.09−0.12 ± 0.14−0.07 ± 0.08−0.15 ± 0.09−0.12 ± 0.05−0.08 ± 0.05
(0.082) #
(0.273)(0.415)(0.374)(0.086)(0.008)
AI (log) *−0.02 ± 0.04−0.08 ± 0.070.08 ± 0.05−0.06 ± 0.04−0.03 ± 0.02−0.03 ± 0.06
(0.599) #
(0.643)(0.223)(0.093)(0.133)(0.208)
Analyses on total studied population were adjusted for waist circumference and/or sex, and/or age, as indicated in variable names and ethnic groups. Total studied population. adjusted models included ethnic groups as fixed effects, in addition to those indicated by the superscript. Variables with non-normal distribution were log transformed (HDL, TG, VLDL and AI). Coefficients for interaction terms are described below and indicated by symbols in the table. * Adjusted for waist circumference, sex, and age. # adjusted for waist circumference and age. ^ adjusted for sex and age. & Adjusted for age. Adjusted for waist circumference. + Interaction polymorphism X age (β = −5.53 ± 2.14, p = 0.010), ++ Interaction polymorphism X age (β = −3.69 ± 2.08, p = 0.076). - interaction polymorphism X sex (β = −0.13 ± 0.06, p = 0.029). ~ Interaction polymorphism X sex (β = −9.71 ± 4.28 p = 0.024), ~~ Interaction polymorphism X sex (β = −8.99 ± 4.27 p = 0.035). / interaction polymorphism X and sex (β = −18.61 ± 7.27 p = 0.012). Significant differences in bold p < 0.05. Abbreviations: SBP (systolic blood pressure), DBP (diastolic blood pressure), TC (total cholesterol), TG (triglycerides), HDL (high density lipoproteins cholesterol), LDL (low density lipoproteins cholesterol), VLDL (very low density lipoproteins cholesterol), and AI (atherogenic index).
Table 8. Association between Pro12Ala PPAR-γ2 or +294T/C PPAR-δ with metabolic traits in the total studied population.
Table 8. Association between Pro12Ala PPAR-γ2 or +294T/C PPAR-δ with metabolic traits in the total studied population.
Metabolic TraitsPro12Ala PPAR-γ2 +294T/C PPAR-δ
CC CG/GGAdjusted OR (95%CI)pTTC/CCAdjusted OR (95%CI)p
(n = 290)(n = 122)(n = 319)(n = 93)
n (%)n (%)n (%)n (%)
BMI (Z-score > 1.04)68 (23)15 (12)0.46 (0.25–0.84)0.01157 (18)26 (28)1.78 (1.04–3.05)0.034
pBP ≥ 90th100 (34)33 (27)0.70 (0.44–1.12)0.142103 (32)30 (32)1 (0.61–1.64)0.996
pWC ≥ 90th9 (3)5 (4)1.33 (0.44–4.07)0.6129 (3)5 (5)1.96 (0.64–5.99)0.239
WHR > 0.45 *108 (37)55 (45)1.79 (1.03–3.08)0.036 131 (41)32 (34)0.74 (0.45–1.21)0.231
Glucose ≥ 110 (mg/dL) **4 (1)1 (1)-- 5 (2)0 (0)--
TG ≥ 110 ^ (mg/dL)76 (26)30 (25)0.91 (0.56–1.50)0.72184 (26)22 (24)0.77 (0.44–1.34)0.350
TC ≥ 150 # (mg/dL)41 (14)9 (7)0.47 (0.22–1)0.050 -39 (12)11 (12)0.87 (0.43–1.81)0.726 ~
HDL ≤ 40 @ (mg/dL)92 (32)65 (53)2.50 (1.60–3.89)<0.01128 (40)29 (31)1 (0.55–1.82)0.995 /
LDL ≥ 110 (mg/dL)13 (4)3 (2)0.54 (0.15–1.96)0.35115 (5)1 (1)0.20 (0.03–1.52)0.118
VLDL ≥ 30 (mg/dL)22 (8)7 (6)0.74 (0.31–1.80)0.51423 (7)6 (6)0.82 (0.32–2.10)0.679
AI (M > 4 F > 3.5) &86 (30)39 (32)1.07 (0.64–1.78)0.79899 (31)26 (28)0.80 (0.45–1.41)0.435
Reference: homozygous group for major allele. * Adjusted by sex,# adjusted by sex and age, ^ adjusted by waist circumference and age, @ adjusted by waist circumference, sex and age, adjusted by waist circumference, & adjusted by sex and WC. Interaction polymorphism X sex (β = 0.46 ± 0.22, p = 0.108), - interaction sex X age (β = 0.39 ± 0.14, p = 0.009), ~ interaction sex X age (β = 0.40 ± 0.14, p = 0.010), / interaction polymorphism X sex (β = 0.27 ± 0.17, p = 0.040). ** Odds ratio was not calculated, due to the low frequency of individuals in the high-glucose category. Abbreviations: pBP (blood pressure percentile), pWC (waist circumference percentile), WHR (waist-height ratio), TC (total cholesterol), TG (triglycerides), HDL (high low density lipoproteins cholesterol), LDL (low density lipoproteins cholesterol), VLDL (very low density lipoproteins cholesterol), and AI (atherogenic index) M = male, F = female.

Share and Cite

MDPI and ACS Style

Carrillo-Venzor, M.A.; Erives-Anchondo, N.R.; Moreno-González, J.G.; Moreno-Brito, V.; Licón-Trillo, A.; González-Rodríguez, E.; Hernández-Rodríguez, P.d.C.; Reza-López, S.A.; Loera-Castañeda, V.; Leal-Berumen, I. Pro12Ala PPAR-γ2 and +294T/C PPAR-δ Polymorphisms and Association with Metabolic Traits in Teenagers from Northern Mexico. Genes 2020, 11, 776. https://doi.org/10.3390/genes11070776

AMA Style

Carrillo-Venzor MA, Erives-Anchondo NR, Moreno-González JG, Moreno-Brito V, Licón-Trillo A, González-Rodríguez E, Hernández-Rodríguez PdC, Reza-López SA, Loera-Castañeda V, Leal-Berumen I. Pro12Ala PPAR-γ2 and +294T/C PPAR-δ Polymorphisms and Association with Metabolic Traits in Teenagers from Northern Mexico. Genes. 2020; 11(7):776. https://doi.org/10.3390/genes11070776

Chicago/Turabian Style

Carrillo-Venzor, Martín A., Nancy R. Erives-Anchondo, Janette G. Moreno-González, Verónica Moreno-Brito, Angel Licón-Trillo, Everardo González-Rodríguez, Pilar del Carmen Hernández-Rodríguez, Sandra A. Reza-López, Verónica Loera-Castañeda, and Irene Leal-Berumen. 2020. "Pro12Ala PPAR-γ2 and +294T/C PPAR-δ Polymorphisms and Association with Metabolic Traits in Teenagers from Northern Mexico" Genes 11, no. 7: 776. https://doi.org/10.3390/genes11070776

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop