Skip to main content
  • Research article
  • Open access
  • Published:

UCP2 and UCP3 variants and gene-environment interaction associated with prediabetes and T2DM in a rural population: a case control study in China

Abstract

Background

There are disparities for the association between uncoupling proteins (UCP) and type 2 diabetes (T2DM). The study was to examine the associations of genetic variants of UCP2 and UCP3 with prediabetes and T2DM in a rural Chinese population.

Methods

A population-based case-control study of 397 adults with T2DM, 394 with prediabetes and 409 with normal glucose tolerance (NGT) was carried out in 2014 in a rural community in eastern China. Three groups were identified through a community survey and the prediabetes and NGT groups were frequently matched by age and gender with the T2DM group and they were not relatives of T2DM subjects. With r2 ≥ 0.8 and minor allele frequency (MAF) ≥0.05 for tag single nucleotide polymorphisms (SNPs) with potential function, three (rs660339, rs45560234 and rs643064) and six (rs7930460, rs15763, rs647126, rs1800849, rs3781907 and rs1685356) SNPs were selected respectively for UCP2 and UCP3 and genotyped in real time using the MassARRAY system (Sequenom; USA). The haplotypes, gene-environmental interaction and association between genetic variants of UCP2 and UCP3 and prediabetes or T2DM were explored.

Results

There were no significant differences in age and sex among three study groups. After the adjustment for possible covariates, the A allele of rs1800849 in UCP3 was significantly associated with prediabetes (aORAA vs GG = 1.68, 95% CI: 1.02–2.78), and the association was also significant under the recessive model (aOR AA vs GA + GG = 1.64, 95% CI: 1.02–2.66). Also, rs15763 was found to be marginally significantly associated with T2DM under dominant model (ORGA + AA vs GG = 0.73, 95% CI: 0.52–1.03, P = 0.072). No haplotype was significantly associated with prediabetes or T2DM. Multiplicative interactions for rs660339-overweight on T2DM were observed. In addition, the AA genotype of rs660339 was associated with an increased risk of T2DM in overweight subjects (OR = 1.48, 95%CI: 0.87–2.52) but with a decreased risk in those with normal weight (OR = 0.54, 95%CI: 0.28–1.05).

Conclusions

Rs1800849 in UCP3 was significantly associated with prediabetes. Overweight might modify the effects of rs660339 of UCP2 on T2DM.

Peer Review reports

Background

Type 2 diabetes (T2DM) is characterized by chronic hyperglycemia caused by a combination of insulin resistance and inadequate compensatory insulin secretion. A recent study found that 11.6% of Chinese adults had diabetes and 50.1% had prediabetes, a precursor to T2DM [1]. T2DM is resulted from both genetic and environmental factors [2].

The uncoupling proteins (UCP) 2 and 3, members of the mitochondrial inner membrane carrier family [3], have been frequently studied for their roles in the development of T2DM. Although UCP2 and UCP3 share structural similarities, they have different tissue expressions in mammals [4]. UCP2 is distributed across a wide range of tissue and cell types, whereas UCP3 is mainly restricted to skeletal muscle and brown adipose tissue of humans [5,6,7]. UCP2 is an important regulator of mitochondrial oxidative stress damage induced by reactive oxygen species, and has been implicated with T2DM [8]. UCP3 is implicated in fatty acid (FA) metabolism, and patients with diabetes have low rather than high UCP3 protein content [9, 10]. But recent studies, including genome wide association studies (GWAS), have shown no evidence for any tagging-SNPs of UCP2-UCP3 tested that are associated with T2DM risk [11, 12], indicating disparities of association between UCP2 or UCP3 and T2DM. Although no associations were found between UCP2/UCP3 and diabetes in GWAS, some studies referred that peroxisome proliferator-activated receptor gamma (PPARG, related to diabetes confirmed in GWAS) interacted with Forkhead box protein O1 (FOXO1) to regulate the expression of UCP2 and rs659366 of UCP2 is unlikely to be identified through a GWAS strategy because it occurs predominantly in already diabetic subjects [12,13,14]. Little is known with regard to the association between UCP2/UCP3 polymorphisms and T2DM or prediabetes in Chinese population.

In recent years, some researchers have focused on the interaction between genetic and environmental factors (such as smoking, alcohol intake, physical activity level, general obesity, abdominal obesity, hypertension, dyslipidemia) for T2DM [15,16,17]. The interrelationship among adiposity, hypertension and UCP is likely complex. Several studies have found that UCP2 polymorphisms are associated with obesity [18,19,20,21] and hypertension [22]. Obesity and hypertension are associated with the development of insulin resistance, linked to the pathogenesis of T2DM [23, 24]. The exploration for the interactions between UCP polymorphisms and environmental factors can help for a better understanding of the association between UCP variants and T2DM.

In the current study, we used bioinformatics analyses to identify and select a set of variants in the UCP2/UCP3 region to determine the association between UCP2/UCP3 and prediabetes or T2DM in a Chinese population. We also examine the interactions of UCP2/UCP3 with adiposity and hypertension for prediabetes and T2DM.

Methods

Study design and subjects

The current study included 397 adults with T2DM, 394 with prediabetes and 409 with normal glucose tolerance (NGT), which was carried out in 2014 in a rural community in eastern China. The three groups of subjects were selected from participants of a cross-sectional study that was conducted during the period from June to December in 2012 in Yuhuan County, an island located in east China. That cross-sectional study was the baseline survey of the Yuhuan Rural Health Cohort Study, with a total 28,155 participants aged 20 years or older (born before 31 December, 1992), and collected information on demography, lifestyle (smoking, alcohol use, regular physical exercise) with a standard questionnaire and performed a free physical examination for those fasting for over 8 h. Venous blood of 5 ml were collected for fast plasma glucose (FPG) test, oral glucose tolerance test (OGTT) and blood lipids tests. Additional blood specimens were preserved under − 20°C. The Institutional Review Board of the Fudan University School of Public Health approved the study and all participants gave a written informed consent. The groups of newly diagnosed prediabetes and NGT subjects were frequently matched in age and sex with the T2DM group by excluding relatives. Subjects with T2DM were defined as those who had FPG ≥ 7.0 mmol/L or 2-h postload plasma glucose (2hPG) ≥ 11.1 mmol/L or having medical treatment history; subjects with prediabetes were defined as those who had 6.1 mmol/L ≤ FPG < 7.0 mmol/L or 7.8 mmol/L ≤ 2hPG < 11.1 mmol/L; NGT was defined as FPG < 6.1 mmol/L and 2hPG < 7.8 mmol/L [25]. Overweight was defined as body mass index (BMI) ≥ 25 kg/m2 [26].

SNP selection and function prediction

A set of SNPs with potential function were selected using SNP Function Prediction of SNPinfo Web Server. According to HapMap dbSNP (CHB/JPT population) and NCBI dbSNP (Asian), a chromosomal from 3000 bp upstream to 3000 bp downstream of UCP2 and UCP3 was searched, and 12 SNPs (rs45560234, rs632862, rs659366, rs660339, rs15763, rs1800006, rs1800849, rs2075577, rs590336, rs647126, rs7926165, rs7930460) were obtained with potential function and minor allele frequency (MAF) > 0.05. Furtherly, linkage disequilibrium (LD) TAG SNP Selection was used to select haplotype-tagging SNPs (tag SNPs). With the same searched chromosomal range and NCBI dbSNP (Asian), r2 ≥ 0.8, MAF ≥ 0.05 for SNPs having potential function, rs660339 and rs1800849 as the mandatory entry site due to function prediction and previous studies [14,15,16], three SNPs of UCP2 (rs660339, rs45560234 and rs643064) and six SNPs of UCP3 (rs7930460, rs15763, rs647126, rs1800849, rs3781907 and rs1685356) were selected for this study (Additional file 1: Table S1).

Genotyping

Genomic DNA used for polymerase chain reaction (PCR) amplification was extracted from peripheral blood using a DNA extraction kit (TIANGEN, China) according to the manufacturer’s instructions. Primers for amplification were designed with AssayDesigner 3.1 software. In the primer extension, ddH2O, 10X Buffer, 25 mM dNTP, 25 mM MgCl2, 0.5 μM Primer, 5 U/ul Hostar Taq and DNA template were used in the PCR reaction system. After purifying the products and transferring to SpectroCHIP, MALDI-TOF mass spectrometry was used for SNP genotyping. Thermocycling was carried out under the following conditions: Initial denaturation, 95 °C for 2 min followed by 45 cycles at 95 °C for 30 s (denaturation), 56 °C for 30 s (annealing), and 72 °C for 1 min (extension), with a final extension step at 72 °C for 5 min. The sequences of primers for amplification of the nine SNPs of UCP2–3 genes were presented in Additional file 2: Table S2. The MAFs of rs660339, rs45560234, rs643064, rs7930460, rs15763, rs647126, rs1800849, rs3781907 and rs1685356 in NGT group were 0.399 (A), 0.048 (A), 0.276 (T), 0.293 (G), 0.222 (A), 0.426 (G), 0.298 (A), 0.426 (G) and 0.392 (T), respectively.

Statistical analysis

Data were described as mean ± standard deviation for continuous variables, and as frequencies (percentage) for non-continuous variables. Comparisons among the three groups were performed using one-way analyses of variance for continuous variables and χ2 test for categorical variables. Deviation from the Hardy-Weinberg equilibrium was assessed by χ2 test for the NGT group. Multinomial logistic regression analysis was performed to calculate crude odds ratios (ORs) and adjusted ORs (aORs) and their 95% confidential intervals (CIs). The genetic models tested in the study included additive, recessive and dominant models. Frequencies of haplotypes were calculated and compared among the three groups. Finally, we employed the general multifactor dimensionality reduction (GMDR) method [27] and logistic regression to estimate the interaction between gene and environmental risk factors.

Multinomial logistic regression, χ2 test and one-way analyses of variance were performed using SPSS version 19 software (Chicago, Illinois, USA). Haploview 4.2 and Hapstat3.0 were used to analyze linkage disequilibrium (LD) and haplotypes. GMDR software were used to estimate the interaction between gene and environmental factors. All tests were two-tailed, and it was considered statistically significant if a p value < 0.05 and marginally significant if a p value between 0.05 and 0.10. The Bonferroni for the number of SNPs was used to adjust the significance level of multiple comparisons when analyzing the association between genotypes and prediabetes or T2DM.

Results

Characteristics of the study population

The present study investigated 397 with newly diagnosed T2DM, 394 with newly diagnosed prediabetes and 409 who had NGT and three groups were matched by age and sex. Table 1 showed that those with T2DM or prediabetes were more likely to be hypertensive and overweight, and had lower income compared with those with NGT. No significant differences among three groups were found for age, sex, regular exercise, smoking and excessive drinking.

Table 1 General characteristics of subjects with normal glucose tolerance (NGT), prediabetes and type 2 diabetes (T2DM)

Associations of UCP2/UCP3 with prediabetes and T2DM

Table 2 and Additional file 3: Table S3 presented the distributions of three genotypes for each SNP and their associations with prediabetes or diabetes. All the genotype frequencies met the assumption of Hardy-Weinberg equilibrium for the NGT group (p > 0.05) except for rs15763. With respect to rs1800849 of UCP3, the genotype of AA (AA vs GG, p = 0.026; ORAA vs GG = 1.75) was found to be a significant risk factor for prediabetes and after adjustment for age, sex, regular exercise, smoking, excessive drinking, hypertension and overweight, AA was still associated with an increased risk of prediabetes (aORAA vs GG = 1.68, 95% CI: 1.02–2.78); The recessive model also revealed a detrimental effect of AA for prediabetes (ORAA vs GA + GG = 1.72, 95% CI: 1.07–2.76; aOR AA vs GA + GG = 1.64, 95% CI: 1.02–2.66). In the recessive model, rs660339 (ORAA vs GA + GG = 1.39, 95% CI: 0.96–2.00) and rs7930460 (ORGG vs AG + AA = 1.52, 95%CI: 0.94–2.47) were marginally significantly associated with prediabetes, and A of rs660339 and G of rs7930460 tended to increase the risk of prediabetes. Additionally, rs15763 was found to be marginally significantly associated with T2DM in the dominant model (ORGA + AA vs GG = 0.73, 95%CI: 0.52–1.03, P = 0.072) after adjusting for covariates, indicating that A could be a protective factor for T2DM (Table 2).

Table 2 Potential prediabetes or T2DM-associated UCP2 and UCP3 polymorphisms

Linkage disequilibrium and haplotype analysis

There was a strong LD between rs643064 and rs660339 (r2 = 0.57) of UCP2 as well as between rs7930460 and rs647126 (r2 = 0.57), and between rs647126 and rs1685356 (r2 = 0.46–0.47). Combining the location information of SNPs, two haplotype blocks, that one included two SNPs (rs643064 and rs660339) of UCP2 and another included three SNPs (rs7930460, rs647126 and rs1685356) of UCP3, were detected. The frequencies of the haplotypes CG and CA formed by rs6413064-rs660339 of UCP2, AAT and AGC formed by rs7930460-rs647126-rs1685356 of UCP3 were lower, while the frequencies of TA formed by rs6413064-rs660339 of UCP2, GGC and AAC formed by rs7930460-rs647126-rs1685356 of UCP3 were higher in the prediabetes and T2DM groups compared to the NGT group. However, these haplotypes were not significantly associated with either prediabetes or T2DM in the regression analysis (Table 3).

Table 3 Associations of haplotypes with prediabetes and type 2 diabetes (T2DM)

Gene-environmental interaction for prediabetes or T2DM

GMDR analysis was conducted to assess the interaction between the 9 SNPs and environment factors (income, regular exercise, smoking, excessive drinking, hypertension and overweight). Our GMDR analysis demonstrated that overweight was associated with prediabetes (testing accuracy = 56.1%, cross-validation consistency = 10/10, P = 0.011) and T2DM (testing accuracy = 60.9%, cross-validation consistency = 10/10, P = 0.001) in the one-factor model. There were potential interactions of overweight-rs1800849 on prediabetes (testing accuracy = 52.7%, cross-validation consistency = 10/10, P = 0.011) and of hypertension-overweight-rs660339 on T2DM (testing accuracy = 61.1%, cross-validation consistency = 10/10, P = 0.011) (Additional file 4: Table S4). We further tested the potential gene-environment interactions by logistic regression to validate the result of GMDR. Because in most cases recessive models showed a better fit for the association, we selected the recessive model in a further interaction analysis. We found that subjects with both overweight and AA of rs1800849 had an increased risk of prediabetes (OR = 3.92, 95% CI: 1.63–9.39), compared with an OR of 1.73 for those only with overweight and of 1.52 for those with only AA genotype, and the subjects with hypertension and overweight and AA of rs660339 also had an increased T2DM risk (OR = 7.78, 95% CI: 2.83–21.40) (Tables 4 and 5). Multiplicative interactions among rs660339, overweight and hypertension in association with T2DM were estimated in logistic regression analysis and only the interaction between rs660339 and overweight was significant (Table 5). Further stratified analysis revealed that the AA genotype of rs660339 was associated with an increased risk of T2DM for overweight subjects (OR = 1.48, 95% CI: 0.87–2.52) but with a decreased risk for those with normal weight (OR = 0.54, 95% CI: 0.28–1.05).

Table 4 Interaction analysis for rs1800849 and overweight on prediabetes
Table 5 Interaction analysis for rs660339 and hypertension and overweight on T2DM

Discussion

In this study, three SNPs of UCP2 and six SNPs of UCP3 were genotyped and the associations between these SNPs and T2DM or prediabetes were evaluated in a rural Chinese population. The findings suggested a significant association between UCP3 rs1800849 and prediabetes in this population. No statistical significance was identified for the association between rs1800849 and T2DM. In previous studies, T allele was associated a lower incidence of T2DM in a French cohort [28], and rs1800849 was associated with elevated HDL-C levels and reduced BMI, independent of modifiable factors [29]. The discrepancies among studies may be due to ethnic differences in the genetic variation.

A number of previous studies have studied rs660339 and rs659366 of UCP2 associated with diabetes and the results are contradictory. G > A variation of rs660339 was found to be significantly associated with an increased risk of T2DM in Koreans [30], but with a decreased risk in Indians [31], while no significant association was found in Iranians and Europeans [32, 33]. de Souza BM et al. performed a comprehensive meta-analysis of 5 studies with a total of 2112 cases and 1841 controls, and demonstrated that rs660339 was significantly associated with an increased risk of T2DM only in Asian population [34]. No statistically significant association between rs660339 and prediabetes or T2DM was found in our study, which was consistent with results from another study of Chinese population [35].

Also, the interactions between hypertension, overweight and rs660339 and between overweight and rs1800849 were observed for T2DM and prediabetes, respectively. Logistic regression confirmed the multiplicative interaction of overweight-rs660339, and rs660339 might play a different role in populations with normal weight and overweight. UCP2 and UCP3 can play important roles in both fatty acid and blood sugar metabolism [8,9,10]. Our findings revealed that the AA genotype of rs660339 was a risk factor of T2DM for overweight subjects (OR = 1.48) but a protective factor for those with normal weight (OR = 0.54). Both UCP2 and overweight/obesity are associated with insulin resistance [21, 36, 37], a fundamental aspect of the etiology of T2DM. It is unknown however, why AA genotype of rs660339 influences T2DM differently in overweight and normal weight people. Previous studies demonstrated that TT of rs660339 increased the risk of overweight [20, 38], indicating that rs660339 could accelerate the development of prediabetes or T2DM through overweight. Another widely studied SNP of UCP2 was rs659366. Previous studies found that rs659366 was associated with decreased insulin sensitivity and obesity [21], and had a strong linkage disequilibrium with rs660339 with D’ = 0.946 and r2 = 0.657 [39], which were not observed in this study. Contradictory results had been reported concerning the association between overweight and rs1800849 [40,41,42].

To our knowledge, no studies have been published about the associations between rs647126 and rs1685356 and prediabetes or T2DM although GG of rs647126 and TT of rs1685356 were reported to be correlated with high BMI [43]. In previous studies of Fins and Caucasians [44, 45], rs15763 met the Hardy-Weinberg equilibrium, but not in the current population.

The current study was based on individuals with newly diagnosed T2DM and prediabetes. Potential gene-environmental interactions were evaluated using GMDR and confirmed by logistic regression. However, there were also some limitations. Firstly, the selection of SNPs was functional tagSNP strategy, and more studies were needed to confirm whether an association was due to SNPs themselves or the variants which were high LD with them. Secondly, genetic variants associated with prediabetes or T2DM were no longer significant after a Bonferroni correction (P < 0.05/9). Although Bonferroni correction was too strict generally, a larger size sample was needed to confirm the reliability of the associations. Additionally, covariates including lifestyle factors, age and sex had been adjusted but known genetic variants such as TCF7L2 and KCNJ11 had not been tested and controlled in this study.

Conclusion

In conclusion, rs1800849 in UCP3 were significantly associated with prediabetes in a rural Chinese population. Overweight modified the effect of rs660339 of UCP2 on T2DM. These findings suggested that rs1800849 in UCP3 and rs660339 in UCP2 might play an important role in the incidence and development of T2DM.

Abbreviations

2hPG:

2-h postload plasma glucose

BMI:

Body mass index

FPG:

Fast plasma glucose

GMDR:

General multifactor dimensionality reduction

GWAS:

Genome wide association study

LD:

Linkage disequilibrium

NGT:

Normal glucose tolerance

OGTT:

Oral glucose tolerance test

PCR:

Polymerase chain reaction

SNPs:

Single nucleotide polymorphisms

T2DM:

Type 2 diabetes mellitus

UCP:

Uncoupling protein

References

  1. Xu Y, Wang L, He J, Bi Y, Li M, Wang T, et al. Prevalence and control of diabetes in Chinese adults. JAMA. 2013;310(9):948–59.

    Article  CAS  PubMed  Google Scholar 

  2. Samsom M, Trivedi T, Orekoya O, Vyas S. Understanding the importance of gene and environment in the etiology and prevention of type 2 diabetes mellitus in high-risk populations. Oral Health Case Rep. 2016;2(1):112.

    PubMed  PubMed Central  Google Scholar 

  3. Diehl AM, Hoek JB. Mitochondrial uncoupling: role of uncoupling protein anion carriers and relationship to thermogenesis and weight control “the benefits of losing control”. J Bioenerg Biomembr. 1999;31(5):493–506.

    Article  CAS  PubMed  Google Scholar 

  4. Erlanson-Albertsson C. Uncoupling proteins--a new family of proteins with unknown function. Nutr Neurosci. 2002;5(1):1–11.

    Article  CAS  PubMed  Google Scholar 

  5. Souza BM, Assmann TS, Kliemann LM, Gross JL, Canani LH, Crispim D. The role of uncoupling protein 2 (UCP2) on the development of type 2 diabetes mellitus and its chronic complications. Arq Bras Endocrinol Metabol. 2011;55(4):239–48.

    Article  PubMed  Google Scholar 

  6. Dalgaard LT, Pedersen O. Uncoupling proteins: functional characteristics and role in the pathogenesis of obesity and type II diabetes. Diabetologia. 2001;44(8):946–65.

    Article  CAS  PubMed  Google Scholar 

  7. Boss O, Samec S, Paoloni-Giacobino A, Rossier C, Dulloo A, Seydoux J, et al. Uncoupling protein-3: a new member of the mitochondrial carrier family with tissue-specific expression. FEBS Lett. 1997;408(1):39–42.

    Article  CAS  PubMed  Google Scholar 

  8. Casteilla L, Rigoulet M, Pénicaud L. Mitochondrial ROS metabolism: modulation by uncoupling proteins. IUBMB Life. 2001;52(3–5):181–8.

    Article  CAS  PubMed  Google Scholar 

  9. Schrauwen P, Hoeks J, Hesselink MK. Putative function and physiological relevance of the mitochondrial uncoupling protein-3: involvement in fatty acid metabolism? Prog Lipid Res. 2006;45(1):17–41.

    Article  CAS  PubMed  Google Scholar 

  10. Schrauwen P, Mensink M, Schaart G, Moonen-Kornips E, Sels JP, Blaak EE, et al. Reduced skeletal muscle uncoupling protein-3 content in prediabetic subjects and type 2 diabetic patients: restoration by rosiglitazone treatment. J Clin Endocrinol Metab. 2006;91(4):1520–5.

    Article  CAS  PubMed  Google Scholar 

  11. Zee RY, Ridker PM, Chasman DI. Mitochondrial uncoupling protein gene cluster variation (UCP2-UCP3) and the risk of incident type 2 diabetes mellitus: the women’s genome health study. Atherosclerosis. 2011;214(1):107–9.

    Article  CAS  PubMed  Google Scholar 

  12. GWAS catalog. http://www.ebi.ac.uk/gwas/. Accessed 2 Mar 2018.

  13. Dalgaard LT. Genetic variance in uncoupling protein 2 in relation to obesity, type 2 diabetes, and related metabolic traits: focus on the functional −866G>a promoter variant (rs659366). J Obes. 2011;2011:340241.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Wu C, Gong Y, Yuan J, Gong H, Zou Y, Ge J. Identification of shared genetic susceptibility locus for coronary artery disease, type 2 diabetes and obesity: a meta-analysis of genome-wide studies. Cardiovasc Diabetol. 2012;11:68.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zhang L, Wang J, Zhang M, Wang G, Shen Y, Wu D, et al. Association of type 2 diabetes mellitus with the interaction between low-density lipoprotein receptor-related protein 5 (LRP5) polymorphisms and overweight and obesity in rural Chinese adults. J Diabetes. 2017;9(11):994–1002.

    Article  CAS  PubMed  Google Scholar 

  16. Perry JR, Voight BF, Yengo L, Amin N, Dupuis J, Ganser M, et al. Stratifying type 2 diabetes cases by BMI identifies genetic risk variants in LAMA1 and enrichment for risk variants in lean compared to obese cases. PLoS Genet. 2012;8(5):e1002741.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Liu L, Chen L, Li Z, Li L, Wang M, Qu J, et al. Association of genetic variants in TOMM7 gene and gene environment interaction with type 2 diabetes in Chinese Dong population. Zhong Nan Da Xue Xue Bao Yi Xue Ban. 2015;40(1):31–8.

    PubMed  Google Scholar 

  18. Gul A, Ateş Ö, Özer S, Kasap T, Ensari E, Demir O, et al. Role of the polymorphisms of uncoupling protein genes in childhood obesity and their association with obesity-related disturbances. Genet Test Mol Biomarkers. 2017;21(9):531–8.

    CAS  PubMed  Google Scholar 

  19. Brondani LA, de Souza BM, Assmann TS, Bouças AP, Bauer AC, Canani LH, et al. Association of the UCP polymorphisms with susceptibility to obesity: case-control study and meta-analysis. Mol Biol Rep. 2014;41(8):5053–67.

    Article  CAS  PubMed  Google Scholar 

  20. Zhang M, Wang M, Zhao ZT. Uncoupling protein 2 gene polymorphisms in association with overweight and obesity susceptibility: a meta-analysis. Meta Gene. 2014;2:143–59.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Andersen G, Dalgaard LT, Justesen JM, Anthonsen S, Nielsen T, Thørner LW, et al. The frequent UCP2-866G>A polymorphism protects against insulin resistance and is associated with obesity: a study of obesity and related metabolic traits among 17 636 Danes. Int J Obes. 2013;37(2):175–81.

    Article  CAS  Google Scholar 

  22. Ji Q, Ikegami H, Fujisawa T, Kawabata Y, Ono M, Nishino M, et al. A common polymorphism of uncoupling protein 2 gene is associated with hypertension. J Hypertens. 2004;22(1):97–102.

    Article  CAS  PubMed  Google Scholar 

  23. Kahn BB, Flier JS. Obesity and insulin resistance. J Clin Invest. 2000;106(4):473–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Reaven GM. Pathophysiology of insulin resistance in human disease. Physiol Rev. 1995;75:473–86.

    Article  CAS  PubMed  Google Scholar 

  25. Alberti KG, Zimmet PZ. Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus. Provisional report of a WHO consultation. Diabet Med. 1998;15(7):539–53.

    Article  CAS  PubMed  Google Scholar 

  26. American Medical Association. Executive summary of the clinical guidelines on the identification, evaluation, and treatment of overweight and obesity in adults. Arch Intern Med. 1998;158(17):1855–67.

    Article  Google Scholar 

  27. Xu HM, Xu LF, Hou TT, Luo LF, Chen GB, Sun XW, et al. GMDR: versatile software for detecting gene-gene and gene-environment interactions underlying complex traits. Curr Genomics. 2016;17(5):396–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Meirhaeghe A, Amouyel P, Helbecque N, Cottel D, Otabe S, Froguel P, et al. An uncoupling protein 3 gene polymorphism associated with a lower risk of developing type II diabetes and with atherogenic lipid profile in a French cohort. Diabetologia. 2000;43(11):1424–8.

    Article  CAS  PubMed  Google Scholar 

  29. Hamada T, Kotani K, Fujiwara S, Sano Y, Domichi M, Tsuzaki K, et al. The common-55 C/T polymorphism in the promoter region of the uncoupling protein 3 gene reduces prevalence of obesity and elevates serum high-density lipoprotein cholesterol levels in the general Japanese population. Metabolism. 2008;57(3):410–5.

    Article  CAS  PubMed  Google Scholar 

  30. Lee HJ, Ryu HJ, Shin HD, Park BL, Kim JY, Cho YM, et al. Associations between polymorphisms in the mitochondrial uncoupling proteins (UCPs) with T2DM. Clin Chim Acta. 2008;398(1–2):27–33.

    Article  CAS  PubMed  Google Scholar 

  31. Vimaleswaran KS, Radha V, Ghosh S, Majumder PP, Sathyanarayana Rao MR, Mohan V, et al. Uncoupling protein 2 and 3 gene polymorphisms and their association with type 2 diabetes in Asian Indians. Diabetes Technol Ther. 2011;13(1):19–25.

    Article  CAS  PubMed  Google Scholar 

  32. Heidari J, Akrami SM, Heshmat R, Amiri P, Fakhrzadeh H, Pajouhi M. Association study of the-866G/A UCP2 gene promoter polymorphism with type 2 diabetes and obesity in a Tehran population: a case control study. Arch Iran Med. 2010;13(5):384–90.

    CAS  PubMed  Google Scholar 

  33. Bielinski SJ, Pankow JS, Boerwinkle E, Bray MS, Kao WH, Folsom AR. Lack of association between uncoupling protein-2 Ala55Val polymorphism and incident diabetes in the atherosclerosis risk in communities study. Acta Diabetol. 2008;45(3):179–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. de Souza BM, Brondani LA, Boucas AP, Sortica DA, Kramer CK, Canani LH, et al. Associations between UCP1-3826A/G, UCP2-866G/A, Ala55Val and Ins/Del, and UCP3-55C/T polymorphisms and susceptibility to type 2 diabetes mellitus: case-control study and meta-analysis. PLoS One. 2013;8(1):e54259.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Shen Y, Wen Z, Wang N, Zheng Z, Liu K, Xia X, et al. Investigation of variants in UCP2 in Chinese type 2 diabetes and diabetic retinopathy. PLoS One. 2014;9(11):e112670.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chan CB, De Leo D, Joseph JW, McQuaid TS, Ha XF, Xu F, et al. Increased uncoupling protein-2 levels in beta-cells are associated with impaired glucose-stimulated insulin secretion mechanism of action. Diabetes. 2001;50(6):1302–10.

    Article  CAS  PubMed  Google Scholar 

  37. Kim JH, Li L, Yun JH, Choi BC, Baek KH. Association study between the -866G/A polymorphism in the promoter of uncoupling protein-2 gene and polycystic ovary syndrome. Mol Med Rep. 2011;4:747–51.

    CAS  PubMed  Google Scholar 

  38. Martinez-Hervas S, Mansego ML, de Marco G, Martinez F, Alonso MP, Morcillo S, et al. Polymorphisms of the UCP2 gene are associated with body fat distribution and risk of abdominal obesity in Spanish population. Eur J Clin Investig. 2012;42(2):171–8.

    Article  CAS  Google Scholar 

  39. Oktavianthi S, Trimarsanto H, Febinia CA, Suastika K, Saraswati MR, Dwipayana P, et al. Uncoupling protein 2 gene polymorphisms are associated with obesity. Cardiovasc Diabetol. 2012;11:41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. de Luis DA, Aller R, Izaola O, Romero E. Effect of -55CT polymorphism of UCP3 on insulin resistance and cardiovascular risk factors after a high protein/low carbohydrate versus a standard hypocaloric diet. Ann Nutr Metab. 2016;68(3):157–63.

    Article  CAS  PubMed  Google Scholar 

  41. Xu YP, Liang L, Wang CL, Fu JF, Liu PN, Lv LQ, et al. Association between UCP3 gene polymorphisms and nonalcoholic fatty liver disease in Chinese children. World J Gastroenterol. 2013;19(35):5897–903.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Zou Z, Mao L, Shi Y, Chen J, Wang L, Cai W. Association between peroxisome proliferator-activated receptor, UCP3 and lipoprotein lipase gene polymorphisms and obesity in Chinese adolescents. Obes Res Clin Pract. 2017;11(1):27–33.

    Article  PubMed  Google Scholar 

  43. van Abeelen AF, de Krom M, Hendriks J, Grobbee DE, Adan RA, van der Schouw YT. Variations in the uncoupling protein-3 gene are associated with specific obesity phenotypes. Eur J Endocrinol. 2008;158(5):669–76.

    Article  CAS  PubMed  Google Scholar 

  44. Salopuro T, Pulkkinen L, Lindström J, Kolehmainen M, Tolppanen AM, Eriksson JG, et al. Variation in the UCP2 and UCP3 genes associates with abdominal obesity and serum lipids: the finnish diabetes prevention study. BMC Med Genet. 2009;10:94.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Acosta A, Camilleri M, Shin A, Vazquez-Roque MI, Iturrino J, Lanza IR, et al. Association of UCP-3 rs1626521 with obesity and stomach functions in humans. Obesity (Silver Spring). 2015;23(4):898–906.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge all the staff who participated in this project, including health workers in Yuhuan County Center for Disease Control and Prevention, and from Fudan School of Public Health, and from local community-based health facilities.

Funding

This work was supported by the National Natural Science Foundation of China (grant number 81473038), Taizhou City Science Program (grant number 1401KY49), Shanghai 3-Year Public Health Action Plan (grant number GWTD2015S04) and Shanghai Leading Academic Discipline Project of Public Health (grant number 15GWZK0801). The sponsors have no role in the study design, survey process, data analysis and manuscript preparation.

Availability of data and materials

The datasets generated and analyzed in the current study are available from the corresponding author on the reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

MS and XC analyzed the data, interpreted the analyses results and wrote the draft of the manuscript. YC, QJ and CF designed the study and revised the manuscript. CW, TX and NW analyzed the data and collected data. MS, ST and XY collected data and supervised the study. All authors had read and approved the final manuscript.

Corresponding author

Correspondence to Chaowei Fu.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Institutional Review Board of the Fudan School of Public Health. Written informed consent was obtained from all participants before the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Main information of tSNPs and their function prediction. (DOCX 17 kb)

Additional file 2:

Table S2. Amplification and extension primers sequences of the nine loci in UCP2–3 genes. (DOCX 16 kb)

Additional file 3:

Table S3. Associations between UCP2, UCP3 and prediabetes or T2DM. (DOCX 31 kb)

Additional file 4:

Table S4. Best gene-environment interaction models identified by GMDR. (DOCX 16 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Su, M., Chen, X., Chen, Y. et al. UCP2 and UCP3 variants and gene-environment interaction associated with prediabetes and T2DM in a rural population: a case control study in China. BMC Med Genet 19, 43 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-018-0554-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-018-0554-4

Keywords