Novel epigenetic determinants of type 2 diabetes in Mexican-American families

Epigenetic regulation of glucose metabolism

PURPOSE OF REVIEW

Glucose metabolism is a central process in mammalian energy homeostasis. Its deregulation is a key factor in development of metabolic disease like diabetes and cancer. In recent decades, our understanding of gene regulation at the signaling, chromatin and posttranscriptional levels has seen dramatic developments.

RECENT FINDINGS

A number of epigenetic mechanisms that do not affect the genetic code can be assessed with new technologies. However, increasing complexity becomes a major challenge for translation into clinical application.

SUMMARY

The current review provides an update of transcriptional control of glucose metabolism, focusing on epigenetic regulators, DNA-methylation, histone modifications and noncoding RNAs. Recent studies heavily support the importance of those mechanisms for future therapeutics and preventive efforts for metabolic diseases.

Role of DNA Methylation in Type 2 Diabetes Etiology: Using Genotype as a Causal Anchor

Several studies have investigated the relationship between genetic variation and DNA methylation with respect to type 2 diabetes, but it is unknown if DNA methylation is a mediator in the disease pathway or if it is altered in response to disease state. This study uses genotypic information as a causal anchor to help decipher the likely role of DNA methylation measured in peripheral blood in the etiology of type 2 diabetes. Illumina HumanMethylation450 BeadChip data were generated on 1,018 young individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. In stage 1, 118 unique associations between published type 2 diabetes single nucleotide polymorphisms (SNPs) and genome-wide methylation (methylation quantitative trait loci [mQTLs]) were identified. In stage 2, a further 226 mQTLs were identified between 202 additional independent non-type 2 diabetes SNPs and CpGs identified in stage 1. Where possible, associations were replicated in independent cohorts of similar age. We discovered that around half of known type 2 diabetes SNPs are associated with variation in DNA methylation and postulated that methylation could either be on a causal pathway to future disease or could be a noncausal biomarker. For one locus (KCNQ1), we were able to provide further evidence that methylation is likely to be on the causal pathway to disease in later life.

Epigenetic Age Acceleration Assessed with Human White-Matter Images

The accurate estimation of age using methylation data has proved a useful and heritable biomarker, with acceleration in epigenetic age predicting a number of age-related phenotypes. Measures of white matter integrity in the brain are also heritable and highly sensitive to both normal and pathological aging processes across adulthood. We consider the phenotypic and genetic interrelationships between epigenetic age acceleration and white matter integrity in humans. Our goal was to investigate processes that underlie interindividual variability in age-related changes in the brain. Using blood taken from a Mexican-American extended pedigree sample (n = 628; age = 23.28-93.11 years), epigenetic age was estimated using the method developed by Horvath (2013). For n = 376 individuals, diffusion tensor imaging scans were also available. The interrelationship between epigenetic age acceleration and global white matter integrity was investigated with variance decomposition methods. To test for neuroanatomical specificity, 16 specific tracts were additionally considered. We observed negative phenotypic correlations between epigenetic age acceleration and global white matter tract integrity (ρ_pheno = -0.119, p = 0.028), with evidence of shared genetic (ρ_gene = -0.463, p = 0.013) but not environmental influences. Negative phenotypic and genetic correlations with age acceleration were also seen for a number of specific white matter tracts, along with additional negative phenotypic correlations between granulocyte abundance and white matter integrity. These findings (i.e., increased acceleration in epigenetic age in peripheral blood correlates with reduced white matter integrity in the brain and shares common genetic influences) provide a window into the neurobiology of aging processes within the brain and a potential biomarker of normal and pathological brain aging.SIGNIFICANCE STATEMENT Epigenetic measures can be used to predict age with a high degree of accuracy and so capture acceleration in biological age, relative to chronological age. The white matter tracts within the brain are also highly sensitive to aging processes. We show that increased biological aging (measured using epigenetic data from blood samples) is correlated with reduced integrity of white matter tracts within the human brain (measured using diffusion tensor imaging) with data from a large sample of Mexican-American families. Given the family design of the sample, we are also able to demonstrate that epigenetic aging and white matter tract integrity also share common genetic influences. Therefore, epigenetic age may be a potential, and accessible, biomarker of brain aging.

Epigenome-wide DNA methylation study of IgE concentration in relation to self-reported allergies

AIM

Epigenetic mechanisms are critical for normal immune development and epigenetic alterations might therefore be possible contributors to immune diseases. To investigate if DNA methylation in whole blood is associated with total and allergen-specific IgE levels.

METHODS

We performed an epigenome-wide association study to investigate the association between DNA methylation and IgE level, allergen-specific IgE and self-reported immune diseases and allergies in 728 individuals.

RESULTS

We identified and replicated 15 CpG sites associated with IgE, mapping to biologically relevant genes, including ACOT7, ILR5A, KCNH2, PRG2 and EPX. A total of 331 loci were associated with allergen-specific IgE, but none of these CpG sites were associated with self-reported allergies and immune diseases.

CONCLUSION

This study shows that IgE levels are associated with DNA methylation levels at numerous CpG sites, which might provide new leads for investigating the links between IgE and allergic inflammation.

Hypermethylation in the promoter of the MTHFR gene is associated with diabetic complications and biochemical indicators

BACKGROUND

DNA methylation is an epigenetic mechanism for regulating the transcription of many genes and has been linked to the development of various diseases. A promising gene to investigate is methylenetetrahydrofolate reductase (MTHFR), since the enzyme methylenetetrahydrofolate reductase (MTHFR) promotes methyl radical synthesis in the homocysteine cycle and can provide methyl groups for DNA methylation. In addition, several studies have correlated gene polymorphisms of this enzyme with a greater risk of diabetes, but little is known regarding the relationship between epigenetic changes in this gene and diabetes and its complications. The aim of this study was to investigate the relationship between methylation profile in the MTHFR gene promoter and biochemical, inflammatory and oxidative stress markers in individuals with type 2 diabetes (T2DM) who have been diagnosed for 5-10 years with or without diabetic retinopathy (DR) and nephropathy (DN).

METHODS

Specific PCR for methylation (MSP) was used to analyze MTHFR methylation profile in leucocytes DNA. Biochemical markers (glycemia, glycated hemoglobin, total cholesterol, LDL, HDL, triglycerides, serum creatinine), inflammatory markers (C-reactive protein and alpha-1 acid glycoprotein) and oxidative stress (total antioxidant and malonaldehyde) were determined in peripheric blood samples and microalbuminuria in 24 h urine samples. The X² and Mann-Whitney statistical tests were performed and p < 0.05 were considered significant.

RESULTS

The hypermethylated profile was most frequently observed in individuals with retinopathy (p < 0.01) and was associated with higher total cholesterol and LDL levels (p = 0.0046, 0.0267, respectively). Individuals with DN and hypermethylated profiles had higher levels of alpha-1 acid glycoprotein (p = 0.0080) and total antioxidant capacity (p = 0.0169) compared to subjects without complications.

CONCLUSIONS

Hypermethylation in the promoter of the MTHFR gene is associated with the occurrence of DR and with biochemical, inflammatory and oxidative stress parameters in the context of chronic complications.

Heritable DNA Methylation in CD4+ Cells among Complex Families Displays Genetic and Non-Genetic Effects

DNA methylation at CpG sites is both heritable and influenced by environment, but the relative contributions of each to DNA methylation levels are unclear. We conducted a heritability analysis of CpG methylation in human CD4+ cells across 975 individuals from 163 families in the Genetics of Lipid-lowering Drugs and Diet Network (GOLDN). Based on a broad-sense heritability (H2) value threshold of 0.4, we identified 20,575 highly heritable CpGs among the 174,445 most variable autosomal CpGs (SD > 0.02). Tests for associations of heritable CpGs with genotype at 2,145,360 SNPs using 717 of 975 individuals showed that ~74% were cis-meQTLs (< 1 Mb away from the CpG), 6% of CpGs exhibited trans-meQTL associations (>1 Mb away from the CpG or located on a different chromosome), and 20% of CpGs showed no strong significant associations with genotype (based on a p-value threshold of 1e-7). Genes proximal to the genotype independent heritable CpGs were enriched for functional terms related to regulation of T cell activation. These CpGs were also among those that distinguished T cells from other blood cell lineages. Compared to genes proximal to meQTL-associated heritable CpGs, genotype independent heritable CpGs were moderately enriched in the same genomic regions that escape erasure during primordial germ cell development and could carry potential for generational transmission.

Identification and validation of seven new loci showing differential DNA methylation related to serum lipid profile: an epigenome-wide approach. The REGICOR study

Lipid traits (total, low-density and high-density lipoprotein cholesterol, and triglycerides) are risk factors for cardiovascular disease. DNA methylation is not only an inherited but also modifiable epigenetic mark that has been related to cardiovascular risk factors. Our aim was to identify loci showing differential DNA methylation related to serum lipid levels. Blood DNA methylation was assessed using the Illumina Human Methylation 450 BeadChip. A two-stage epigenome-wide association study was performed, with a discovery sample in the REGICOR study (n = 645) and validation in the Framingham Offspring Study (n = 2,542). Fourteen CpG sites located in nine genes (SREBF1, SREBF2, PHOSPHO1, SYNGAP1, ABCG1, CPT1A, MYLIP, TXNIP and SLC7A11) and 2 intergenic regions showed differential methylation in association with lipid traits. Six of these genes and 1 intergenic region were new discoveries showing differential methylation related to total cholesterol (SREBF2), HDL-cholesterol (PHOSPHO1, SYNGAP1 and an intergenic region in chromosome 2) and triglycerides (MYLIP, TXNIP and SLC7A11). These CpGs explained 0.7%, 9.5% and 18.9% of the variability of total cholesterol, HDL cholesterol and triglycerides in the Framingham Offspring Study, respectively. The expression of the genes SREBF2 and SREBF1 was inversely associated with methylation of their corresponding CpGs (P-value = 0.0042 and 0.0045, respectively) in participants of the GOLDN study (n = 98). In turn, SREBF1 expression was directly associated with HDL cholesterol (P-value = 0.0429). Genetic variants in SREBF1, PHOSPHO1, ABCG1 and CPT1A were also associated with lipid profile. Further research is warranted to functionally validate these new loci and assess the causality of new and established associations between these differentially methylated loci and lipid metabolism.

Associations of polymorphisms in TXNIP and gene-environment interactions with the risk of coronary artery disease in a Chinese Han population

Single nucleotide polymorphisms (SNPs) in thioredoxin-interacting protein (TXNIP) gene may modulate TXNIP expression, then increase the risk of coronary artery disease (CAD). In a two-stage case-control study with a total of 1818 CAD patients and 1963 controls, we genotyped three SNPs in TXNIP and found that the variant genotypes of SNPs rs7212 [odds ratio (OR) = 1.26, P = 0.001] and rs7211 (OR = 1.23, P = 0.005) were significantly associated with increased CAD risk under a dominant model. In haplotype analyses, compared with the reference haplotype, haplotype 'G-T' had a 1.22-fold increased risk of CAD (P = 0.003). We also observed the cumulative effects of SNPs rs7212 and rs7211 on CAD risk and the severity of coronary atherosclerosis. Moreover, the gene-environment interactions among the variant genotypes of SNP rs7212, smoking habit, alcohol drinking habit and history of type 2 diabetes were associated with a 3.70-fold increased risk of CAD (P < 0.001). Subsequent genotype-phenotype correlation analyses further observed the significant effects of SNP rs7212 on TXNIP mRNA expression, plasma TXNIP and malondialdehyde levels. Taken together, our data suggest that TXNIP SNPs may individually and cumulatively affect CAD risk through a possible mechanism for regulating TXNIP expression and gene-environment interactions.

The role of global and regional DNA methylation and histone modifications in glycemic traits and type 2 diabetes: A systematic review

BACKGROUND

New evidence suggests the potential involvement of epigenetic mechanisms in type 2 diabetes (T2D) as a crucial interface between the effects of genetic predisposition and environmental influences.

AIM

To systematically review studies investigating the association between epigenetic marks (DNA methylation and histone modifications) with T2D and glycemic traits (glucose and insulin levels, insulin resistance measured by HOMA-IR).

METHOD AND RESULTS

Six bibliographic databases (Embase.com, Medline (Ovid), Web-of-Science, PubMed, Cochrane Central and Google Scholar) were screened until 28th August 2015. We included randomized controlled trials, cohort, case-control and cross-sectional studies in humans that examined the association between epigenetic marks (global, candidate or genome-wide methylation of DNA and histone modifications) with T2D, glucose and insulin levels and insulin metabolism. Of the initially identified 3879 references, 53 articles, based on 47 unique studies met our inclusion criteria. Overall, data were available on 10,823 participants, with a total of 3358 T2D cases. There was no consistent evidence for an association between global DNA-methylation with T2D, glucose, insulin and insulin resistance. The studies reported epigenetic regulation of several candidate genes for diabetes susceptibility in blood cells, muscle, adipose tissue and placenta to be related with T2D without any general overlap between them. Histone modifications in relation to T2D were reported only in 3 observational studies.

CONCLUSIONS AND RELEVANCE

Current evidence supports an association between epigenetic marks and T2D. However, overall evidence is limited, highlighting the need for further larger-scale and prospective investigations to establish whether epigenetic marks may influence the risk of developing T2D.

Epigenomic profiling reveals an association between persistence of DNA methylation and metabolic memory in the DCCT/EDIC type 1 diabetes cohort

We examined whether persistence of epigenetic DNA methylation (DNA-me) alterations at specific loci over two different time points in people with diabetes are associated with metabolic memory, the prolonged beneficial effects of intensive vs. conventional therapy during the Diabetes Control and Complications Trial (DCCT) on the progression of microvascular outcomes in the long-term follow-up Epidemiology of Diabetes Interventions and Complications (EDIC) Study. We compared DNA-me profiles in genomic DNA of whole blood (WB) isolated at EDIC Study baseline from 32 cases (DCCT conventional therapy group subjects showing retinopathy or albuminuria progression by EDIC Study year 10) vs. 31 controls (DCCT intensive therapy group subjects without complication progression by EDIC year 10). DNA-me was also profiled in blood monocytes (Monos) of the same patients obtained during EDIC Study years 16-17. In WB, 153 loci depicted hypomethylation, and 225 depicted hypermethylation, whereas in Monos, 155 hypomethylated loci and 247 hypermethylated loci were found (fold change ≥1.3; P < 0.005; cases vs. controls). Twelve annotated differentially methylated loci were common in both WB and Monos, including thioredoxin-interacting protein (TXNIP), known to be associated with hyperglycemia and related complications. A set of differentially methylated loci depicted similar trends of associations with prior HbA1c in both WB and Monos. In vitro, high glucose induced similar persistent hypomethylation at TXNIP in cultured THP1 Monos. These results show that DNA-me differences during the DCCT persist at certain loci associated with glycemia for several years during the EDIC Study and support an epigenetic explanation for metabolic memory.

DNA methylation of loci within ABCG1 and PHOSPHO1 in blood DNA is associated with future type 2 diabetes risk

Identification of subjects with a high risk of developing type 2 diabetes (T2D) is fundamental for prevention of the disease. Consequently, it is essential to search for new biomarkers that can improve the prediction of T2D. The aim of this study was to examine whether 5 DNA methylation loci in blood DNA (ABCG1, PHOSPHO1, SOCS3, SREBF1, and TXNIP), recently reported to be associated with T2D, might predict future T2D in subjects from the Botnia prospective study. We also tested if these CpG sites exhibit altered DNA methylation in human pancreatic islets, liver, adipose tissue, and skeletal muscle from diabetic vs. non-diabetic subjects. DNA methylation at the ABCG1 locus cg06500161 in blood DNA was associated with an increased risk for future T2D (OR = 1.09, 95% CI = 1.02–1.16, P-value = 0.007, Q-value = 0.018), while DNA methylation at the PHOSPHO1 locus cg02650017 in blood DNA was associated with a decreased risk for future T2D (OR = 0.85, 95% CI = 0.75–0.95, P-value = 0.006, Q-value = 0.018) after adjustment for age, gender, fasting glucose, and family relation. Furthermore, the level of DNA methylation at the ABCG1 locus cg06500161 in blood DNA correlated positively with BMI, HbA1c, fasting insulin, and triglyceride levels, and was increased in adipose tissue and blood from the diabetic twin among monozygotic twin pairs discordant for T2D. DNA methylation at the PHOSPHO1 locus cg02650017 in blood correlated positively with HDL levels, and was decreased in skeletal muscle from diabetic vs. non-diabetic monozygotic twins. DNA methylation of cg18181703 (SOCS3), cg11024682 (SREBF1), and cg19693031 (TXNIP) was not associated with future T2D risk in subjects from the Botnia prospective study.

Recent progress in genetic and epigenetic research on type 2 diabetes

Type 2 diabetes (T2DM) is a common complex metabolic disorder that has a strong genetic predisposition. During the past decade, progress in genetic association studies has enabled the identification of at least 75 independent genetic loci for T2DM, thus allowing a better understanding of the genetic architecture of T2DM. International collaborations and large-scale meta-analyses of genome-wide association studies have made these achievements possible. However, whether the identified common variants are causal is largely unknown. In addition, the detailed mechanism of how these genetic variants exert their effect on the pathogenesis of T2DM requires further investigation. Currently, there are ongoing large-scale sequencing studies to identify rare, functional variants for T2DM. Environmental factors also have a crucial role in the development of T2DM. These could modulate gene expression via epigenetic mechanisms, including DNA methylation, histone modification and microRNA regulation. There is evidence that epigenetic changes are important in the development of T2DM. Recent studies have identified several DNA methylation markers of T2DM from peripheral blood and pancreatic islets. In this review, we will briefly summarize the recent progress in the genetic and epigenetic research on T2DM and discuss how environmental factors, genetics and epigenetics can interact in the pathogenesis of T2DM.

Epigenetic associations of type 2 diabetes and BMI in an Arab population

Background

The prevalence of type 2 diabetes (T2D) and obesity has dramatically increased within a few generations, reaching epidemic levels. In addition to genetic risk factors, epigenetic mechanisms triggered by changing environment are investigated for their role in the pathogenesis of these complex diseases. Epigenome-wide association studies (EWASs) have revealed significant associations of T2D, obesity, and BMI with DNA methylation. However, populations from the Middle East, where T2D and obesity rates are highest worldwide, have not been investigated so far.

Methods

We performed the first EWAS in an Arab population with T2D and BMI and attempted to replicate 47 EWAS associations previously reported in Caucasians. We used the Illumina Infinium HumanMethylation450 BeadChip to quantify DNA methylation in whole blood DNA from 123 subjects of 15 multigenerational families from Qatar. To investigate the effect of differing genetic background and environment on the epigenetic associations, we further assessed the effect of replicated loci in 810 twins from UK.

Results

Our EWAS suggested a novel association between T2D and cg06721411 (DQX1; p value = 1.18 × 10⁻⁹). We replicated in the Qatari population seven CpG associations with BMI (SOCS3, p value = 3.99 × 10⁻⁶; SREBF1, p value = 4.33 × 10⁻⁵; SBNO2, p value = 5.87 × 10⁻⁵; CPT1A, p value = 7.99 × 10⁻⁵; PRR5L, p value = 1.85 × 10⁻⁴; cg03078551, intergenic region on chromosome 17; p value = 1.00 × 10⁻³; LY6G6E, p value = 1.10 × 10⁻³) and one with T2D (TXNIP, p value = 2.46 × 10⁻⁵). All the associations were further confirmed in the UK cohort for both BMI and T2D. Meta-analysis increased the significance of the observed associations and revealed strong heterogeneity of the effect sizes (apart from CPT1A), although associations at these loci showed concordant direction in the two populations.

Conclusions

Our study replicated eight known CpG associations with T2D or BMI in an Arab population. Heterogeneity of the effects at all loci except CPT1A between the Qatari and UK studies suggests that the underlying mechanisms might depend on genetic background and environmental pressure. Our EWAS results provide a basis for comparison with other ethnicities.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0177-6) contains supplementary material, which is available to authorized users.

Genome- and epigenome-wide association study of hypertriglyceridemic waist in Mexican American families

BACKGROUND

There is growing interest in the hypertriglyceridemic waist (HTGW) phenotype, defined as high waist circumference (≥95 cm in males and ≥80 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) as a marker of type 2 diabetes (T2D) and cardiovascular disease. However, the prevalence of this phenotype in high-risk populations, its association with T2D, and the genetic or epigenetic influences on HTGW are not well explored. Using data from large, extended families of Mexican Americans (a high-risk minority population in the USA) we aimed to: (1) estimate the prevalence of this phenotype, (2) test its association with T2D and related traits, and (3) dissect out the genetic and epigenetic associations with this phenotype using genome-wide and epigenome-wide studies, respectively.

RESULTS

Data for this study was from 850 Mexican American participants (representing 39 families) recruited under the ongoing San Antonio Family Heart Study, 26 % of these individuals had HTGW. This phenotype was significantly heritable (h (2) r = 0.52, p = 1.1 × 10(-5)) and independently associated with T2D as well as fasting glucose levels and insulin resistance. We conducted genome-wide association analyses using 759,809 single nucleotide polymorphisms (SNPs) and epigenome-wide association analyses using 457,331 CpG sites. There was no evidence of any SNP associated with HTGW at the genome-wide level but two CpG sites (cg00574958 and cg17058475) in CPT1A and one CpG site (cg06500161) in ABCG1 were significantly associated with HTGW and remained significant after adjusting for the closely related components of metabolic syndrome. CPT1A holds a cardinal position in the metabolism of long-chain fatty acids while ABCG1 plays a role in triglyceride metabolism.

CONCLUSIONS

Our results reemphasize the value of HTGW as a marker of T2D. This phenotype shows association with DNA methylation within CPT1A and ABCG1, genes involved in fatty acid and triglyceride metabolism. Our results underscore the importance of epigenetics in a clinically informative phenotype.

Type 2 diabetes and leucocyte DNA methylation: an epigenome-wide association study in over 1,500 older adults

AIMS/HYPOTHESIS

Development of type 2 diabetes depends on environmental and genetic factors. We investigated the epigenome-wide association of prevalent diabetes with DNA methylation (DNAm) in peripheral blood.

METHODS

DNAm was measured in whole blood with the Illumina Infinium HumanMethylation450 BeadChip in two subsamples of participants from the ESTHER cohort study. Cohort 1 included 988 participants, who were consecutively recruited between July and October 2000 and cohort 2 included 527 randomly selected participants. The association of DNAm with prevalent type 2 diabetes at recruitment was estimated using median regression analysis adjusting for sex, age, BMI, smoking behaviour, cell composition and batch at 361,922 CpG sites.

RESULTS

Type 2 diabetes was prevalent in 16% of the participants, and diabetes was poorly controlled in 45% of the diabetic patients. In cohort 1 (discovery) DNAm at 39 CpGs was significantly associated with prevalent diabetes after correction for multiple testing. In cohort 2 (replication) at one of these CpGs, DNAm was still significantly associated. Decreasing methylation levels at cg19693031 with increasing fasting glucose and HbA1c concentrations were observed using restricted cubic spline analysis. In diabetic patients with poorly controlled diabetes, the decrease in estimated DNAm levels was approximately 5% in comparison with participants free of diagnosed diabetes.

CONCLUSIONS/INTERPRETATION

Cg19693031, which is located within the 3'-untranslated region of TXNIP, might play a role in the pathophysiology of type 2 diabetes. This result appears biologically plausible given that thioredoxin-interacting protein is overexpressed in diabetic animals and humans and 3'-untranslated regions are known to play a regulatory role in gene expression.

Epigenome-wide association study identifies TXNIP gene associated with type 2 diabetes mellitus and sustained hyperglycemia

Type 2 diabetes mellitus (DM) is an established risk factor for a wide range of vascular diseases, including ischemic stroke (IS). Glycated hemoglobin A1c (HbA1c), a marker for average blood glucose levels over the previous 12 weeks, is used as a measure of glycemic control and also as a diagnostic criterion for diabetes (HbA1c levels ≥ 6.5%). Epigenetic mechanisms, such as DNA methylation, may be associated with aging processes and with modulation of the risk of various pathologies, such as DM. Specifically, DNA methylation could be one of the mechanisms mediating the relation between DM and environmental exposures. Our goal was to identify new CpG methylation sites associated with DM. We performed a genome-wide methylation study in whole-blood DNA from an IS patient cohorts. Illumina HumanMethylation450 BeadChip array was used to measure DNA methylation in CpG sites. All statistical analyses were adjusted for sex, age, hyperlipidemia, body mass index (BMI), smoking habit and cell count. Findings were replicated in two independent cohorts, an IS cohort and a population-based cohort, using the same array. In the discovery phase (N = 355), we identified a CpG site, cg19693031 (located in the TXNIP gene) that was associated with DM (P = 1.17 × 10(-12)); this CpG was replicated in two independent cohorts (N = 167 and N = 645). Methylation of TXNIP was inversely and intensely associated with HbA1c levels (P = 7.3 × 10(-16)), specifically related to diabetic patients with poor control of glucose levels. We identified an association between the TXNIP gene and DM through epigenetic mechanisms, related to sustained hyperglycemia levels (HbA1c ≥ 7%).