Interaction between a diabetes-related methylation site (TXNIP cg19693031) and variant (GLUT1 rs841853) on fasting blood glucose levels among non-diabetics

BACKGROUND: Type 2 diabetes mellitus (T2DM) is caused by a combination of environmental, genetic, and epigenetic factors including, fasting blood glucose (FBG), genetic variant rs841853, and cg19693031 methylation. We evaluated the interaction between rs841853 and cg19693031 on the FBG levels of non-diabetic Taiwanese adults. METHODS: We used Taiwan Biobank (TWB) data collected between 2008 and 2016. The TWB data source contains information on basic demographics, personal lifestyles, medical history, methylation, and genotype. The study participants included 1300 people with DNA methylation data. The association of cg19693031 methylation (stratified into quartiles) with rs841853 and FBG was determined using multiple linear regression analysis. The beta-coefficients (β) and p-values were estimated. RESULTS: The mean ± standard deviation (SD) of FBG in rs841853-CC individuals (92.07 ± 7.78) did not differ significantly from that in the CA + AA individuals (91.62 ± 7.14). However, the cg19693031 methylation levels were significantly different in the two groups (0.7716 ± 0.05 in CC individuals and 0.7631 ± 0.05 in CA + AA individuals (p = 0.002). The cg19693031 methylation levels according to quartiles were β < 0.738592 (< Q1), 0.738592 ≤ 0.769992 (Q1–Q2), 0.769992 ≤ 0.800918 (Q2–Q3), and β ≥ 0.800918 (≥ Q3). FBG increased with decreasing cg19693031 methylation levels in a dose–response manner (p(trend) = 0.005). The β-coefficient was − 0.0236 (p = 0.965) for Q2–Q3, 1.0317 (p = 0.058) for Q1–Q2, and 1.3336 (p = 0.019 for < Q1 compared to the reference quartile (≥ Q3). The genetic variant rs841853 was not significantly associated with FBG. However, its interaction with cg19693031 methylation was significant (p-value = 0.036). Based on stratification by rs841853 genotypes, only the CC group retained the inverse and dose–response association between FBG and cg19693031 methylation. The β (p-value) was 0.8082 (0.255) for Q2–Q3, 1.6930 (0.022) for Q1–Q2, and 2.2190 (0.004) for < Q1 compared to the reference quartile (≥ Q3). The p(trend) was 0.002. CONCLUSION: Summarily, methylation at cg19693031 was inversely associated with fasting blood glucose in a dose-dependent manner. The inverse association was more prominent in rs841853-CC individuals, suggesting that rs841853 could modulate the association between cg19693031 methylation and FBG. Our results suggest that genetic variants may be involved in epigenetic mechanisms associated with FBG, a hallmark of diabetes. Therefore, integrating genetic and epigenetic data may provide more insight into the early-onset of diabetes. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12967-022-03269-y.