Next Article in Journal
Improvement in Clinical Features of L-NAME-Induced Preeclampsia-like Rats through Reduced SERPINA5 Expression
Next Article in Special Issue
PPARγ Modulators in Lung Cancer: Molecular Mechanisms, Clinical Prospects, and Challenges
Previous Article in Journal
Epileptic Encephalopathy GABRB Structural Variants Share Common Gating and Trafficking Defects
Previous Article in Special Issue
Current Clinical Trial Status and Future Prospects of PPAR-Targeted Drugs for Treating Nonalcoholic Fatty Liver Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolic Biomarkers in Adults with Type 2 Diabetes: The Role of PPAR-γ2 and PPAR-β/δ Polymorphisms

by
Sandra A. Reza-López
1,
Susana González-Gurrola
2,
Oscar O. Morales-Morales
1,
Janette G. Moreno-González
1,
Ana M. Rivas-Gómez
2,
Everardo González-Rodríguez
1,
Verónica Moreno-Brito
1,
Angel Licón-Trillo
1 and
Irene Leal-Berumen
1,*
1
Facultad de Medicina y Ciencias Biomédicas, Universidad Autónoma de Chihuahua, Campus II. Circuito Universitario S/N, Chihuahua 31109, CP, Mexico
2
Instituto Mexicano del Seguro Social UMF 33, Avenida Melchor Ocampo y Arroyo de los Perros S/N, Col. El Palomar, Chihuahua 31204, CP, Mexico
*
Author to whom correspondence should be addressed.
Biomolecules 2023, 13(12), 1791; https://doi.org/10.3390/biom13121791
Submission received: 9 November 2023 / Revised: 27 November 2023 / Accepted: 2 December 2023 / Published: 14 December 2023
(This article belongs to the Special Issue PPARs as Key Regulators in Different Diseases)

Abstract

:
Glucose and lipid metabolism regulation by the peroxisome proliferator-activated receptors (PPARs) has been extensively reported. However, the role of their polymorphisms remains unclear. Objective: To determine the relation between PPAR-γ2 rs1801282 (Pro12Ala) and PPAR-β/δ rs2016520 (+294T/C) polymorphisms and metabolic biomarkers in adults with type 2 diabetes (T2D). Materials and Methods: We included 314 patients with T2D. Information on anthropometric, fasting plasma glucose (FPG), HbA1c and lipid profile measurements was taken from clinical records. Genomic DNA was obtained from peripheral blood. End-point PCR was used for PPAR-γ2 rs1801282, while for PPAR-β/δ rs2016520 the PCR product was digested with Bsl-I enzyme. Data were compared with parametric or non-parametric tests. Multivariate models were used to adjust for covariates and interaction effects. Results: minor allele frequency was 12.42% for PPAR-γ2 rs1801282-G and 13.85% for PPAR-β/δ rs2016520-C. Both polymorphisms were related to waist circumference; they showed independent effects on HbA1c, while they interacted for FPG; carriers of both PPAR minor alleles had the highest values. Interactions between FPG and polymorphisms were identified in their relation to triglyceride level. Conclusions: PPAR-γ2 rs1801282 and PPAR-β/δ rs2016520 polymorphisms are associated with anthropometric, glucose, and lipid metabolism biomarkers in T2D patients. Further research is required on the molecular mechanisms involved.

Graphical Abstract

1. Introduction

Diabetes is a serious public health issue with a 10.5% worldwide prevalence [1]. In Mexico, it reaches about 18.3% in the adult population [2]. Type 2 diabetes (T2D) is the most common form of diabetes, and accounts for over 90% of cases [1]. Numerous and diverse environmental and genetic factors are involved in T2D development. Medical history, serum biomarkers, psychosocial and lifestyle factors, unhealthy diet, adiposity, family history, age, ethnicity, and genetics have all been related to T2D and glycemic dysregulation [1,3,4,5].
Diabetes can lead to complications such as cardiovascular disease, nephropathy, or retinopathy [6], which have in turn been related to hyperglycemia, hypertension, and hypercholesterolemia. Decreasing HbA1c may reduce the risk of microvascular complications by 40% [7]. Therefore, glycemic control is a common treatment target for T2D patients. However, several reports indicate that this goal is achieved by less than 50% of patients. Lifestyle variables, medication, and individual characteristics [8,9], as well as genetics have been related to glycemic control. Among the latter, several gene polymorphisms have been associated with the response to medications [10]. However, few studies have reported the relation of polymorphisms and glycemic control by other pathways [11,12]. For instance, poor glycemic control has been observed in carriers of the T allele of rs2241766 SNP in ADIPOQ, a gene that encodes for adiponectin, which is an adipokine related to insulin sensitivity [11].
Other genes that play important roles in glucose and lipid metabolic networks are those encoding for the peroxisome-proliferator-activated-receptors (PPAR) [13,14,15]. There are three main PPAR nuclear receptors, each isoform with different tissue distribution and activated by specific ligands with wide physiological functions [16,17]. PPAR-β/δ has been shown to promote glycolysis, glucose absorption, glycogen storage, and gluconeogenesis downregulation in animal models [18,19]; similar to PPAR-α, it is related to fatty acid oxidation, increase serum high-density lipoproteins (HDL) levels with lower plasma triglycerides (TG), and very-low-density lipoproteins (VLDL) [15]. PPAR-γ has three isoforms expressed in different tissues. PPAR-γ contributes to lipid synthesis and storage in white adipose tissue, and plays a role in insulin sensitivity and glucose homeostasis [20,21].
Several polymorphisms of PPARG and PPARD genes have been described. PPAR-γ2 rs1801282 is a missense single-nucleotide-polymorphism (SNP) in the B exon on chromosome 3p25.3 (OMIM number 601487), where a cytosine/guanine (C/G) transversion promotes a polypeptide substitution of proline to alanine in the 12th amino acid (Pro12Ala) [22]. This variant has been related to anthropometric measures, insulin resistance and serum lipids levels [23,24]. PPAR-β/δ rs2016520 (+294T/C) is a transition SNP located in the 5’UTR exon 4 region on chromosome 6p21.31 (OMIM number 600409). While it is not associated with T2D risk [25], studies in different populations have reported the association of this PPAR-β/δ SNP with glucose and lipid biomarkers [26,27]; however, studies in subjects with T2D have shown diverse results [25,28]. In this study, we aimed to determine the frequencies of PPAR-γ2 rs1801282 (Pro12Ala) and PPAR-β/δ rs2016520 (+294T/C) polymorphisms and their relation to anthropometric and metabolic biomarkers in adults with T2D from Chihuahua, Northern Mexico.

2. Materials and Methods

2.1. Study Design and Population

This cross-sectional study included 314 patients diagnosed with T2D by their treating physicians at the Unidad de Medicina Familiar, Clínica 33 of the Instituto Mexicano del Seguro Social (UMF33-IMSS) Chihuahua, Chihuahua, Mexico. Patients were older than 18 years of age, affiliated with UMF33-IMSS Chihuahua, they had a clinical diagnosis for T2D in their file, and were willing to donate a blood sample. Patients who had interest in participating and who signed a written informed consent were included in the study. The study was approved by the institutional ethics research committees with the number R-2020-805-028 and CI-020-19 for the IMSS and the Facultad de Medicina y Ciencias Biomédicas, Universidad Autónoma de Chihuahua, respectively.

2.2. Phenotyping and Biochemical Data

Patients’ anthropometric and biochemical data were obtained from their medical records. Anthropometric measures were: height (m), weight (kg), and waist circumference (cm). Body mass index (BMI) was calculated as the weight in kilograms divided by the square of the height in meters (kg/m2) and classified as follows: ≤24.9 as normal weight, 25–29.9 overweight, 30–34.9 obesity class I, 35–39.9 obesity class II, and ≥40 obesity class III. Glycemic control was considered with a glycosylated hemoglobin A1c (HbA1c) < 7% according to the Mexican guidelines [29]. The TC/HDL-C index was calculated [30].

2.3. Genotyping

The genomic DNA was isolated from peripheral blood buffy coats using a Master Pure Epicenter DNA extraction Kit (Madison, Wisconsin, USA) following the manufacturer’s recommendations. Genotyping of PPAR-γ2 rs1801282 (Pro12Ala) gene and the PPAR-β/δ rs2016520 (+294T/C) were performed as previously described [27].
The polymerase chain reaction (PCR) conditions used for both SNPs were: an initial 3 min denaturation 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 consisted of 5 min at 72 °C (Agilent SureCycler 8800, Santa Clara, CA, USA).
Electrophoresis in 3.5% agarose gels was used to identify PPAR-γ2 rs1801282 genotypes. The PPAR-β/δ rs2016520 -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 resulting fragments and the primers used for genotyping are presented in Table 1.

2.4. Statistical Analysis

The Χ2 test verified Hardy–Weinberg equilibrium (H-WE). Quantitative measures were described as means and standard deviation (SD) or median and interquartile range (IQR) for variables with normal or non-normal distribution, respectively. Absolute and relative frequencies are presented for categorical variables. Three inheritance models were analyzed: codominant, dominant, and over-dominant. Anthropometric and biochemical measures were compared by inheritance model with Student’s t-test or Wilcoxon’s rank sum test (Mann–Whitney U test); analysis of variance or the Kruskal–Wallis test, followed by Bonferroni’s post hoc test or Dunn’s test, as appropriate. Categorical variables were compared by Pearson’s Χ2 test. A p < 0.05 value was considered significant.
Multiple linear regression models were used to adjust the relation between biochemical measures as continuous variables (log-transformed when needed to reach normality) and PPAR for variables such as age, sex, age at T2D diagnosis, BMI, and waist circumference as potential confounders and to test interaction effects. Variables with p < 0.05 were kept in multivariate analyses. A stratified analysis was conducted when an interaction effect was identified (p ≤ 0.10). The normality of residuals, homoscedasticity, and variance inflation factors were tested. Logistic regression was used to analyze the relation between HbA1c as categorical variable (<7% vs. ≥7%) and PPAR SNPs, adjusting for potential confounders as covariates. Stata v. 14.2 (2015, Stata Corp., College Station, TX, USA) was used for statistical analysis.

3. Results

3.1. Population Characteristics

A total of 314 participants with the following characteristics were included: age of 57.9 ± 10.6 years (mean ± SD); age at the time of T2D diagnosis of 48.9 ± 9.9 years and time from T2D diagnosis of 9 (4–13) years (median and IQR); height of 1.62 ± 0.09 m and weight 81 (72–93). According to BMI, 11.2% (n = 35) showed normal weight, 31.2% (n = 98) were overweight, and 57.6% (n = 181) had obesity (33.1% grade I, 15.9% grade II, and 8.6% grade III).
Almost 60% (59.9%) of T2D patients had uncontrolled glycaemia, defined as HbA1c ≥ 7%: 70% in those who were diagnosed ≤45 years old, 62% in the age group 45–53, and 47% in patients older than 53 years. Ninety percent of participants’ treatment included biguanides, 33% were treated with insulin, and 32% received dipeptidyl peptidase 4 inhibitors (DPP4); moreover, 19% received sulfonylureas, 7% thiazolidinediones, 6% sodium–glucose cotransporter inhibitors, 3.5% alpha-glucosidase inhibitors, and 1% glucagon-like peptide I (GLP-I) agonists. Thirty-four percent had one glucose-lowering drug prescribed, 44% had two of them, and 22% were using three or more. For lipid control, 73% of participants had one medication; 49% were taking statins, 26% used fibrates, and less than 1% used sterol absorption inhibitors.

3.2. PPAR-γ2 rs18 01282 and PPAR-β/δ rs2016520

The genotype frequencies for both PPAR polymorphisms showed H-WE equilibrium (Table 2).
Anthropometric and biochemical measures are presented in Table 3 by PPAR-γ2 rs1801282 and PPAR-β/δ rs2016520 (dominant model). Carriers of at least one minor allele of PPAR-γ2 rs1801282-G had higher values of waist circumference (p = 0.029), FPG (p = 0.026), and marginally higher HbA1c (p = 0.057), while carriers of at least one minor allele of PPAR-β/δ rs2016520-C had higher TG (p = 0.026) and VLDL (p = 0.026) and tended to have increased waist circumference (p = 0.052), FPG (p = 0.080) and HbA1c values (p = 0.089). According to genotype (codominant model), FPG was higher in rs1801282-CG compared to the CC homozygous for the major allele (p = 0.040); TG was greater in rs2016520-TC compared to the homozygous major allele genotype (p = 0.046). Similarly, rs2016520-TC tended to higher HbA1c values (p = 0.063), and a trend was observed in waist circumference for the heterozygous of both PPARs (p = 0.072 and p = 0.103 for the rs1801282-CG and rs2016520-TC, respectively) [see Supplementary Table S1].

3.3. PPAR SNPs and Anthropometric Measures

The relation between waist circumference and the studied SNPs were statistically significant after adjusting for the T2D duration (current age–age at the time of diagnosis). Having at least one minor allele of PPAR-γ2 rs1801282-CG + GG was related to higher value of waist circumference (log cm, β = 0.04, p = 0.026), while in those with at least one minor allele of PPAR-β/δ rs2016520-TC + CC the relation with waist circumference was significant (log cm, β = 0.04, p = 0.029). BMI (log-transformed) was related to PPAR-β/δ rs2016520-TC + CC (β = 0.07, p = 0.015) after adjusting for HbA1c and age. PPAR-γ2 rs1801282-G was not significantly related to BMI (p = 0.802).

3.4. PPAR SNPs and Glycemic Biomarkers

For FPG as dependent variable, an interaction effect between both polymorphisms was observed (p = 0.013). Therefore, the geometric means after adjusting for age at the time of diagnosis are presented for rs2016520 separately (TT vs. TC + CC and TT + CC vs. TC) [Figure 1]. The subgroup carrying at least one minor allele of both PPAR SNPs had the highest FPG values after adjusting for age at the time of T2D diagnosis. Similar results were observed in the over-dominant inheritance model. Including the glucose-lowering drugs in the multivariate model did not modify these results. Among patients with insulin treatment (alone or in combination), those with at least a minor allele of the PPAR-γ2 Pro12Ala had higher FPG than the major allele homozygous (p for interaction = 0.063).
Both PPAR SNPs were independently associated with HbA1c values (log transformed) in multivariate analyses after adjusting for BMI category, age at the time of T2D diagnosis, and HDL-C. Carriers of at least one minor allele of PPAR-γ2 rs1801282 had higher values of HbA1c. Similarly, carriers of at least one minor allele of PPAR-β/δ rs2016520 had increased HbA1c (Table 4). Treatment with sulfonylureas, insulin, or IDPP4 was related to HbA1c values (p < 0.05); however, including pharmacological treatment did not change the direction of the results.
HbA1c as categorical variable (<7 vs. ≥7%) was associated with the same variables by logistic regression analysis (PPAR-γ2 rs1801282 CC vs. CG + GG: OR =2.14, 95%CI 1.01, 4.58, p = 0.048); PPAR-β/δ rs2016520 TT vs. TC + CC: OR =2.08, 95%CI 1.03, 4.21, p = 0.041), after adjusting for the same variables as the linear regression model. Among glucose lowering drugs, insulin, IDPP4, and sulfonylureas were associated with uncontrolled glycaemia (OR =5.46, OR =4.40, and OR =2.95, p < 0.05, respectively) in the fully adjusted model. No interaction effects between glucose-lowering medications and PPAR SNPs were detected. Similar results were observed in the over-dominant model.

3.5. PPAR SNPs and Lipid Biomarkers

For TG, we observed an interaction effect between FPG and PPAR-γ2 and PPAR-β/δ polymorphisms; thus, we present them separately (Table 5). FPG was directly related to TG levels in all groups except those carrying both major alleles of the PPAR-β/δ rs2016520-TT and at least one minor allele of PPAR-γ2 rs1801282-CG + GG. For those carrying at least one minor allele of both SNPs, FPG and BMI explained more than 60% of the variability in TG levels (Adj R2 = 0.61).

4. Discussion

This study analyzed the frequencies of PPAR-γ2 rs1801282 (Pro12Ala) and PPAR-β/δ rs2016520 (+294T/C) polymorphisms and their association with anthropometric and metabolic biomarkers in adults with T2D from a public health institution in Northern Mexico. An interaction effect between these polymorphisms resulted in increased values of FPG in minor allele carriers. Both SNPs were independently associated with increased HbA1c and modified the relation between FPG and TG values.
The prevalence of PPAR-γ2 rs1801282 (Pro12Ala) polymorphism varies across populations. According to genotype frequency, several studies have reported 0% homozygous GG in Asian Indians [31], original South African [32], while in Romanian children and Caucasian adults from Poland the frequency has ranged between 2.7 and 3.7% [33,34]. In teenagers from our region’s main ethnic groups, it ranges from 3 to 5%, [27]. In Mexican adults with diabetes, the frequency of the G allele was identified as 13.41% [35], similar to that observed in this study (12.42%). PPAR-β/δ rs2016520 (+294T/C) frequency is variable among populations as well. Reported frequencies range from 28.8 to 44.3% of heterozygous (TC) and 3.09 to 12.9% of CC homozygous [36,37,38] in Russian, Chinese, and Mexican populations. Consistent with previous studies in Mexican teenagers from Northern Mexico [27] reporting 1% of CC and 22% of TC, in this study we identified 1.27% of CC homozygous and 25.16% heterozygous in adults with T2D.
Gene polymorphisms involved in glucose and lipid regulation influence glycemic control in T2D. A recent systematic review revealed a high prevalence of poor glycemic control, ranging between 45.2% to 93% in different countries [39]. In Mexico, glycemic control prevalence was 39% (i.e., uncontrolled 61%) [40]. Similarly, the prevalence of patients with uncontrolled glycaemia found in the present study was 59.9%. In addition, we observed that carriers of minor alleles of the studied PPAR SNPs were independently associated with higher HbA1c values and increased odds for uncontrolled glycaemia (adjOR = 2.14 and adjOR = 2.08 for PPAR-γ2 rs1801282 and PPAR-β/δ rs2016520, respectively; p < 0.05). In multivariate models, HDL-C levels were inversely associated with glycemic control, consistent with other reports in adults with diabetes from the National Health and Nutrition Examination Survey of the United States [41]. In addition, receiving pharmacological treatments such as insulin, IDPP4, and sulfonylureas was associated with poor control. Similarly, lower HDL-C, longer duration of T2D, and treatment with insulin either alone or in combination with oral hypoglycemic agents have all been associated with poor glycemic control in other studies [8,42]. The unexpected associations may be partially explained by the prescription of insulin and other glucose-lowering treatments being recommended to individuals with uncontrolled glycaemia. According to clinical guidelines, dual therapy would be indicated for those patients with uncontrolled glycaemia receiving monotherapy (usually metformin) and triple or insulin schemes when therapeutic goals are not reached [43]. Because we examined prevalent cases with different disease duration, we speculate that cases with poor control are prescribed with more medication. In this way, patients’ glycemic control status could have influenced the number and type of prescribed medications.
HbA1c is the main biochemical measure recommended to evaluate glycemic control, reflecting average glucose levels. However, other metabolic biomarkers are assessed in T2D patients as well, such as fasting plasma glucose levels and lipid profile [44]. In our study, PPAR-γ2 rs1801282-CG + GG and PPAR-β/δ rs2016520-TC + CC were associated with FPG, HbA1c, and TG. Similar to our findings, in a Bangladeshi T2D population, carriers of PPAR γ2 rs1801282-CG + GG had higher values of FPG and HbA1c [45], while in a T2D study from a Han Chinese population, heterozygous tended to have higher baseline values of FPG, although the difference did not reach statistical significance. However, the decrease in postprandial glucose and TG in response to pioglitazone was greater than in major allele homozygous [46]. PPAR-γ2 rs1801282-CG + GG minor allele carriers had higher values of TG in subjects with obesity, whereas they were lower in T2D patients from Australian and Russian populations, respectively [47,48]. In contrast, FPG, cholesterol, and TG were lower in patients with the rs1801282-CG genotype compared to the homozygous CC in an Iranian population with coronary artery disease (CAD) and T2D [49]. Of note, the levels of these biomarkers were lower than those found in our study, and given the antecedent of CAD, probably received different pharmacological treatment.
There are few studies of PPAR-β/δ rs2016520 in T2D. In a study investigating the role of PPARD polymorphism in the response to exenatide treatment in newly diagnosed T2D patients from China, baseline values of BMI were higher in patients with TC and CC than in those with TT genotype, while TG values were higher in minor allele homozygous CC of the PPAR-β/δ rs2016520. In their study, fasting insulin and HOMA-IR in carriers of at least one C allele were lower after 6 months of exenatide treatment. However, in multivariate models the rs2016520 SNP was not significantly associated with HbA1c improvement [50]. In another study, TT major allele homozygous responded to nateglinide with a greater decrease in postprandial glucose and increase in HOMA-B than in subjects with the TC/CC genotype [28].
The molecular and physiological impacts of these polymorphisms on glucose metabolism are not fully understood. One of the most studied is the PPAR-γ2 Pro12Ala. Carrying the G allele is associated with higher BMI, waist circumference, and TC than in CC homozygous individuals [51]. The PPAR-γ2 rs1801282-G variant has shown decreased transactivation because of its lower affinity to the responsive element and improved insulin sensitivity [52]. By forming a heterodimer with the Retinoid X receptor, PPARγ may induce the adiponectin gene, binding to PPAR-responsive element in human adiponectin promoter (the PPAR-γ isoform was not specified) [53]. Thus, it is possible that PPAR-γ2 polymorphism could modulate adiponectin expression. In an animal study, Heikkinen et al. (2009) analyzed the phenotype and gene expression of homozygous Ala/Ala compared to Pro/Pro under a normal and high fat diet. Ala/Ala mice were leaner, had better insulin sensitivity, and lived longer than Pro/Pro mice. Gene expression analysis of white adipose tissue, muscle, and liver revealed an upregulation of adiponectin receptor 2 expression in adipose tissue and muscle in Ala/Ala mice and in adiponectin expression in Ala/Ala muscle of mice fed a high fat diet. According to the authors, these findings suggest that adiponectin signaling could be involved in the observed Ala/Ala mice characteristics [54]. Studies in Finnish servicemen showed that the Ala allele was associated with increased adiponectin levels in those with more than 10% weight loss [55]. These studies show the complexity of glucose metabolism regulation and the participation of environmental factors.
Less is known about the effect of the PPARD +294T/C polymorphism. PPAR-β/δ plays a role in fatty acid catabolism and energy homeostasis through cellular-level metabolic pathways [56]. Because PPAR-β/δ is an important regulator of lipid metabolism, most studies have focused on investigating lipid profile according to this SNP. The minor allele C had been shown to increase its transcriptional activity compared to the major allele T; in healthy men, carriers of the CC allele had higher LDL-C [57]. In our study in T2D patients, carriers of the C allele had higher TG values (p < 0.05) in the dominant inheritance model. In the same line, recent studies in a Chinese population with diabetes found higher TG levels in CC homozygous at the baseline of a study investigating the effects of exenatide. Regarding glycemic control, the same authors reported lower expression of PPAR-β/δ in liver tissue of db/db mice and in an insulin resistance model in HepG2 cells. Their results suggest a role of PPAR-β/δ activation in GLP-1R expression regulation, explaining the response to exenatide [50]. Thus, while the C allele might increase transcriptional activity, this could be modified in part by an insulin resistance or dysregulation in T2D patients.
Outcomes related to the studied PPAR SNPs vary across populations, by ethnic background, individual characteristics, and by conditions such as obesity or T2D. This can be illustrated by the results of a recent meta-analysis by Li et al. (2022) on PPARG rs1801282. Carriers of the G allele had a higher BMI, waist circumference, or differences in lipid profile in several Asian or African populations, while no association with obesity indexes or lipid profile was found in Caucasians, either European or American. The same authors reported interactions of this polymorphism with gender; whereas male G allele carriers had higher BMI, female carriers had higher values of waist circumference [51]. In the same line, some target genes do not show differences between heterozygous and major allele carriers in subjects with morbid obesity, suggesting an interaction with other characteristics [58]. Lifestyle factors contribute to the high variation observed in glycemic responses among individuals as well as interactions between SNPs. Factors such as diet [54], weight loss [55], exercise [59], and fatty acid levels [60] have been documented to influence or interact with SNPs. Furthermore, gene–gene interactions should be taken into account when assessing the influence of SNPs in metabolic biomarkers. Examining polymorphisms, particularly in patients with poor glycemic control, may be important in clinical practice, as genetic variants may influence the response to pharmacological treatment. Furthermore, the effects of genetic variants on metabolic pathways, independently of those for pharmacological treatment actions, remain to be investigated.
Taken together, our results suggest that PPAR polymorphisms are related to metabolic measures in T2D patients. We have to consider several limitations of this study: first, data on lifestyle variables such as nutrient intake and physical activity were not collected; thus, we were unable to assess their role in conjunction with PPAR polymorphisms. Having information on treatment adherence would also have been desirable. Adherence in T2D patients has been reported in about 47% of T2D patients from a similar population in North Mexico [61]. Although treatment adherence could have influenced metabolic biomarkers, we are not aware of differences by the studied SNP genotypes. Second, we included subjects with variable T2D duration, which may influence biomarker values; however, we adjusted models for either age at the time of diagnosis or duration of the disease. As this was a cross-sectional study, only one-time point was selected. Third, we had sample size limitations when fully evaluating biomarkers in carriers of both minor alleles or with haplotype combinations. For instance, heterozygous PPAR-β/δ rs2016520-TC had the highest values of glucose, HbA1c, and TG, whereas minor allele homozygous had the lowest. However, in the dominant model (TC + CC), the overall effect indicated an increase in such biomarkers because of the small number of subjects with the CC genotype (n = 4). Thus, the sample size for the minor allele homozygous did not allow us to analyze them separately. in addition, the sample size limited certain medication–gene interaction analyses for less commonly prescribed medications. Our study included patients from a public health institution with which about 50% of the population is affiliated [62]; thus, our results can be generalized to populations with similar characteristics.

5. Conclusions

In conclusion, PPAR-γ2 rs1801282 (Pro12Ala) and PPAR-β/δ rs2016520 (+294T/C) polymorphisms were associated with anthropometric, glucose and lipid metabolism biomarkers. While there was an interaction effect between both studied polymorphisms for FPG, they showed independent effects for HbA1c. Furthermore, these polymorphisms modified the relation between FPG and TG in adults with T2D. For the best of our knowledge, this is the first study to report a gene–gene interaction between PPAR-β/δ rs2016520 and PPAR-γ rs1801282 in T2D metabolic biomarkers. Further research is needed to gain insight about the mechanisms involved in SNP-related differences, their impact on overall health, and the interaction between lifestyle components and genetic background. The high prevalence of overweight and diabetes patients in Mexico highlights the importance of strategies for primary prevention as well as for preventing disease complications.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biom13121791/s1, Table S1: Anthropometric and biochemical measures by PPAR genotypes (codominant model); Table S2: Multivariate analysis of the relation between PPARs polymorphisms and fasting plasma glucose (log mg/dL).

Author Contributions

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

Funding

This study was partially supported by the Universidad Autónoma de Chihuahua (Facultad de Medicina y Ciencias Biomédicas).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved 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” with registration number CI-020-19 and R-2020-805-028 number for IMSS.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are publicly unavailable due to privacy and ethical restrictions. However, they could be made available through a request letter to the corresponding author previous to the ethics committee authorization.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Magliano, D.; Boyko, E.J. IDF Diabetes Atlas, 10th ed.; International Diabetes Federation: Brussels, Belgium, 2021; ISBN 978-2-930229-98-0. [Google Scholar]
  2. Basto-Abreu, A.; López-Olmedo, N.; Rojas-Martínez, R.; Aguilar-Salinas, C.A.; Moreno-Banda, G.L.; Carnalla, M.; Rivera, J.A.; Romero-Martinez, M.; Barquera, S.; Barrientos-Gutiérrez, T. Prevalencia de prediabetes y diabetes en México: Ensanut 2022. Salud Pública México 2023, 65, s163–s168. [Google Scholar] [CrossRef]
  3. Bellou, V.; Belbasis, L.; Tzoulaki, I.; Evangelou, E. Risk Factors for Type 2 Diabetes Mellitus: An Exposure-Wide Umbrella Review of Meta-Analyses. PLoS ONE 2018, 13, e0194127. [Google Scholar] [CrossRef]
  4. Goyal, S.; Rani, J.; Bhat, M.A.; Vanita, V. Genetics of Diabetes. World J. Diabetes 2023, 14, 656–679. [Google Scholar] [CrossRef]
  5. Liu, X.; Collister, J.A.; Clifton, L.; Hunter, D.J.; Littlejohns, T.J. Polygenic Risk of Prediabetes, Undiagnosed Diabetes, and Incident Type 2 Diabetes Stratified by Diabetes Risk Factors. J. Endocr. Soc. 2023, 7, bvad020. [Google Scholar] [CrossRef]
  6. Bancks, M.P.; Carnethon, M.; Chen, H.; Cotch, M.F.; Klein, B.; Klein, R.; Szklo, M.; Bertoni, A. Diabetes Subgroups and Risk for Complications: The Multi-Ethnic Study of Atherosclerosis (MESA). J. Diabetes Complicat. 2021, 35, 107915. [Google Scholar] [CrossRef]
  7. Deshpande, A.D.; Harris-Hayes, M.; Schootman, M. Epidemiology of Diabetes and Diabetes-Related Complications. Phys. Ther. 2008, 88, 1254–1264. [Google Scholar] [CrossRef]
  8. Haghighatpanah, M.; Nejad, A.S.M.; Haghighatpanah, M.; Thunga, G.; Mallayasamy, S. Factors That Correlate with Poor Glycemic Control in Type 2 Diabetes Mellitus Patients with Complications. Osong Public Health Res. Perspect. 2018, 9, 167–174. [Google Scholar] [CrossRef]
  9. Souliotis, K.; Koutsovasilis, A.; Vatheia, G.; Golna, C.; Nikolaidi, S.; Hatziagelaki, E.; Kotsa, K.; Koufakis, T.; Melidonis, A.; Papazafiropoulou, A.; et al. Profile and Factors Associated with Glycaemic Control of Patients with Type 2 Diabetes in Greece: Results from the Diabetes Registry. BMC Endocr. Disord. 2020, 20, 16. [Google Scholar] [CrossRef]
  10. Mofo Mato, E.P.; Guewo-Fokeng, M.; Essop, M.F.; Owira, P.M.O. Genetic Polymorphisms of Organic Cation Transporter 1 (OCT1) and Responses to Metformin Therapy in Individuals with Type 2 Diabetes: A Systematic Review. Medicine 2018, 97, e11349. [Google Scholar] [CrossRef]
  11. Alfaqih, M.A.; Al-hawamdeh, A.; Amarin, Z.O.; Khader, Y.S.; Mhedat, K.; Allouh, M.Z. Single Nucleotide Polymorphism in the ADIPOQ Gene Modifies Adiponectin Levels and Glycemic Control in Type Two Diabetes Mellitus Patients. BioMed Research International 2022, 2022, 6632442. [Google Scholar] [CrossRef]
  12. Zhang, X.; Sun, J.; Han, W.; Jiang, Y.; Peng, S.; Shan, Z.; Teng, W. The Type 2 Deiodinase Thr92Ala Polymorphism Is Associated with Worse Glycemic Control in Patients with Type 2 Diabetes Mellitus: A Systematic Review and Meta-Analysis. J. Diabetes Res. 2016, 2016, 5928726. [Google Scholar] [CrossRef]
  13. Lamichane, S.; Dahal Lamichane, B.; Kwon, S.-M. Pivotal Roles of Peroxisome Proliferator-Activated Receptors (PPARs) and Their Signal Cascade for Cellular and Whole-Body Energy Homeostasis. Int. J. Mol. Sci. 2018, 19, 949. [Google Scholar] [CrossRef]
  14. Mirza, A.Z.; Althagafi, I.I.; Shamshad, H. Role of PPAR Receptor in Different Diseases and Their Ligands: Physiological Importance and Clinical Implications. Eur. J. Med. Chem. 2019, 166, 502–513. [Google Scholar] [CrossRef]
  15. Hong, F.; Pan, S.; Guo, Y.; Xu, P.; Zhai, Y. PPARs as Nuclear Receptors for Nutrient and Energy Metabolism. Molecules 2019, 24, 2545. [Google Scholar] [CrossRef]
  16. Levine, T.B.; Levine, A.B. Metabolic Syndrome and Cardiovascular Disease; Saunders/Elsevier: Philadelphia, PA, USA, 2006; ISBN 978-1-4160-2545-0. [Google Scholar]
  17. Frkic, R.L.; Richter, K.; Bruning, J.B. The Therapeutic Potential of Inhibiting PPARγ Phosphorylation to Treat Type 2 Diabetes. J. Biol. Chem. 2021, 297, 101030. [Google Scholar] [CrossRef]
  18. Lee, C.-H.; Olson, P.; Hevener, A.; Mehl, I.; Chong, L.-W.; Olefsky, J.M.; Gonzalez, F.J.; Ham, J.; Kang, H.; Peters, J.M.; et al. PPARδ Regulates Glucose Metabolism and Insulin Sensitivity. Proc. Natl. Acad. Sci. USA 2006, 103, 3444–3449. [Google Scholar] [CrossRef]
  19. Liu, S.; Hatano, B.; Zhao, M.; Yen, C.-C.; Kang, K.; Reilly, S.M.; Gangl, M.R.; Gorgun, C.; Balschi, J.A.; Ntambi, J.M.; et al. Role of Peroxisome Proliferator-Activated Receptor δ/β in Hepatic Metabolic Regulation. J. Biol. Chem. 2011, 286, 1237–1247. [Google Scholar] [CrossRef]
  20. Dubois, V.; Eeckhoute, J.; Lefebvre, P.; Staels, B. Distinct but Complementary Contributions of PPAR Isotypes to Energy Homeostasis. J. Clin. Investig. 2017, 127, 1202–1214. [Google Scholar] [CrossRef]
  21. Ahmadian, M.; Suh, J.M.; Hah, N.; Liddle, C.; Atkins, A.R.; Downes, M.; Evans, R.M. PPARγ Signaling and Metabolism: The Good, the Bad and the Future. Nat. Med. 2013, 19, 557–566. [Google Scholar] [CrossRef]
  22. Yen, C.-J.; Beamer, B.A.; Negri, C.; Silver, K.; Brown, K.A.; Yarnall, D.P.; Burns, D.K.; Roth, J.; Shuldiner, A.R. Molecular Scanning of the Human Peroxisome Proliferator Activated Receptor g (hPPARg) Gene in Diabetic Caucasians: Identification of a Pro12Ala PPARg2 Missense Mutation. Biochem. Biophys. Res. Commun. 1997, 241, 270–274. [Google Scholar] [CrossRef]
  23. Lagou, V.; Scott, R.A.; Manios, Y.; Chen, T.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]
  24. Stryjecki, C.; Peralta-Romero, J.; Alyass, A.; Karam-Araujo, R.; Suarez, F.; Gomez-Zamudio, J.; Burguete-Garcia, A.; Cruz, M.; Meyre, D. Association between PPAR-Γ2 Pro12Ala Genotype and Insulin Resistance Is Modified by Circulating Lipids in Mexican Children. Sci. Rep. 2016, 6, 24472. [Google Scholar] [CrossRef]
  25. 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]
  26. Attianese, G.M.G.; Desvergne, B. Integrative and Systemic Approaches for Evaluating PPARβ/δ (PPARD) Function. Nucl. Recept. Signal. 2015, 13, e001. [Google Scholar] [CrossRef]
  27. 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. [Google Scholar] [CrossRef]
  28. Wang, T.; Song, J.-F.; Zhou, X.-Y.; Li, C.-L.; Yin, X.-X.; Lu, Q. PPARD Rs2016520 (T/C) and NOS1AP Rs12742393 (A/C) Polymorphisms Affect Therapeutic Efficacy of Nateglinide in Chinese Patients with Type 2 Diabetes Mellitus. BMC Med. Genom. 2021, 14, 267. [Google Scholar] [CrossRef] [PubMed]
  29. Secretaría de Salud. NORMA Oficial Mexicana NOM-015-SSA2-2010, Para La Prevención, Tratamiento y Control de La Diabetes Mellitus. Diario Oficial de La Federación 23 noviembre 2010. Estados Unidos Mexicanos. Available online: https://dof.gob.mx/nota_detalle.php?codigo=5168074&fecha=23/11/2010#gsc.tab=0 (accessed on 20 November 2023).
  30. Millán, J.; Pintó, X.; Muñoz, A.; Zúñiga, M.; Rubiés-Prat, J.; Pallardo, L.F.; Masana, L.; Mangas, A.; Mijares, A.H.; Santos, P.G.; et al. Cocientes lipoproteicos: Significado fisiológico y utilidad clínica de los índices aterogénicos en prevención cardiovascular. Clínica E Investig. Arterioscler. 2010, 22, 25–32. [Google Scholar] [CrossRef]
  31. Syed, R.; Jamil, K.; Asimuddin, M.; Alqahtani, M.S.; Alshehri, M.; Mateen, A.; Wahab Ali Aduderman, A.; Shamsul Ola, M.; Malik, A. Molecular & Biochemical Analysis of Pro12Ala Variant of PPAR-Γ2 Gene in Type 2 Diabetes Mellitus. Saudi J. Biol. Sci. 2020, 27, 2439–2443. [Google Scholar] [CrossRef]
  32. Abrahams-October, Z.; Xhakaza, L.; Pearce, B.; Mandisa Masilela, C.; Benjeddou, M.; Vincent Adeniyi, O.; Johnson, R.; Jebio Ongole, J. Genetic Association of Solute Carrier Transporter Gene Variants with Metformin Response. Balk. J. Med. Genet. 2021, 24, 47–56. [Google Scholar] [CrossRef]
  33. Grygiel-Górniak, B.; Mosor, M.; Marcinkowska, J.; Przysławski, J.; Nowak, J. Impact of the PPAR Gamma-2 Gene Polymorphisms on the Metabolic State of Postmenopausal Women. J. Biosci. 2016, 41, 427–437. [Google Scholar] [CrossRef]
  34. Muntean, C.; Sasaran, M.O.; Crisan, A.; Banescu, C. Effects of PPARG and PPARGC1A Gene Polymorphisms on Obesity Markers. Front. Public Health 2022, 10, 962852. [Google Scholar] [CrossRef]
  35. Martínez-Gómez, L.E.; Cruz, M.; Martínez-Nava, G.A.; Madrid-Marina, V.; Parra, E.; García-Mena, J.; Espinoza-Rojo, M.; Estrada-Velasco, B.I.; Piza-Roman, L.F.; Aguilera, P.; et al. A Replication Study of the IRS1, CAPN10, TCF7L2, and PPARG Gene Polymorphisms Associated with Type 2 Diabetes in Two Different Populations of Mexico: IRS1, CAPN10, TCF7L2, and PPARG Polymorphisms and T2D in Mexico Population. Ann. Hum. Genet. 2011, 75, 612–620. [Google Scholar] [CrossRef]
  36. Nikitin, A.G.; Chistiakov, D.A.; Minushkina, L.O.; Zateyshchikov, D.A.; Nosikov, V.V. Association of the CYBA, PPARGC1A, PPARG3, and PPARD Gene Variants with Coronary Artery Disease and Metabolic Risk Factors of Coronary Atherosclerosis in a Russian Population. Heart Vessels 2010, 25, 229–236. [Google Scholar] [CrossRef]
  37. 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]
  38. Wei, X.-L.; Yin, R.-X.; Miao, L.; Wu, D.-F. The Peroxisome Proliferator-Activated Receptor Delta +294T > C Polymorphism and Alcohol Consumption on Serum Lipid Levels. Lipids Health Dis. 2011, 10, 242. [Google Scholar] [CrossRef]
  39. Bin Rakhis, S.A.; AlDuwayhis, N.M.; Aleid, N.; AlBarrak, A.N.; Aloraini, A.A. Glycemic Control for Type 2 Diabetes Mellitus Patients: A Systematic Review. Cureus 2022, 14, e26180. [Google Scholar] [CrossRef]
  40. Basto-Abreu, A.C.; López-Olmedo, N.; Rojas-Martínez, R.; Aguilar-Salinas, C.A.; De La Cruz-Góngora, V.V.; Rivera-Dommarco, J.; Shamah-Levy, T.; Romero-Martínez, M.; Barquera, S.; Villalpando, S.; et al. Prevalence of Diabetes and Glycemic Control in Mexico: National Results from 2018 and 2020. Salud Pública México 2021, 63, 725–733. [Google Scholar] [CrossRef]
  41. Huang, R.; Yan, L.; Lei, Y. The Relationship between High-Density Lipoprotein Cholesterol (HDL-C) and Glycosylated Hemoglobin in Diabetic Patients Aged 20 or above: A Cross-Sectional Study. BMC Endocr. Disord. 2021, 21, 198. [Google Scholar] [CrossRef]
  42. Suprapti, B.; Izzah, Z.; Anjani, A.G.; Andarsari, M.R.; Nilamsari, W.P.; Nugroho, C.W. Prevalence of Medication Adherence and Glycemic Control among Patients with Type 2 Diabetes and Influencing Factors: A Cross-Sectional Study. Glob. Epidemiol. 2023, 5, 100113. [Google Scholar] [CrossRef]
  43. Instituto Mexicano del Seguro Social. Diagnóstico y Tratamiento Farmacológico de la Diabetes Mellitus Tipo 2 en el Primer Nivel de Atención. Guía de Evidencias y Recomendaciones: Guía de Práctica Clínica; Instituto Mexicano del Seguro Social: Cd. de Mexico, Mexico, 2018; pp. 1–55.
  44. World Health Organization Global Report on Diabetes; World Health Organization: Geneva, Switzerland, 2016; ISBN 978-92-4-156525-7.
  45. Saleh, R.; Zahid, Z.I.; Rahman, M.A.; Jain, P.; Alam, A.; Kawaichi, M.; Reza, H.M. Prevalence of PPAR-Γ2 (Rs1801282), RETN (Rs3745367) and ADIPOQ (Rs2241766) SNP Markers in the Bangladeshi Type 2 Diabetic Population. Meta Gene 2016, 10, 100–107. [Google Scholar] [CrossRef]
  46. Pei, Q.; Huang, Q.; Yang, G.; Zhao, Y.; Yin, J.; Song, M.; Zheng, Y.; Mo, Z.; Zhou, H.; Liu, Z. PPAR-Γ2 and PTPRD Gene Polymorphisms Influence Type 2 Diabetes Patients’ Response to Pioglitazone in China. Acta Pharmacol. Sin. 2013, 34, 255–261. [Google Scholar] [CrossRef]
  47. Swarbrick, M.M.; Chapman, C.M.; McQuillan, B.M.; Hung, J.; Thompson, P.L.; Beilby, J.P. A Pro12Ala Polymorphism in the Human Peroxisome Proliferator-Activated Receptor-Gamma 2 Is Associated with Combined Hyperlipidaemia in Obesity. Eur. J. Endocrinol. 2001, 144, 277–282. [Google Scholar] [CrossRef] [PubMed]
  48. Chistiakov, D.A.; Potapov, V.A.; Khodirev, D.S.; Shamkhalova, M.S.; Shestakova, M.V.; Nosikov, V.V. The PPARγ Pro12Ala Variant Is Associated with Insulin Sensitivity in Russian Normoglycaemic and Type 2 Diabetic Subjects. Diab. Vasc. Dis. Res. 2010, 7, 56–62. [Google Scholar] [CrossRef]
  49. Saremi, L. Association between PRO12ALa Polymorphism of PARγ2 Gene and Coronary Artery Disease in Iranian Population with Type Two Diabetes Mellitus. Acta Endocrinol. Buchar. 2022, 18, 139–144. [Google Scholar] [CrossRef] [PubMed]
  50. Song, J.; Li, N.; Hu, R.; Yu, Y.; Xu, K.; Ling, H.; Lu, Q.; Yang, T.; Wang, T.; Yin, X. Effects of PPARD Gene Variants on the Therapeutic Responses to Exenatide in Chinese Patients with Type 2 Diabetes Mellitus. Front. Endocrinol. 2022, 13, 949990. [Google Scholar] [CrossRef]
  51. Li, S.; He, C.; Nie, H.; Pang, Q.; Wang, R.; Zeng, Z.; Song, Y. G Allele of the Rs1801282 Polymorphism in PPARγ Gene Confers an Increased Risk of Obesity and Hypercholesterolemia, While T Allele of the Rs3856806 Polymorphism Displays a Protective Role Against Dyslipidemia: A Systematic Review and Meta-Analysis. Front. Endocrinol. 2022, 13, 919087. [Google Scholar] [CrossRef]
  52. 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 PPARgamma2 Associated with Decreased Receptor Activity, Lower Body Mass Index and Improved Insulin Sensitivity. Nat. Genet. 1998, 20, 284–287. [Google Scholar] [CrossRef]
  53. Iwaki, M.; Matsuda, M.; Maeda, N.; Funahashi, T.; Matsuzawa, Y.; Makishima, M.; Shimomura, I. Induction of Adiponectin, a Fat-Derived Antidiabetic and Antiatherogenic Factor, by Nuclear Receptors. Diabetes 2003, 52, 1655–1663. [Google Scholar] [CrossRef]
  54. Heikkinen, S.; Argmann, C.; Feige, J.N.; Koutnikova, H.; Champy, M.-F.; Dali-Youcef, N.; Schadt, E.E.; Laakso, M.; Auwerx, J. The Pro12Ala PPARγ2 Variant Determines Metabolism at the Gene-Environment Interface. Cell Metab. 2009, 9, 88–98. [Google Scholar] [CrossRef]
  55. Mousavinasab, F.; Tähtinen, T.; Jokelainen, J.; Koskela, P.; Vanhala, M.; Oikarinen, J.; Keinänen-Kiukaanniemi, S.; Laakso, M. Effect of the Pro12Ala Polymorphism of the PPARγ2 Gene on Serum Adiponectin Changes. Endocrine 2005, 27, 307–309. [Google Scholar] [CrossRef]
  56. Maciejewska-Skrendo, A.; Massidda, M.; Tocco, F.; Leźnicka, K. The Influence of the Differentiation of Genes Encoding Peroxisome Proliferator-Activated Receptors and Their Coactivators on Nutrient and Energy Metabolism. Nutrients 2022, 14, 5378. [Google Scholar] [CrossRef] [PubMed]
  57. Skogsberg, J.; Kannisto, K.; Cassel, T.N.; Hamsten, A.; Eriksson, P.; Ehrenborg, E. Evidence That Peroxisome Proliferator–Activated Receptor Delta Influences Cholesterol Metabolism in Men. Arterioscler. Thromb. Vasc. Biol. 2003, 23, 637–643. [Google Scholar] [CrossRef] [PubMed]
  58. Kolehmainen, M.; Uusitupa, M.I.J.; Alhava, E.; Laakso, M.; Vidal, H. Effect of the Pro12Ala Polymorphism in the Peroxisome Proliferator-Activated Receptor (PPAR) Γ2 Gene on the Expression of PPARγ Target Genes in Adipose Tissue of Massively Obese Subjects. J. Clin. Endocrinol. Metab. 2003, 88, 1717–1722. [Google Scholar] [CrossRef] [PubMed]
  59. Leońska-Duniec, A.; Cieszczyk, P.; Jastrzębski, Z.; Jażdżewska, A.; Lulińska-Kuklik, E.; Moska, W.; Ficek, K.; Niewczas, M.; Maciejewska-Skrendo, A. The Polymorphisms of the PPARD Gene Modify Post-Training Body Mass and Biochemical Parameter Changes in Women. PLoS ONE 2018, 13, e0202557. [Google Scholar] [CrossRef]
  60. Damsgaard, C.T.; Vuholm, S.; Teisen, M.N.; Stark, K.D.; Lauritzen, L. Does Polymorphisms in PPAR and APOE Genes Modify Associations between Fatty Acid Desaturase (FADS), n-3 Long-Chain PUFA and Cardiometabolic Markers in 8–11-Year-Old Danish Children? Br. J. Nutr. 2021, 125, 369–376. [Google Scholar] [CrossRef]
  61. Loera Morales, J., III; Martínez Bautista, H.; Torres Yañez, E.E.; Almazán Orta, L.; Vázquez Martínez, V.H. Prevalencia y factores asociados a adherencia terapéutica de pacientes con diabetes mellitus 2: Post pandemia por COVID 19. Arch. En Med. Fam. 2023, 25, 235–245. [Google Scholar]
  62. Instituto Nacional de Estadística y Geografía, (INEGI) Censo de Población y Vivienda 2020. Derechohabiencia. Available online: https://www.inegi.org.mx/temas/derechohabiencia/ (accessed on 20 November 2023).
Figure 1. Fasting plasma glucose (geometric means ± SD) by PPAR polymorphisms (inheritance dominant model) adjusted for age at the time of T2D diagnosis. * p < 0.05; p = 0.013 for the interaction effect. Full data in Table S2 including over-dominant model results.
Figure 1. Fasting plasma glucose (geometric means ± SD) by PPAR polymorphisms (inheritance dominant model) adjusted for age at the time of T2D diagnosis. * p < 0.05; p = 0.013 for the interaction effect. Full data in Table S2 including over-dominant model results.
Biomolecules 13 01791 g001
Table 1. Genotyping primers for PPAR-γ2 and PPAR-β/δ.
Table 1. Genotyping primers for PPAR-γ2 and PPAR-β/δ.
SNPSequenceMajor Alelle
Genotype
Heterozygous
Genotype
Minor
Alelle Genotype
rs1801282P1: 5′-GTGTATCAGTGAAGGAATCGCTTTCTTG-3
P2: 5′-TTGTGATATGTTTGCAGACAAGGTATCAG
TGAAGGAATCGCTTTGTGC-3′
P3: 5′-TTTCTGTGTTTATTCCCATCTCTCCC-3′
230 bp230 and 250 bp250 bp
rs2016520F: 5′-CATGGTATAGCACTGCAGGAA-3′
R: 5′-CTTCCTCCTGTGGCTGCTC-3′
269 bp269, 167, and 102 bp167 and 102 bp
PPAR-γ2 rs1801282: Major homozygous genotype CC (Pro/Pro alleles), heterozygous CG (Pro/Ala alleles), and minor homozygous genotype GG (Ala/Ala alleles). The PPAR-β/δ rs2016520: Major homozygous genotype TT alleles, heterozygous genotype TC alleles, and minor homozygous genotype CC alleles.
Table 2. Genotype, allele frequencies, and H-WE of PPAR gene polymorphisms.
Table 2. Genotype, allele frequencies, and H-WE of PPAR gene polymorphisms.
PPAR-γ2 rs1801282PPAR-β/δ rs2016520

Genotype
n
(%)

CC
238
(76.77)

CG
67
(21.61)

GG
5
(1.61)
H-WE
p

TT
231
(73.57)

TC
79
(25.16)

CC
4
(1.27)
H-WE
p
0.990.63
AlleleCGTC
n5437754187
(%)(87.58)(12.42)(86.15)(13.85)
Χ2 test.rs2801282-G = minor allele, rs2016520-C = minor allele. Note: four samples did not amplify for PPAR-γ2 rs2801282.
Table 3. Anthropometric and biochemical measures by PPAR genotypes (dominant model).
Table 3. Anthropometric and biochemical measures by PPAR genotypes (dominant model).
VariablePPAR-γ2 rs1801282PPAR-β/δ rs2016520
CC (n = 238) Md (IQR)CG + GG (n = 73)
Md (IQR)
pTT (n = 231)
Md (IQR)
TC + CC (n = 83) Md (IQR)p
BMI (kg/m2)31
(27–35)
31
(28–33)
0.76831
(27–35)
32
(28–35)
0.238
Waist circumference (cm)101
(95–110)
106
(99–115)
0.029100
(95–110)
103
(100–115)
0.052
FPG (mg/dL)144
(114–179)
157
(127–207)
0.026145
(117–177)
158
(125–195)
0.089
HbA1c (%)7
(6–9)
8
(6–10)
0.0577
(6–9)
8
(6–10)
0.080
TG (mg/dL)177
(126–235)
181
(138–254)
0.250174
(124–227)
190
(155–282)
0.026
TC (mg/dL)194
(164–220)
192
(168–218)
0.838193
(165–217)
189
(163–226)
0.526
HDL-C (mg/dL)44
(36–54)
47
(38–56)
0.26044
(37–54)
46
(37–56)
0.457
LDL-C (mg/d L)110
(80–134)
99
(73–132)
0.224109
(79–132)
103
(65–140)
0.671
VLDL (mg/dL)35
(25–47)
36
(28–51)
0.25035
(25–45)
38
(31–56)
0.026
TC/HDL-C index4
(3–5)
4
(3–5)
0.3034
(3–5)
4
(3–5)
0.964
BMI: body mass index; FPG: fasting plasma glucose; TG: triglycerides; TC: total cholesterol; HDL-C: high-density lipoprotein; LDL-C: low-density lipoprotein; Md: median; IQR: interquartile range.
Table 4. Multivariate analysis of the relation between PPAR polymorphisms and HbA1c (log-transformed).
Table 4. Multivariate analysis of the relation between PPAR polymorphisms and HbA1c (log-transformed).
VariableDominant Model
n = 213
Over-Dominant Model
n = 213
β (95%CI)
p
β (95%CI)
p
PPAR-γ2 rs18012820.09 (0.02, 0.16) a
0.011
0.08 (0.01, 0.16) c
0.020
PPAR-β/δ rs20165200.10 (0.03, 0.16) b
0.005
0.10 (0.04, 0.17) d
0.003
Body mass index e
Overweight−0.11 (−0.21, −0.006)
0.038
−0.11 (−0.21, −0.005)
0.041
Obesity grade I−0.15 (−0.25, −0.05)
0.004
−0.15 (−0.25, −0.04)
0.005
Obesity grade II−0.11 (−0.22, 0.01)
0.064
−0.11 (−0.22, 0.01)
0.066
Obesity grade III−0.18 (−0.31, −0.05)
0.005
−0.18 (−0.31, −0.05)
0.006
Age at the time ofT2D diagnosis (years)−0.01 (−0.01, −0.004)
<0.001
−0.01 (−0.01, −0.004)
<0.001
HDL-C (mg/dL)−0.003 (−0.005, −0.0004)
0.020
−0.003 (−0.005, −0.0003)
0.024
a reference CC, b reference TT, c reference CC + GG, d reference TT + CC, e reference: normal weight. 95%CI; 95% confidence interval. Adj R2 = 0.16 in both models.
Table 5. Relation between triglycerides (log mg/dL) and fasting plasma glucose by PPAR polymorphisms (inheritance dominant model).
Table 5. Relation between triglycerides (log mg/dL) and fasting plasma glucose by PPAR polymorphisms (inheritance dominant model).
VariablePPAR-γ2 rs1801282-CC
&
PPAR-β/δ rs2016520-TT
n = 168
PPAR-γ2 rs1801282-CG + GG &
PPAR-β/δ rs2016520-TT
n = 61
PPAR-γ2 rs1801282-CC
&
PPAR-β/δ rs2016520-TC + CC
n = 50
PPAR-γ2 rs1801282-CG + GG
&
PPAR-β/δ rs2016520-TC + CC
n = 16
β
(95%CI)
p
β
(95%CI)
p
β
(95%CI)
p
β
(95%CI)
p
Body mass index
Overweight0.23
(0.013, 0.454)
0.038
0.31
(−0.457, 1.072)
0.422
0.44
(−0.008, 0.890)
0.054
0.49
(−0.257, 1.242)
0.174
Obesity class I0.26
(0.038, 0.474)
0.022
0.27
(−0.500, 1.048)
0.479
0.24
(−0.205, 0.690)
0.282
0.93
(0.263, 1.598)
0.011
Obesity class II0.28
(0.030, 0.521)
0.028
0.17
(−0.658, 1.004)
0.677
0.18
(−0.307, 0.672)
0.458
0.99
(0.247, 1.726)
0.014
Obesity class III0.19
(−0.093, 0.470)
0.187
−0.03
(−0.970, 0.914)
0.952
0.52
(−0.108, 1.144)
0.103
0.18
(−0.577, 0.946)
0.600
FPG (mg/dL)0.002
(0.00, 0.003)
0.005
0.00
(−0.002, 0.002)
0.983
0.004
(0.002, 0.005)
<0.001
0.002
(0.00, 0.004)
0.033
Adj R20.060.000.250.61
FPG: fasting plasma glucose; 95%CI: 95% confidence interval.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Reza-López, S.A.; González-Gurrola, S.; Morales-Morales, O.O.; Moreno-González, J.G.; Rivas-Gómez, A.M.; González-Rodríguez, E.; Moreno-Brito, V.; Licón-Trillo, A.; Leal-Berumen, I. Metabolic Biomarkers in Adults with Type 2 Diabetes: The Role of PPAR-γ2 and PPAR-β/δ Polymorphisms. Biomolecules 2023, 13, 1791. https://doi.org/10.3390/biom13121791

AMA Style

Reza-López SA, González-Gurrola S, Morales-Morales OO, Moreno-González JG, Rivas-Gómez AM, González-Rodríguez E, Moreno-Brito V, Licón-Trillo A, Leal-Berumen I. Metabolic Biomarkers in Adults with Type 2 Diabetes: The Role of PPAR-γ2 and PPAR-β/δ Polymorphisms. Biomolecules. 2023; 13(12):1791. https://doi.org/10.3390/biom13121791

Chicago/Turabian Style

Reza-López, Sandra A., Susana González-Gurrola, Oscar O. Morales-Morales, Janette G. Moreno-González, Ana M. Rivas-Gómez, Everardo González-Rodríguez, Verónica Moreno-Brito, Angel Licón-Trillo, and Irene Leal-Berumen. 2023. "Metabolic Biomarkers in Adults with Type 2 Diabetes: The Role of PPAR-γ2 and PPAR-β/δ Polymorphisms" Biomolecules 13, no. 12: 1791. https://doi.org/10.3390/biom13121791

APA Style

Reza-López, S. A., González-Gurrola, S., Morales-Morales, O. O., Moreno-González, J. G., Rivas-Gómez, A. M., González-Rodríguez, E., Moreno-Brito, V., Licón-Trillo, A., & Leal-Berumen, I. (2023). Metabolic Biomarkers in Adults with Type 2 Diabetes: The Role of PPAR-γ2 and PPAR-β/δ Polymorphisms. Biomolecules, 13(12), 1791. https://doi.org/10.3390/biom13121791

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