Oxidative Medicine and Cellular Longevity
Volume 2016 (2016), Article ID 7367641, 8 pages
http://dx.doi.org/10.1155/2016/7367641

Research Article

Association of Nuclear Factor-Erythroid 2-Related Factor 2, Thioredoxin Interacting Protein, and Heme Oxygenase-1 Gene Polymorphisms with Diabetes and Obesity in Mexican Patients

Angélica Saraí Jiménez-Osorio,¹ Susana González-Reyes,¹ Wylly Ramsés García-Niño,¹ Hortensia Moreno-Macías,² Martha Eunice Rodríguez-Arellano,³ Gilberto Vargas-Alarcón,⁴ Joaquín Zúñiga,⁵ Rodrigo Barquera,⁶ and José Pedraza-Chaverri¹

¹Department of Biology, National Autonomous University of Mexico (UNAM), 04510 Mexico City DF, Mexico
²Economy Department, Autonomous Metropolitan University-Iztapalapa, 09340 Mexico City, DF, Mexico
³Research Department, Regional Hospital “Lic. Adolfo López Mateos”, ISSSTE, 01030 Mexico City, DF, Mexico
⁴Molecular Biology Department, National Institute of Cardiology “Ignacio Chávez”, 14080 Mexico City, DF, Mexico
⁵Department of Immunology, National Institute of Respiratory Diseases “Ismael Cosío Villegas”, 14080 Mexico City, DF, Mexico
⁶Molecular Genetics Laboratory, National School of Anthropology and History, 14030 Mexico City, DF, Mexico

Received 4 January 2016; Revised 29 March 2016; Accepted 5 April 2016

Academic Editor: Silvana Hrelia

Copyright © 2016 Angélica Saraí Jiménez-Osorio et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The nuclear factor-erythroid 2- (NF-E2-) related factor 2 (Nrf2) is abated and its ability to reduce oxidative stress is impaired in type 2 diabetes and obesity. Thus, the aim of this study was to explore if polymorphisms in Nrf2 and target genes are associated with diabetes and obesity in Mexican mestizo subjects. The rs1800566 of NAD(P)H:quinone oxidoreductase 1 (NQO1) gene, rs7211 of thioredoxin interacting protein (TXNIP) gene, rs2071749 of heme oxygenase-1 (HMOX1) gene, and the rs6721961 and the rs2364723 from Nrf2 gene were genotyped in 627 diabetic subjects and 1020 controls. The results showed that the rs7211 polymorphism is a protective factor against obesity in nondiabetic subjects (CC + CT versus TT, OR = 0.40, ) and in women (CC versus CT + TT, OR = 0.7, ). TT carriers had lower high-density lipoprotein cholesterol levels and lower body mass index. The rs2071749 was positively associated with obesity (AA versus AG + GG, OR = 1.25, ). Finally, the rs6721961 was negatively associated with diabetes in men (CC versus CA + AA, OR = 0.62, ). AA carriers showed lower glucose concentrations. No association was found for rs1800566 and rs2364723 polymorphisms. In conclusion, the presence of Nrf2 and related genes polymorphisms are associated with diabetes and obesity in Mexican patients.

1. Introduction

The International Diabetes Federation has reported that there are 382 million people living with diabetes worldwide [1]. Type 2 diabetes (T2DM) is a public health problem in developing countries with a direct impact in the economic and social sectors. Being overweight and obesity are risk factors to develop T2DM and could explain the dramatic increase in the incidence and prevalence of T2DM [2, 3]. Both, obesity and T2DM, feature insulin resistance and atherogenic lipid profiles such as increased cholesterol and triglycerides and decreased high-density lipoprotein cholesterol (HDL-C). It is widely accepted that oxidative stress is a common mechanism in the development and progression of these pathologies, with an increased free radical production and reduced antioxidant capacity [4, 5].

The master antioxidant regulator, the nuclear factor-erythroid 2- (NF-E2-) related factor 2 (Nrf2) is a member of the cap’n’collar family of basic leucine zipper transcription factors that regulates the expression of many antioxidant genes including NAD(P)H:quinone oxidoreductase 1 (NQO1) and heme oxygenase-1 (HMOX1) to avoid oxidative damage [6]. Also, it has been demonstrated that Nrf2 activity is abated in diabetes and factors like age, body weight, and blood glucose could modify its activity [7–9], but genetic factors have been poorly studied. Many single nucleotide polymorphisms (SNPs) have been identified in the Nrf2 gene [10]. In particular, the rs6721961 (C-617C) polymorphism has been associated with oxidative stress and risk of newly diagnosed T2DM [7] and increased blood pressure [11]. The rs2364723 (C107G) polymorphism was associated with reduced risk of cardiovascular mortality [12]. The NQO1 gene polymorphism rs1800566 (C609T) is characterized by the replacement of proline by serine at position number 187 of the functional protein, which causes destabilization and inactivation of the enzyme [13, 14]. The rs1800566 polymorphism has been associated with lower levels of blood coagulation factors [15], higher risk of coronary artery disease in T2DM patients [16, 17], and increased triglycerides levels and decreased HDL-C levels in individuals with metabolic syndrome [18]. Polymorphisms in the promoter region of the HMOX1 gene such as rs2071746 (T-413A) and microsatellite have been associated with various human diseases [19]; however, the rs2071749 polymorphism (A/G) is not associated with hypertension and/or blood pressure in hypertensive patients [20]. It has been reported that the endogenous inhibitor of the thioredoxin (TXN) system, the thioredoxin interacting protein (TXNIP), can be suppressed by Nrf2 [21, 22]. The rs7211 polymorphism in TXNIP was associated with inhibition of TNX, glucose homeostasis, diabetes, and hypertension [23]. Individuals carrying the T allele of the rs7211 polymorphism in the TXNIP gene showed higher plasma triglycerides levels [24].

However, despite the recent efforts to understand the association of antioxidant gene polymorphism with T2DM, the number of studies is still very limited in Mexican population. Thus, the aim of this study was to investigate potential associations between Nrf2 (rs6721961 and rs2364723), NQO1 (rs1800566), HMOX1 (rs2071749), and TXNIP (rs7211) polymorphisms in a T2DM population.

2. Materials and Methods

2.1. Study Design and Subjects

A case control study was performed which included 1647 Mexican mestizo subjects. This study was approved by the Ethics Committee on Human Studies from the Committee on Research, Ethics and Safety of the Hospital Regional “Lic. Adolfo López Mateos” with registration number 254.2013 and conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all subjects who were recruited in Family Medical Clinics of the Instituto de Seguridad y Servicios Sociales de los Trabajadores del Estado (ISSSTE) in Mexico City.

A standard questionnaire was applied to obtain demographic information, family and personal health history (diabetes and/or further diseases), and information about physical activity, alcohol consumption, smoking, and drugs consumption. Their medical history data, weight (kg), and height (m) were obtained. The body mass index (BMI) was estimated by dividing weight by the square of height. Criteria for inclusion/exclusion of participants were as follows: the study included men and women aged 35 and older, nonpregnant and nonlactating women, and subjects without excessive alcohol consumption. The criteria for classification and diagnosis of diabetes were according to the standards in medical care on diabetes from American Diabetes Association (ADA). Subjects of the control group () had fasting glucose <5.5 mmol/L and glycated hemoglobin (HbA1c) <5.7%. Subjects with T2DM () had plasma glucose >7.0 mmol/L and HbA1c >6.5%. Subjects with mental health problems (senile dementia and Alzheimer and Parkinson’s disease) and cancer were also excluded. The recruitment period was from May 2012 to October 2014. A total of 1627 individuals were recruited that self-reported Mexican mestizo ancestry (three generations). To each individual the information was given about the protocol and those who decided to enter gave written informed consent to participate. Anthropometric data (waist circumference, height, and weight) and arterial pressure measurements were obtained before blood collection.

2.2. Blood Collection and Biochemical Analyses

Whole blood samples (20 mL) were collected from patients and controls (with fasting of 8–12 h) in order to determine the levels of glucose, triglycerides, total cholesterol, HDL-C, low-density lipoprotein cholesterol (LDL-C), and creatinine in serum using an automatized analyzer (Miura 200, ISE, Rome, Italy). Moreover, total blood with ethylenediaminetetraacetic (EDTA) acid as anticoagulant was used to obtain genomic DNA.

2.3. Genotyping

Genomic DNA from whole blood containing EDTA was isolated by standard techniques [25]. The rs2364723 (Nrf2), rs1800566 (NQO1), rs7211 (TXNIP), and rs2071749 (HMOX1) SNPs were genotyped using predesigned 5′ exonuclease TaqMan genotyping assays on a 7500 series Real-Time PCR system, according to manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA). The rs6721961 (Nrf2) TaqMan probe was designed according to manufacturer’s directions.

2.4. Short Tandem Repeats (STRs) Genotyping and Admixture Estimations

Fifteen autosomal STR markers (CSF1PO, FGA, THO1, TPOX, VWA, D3S11358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D21S11, D19S433, and D2S1338) along with amelogenin were genotyped in 200 controls and 200 cases, using the AmpFlSTR Identifiler Kit (Applied Biosystems, Foster City, CA, USA) as previously described [26].

We performed admixture estimations using the STR’s allele distribution by a model-based clustering method with the Structure software v. 2.3.4, assuming populations and dememorisation steps. The estimations of Amerindian, European, and African components in our Mexican studied groups were performed using the distribution of STRs alleles in different populations including Spaniards [27], Fang Africans [28], and a Native American pool of Huastecos [29] and Tepehuas [30] from the central region of Mexico, who were considered as parental populations.

2.5. Statistical Analyses

Continuous variables were analyzed using a -test and were presented as mean ± standard deviation (SD). For categorical variables, a Chi-square test or Fischer exact test was applied and data was presented as percentage.

The Hardy Weinberg equilibrium was calculated in controls using StatCalc software (Epi Info 2005 v3.3.2; Centers of Disease Control and Prevention, Atlanta, GA, USA). Multivariable linear regression models were carried out for adjustment of glucose for potential confounders like age, gender, BMI, and tobacco. In order to assess the genetic risk factor for diabetes and obesity, logistic regression analyses were applied to estimate the OR for each polymorphism. Because the associations between the SNPs and the outcomes have been previously reported, it is unlikely to detect effects due to statistical fluctuations only. Therefore, correction by multiple comparisons was not applied. Associations were considered statistically significant at a nominal value ≤0.05. Haplotype analysis was performed using PLINK and Haploview software. The -test, Chi-square test, Fischer exact test, and multivariable regressions were performed using the statistical software Stata 12.0 (StataCorp LP, College Station, TX, USA).

3. Results

3.1. Clinical and Anthropometrical Measures

General characteristics of the study population are shown in Table 1. The mean of age, BMI, glucose, HbA1c, triglycerides, LDL-C, and systolic blood pressure were significantly higher in subjects with diabetes as compared with control group. On the other hand, the HDL-C levels were found lower in diabetic subjects (), which together show some of the metabolic abnormalities associated with T2DM. Taking into account that obesity frequency was found higher in the control group (), it was considered as a potential confounder in multivariate analyses.

Table 1: Clinical and anthropometric characteristics of the study groups.

Characteristic
Control
Diabetes

1020
627

Age (years)
42.1 ± 7.3
53.2 ± 9.4
<0.001

Males, (%)
474 (46)
302 (48)
0.759

BMI (kg⋅m⁻²)
28.7 ± 4.9
29.7 ± 6.9
0.0016

Obesity, (%)
350 (34)
289 (46)
<0.001

Glucose (mg⋅dL⁻¹)
95.5 ± 9.6
154 ± 73
<0.001

HbA1c
4.45 ± 0.8
7.27 ± 2.1
<0.001

Triglycerides (mg⋅dL⁻¹)
183 ± 122
214 ± 140
<0.001

Total cholesterol (mg⋅dL⁻¹)
201 ± 43
198 ± 52
0.2257

HDL-C (mg⋅dL⁻¹)
49.5 ± 12.5
45.4 ± 13
<0.001

LDL-C (mg⋅dL⁻¹)
129.2 ± 36
152.5 ± 56
<0.001

Creatinine (mg⋅dL⁻¹)
1.17 ± 0.22
1.2 ± 0.92
0.4225

SBP
116 ± 14
121 ± 16
<0.001

DBP
77 ± 10
77.4 ± 10.6
0.102

Smoking, (%)
256 (25)
159 (25.3)
0.918

Continuous variables are presented as means ± SD and categorical variables as numbers (percentage). BMI, body mass index; HbA1c, glycated hemoglobin; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; SBP, systolic blood pressure; DBP: diastolic blood pressure.

3.2. Genotype Frequency in Diabetic and Obese Subjects

The distribution of the Native American (NAM), the European (EUR), and the African (AFR) individual admixture proportions was comparable between the diabetic subjects and controls (NAM: ; EUR: ; and AFR: ). The analysis did not include the individual ancestry proportions, which were not available for all the subjects.

Table 2 shows genotype and allelic frequency of the polymorphisms studied. Hardy Weinberg equilibrium (HWE) test showed no deviation inside the population. There was no association with diabetes and genotype frequency found in any of the polymorphisms studied (). However, taking into account that the polymorphisms studied have been associated with obesity, the genotype frequencies were analyzed according to obesity (Table 3). Lower frequency of the TT genotype in obese people was observed (4.7% versus 7.7%) which provides a significant protection against obesity (OR: 0.54, ). A marginal association, but not significant (), was found with rs2071749 polymorphism and obesity. Therefore, subanalyses were carried out according to genetic inheritance models and diabetes and obesity status, and stratification by gender.

Table 2: Genotype frequencies of the polymorphisms studied in diabetic patients and controls.

Gene/polymorphism
Genotypes alleles
Diabetes
Controls
OR (95% CI)

of HWE

TXNIP
rs7211
CC, (%)
345 (55.4)
528 (54.5)
Reference

0.880

CT, (%)
239 (38.3)
376 (38.8)
0.88 (0.6–1.3)
0.577

TT, (%)
39 (6.3)
65 (6.7)
0.9 (0.8–1.2)
0.776

NQO1
rs1800566
CC, (%)
216 (34.7)
327 (33)
Reference

0.406

CT, (%)
288 (46.2)
483 (48.6)
0.90 (0.7–1.1)
0.373

TT, (%)
119 (19.1)
183 (18.4)
0.98 (0.7–1.3)
0.915

HMOX1
rs2071749
AA, (%)
269 (43.8)
413 (43.2)
Reference

0.625

AG, (%)
267 (43.5)
417 (43.6)
1.05 (0.76–1.4)
0.757

GG, (%)
78 (12.7)
126 (13.2)
1.03 (0.74–1.4)
0.837

NRF2
rs2364723
CC, (%)
210 (33.6)
301 (30.3)
Reference

0.092

CG, (%)
286 (56.8)
471 (47.5)
0.87 (0.7–1.1)
0.236

GG, (%)
129 (20.6)
220 (22.2)
0.84 (0.63–1.1)
0.223

NRF2
rs6721961
CC, (%)
407 (65.3)
618 (62.5)
Reference

0.281

CA, (%)
189 (30.4)
317 (32)
0.9 (0.7–1.1)
0.374

AA, (%)
27 (4.3)
54 (5.5)
0.76 (0.5–1.2)
0.259

CI, confidence interval; HWE, Hardy-Weinberg equilibrium; HMOX1, heme oxygenase-1; NQO1, NAD(P)H quinone oxidoreductase 1; NRF2, Nuclear factor-erythroid 2- (NF-E2-) related factor 2; OR, odds ratio; and TXNIP, thioredoxin interacting protein.

Table 3: Genotype and allele frequencies of the polymorphisms studied in obese and nonobese subjects.

Gene/polymorphism
Genotype
Obesity
No obesity
OR (95% CI)

TRXNIP
rs7211
CC, n (%)
350 (56.6)
523 (53.7)
Reference

CT, n (%)
239 (38.7)
376 (38.6)
0.95 (0.7–1.2)
0.627

TT, n (%)
29 (4.7)
75 (7.7)
0.56 (0.35–0.87)
0.012

NQO1
rs1800566
CC, n (%)
212 (34.2)
331 (33.3)
Reference

CT, n (%)
302 (48.8)
469 (47)
1 (0.8–1.25)
0.963

TT, n (%)
106 (17)
196 (19.7)
0.84 (0.6–1.1)
0.257

HMOX1
rs2071749
AA, n (%)
261 (40)
419 (45.7)
Reference

AG, n (%)
306 (46.9)
378 (41.3)
1.3 (1–1.6)
0.019

GG, n (%)
85 (13.1)
119 (13)
1.1 (0.8–1.5)
1.14

NRF2
rs2364723
CC, n (%)
194 (31)
317 (32)
Reference

CG, n (%)
300 (48)
457 (46)
1.1 (0.85–1.35)
0.551

GG, n (%)
131 (21)
210 (22)
0.98 (0.7–1.3)
0.899

NRF2
rs6721961
CC, n (%)
390 (63.1)
635 (63.9)
Reference

CA, n (%)
195 (31.6)
311 (31.3)
1 (0.8–1.3)
0.853

AA, n (%)
33 (5.3)
48 (4.8)
1.1 (0.7–1.7)
0.631

CI, confidence interval; HMOX1, heme oxygenase-1; NQO1, NAD(P)H quinone oxidoreductase 1; NRF2, Nuclear factor-erythroid 2- (NF-E2-) related factor 2; OR, odds ratio; and TXNIP, thioredoxin interacting protein.

When CC + CT genotype of the rs7211 polymorphism was compared with TT, lower frequency of TT genotype associated with lower BMI (TT: 27.6 ± 4 versus CC + CT: 29 ± 5, ) and higher HDL-C levels (TT: 58 ± 13 ± 4 versus CC + CT: 49 ± 12, ) in obese subjects without diabetes was found. The OR revealed that TT genotype is a protective factor against obesity in nondiabetic individuals (Table 4). In the multivariable analysis, genotype (as recessive model) was negatively associated with BMI, as well as HDL-C levels in nondiabetic subjects (Table 5). However, these associations were not found in diabetic or diabetic obese subjects.

Table 4: Genotype frequency of the rs7211 polymorphism in subjects without diabetes and women.

Obese
Nonobese
Crude OR (95% CI)

Adjusted OR (95% CI)

Nondiabetic

    CC
189 (56.6)
339 (53.4)
Reference

Reference

    CT
133 (40)
241 (38)
0.98 (0.74–1.3)
0.903
1 (0.76–1.4)
0.863

    TT
55 (8.6)
12 (3.6)
0.3 (0.2–0.7)
0.004
0.3 (0.15–0.7)
0.003

    CC + CT
322 (96.4)
580 (91.5)
0.4 (0.2–0.76)
0.005
0.3 (0.15–0.7)
0.003

Women

    CC
197 (59)
259 (51)
Reference

Reference

    CT
118 (36)
203 (40)
0.7 (0.6–1)
0.072
0.9 (0.6–1.2)
0.418

    TT
17 (5)
47 (9)
0.5 (0.26–0.85)
0.013
0.5 (0.25–0.96)
0.04

    CT + TT
135 (41)
250 (49)
0.70 (0.5–0.9)
0.016
0.7 (0.5–0.96)
0.028

CI, confidence interval; OR, odds ratio.

Obesity in logistic regression was adjusted by age, gender (except in women model), glucose, triglycerides, LDL-C, and HDL-C levels.

Table 5: Multiple linear regression of BMI as dependent variable in nondiabetic subjects and women.

Nondiabetic subjects

    HDL-C

<0.001
0.11

    Glucose

<0.001

    rs7211 (CT + TT: 0, TT: 1)

0.033

Women

    HDL

0.003
0.048

    Glucose

0.001

    rs7211 (CC: 0, CT + TT: 1)

0.037

Models also include age, gender (in nondiabetic subjects model), LDL-C, and triglycerides. Significant variables were presented. BMI, body mass index; HDL-C, high-density lipoprotein cholesterol.

After gender stratification, CC genotype was higher in obese women compared with nonobese women (Table 4). TT genotype was also lower in obese women (TT frequency in obese women = 5.12 versus nonobese women = 9.23, ) but significance was achieved when multivariate analysis and logistic regression were carried out including genotype as dominant model. Results showed that carrying one copy of T allele decreases the risk of obesity in women (Table 4). However the association of BMI with genotype was weak and explains just 0.48% of BMI variability (Table 5).

The rs2071749 polymorphism in the HMOX1 gene was significantly associated with obesity when it was analyzed under a dominant model (Table 6). AA carriers showed lower BMI than AG + GG carriers and had lower risk of obesity. The association of the rs2071749 polymorphism persists with obesity, after being adjusted by age, gender, HDL-C levels, LDL-C levels, and triglycerides concentrations.

Table 6: Association of the rs2071749 polymorphism of the HMOX1 gene with obesity.

Variable
Obese, (%)
Nonobese, (%)

BMI

AA
261 (40)
410 (46)
0.026
28.9 ± 4.9
0.027

AG + GG
390 (60)
497 (54)

29.5 ± 5.3

Crude OR (95% CI)
1.25 (1.02–1.54)

0.026

OR adjusted (95% CI)
1.34 (1.06–1.7)

0.013

Beta coefficient/
0.063/0.03

0.023

BMI values are presented as mean ± SD. BMI, body mass index, CI: confidence interval, and OR: odds ratio.

Obesity in logistic regression was adjusted by age, gender, glucose, triglycerides, and LDL-C and HDL-C levels. Genotype was included as dominant model in logistic and multivariate analysis.

Linear regression of BMI as dependent variable was adjusted by age, gender, glucose, triglycerides, and LDL-C and HDL-C levels.

The rs6721961 polymorphism of the Nrf2 gene was associated with diabetes in men but not in women, after stratification by gender (Table 7). CC carriers had higher glucose levels in comparison with CA + CC carriers when genotype was compared as dominant model. Thus, the presence of the A allele represents a protective factor against diabetes in men. The pairwise linkage disequilibrium (LD) analysis of the rs6721961 and rs2364723 in NRF2 gene revealed a low LD between them (). Additionally, the analyses of potential haplotype effects showed no association with diabetes () or obesity ().

Table 7: Association of the rs6721961 polymorphism of the Nrf2 gene with diabetes in men.

Variable
Diabetic, (%)
Nondiabetic, (%)

Glucose

CC
203 (69)
271 (58)
0.011
126.6 ± 59
0.056

CA
79 (27)
171 (37)

116 ± 56

AA
12 (4)
23 (5)

105 ± 32

CA + AA
91 (31)
194 (42)
0.003
114.4 ± 53
0.031

Crude OR (95% CI)
0.62 (0.46–0.85)

0.003

OR adjusted (95% CI)
0.56 (0.38–0.82)

0.003

Crude beta coefficient
−0.078

0.031

Glucose values (mg⋅dL⁻¹) are presented as mean ± SD. CI: confidence interval; OR: odds ratio.

Diabetes in logistic regression was adjusted by age, triglycerides, and LDL-C and HDL-C levels. Genotype was included as dominant model in logistic and linear regression.

⁢ value when comparing each genotype.

⁢ value in comparison of CC with CA + AA.

4. Discussion

Oxidative stress is one of the major metabolic factors that lead to the onset of chronic disease like insulin resistance, hypertension, metabolic syndrome, prediabetes, diabetes, and its complications. It has been recognized that the antioxidant system is abated in diabetes and obesity, accompanied by an increased production of inflammatory cytokines. A defect in Nrf2 activation in many organs has been documented widely in experimental models of insulin resistance and diabetes leading to decreased expression of its target genes [31–34]. In this study, it was found that the rs6721961 (−617C/A) polymorphism of Nrf2 gene was associated with diabetes in Mexican mestizo men, while the rs7211 polymorphism of TXNIP gene and the rs2071749 polymorphism of HMOX1 gene were associated with obesity. Nevertheless, the polymorphisms rs1800566 (NQO1) and rs2364723 (Nrf2) were not found to be associated with diabetes or obesity.

Wang et al. [35] did not found association between the rs1800566 polymorphism in NQO1 gene and the risk of T2DM in Chinese population. On the other hand, Kim [36] reported that there were no associations between the rs1800566 polymorphism and BMI, blood pressure, lipid profile, HbA1c, postprandial glucose, and homeostasis model assessment-insulin resistance (HOMA-IR). Nevertheless, Martínez-Hernández et al. [18] showed that the T allele of the rs1800566 polymorphism was associated with increased triglycerides levels and decreased HDL-C levels in Mexican mestizo individuals with metabolic syndrome. However, in this study this association was not found and we suggest that this difference is due to the fact that 57% and 47% of the subjects recruited in our study were under medical treatment with fibrates and statins, respectively, and attend to clinical detection annually.

The rs7211 has been poorly studied in diabetic and obese subjects. First, van Greevenbroek et al. [23] found that triglycerides levels were higher in diabetic subjects with the T allele than the C allele carriers, which was associated with higher glucose concentrations. Subsequently, Ferreira et al. [24] found that the rs7211 was associated in Brazilian subjects with diabetes and hypertension, and T carriers showed higher concentrations of blood glucose and systolic blood pressure. Das et al. [37] reported that, in European-American or African-American subjects recruited in USA, the TXNIP gene expression was negatively correlated with obesity. In this study, the T allele was found as a protection factor against metabolic traits like BMI and HDL-C. The T allele carriers showed that lower BMI values and HDL-C concentrations in nondiabetic subjects and in women are still associated with obesity after adjustment by confounding factors. Some authors reported that this polymorphism can increase TXNIP expression [22], but this still remains unclear. However, this is the first report in which it was found that T allele is a protective factor against obesity and we suggest that TXNIP expression or function would be affected allowing TXN to exert its antioxidant action.

The SNP rs2071746 (T-413A) and microsatellite are polymorphisms in the promoter region of HMOX1 gene which can modulate its transcriptional activity and have been associated with various human diseases [19]. Lin et al. [20] showed that the repeat in the HMOX1 promoter was significantly associated with essential hypertension, systolic blood pressure, and diastolic blood pressure, whereas the other two SNPs rs2071746 and rs2071749 were not associated. Although the rs2071749 has been poorly studied with metabolic traits, in this study it was found that AA carriers showed lower BMI values which were associated with the risk of presenting obesity, whereas GG carriers showed higher BMI values. This is the first report showing the association of the rs2071749 polymorphism and obesity. This polymorphism is an intronic tag SNP and it has reported that it is in LD with the rs3761439, which is located in the promoter region of the HMOX1 gene. The consensus sequence of a nuclear factor κB (NF-κB) binding site to HMOX1 is altered by rs3761439 polymorphism that increases the risk of developing impaired lung function [38]. This suggests that the mechanisms by which this polymorphism could be involved in obesity may be secondary to increased oxidative stress and inflammation. The results in this work open the gate to new investigations about the role of the HMOX1 gene in Mexican obese patients.

NRF2 gene regulates the enzymes studied here. Many SNPs have been identified in this gene. The rs35652124 and rs6721961 SNPs are predicted to affect Nrf2 myeloid zinc finger 1 (MZF1) and antioxidant response elements- (ARE-) like promoter binding sites, respectively. These SNPs affect the efficient binding of proteins such as Nrf2 to the MZF1 and ARE-like promoter binding sites [39].

In metabolic diseases, the rs6721961 polymorphism has been associated with blood pressure in Japanese subjects [39]. Wang et al. [7] showed that the rs6721961 polymorphism in the NRF2 gene was significantly associated with oxidative stress, antioxidant status, and risk of newly diagnosed T2DM, as well as with impaired insulin secretory capacity and increased insulin resistance in T2DM patients of a Chinese population. Individuals with the CC genotype had lower total antioxidant capacity, glutathione levels and superoxide dismutase, catalase, and glutathione peroxidase activities as well as lower homeostasis model assessment of β-cell function index (HOMA-β) in comparison with individuals with the CC genotype. Those with the AA genotype also had a higher malondialdehyde concentration and HOMA-IR index values. The frequency of allele A was significantly higher in T2DM subjects (29.4%). Individuals with the AA genotype had a significantly higher risk of developing T2DM, relative to those with the CC genotype, even after adjusting for known T2DM risk factors. However, this investigation showed that the AA carriers had lower glucose concentrations and the OR revealed that A carriers had lower risk of developing diabetes in men after stratification by gender. We had the limitation that Nrf2 gene expression was not determined in these patients and it is essential to elucidate the mechanisms by which the rs6721961 polymorphism can exert its protective effect in our population. Larger sample size could support this.

We know that this research has some limitations. The gene expressions and activities were not determined and just candidate genes were included. Dietary information was not available in all patients. However, this study opens the gate to new researches about the role of HMOX1, TXNIP, and Nrf2 in obesity, diabetes, and worse metabolic traits.

5. Conclusions

This study shows that the rs7211 polymorphism in the TXNIP gene and the rs2071749 of the HMOX1 gene are associated with obesity in Mexican mestizo people. In this sense, rs2071749 increases the risk of having obesity and the rs7211 polymorphism may be a protective factor to develop obesity and to have lower concentrations of HDL-C in Mexican mestizo women. In addition, rs6721961 gen polymorphism in Nrf2 was negatively associated with diabetes in men. All of these results open the doors for new research, taking into account the effect on metabolic traits, and may be useful tools to design new therapeutic strategies in obesity and diabetes type 2.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to acknowledge the Medical Family Clinics from ISSSTE: Ignacio Chávez, Del Valle, Xochimilco, and Fuentes Brotantes. This work was financially supported by CONACyT-Mexico (Grant no. SALUD-2013-01-201519).

References

International Diabetes Federation, Diabetes Atlas, International Diabetes Federation, Brussels, Belgium, 7th edition, 2015, http://www.diabetesatlas.org.

R. H. Eckel, S. E. Kahn, E. Ferrannini et al., “Obesity and type 2 diabetes: what can be unified and what needs to be individualized?” Diabetes Care, vol. 34, no. 6, pp. 1424–1430, 2011.

S. Tyrovolas, A. Koyanagi, N. Garin et al., “Diabetes mellitus and its association with central obesity and disability among older adults: a global perspective,” Experimental Gerontology, vol. 64, pp. 70–77, 2015.

M. Serrano-Rios, “Relationship between obesity and the increased risk of major complications in non-insulin-dependent diabetes mellitus,” European Journal of Clinical Investigation, vol. 28, supplement 2, pp. 14–18, 1998.

S.-I. Yamagishi, S. Maeda, T. Matsui, S. Ueda, K. Fukami, and S. Okuda, “Role of advanced glycation end products (AGEs) and oxidative stress in vascular complications in diabetes,” Biochimica et Biophysica Acta—General Subjects, vol. 1820, no. 5, pp. 663–671, 2012.

Y.-J. Surh, J. K. Kundu, and H.-K. Na, “Nrf2 as a master redox switch in turning on the cellular signaling involved in the induction of cytoprotective genes by some chemopreventive phytochemicals,” Planta Medica, vol. 74, no. 13, pp. 1526–1539, 2008.

X. Wang, H. Chen, J. Liu et al., “Association between the NF-E2 related factor 2 gene polymorphism and oxidative stress, anti-oxidative status, and newly-diagnosed type 2 diabetes mellitus in a Chinese population,” International Journal of Molecular Sciences, vol. 16, no. 7, pp. 16483–16496, 2015.

A. S. Jiménez-Osorio, A. Picazo, S. González-Reyes, D. Barrera-Oviedo, M. E. Rodríguez-Arellano, and J. Pedraza-Chaverri, “Nrf2 and redox status in prediabetic and diabetic patients,” International Journal of Molecular Sciences, vol. 15, no. 11, pp. 20290–20305, 2014.

J. Pi, L. Leung, P. Xue et al., “Deficiency in the nuclear factor E2-related factor-2 transcription factor results in impaired adipogenesis and protects against diet-induced obesity,” The Journal of Biological Chemistry, vol. 285, no. 12, pp. 9292–9300, 2010.

H.-Y. Cho, “Genomic structure and variation of nuclear factor (erythroid-derived 2)-like 2,” Oxidative Medicine and Cellular Longevity, vol. 2013, Article ID 286524, 24 pages, 2013.

Y. Shimoyama, Y. Mitsuda, N. Hamajima, and T. Niwa, “Polymorphisms of Nrf2, an antioxidative gene, are associated with blood pressure in Japanese,” Nagoya Journal of Medical Science, vol. 76, no. 1-2, pp. 113–120, 2014.

S. M. Figarska, J. M. Vonk, and H. M. Boezen, “NFE2L2 polymorphisms, mortality, and metabolism in the general population,” Physiological Genomics, vol. 46, no. 12, pp. 411–417, 2014.

D. Ross and D. Siegel, “NAD(P)H:quinone oxidoreductase 1 (NQO1, DT-diaphorase), functions and pharmacogenetics,” Methods in Enzymology, vol. 382, pp. 115–144, 2004.

J.-E. Chung, B. C. Chang, K. E. Lee, J. H. Kim, and H. S. Gwak, “Effects of NAD(P)H quinone oxidoreductase 1 polymorphisms on stable warfarin doses in Korean patients with mechanical cardiac valves,” European Journal of Clinical Pharmacology, vol. 71, no. 10, pp. 1229–1236, 2015.

H.-Y. Shyu, C.-S. Fong, Y.-P. Fu et al., “Genotype polymorphisms of GGCX, NQO1, and VKORC1 genes associated with risk susceptibility in patients with large-artery atherosclerotic stroke,” Clinica Chimica Acta, vol. 411, no. 11-12, pp. 840–845, 2010.

S. J. Han, E. S. Kang, H. J. Kim et al., “The C609T variant of NQO1 is associated with carotid artery plaques in patients with type 2 diabetes,” Molecular Genetics and Metabolism, vol. 97, no. 1, pp. 85–90, 2009.

T. Ramprasath, P. S. Murugan, E. Kalaiarasan, P. Gomathi, A. Rathinavel, and G. S. Selvam, “Genetic association of Glutathione peroxidase-1 (GPx-1) and NAD(P)H:quinone oxidoreductase 1(NQO1) variants and their association of CAD in patients with type-2 diabetes,” Molecular and Cellular Biochemistry, vol. 361, no. 1-2, pp. 143–150, 2012.

A. Martínez-Hernández, E. J. Córdova, O. Rosillo-Salazar et al., “Association of HMOX1 and NQO1 polymorphisms with metabolic syndrome components,” PLoS ONE, vol. 10, no. 5, Article ID e0123313, 2015.

E. Y. Lee, Y. H. Lee, S. H. Kim et al., “Association between heme oxygenase-1 promoter polymorphisms and the development of albuminuria in type 2 diabetes: a case-control study,” Medicine, vol. 94, no. 43, article e1825, 2015.

R. Lin, W. Fu, W. Zhou et al., “Association of heme oxygenase-1 gene polymorphisms with essential hypertension and blood pressure in the Chinese Han population,” Genetic Testing and Molecular Biomarkers, vol. 15, no. 1-2, pp. 23–28, 2011.

E. Junn, S. H. Han, J. Y. Im et al., “Vitamin D3 up-regulated protein 1 mediates oxidative stress via suppressing the thioredoxin function,” The Journal of Immunology, vol. 164, no. 12, pp. 6287–6295, 2000.

X. He and Q. Ma, “Redox regulation by nuclear factor erythroid 2-related factor 2: gatekeeping for the basal and diabetes-induced expression of thioredoxin-interacting protein,” Molecular Pharmacology, vol. 82, no. 5, pp. 887–897, 2012.

M. M. J. van Greevenbroek, V. M. M.-J. Vermeulen, E. J. M. Feskens et al., “Genetic variation in thioredoxin interacting protein (TXNIP) is associated with hypertriglyceridaemia and blood pressure in diabetes mellitus,” Diabetic Medicine, vol. 24, no. 5, pp. 498–504, 2007.

N. E. Ferreira, S. Omae, A. Pereira et al., “Thioredoxin interacting protein genetic variation is associated with diabetes and hypertension in the Brazilian general population,” Atherosclerosis, vol. 221, no. 1, pp. 131–136, 2012.

J. Grimberg, S. Nawoschik, L. Belluscio, R. McKee, A. Turck, and A. Eisenberg, “A simple and efficient non-organic procedure for the isolation of genomic DNA from blood,” Nucleic Acids Research, vol. 17, no. 20, p. 8390, 1989.

J. Zúñiga, N. Yu, R. Barquera et al., “HLA class I and class II conserved extended haplotypes and their fragments or blocks in mexicans: implications for the study of genetic diversity in admixed populations,” PLoS ONE, vol. 8, no. 9, Article ID e74442, 2013.

P. Sanz, V. Prieto, I. Flores, Y. Torres, M. López-Soto, and M. J. Farfán, “Population data of 13 STRS in southern Spain (Andalusia),” Forensic Science International, vol. 119, no. 1, pp. 113–115, 2001.

P. Calzada, I. Suárez, S. García et al., “The Fang population of Equatorial Guinea characterised by 15 STR-PCR polymorphisms,” International Journal of Legal Medicine, vol. 119, no. 2, pp. 107–110, 2005.

C. Barrot, C. Sánchez, M. Ortega et al., “Characterisation of three Amerindian populations from Hidalgo State (Mexico) by 15 STR-PCR polymorphisms,” International Journal of Legal Medicine, vol. 119, no. 2, pp. 111–115, 2005.

A. González-Martín, A. Gorostiza, H. Rangel-Villalobos et al., “Analyzing the genetic structure of the Tepehua in relation to other neighbouring Mesoamerican populations. A study based on allele frequencies of STR markers,” American Journal of Human Biology, vol. 20, no. 5, pp. 605–613, 2008.

Y. Tan, T. Ichikawa, J. Li et al., “Diabetic downregulation of Nrf2 activity via ERK contributes to oxidative stress-induced insulin resistance in cardiac cells in vitro and in vivo,” Diabetes, vol. 60, no. 2, pp. 625–633, 2011.

H.-A. Seo and I.-K. Lee, “The role of NRF2: adipocyte differentiation, obesity, and insulin resistance,” Oxidative Medicine and Cellular Longevity, vol. 2013, Article ID 184598, 7 pages, 2013.

J. Xu, A. C. Donepudi, V. R. More et al., “Deficiency in Nrf2 transcription factor decreases adipose tissue mass and hepatic lipid accumulation in leptin-deficient mice,” Obesity, vol. 23, no. 2, pp. 335–344, 2015.

A. Uruno, Y. Yagishita, and M. Yamamoto, “The Keap1-Nrf2 system and diabetes mellitus,” Archives of Biochemistry and Biophysics, vol. 566, pp. 76–84, 2015.

G. Wang, L. Zhang, and Q. Li, “Genetic polymorphisms of GSTT1, GSTM1, and NQO1 genes and diabetes mellitus risk in Chinese population,” Biochemical and Biophysical Research Communications, vol. 341, no. 2, pp. 310–313, 2006.

D. Kim, “An association between 609 C→T Polymorphism in NAD(P)H:quinone oxidoreductase 1 (NQO1) gene and blood glucose levels in Korean population,” Korean Diabetes Journal, vol. 33, no. 1, pp. 24–30, 2009.

S. K. Das, N. K. Sharma, S. J. Hasstedt et al., “An integrative genomics approach identifies activation of thioredoxin/thioredoxin reductase-1-mediated oxidative stress defense pathway and inhibition of angiogenesis in obese nondiabetic human subjects,” The Journal of Clinical Endocrinology & Metabolism, vol. 96, no. 8, pp. E1308–E1313, 2011.

G. Tanaka, F. Aminuddin, L. Akhabir et al., “Effect of heme oxygenase-1 polymorphisms on lung function and gene expression,” BMC Medical Genetics, vol. 12, article 117, 2011.

Y. Shimoyama, Y. Mitsuda, Y. Tsuruta, N. Hamajima, and T. Niwa, “Polymorphism of Nrf2, an antioxidative gene, is associated with blood pressure and cardiovascular mortality in hemodialysis patients,” International Journal of Medical Sciences, vol. 11, no. 7, pp. 726–731, 2014.