IRS1, TCF7L2, ADRB1, PPARG, and HHEX Polymorphisms Associated with Atherogenic Risk in Mexican Population
Objective. We aimed to explore the association between polymorphisms of IRS1 (rs1801278), TCF7L2 (rs7903146 and rs12255372), ADRB1 (rs1801253), PPARG (rs1801282), and HHEX (rs5015480) genes with atherogenic risk (AI = Total cholesterol/HDL) in MetS, T2D, and healthy populations from the Mexican Social Security Institute. Methodology and Results. Four hundred thirty-five MetS, 517 T2D, and 547 healthy individuals were selected. The association between the SNPs and the atherogenic index was evaluated by multiple linear regression and multinomial logistic regression models. The ADRB1 gene showed a statistically significant association with high-risk atherogenic index, (IC 95% 1.64–5.24; ) for the Arg/Gly variant, under the dominant model an (IC 95% 1.67–5.25; ), and under the Log additive model an (IC 95% 1.54–4.15; ). Conclusions. The Arg389Gly polymorphism of the ADRB1 gene may be a worthy biological marker to predict the risk of developing cardiovascular diseases given a high-risk atherogenic index.
The metabolic syndrome (MetS) is a group of metabolic disorders that directly promote the development of cardiovascular diseases (CVD) [1, 2]. One of these disorders, insulin resistance (IR) , is at the core of the physiopathology of the MetS and type 2 diabetes (T2D). Lipid metabolism is affected by the consequences of IR, culminating in the MetS characteristic lipid disorders: IR leads to high levels of very low density lipoproteins (VLDL) that cause high triglyceride levels and low HDL-c levels, a profile known as atherogenic lipoprotein phenotype . IR can, thus, be associated with fatty acid metabolism disequilibrium .
The last two decades have seen a marked increase in the number of people with MetS worldwide . In Mexico, the prevalence of MetS under the ATP III, AHA/NHLBI, and IDF criteria was 36.8%, 41.6%, and 49.8%, respectively .
Several susceptibility loci for MetS component traits have been soundly identified to play a role in the risk of IR , obesity [9, 10], hypertriglyceridemia, and low HDL-c levels [11, 12]. The IRS1 gene Gly972Arg polymorphism affects the PI3-k/AKT/GSK3 signaling pathway, producing IR due to a decrease in the ability of insulin to transport glucose, translocation of the glucose transporter, and glycogen synthesis . This polymorphism has been associated with elevated serum triglyceride, low HDL-c levels, high free fatty acid levels, and elevated systolic blood pressure . Variants in the TCF7L2, SLC30A8, HHEX, FTO, PPARG, and KCNJ11 genes are associated with T2D and cardiovascular disease risk [15–17]. Polymorphisms in the beta-adrenergic receptors genes are associated with cardiovascular and metabolic phenotypes. The adrenergic receptors play a major role in energy expenditure through thermogenesis stimulation . The aim of this study was to analyze the association between the IRS1 (rs1801278), TCF7L2 (rs7903146 and rs12255372), ADRB1 (rs1801253), PPARG (rs1801282), and HHEX (rs5015480) polymorphisms with IR, lipid profile, and atherogenic risk.
2.1. Design and Study Population
In a cross-sectional study, we recruited 431 patients with MetS and 517 patients with T2D defined according to the criteria of the American Heart Association/National Heart, Lung, and Blood Institute Scientific Statement (AHA/NHLBI)  and the American Diabetes Association (ADA), respectively, and 547 health individuals (without MetS or T2D) between 35 and 65 years old. The AHA/NHLBI criteria are: waist circumference ≥102 cm for men and ≥88 cm for women, fasting glucose ≥100 mg/dL, systolic blood pressure (SBP) ≥130 mm Hg, diastolic blood pressure (DBP) ≥85 mm Hg, HDL-c <40 mg/dL for men and <50 mg/dL for women, and triglycerides ≥150 mg/dL. The study population was selected at the IMSS Medical Center of Mexico City. The inclusion criterion for healthy individuals was absence of a family history of T2D among parents, siblings, or offspring. Written consent was obtained from the participants, and the protocol was approved by the Institutional Review Board (IRB) of the National Public Health Institute and the National Committee of the IMSS.
2.2. Biochemical and Anthropometric Measurements
All participants answered a questionnaire used for collecting socioeconomic data, as well as family history of diabetes and other diseases. The following anthropometric measurements were performed: weight (kg), height (cm), body mass index (BMI) (kg/m2), waist (cm) and hip (cm) circumference, and the waist/hip ratio (WHR). SBP and DBP were measured with a mercury column sphygmomanometer (American Diagnostic Inc., New York), after the patients had rested and were seated for at least 5 minutes. A blood sample was obtained from all participants, after a 12-hour fast, for the determination of biochemical profiles of glucose (mg/dL), total cholesterol, LDL cholesterol (LDL-C, mg/dL), HDL cholesterol (HDL-C, mg/dL), and triglycerides (mg/dL). These parameters were determined by the ILab 30 Clinical Chemistry System instrument (Instrumentation Laboratory, Barcelona, Spain). IR was estimated by means of the Homeostasis Model Assessment-Insulin Resistance, mmol/pmol (HOMA-IR), and the atherogenic index (AI) was estimated with the following formula: AI = Total cholesterol/HDL-c, with a cut-off point for cardiovascular risk of >4.5 for women and >5 for men .
DNA extraction from peripheral vein blood samples was carried out by a semiautomated method with the commercial QIAamp DNA Blood Mini/Kit column purification kits (Qiagen, Germany). DNA integrity was determined by 0.8% agarose gel electrophoresis, stained with ethidium bromide, and visualized through the Gel Doc 2000 equipment (BIORADCA, USA). DNA concentration was determined by spectrophotometry, using the VICTOR3 1420 equipment (Perkin-Elmer, Germany). Single nucleotide polymorphism (SNP) analysis was carried out in duplicate with real-time PCR TaqMan technology (Applied Biosystems 7900HT equipment, Foster City, California, USA), following standardized techniques. Studied SNPs were: IRS1 (rs1801278), TCF7L2 (rs7903146 and rs12255372), ADRB1 (rs1801253), PPARG (rs1801282), and HHEX (rs5015480) (Table 1). For quality control (QC) of the genotyping, we used a call rate of 0.99 for cases and controls; all the samples with 0.99 call rate were analyzed by two independent researchers who scored genotype results using the graphical view (SD 2.1.1, Applied Biosystems), and all differences were resolved by discussions between the independent researchers and by rerunning the samples if necessary. If the concordance was less than 99.99%, this was considered as a missing value. Only 3.35% of the samples had less than 99.99% concordance.
In order to control the possible effect of population stratification, a panel of ancestry informative markers (AIM) was also genotyped. The panel was integrated by 27 AIMS (rs2814778, rs723822, rs1008984, rs1435090, rs17203, rs768324, rs719776, rs1112828, rs1403454, rs3340, rs2077681, rs1935946, rs1320892, rs1373302, rs2695, rs1980888, rs1327805, rs2207782, rs1487214, rs2078588, rs724729, rs292932, rs1369290, rs386569, rs718092, rs878825, and rs16383) taken from a previous assessed ancestry panel . These markers showed differences in allele frequency between European 30%, African 5%, and Amerindian populations 65%, in Mexican population in a previous report of our research group .
2.4. Statistical Analysis
In order to compare the relevant variants between cases and controls, the Kruskal-Wallis test for group comparison was used for continuous variables, and the chi square test () was used for categorical data. Hardy-Weinberg equilibrium was evaluated for each variant under study. A multiple linear regression analysis was performed to evaluate the effects of the selected variants on the IR biochemical markers (glucose, insulin, and HOMA-IR) and the atherogenic profile biomarkers (total cholesterol, LDL-c, HDL-c, and triglycerides). A multinomial logistic regression was also carried out to evaluate the effect of each and every variant on AI risk (AI = Total cholesterol/HDL; with the following cut-off points: low risk <4.5 for women and <5 for men, medium risk 4.5–7 for women and 5–9 for men, and high risk >7 for women and >9 for men)  for the different genotypes in the codominant, dominant, recessive, and log-additive heritability models. All models were adjusted for age, sex, disease status, BMI, and ancestry percentage. The study’s power average was 80%. The possibility of multiple testing burden was avoided by using the Bonferroni method (0.05/42 comparisons = for linear regression and 0.05/12 comparisons = for multinomial logistic regression) and false-positive report probability (FPRP) for correction of multiple comparisons . The STATA SE version 12.0 statistical package software (College Station, TX) was used for the analyses.
The general characteristics of the groups studied are shown in Table 2. The Low risk AI group had lower levels in mg/dL to total cholesterol (190.93 ± 39.92, ), LDL (119.21 ± 30.88, ), triglycerides (158.95 ± 77.54, ), glucose (118.41 ± 56.90, ), and insulin (11.73 ± 8.89, ), compared to medium and high risk AI groups. There was no difference in BMI between the three groups. The ancestry percentage was similar in the three groups, with the Amerindian group having the highest percentage, followed by the European group, whereas the African group had a minimum percentage.
Table 3 shows that all studied polymorphisms are in Hardy-Weinberg equilibrium. For all polymorphisms, genotype and allele frequencies were similar in the three groups of risk AI. All polymorphisms were in HW equilibrium in the control group of low risk AI (IRS1 , TCF7L2 rs12255372 , rs7903146 , PPARG , ADRB1 , and HHEX ).
In order to evaluate the variant’s effect on the IR biochemical and on atherogenic profile biomarkers, measured by total cholesterol, LDL-c, HDL-c, triglycerides, glucose, insulin, and HOMA-IR, a multiple linear regression analysis was carried out according to the codominant, dominant, and recessive models, adjusted for age, sex, BMI, ancestry percentage, and disease status. No statistically significant associations were found (data not shown).
A possible association between the polymorphic variants and other characteristics such as SBP, DBP, and WHR, was assessed through a multiple linear regression analysis adjusted for age, sex, BMI, ancestry percentage, and disease status. However, the results of this analysis did not show any potential association (data not shown).
Table 4 shows multinomial logistic regression results that evaluates the effect of each and every variant on the risk of presenting an AI category (medium or high risk; low risk was considered the reference category) for the log-additive model, adjusted for age, sex, BMI, disease status, and ancestry percentage; only the ADRB1 polymorphism (Arg389Gly) showed a statistically significant association with high-risk AI (OR= 2.52, IC 95% 1.54–4.15; ).
Finally, the multinomial logistic regression results that evaluate the effect ADRB1 variant on the risk of presenting an AI category (medium or high risk; low risk was considered the reference category), for the codominant, dominant, and recessive models, adjusted for age, sex, BMI, disease status, and ancestry percentage, showed a statistically significant association of the ADRB1 gene with high-risk AI, (OR = 2.94, IC 95% 1.64–5.24; ) for the Arg/Gly variant, for the dominant model an OR = 2.96 (IC 95% 1.67–5.25; ) (Table 5).
4. Discussion and Conclusions
Identification of the genes involved in the metabolic alterations produced in complex diseases like MetS should be a priority to study. We include in the study SNPs of genes related to glucose metabolism [22, 23], fatty acid metabolism, and development of an atherogenic profile in MetS patients.
The main findings of this study are that the ADRB1 gene heterozygote and dominant model showed an association with high risk AI, with an OR of 2.94 (IC 95% 1.64–5.24; ) for the Arg/Gly variant, for the dominant model an OR of 2.96 (IC 95% 1.67–5.25; ), and for the Log additive model an OR of 2.52 (IC 95% 1.54–4.15; ).
The ADRB1 gene association is in agreement with reports that associate this gene variant with fat accumulation, BMI, and obesity . The significant association of Arg389Gly variant means that carriers of these polymorphisms are almost three times more likely to present a high-risk AI than noncarriers. These associations are independent of age, sex, BMI, and ancestry percentage. These observations are supported by studies that associated the variant with a poor coupling of the receptor with the G-protein, which translates to functional consequences . Associations between the ADRB1 gene Arg389Gly polymorphism and heart diseases and an increase in fat accumulation have also been found .
The HHEX T/C heterozygote and the dominant model showed an association although not significant with an increase in total cholesterol levels and in LDL-c levels, suggesting that carriers of this variant may be prone to develop lipid metabolism alterations. This is consistent with recent results for the Mexican population . The HHEX gene regulates pancreatic development by controlling endodermic cell proliferation and the positioning of the cells that will later turn into pancreatic tissue. This probably occurs through the effects of HHEX on the transcription and translation of the cell proliferator. HHEX expression remains activated in the adult pancreas, which suggests that HHEX may have additional functions in differentiated pancreatic cells. Because the HHEX gene plays a role in pancreas development, as described above, it may be involved in glucose metabolism . Some HHEX gene variants have been strongly associated with T2D and insulin-stimulated glucose release .
The results of the evaluation of the effects of the polymorphic variants on the IR biochemical markers (glucose, insulin, and HOMA-IR) showed that the IRS1 Gly/Arg variant is associated although not significant with an increase in HOMA-IR, consistentl with previous reports of an association between the IRS1 gene Gly972Arg polymorphism and IR, and with a higher risk of T2D and CVD in various populations, including the Mexican group, where this association was replicated [27–29]. Some studies suggest that this polymorphism may contribute to IR development due to impairment in the insulin ability to activate the IRS1/PI3-kinase/Akt/GSK-3 signaling pathway, causing deficient glucose transportation, translocation of glucose transporters, and glycogen synthesis .
We found no significant association with PPARγ Pro12Ala polymorphism; this is contrary to previous studies which report that Ala12 homozygotes have an increased risk of cardiovascular disease (CAD), but a meta-analysis conducted to investigate the relationship between the PPARγ polymorphism and the potential risk of CAD reported that this risk is higher in Caucasian populations compared to a Asiatic population, suggesting that the risk for cardiovascular disease conferred by the Pro12Ala polymorphism differs across ethnic populations .
Even TCF7L2 may contribute not only to the development of T2DM but also to the development of CAD because its important role in vascular remodeling through the regulation of smooth muscle cell proliferation and endothelial cell growth [31, 32]; in our study, we did not find significant association. Our result is similar to a recent report by the Atherosclerosis Risk in Communities (ARIC) study where they did not find significant associations between TCF7L2 SNPs and incident vascular events .
The results of this assessment of the effects of the polymorphic variants on the medium or high risk of atherogenesis confirm the association between the ADRB1 gene Arg389Gly . Our results suggest that the ADRB1 gene Arg389Gly polymorphism can be a good biological marker to predict the risk of cardiovascular diseases development given a high-risk AI. This is in agreement with the notion that the genetic susceptibility to atherosclerosis is due in part to known metabolic risk factors (in this case dyslipidemia) . Further studies of this genetic effect are required to evaluate its potential interaction with other risk factors involved in the same metabolic pathway, since the known genetic risk factors do not fully explain the hereditary propensity to atherosclerosis .
|AIM:||Ancestry informative markers|
|BMI:||Body mass index|
|DBP:||Diastolic blood pressure|
|HDL:||High density lipoproteins|
|HOMA-IR:||Homeostasis model assessment-insulin resistance|
|LDL:||Low density lipoproteins|
|SBP:||Systolic blood pressure|
|T2D:||Type 2 diabetes|
|VLDL:||Very low density lipoproteins.|
Conflict of Interests
There is no conflict of interests that could be perceived as prejudicing the impartiality of the research reported.
B. I. Estrada-Velasco and A. I. Burguete-García have equal contribution.
This work was supported by a Grant from CONACYT SALUD (2005-C02-14412). Proyectos Estratégicos; IMSS-2004/014, 2005/2/I/363, and 2006/1B/1/002, Fundación IMSS A. C, Fundación Gonzalo Rio Arronte I. A. P. Mexico, and the “Vinculación del Cinvestav con el Sector Salud” Support Program, 2006. M. Cruz is a recipient of a Fundación IMSS Scholarship, Mexico.
B. Isomaa, P. Almgren, T. Tuomi et al., “Cardiovascular morbidity and mortality associated with the metabolic syndrome,” Diabetes Care, vol. 24, no. 4, pp. 683–689, 2001.View at: Google Scholar
R. Rojas, C. A. Aguilar-Salinas, A. Jiménez-Corona et al., “Metabolic syndrome in Mexican adults. Results from the National Health and Nutrition Survey 2006,” Salud Publica de Mexico, vol. 52, supplement 1, pp. S11–S18, 2010.View at: Google Scholar
S. Canizales-Quinteros, C. A. Aguilar-Salinas, M. G. Ortiz-López et al., “Association of PPARG2 Pro12Ala variant with larger body mass index in Mestizo and Amerindian populations of Mexico,” Human Biology, vol. 79, no. 1, pp. 111–119, 2007.View at: Google Scholar
K. Kristiansson, M. Perola, E. Tikkanen et al., “Genome-wide screen for metabolic syndrome susceptibility loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits,” Circulation Cardiovascular Genetics, vol. 5, no. 2, pp. 242–249, 2012.View at: Publisher Site | Google Scholar
M. A. Marini, S. Frontoni, D. Mineo et al., “The Arg972 variant in insulin receptor substrate-1 is associated with an atherogenic profile in offspring of type 2 diabetic patients,” Journal of Clinical Endocrinology and Metabolism, vol. 88, no. 7, pp. 3368–3371, 2003.View at: Publisher Site | Google Scholar
Diabetes Genetics Initiative of Broad Institute of Harvard and MIT, Lund University, and Novartis Institutes for BioMedical Research, “Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels,” Science, vol. 316, article 5829, pp. 1331–1336, 2007.View at: Publisher Site | Google Scholar
S. Wacholder, S. Chanock, M. Garcia-Closas, L. El Ghormli, and N. Rothman, “Assessing the probability that a positive report is false: an approach for Molecular Epidemiology Studies,” Journal of the National Cancer Institute, vol. 96, no. 6, pp. 434–442, 2004.View at: Google Scholar
M. Cruz, J. García-Mena, E. López et al., “Genes candidatos como posibles marcadores de susceptibilidad a Diabetes tipo 2,” REB, vol. 24, pp. 81–86, 2005.View at: Google Scholar
M. Cruz, A. Valladares, J. Garcia-Mena et al., “Candidate Gene Association Study conditioning on individual ancestry in type 2 diabetes and metabolic syndrome patients from Mexico City,” Diabetes/Metabolism Research and Reviews, vol. 26, no. 4, pp. 261–270, 2010.View at: Publisher Site | Google Scholar
S. Podlowski, K. Wenzel, H. P. Luther et al., “β1-adrenoceptor gene variations: a role in idiopathic dilated cardiomyopathy?” Journal of Molecular Medicine, vol. 78, no. 2, pp. 87–93, 2000.View at: Google Scholar
H. Staiger, A. Stančáková, J. Zilinskaite et al., “A candidate type 2 diabetes polymorphism near the HHEX locus affects acute glucose-stimulated insulin release in European populations results from the EUGENE2 Study,” Diabetes, vol. 57, no. 2, pp. 514–517, 2008.View at: Publisher Site | Google Scholar
A. I. Burguete-Garcia, M. Cruz-Lopez, V. Madrid-Marina et al., “Association of Gly972Arg polymorphism of IRS1 gene with type 2 diabetes mellitus in lean participants of a national health survey in Mexico: a Candidate Gene Study,” Metabolism Clinical and Experimental, vol. 59, no. 1, pp. 38–45, 2010.View at: Publisher Site | Google Scholar
L. E. Martínez-Gómez, M. Cruz, G. A. Martínez-Nava 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,” Annals of Human Genetics, vol. 75, no. 5, pp. 612–620, 2011.View at: Publisher Site | Google Scholar
Z. Wu, Y. Lou, W. Jin, Y. Liu, L. Lu, and G. lu, “The Pro12Ala polymorphism in the peroxisome proliferator-activated receptor gamma-2 gene (PPARc2) is associated with increased risk of coronary artery disease: a meta-analysis,” PLoS ONE, vol. 7, no. 12, Article ID e53105, 2012.View at: Google Scholar
S. J. Bielinski, J. S. Pankow, A. R. Folsom, K. E. North, and E. Boerwinkle, “TCF7L2 single nucleotide polymorphisms, cardiovascular disease and all-cause mortality: the Atherosclerosis Risk in Communities (ARIC) Study,” Diabetologia, vol. 51, no. 6, pp. 968–970, 2008.View at: Publisher Site | Google Scholar