Email Alert | RSS

Journal of Peking University (Health Sciences) ›› 2020, Vol. 52 ›› Issue (3): 425-431. doi: 10.19723/j.issn.1671-167X.2020.03.005

Previous Articles     Next Articles

Correlation between fasting plasma glucose, HbA1c and DNA methylation in adult twins

Zhao-nian WANG1,Wen-jing GAO1,△(),Bi-qi WANG1,Wei-hua CAO1,Jun LV1,Can-qing YU1,Zeng-chang PANG2,Li-ming CONG3,Hua WANG4,Xian-ping WU5,Yu LIU6,Li-ming LI1   

  1. 1. Department of Epidemiology and Biostatistics, Peking University School of Public Health, Beijing 100191, China
    2. Qingdao Municipal Center for Disease Control and Prevention, Qingdao 266033, Shandong, China
    3. Zhejiang Center for Disease Control and Prevention, Hangzhou 310051, China
    4. Jiangsu Center for Disease Control and Prevention, Nanjing 210009, China
    5. Sichuan Center for Disease Control and Prevention, Chengdu 610041, China
    6. Center for Disease Control and prevention, Heilongjiang Agricultural Reclamation Bureau, Harbin 150090, China
  • Received:2020-02-15 Online:2020-06-18 Published:2020-06-30
  • Contact: Wen-jing GAO E-mail:pkuepigwj@126.com
  • Supported by:
    Special Found for Health Scienti-fic Research in the Public Welfare(201502006);Special Found for Health Scienti-fic Research in the Public Welfare(201002007);National Natural Science Foundation of China(81573223)

RICH HTML

 20   

PDF (PC)

176

Abstract:

Objective: To explore the cytidine-phosphate-guanosine (CPG) sites associated with fas-ting plasma glucose (FPG) and glycated haemoglobin (HbA1c) in twins.Methods: In the study, 169 pairs of monozygotic twins were recruited in Qingdao, Zhejiang, Jiangsu, Sichuan and Heilongjiang in June to December of 2013 and June 2017 to October 2018. The methylation was detected by Illumina Infinium HumanMethylation450 BeadChip and Illumina Infinium MethylationEPIC BeadChip. According to the Linear Mixed Effect model (LME model), fasting plasma glucose and HbA1c were taken as the main effects, the methylation level (β value) was taken as the dependent variable, continuous variables, such as age, body mass index (BMI), blood pressure, components of blood cells, surrogate variables generated by SVA, and categorical variables, such as gender, smoking and drinking status, hypoglycemic drugs taking, were included in the fixed effect model as covariates, and the identity numbers (ID) of the twins was included in the random effect model. The intercept was set as a random. Regression analysis was carried out to find out the CpG sites related to fasting blood glucose or HbA1c, respectively.Results: In this study, 338 monozygotic twins (169 pairs) were included, with 412 459 CpG loci. Among them, 114 pairs were male, and 55 pairs were female, with an average age of (48.2±11.9) years. After adjustment of age, gender, BMI, blood pressure, smoking, drinking, blood cell composition, and other covariates, and multiple comparison test, 7 CpG sites (cg19693031, cg01538969, cg08501915, cg04816311, ch.8.1820050F, cg06721411, cg26608667) were found related to fasting blood glucose, 3 of which (cg08501915, ch.8.1820050f, cg26608667) were the newly found sites in this study; whereas 10 CpG sites (cg19693031, cg04816311, cg01538969, cg01339781, cg01676795, cg24667115, cg09029192, cg20697417, ch.4.1528651F, cg16097041) were found related to HbA1c, and 4 of which(cg01339781, cg24667115, cg20697417, and ch.4.1528651f) were new. We found that cg19693031 in TXNIP gene was the lowest P-value site in the association analysis between DNA methylation and fas-ting plasma glucose and HbA1c (PFPG=2.42×10-19, FDRFPG<0.001; PHbA1c=1.72×10-19, FDRHbA1c<0.001).Conclusion: In this twin study, we found new CpG sites related to fasting blood glucose and HbA1c, and provided some clues that partly revealed the potential mechanism of blood glucose metabolism in terms of DNA methylation, but it needed further verification in external larger samples.

Key words: Fasting plasma glucose, Glycated haemoglobin, DNA methylation, Twin

CLC Number: 

  • R179

Table 1

Demographic characteristics and basic information of related variables among twins"

Variables Values
Pair of twins, n 169
Age/years 48.2 ± 11.9
FPG/(mmol/L) 5.9 ± 1.9
HbA1c/% 5.9 ± 1.0
BMI/(kg/m2) 25.0 ± 3.5
SBP/mmHg 135.0 ± 18.7
DBP/mmHg 81.3 ± 11.2
Gender, n(%)
Male 228 (67.5)
Female 110 (32.5)
Type 2 diabetes mellitus, n(%) 35 (10.4)
Hypertension, n(%) 6 (1.8)
Hypoglycemic drug use, n(%) 29 (8.6)
Drinking, n(%)
Never 160 (47.3)
Used to 4 (1.2)
Now 174 (51.5)
Smoking, n(%)
Never 184 (54.4)
Used to 48 (14.2)
Now 106 (31.4)

Figure 1

Manhattan plot of correlation analysis between blood glucose related indicators and DNA methylation of the whole genome The abscissa represents the chromosome, the ordinate represents the negative logarithm with 10 as the base of the original P value, the points above the red line represent P < 5 × 10-8, and the points in green are those which still significant after Bonferroni test. A, fasting blood glucose; B, HbA1c."

Table 2

Association analysis between DNA methylation and indicators of plasma glucose levels among the whole genome"

G site Chromosome Position on
chromosome		 Gene Position
on gene		 Relation to
CpG island		 Slope SE P value FDR
FPG
cg19693031 1 145441552 TXNIP 3'UTR OpenSea -1.44×10-2 1.32×10-3 2.42×10-19 <0.001
cg01538969 6 30624636 DHX16 Body OpenSea 4.76×10-3 6.87×10-4 2.68×10-10 <0.001
cg08501915 4 129208894 PGRMC2 1st exon Island 5.80×10-3 8.36×10-4 2.70×10-10 <0.001
cg04816311 7 1066650 C7orf50 Body N_Shore 7.16×10-3 1.04×10-3 3.52×10-10 <0.001
ch.8.1820050F 8 89645708 - - OpenSea 4.12×10-3 6.04×10-4 4.75×10-10 <0.001
cg06721411 2 74753759 DQX1 TSS1500 N_Shelf 4.53×10-3 7.73×10-4 4.73×10-8 0.019
cg26608667 7 1196370 ZFAND2A Body N_Shelf 4.52×10-3 7.96×10-4 1.04×10-7 0.043
HbA1c
cg19693031 1 145441552 TXNIP 3'UTR OpenSea -2.91×10-2 2.66×10-3 1.72×10-19 <0.001
cg04816311 7 1066650 C7orf50 Body N_Shore 1.66×10-2 2.04×10-3 6.00×10-13 <0.001
cg01538969 6 30624636 DHX16 Body OpenSea 9.55×10-3 1.39×01-3 3.31×10-10 <0.001
cg01339781 6 116989657 ZUFSP 5'UTR OpenSea 2.52×10-3 4.03×10-4 7.05×10-9 0.003
cg01676795 7 75586348 POR Body OpenSea 1.21×10-02 1.97×10-3 1.24×10-8 0.005
cg24667115 6 91004482 BACH2 5'UTR N_Shore 1.36×10-2 2.24×10-3 1.90×10-8 0.008
cg09029192 17 76015204 TNRC6C 5'UTR OpenSea 8.13×10-3 1.39×10-3 4.75×10-8 0.020
cg20697417 1 41786797 - - OpenSea 8.82×10-3 1.51×10-3 5.22×10-8 0.021
ch.4.1528651F 4 79200030 FRAS1 Body OpenSea 6.73×10-3 1.17×10-3 6.99×10-8 0.029
cg16097041 1 154965544 FLAD1 3'UTR OpenSea 8.30×10-3 1.44×10-3 7.77×10-8 0.032

Table 3

Overlapping sites in association analysis between DNA methylation and indictors of plasma glucose level"

CpG site Chromosome Position on
chromosome		 Gene Position on gene HbA1c FPG
P value FDR P value FDR
cg19693031 1 145441552 TXNIP 3'UTR 1.72×10-19 <0.001 2.42×10-19 <0.001
cg01538969 6 30624636 DHX16 Body 3.31×10-10 <0.001 2.68×10-10 <0.001
cg04816311 7 1066650 C7orf50 Body 6.00×10-13 <0.001 3.52×10-10 <0.001

Figure 2

Q-Q chart of correlation analysis between blood glucose related indexes and DNA methylation of the whole genome The abscissa represents the negative logarithm based on 10 of the predicted P values, and the ordinate represents the negative logarithm based on 10 of the original P values. A, Fasting blood glucose; B, HbA1c."

[1] Cho NH, Shaw JE, Karuranga S, et al. IDF diabetes atlas: Glo-bal estimates of diabetes prevalence for 2017 and projections for 2045[J]. Diabetes Res Clin Pract, 2018,138:271-281.
[2] Wang L, Gao P, Zhang M, et al. Prevalence and ethnic pattern of diabetes and prediabetes in China in 2013[J]. JAMA, 2017,317(24):2515-2523.
doi: 10.1001/jama.2017.7596 pmid: 28655017
[3] Zheng Y, Ley SH, Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications[J]. Nat Rev Endocrinol, 2018,14(2):88-98.
[4] American Diabetes Association. Classification and diagnosis of diabetes: Standards of medical care in diabetes-2018[J]. Diabetes Care, 2018,41(Suppl 1):13-27.
[5] World Health Organization. Use of glycated haemoglobin (HbA1c) in the diagnosis of diabetes mellitus: Abbreviated report of a who consultation[R/OL]. (2011-01-12)[2020-01-05]. https://www.who.int/diabetes/publications/report-hba1c_2011.pdf.
[6] Singh GM, Danaei G, Farzadfar F, et al. The age-specific quantitative effects of metabolic risk factors on cardiovascular diseases and diabetes: A pooled analysis[J]. PLoS One, 2013,8(7):e65174.
pmid: 23935815
[7] 高文静, 李立明. 以分开抚养双生子为基础的研究进展[J]. 中华医学遗传学杂志, 2014,31(3):327-329.
[8] Tan Q, Christiansen L, von Bornemann Hjelmborg J, et al. Twin methodology in epigenetic studies[J]. J Exp Biol, 2015,218(Pt 1):134-139.
pmid: 25568460
[9] Hwang JY, Lee HJ, Go MJ, et al. Genome-wide methylation ana-lysis identifies ELOVL5 as an epigenetic biomarker for the risk of type 2 diabetes mellitus[J]. Sci Rep, 2018,8(1):14862.
pmid: 30291282
[10] Liu F, Sun Q, Wang L, et al. Bioinformatics analysis of abnormal DNA methylation in muscle samples from monozygotic twins discordant for type 2 diabetes[J]. Mol Med Rep, 2015,12(1):351-356.
pmid: 25760736
[11] Yuan W, Xia Y, Bell CG, et al. An integrated epigenomic analysis for type 2 diabetes susceptibility loci in monozygotic twins[J]. Nat Commun, 2014,5:5719.
[12] Fortin JP, Triche TJ Jr., Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi[J]. Bioinformatics, 2017,33(4):558-560.
pmid: 28035024
[13] Pidsley R, Cc YW, Volta M, et al. A data-driven approach to preprocessing Illumina 450K methylation array data[J]. BMC Genomics, 2013,14:293.
[14] Wang B, Gao W, Yu C, et al. Determination of zygosity in adult Chinese twins using the 450K methylation array versus questionnaire data[J]. PLoS One, 2015,10(4):e0123992.
pmid: 25927701
[15] Houseman EA, Accomando WP, Koestler DC, et al. DNA methylation arrays as surrogate measures of cell mixture distribution[J]. BMC Bioinformatics, 2012,13:86.
doi: 10.1186/1471-2105-13-86 pmid: 22568884
[16] Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis[J]. PLoS Genet, 2007,3(9):1724-1735.
pmid: 17907809
[17] Michels KB, Binder AM. Considerations for design and analysis of DNA methylation studies[J]. Methods Mol Biol, 2018,1708:31-46.
[18] Carlin JB, Gurrin LC, Sterne JA, et al. Regression models for twin studies: A critical review[J]. Int J Epidemiol, 2005,34(5):1089-1099.
pmid: 16087687
[19] Soriano-Tarraga C, Jimenez-Conde J, Giralt-Steinhauer E, et al. Epigenome-wide association study identifies TXNIP gene associated with type 2 diabetes mellitus and sustained hyperglycemia[J]. Hum Mol Genet, 2016,25(3):609-619.
[20] Meeks KAC, Henneman P, Venema A, et al. Epigenome-wide association study in whole blood on type 2 diabetes among sub-Saharan African individuals: Findings from the RODAM study[J]. Int J Epidemiol, 2019,48(1):58-70.
[21] Pena GG, Dutra MS, Gazzinelli A, et al. Heritability of phenotypes associated with glucose homeostasis and adiposity in a rural area of Brazil[J]. Ann Hum Genet, 2014,78(1):40-49.
[22] Cardona A, Day FR, Perry JRB, et al. Epigenome-wide association study of incident type 2 diabetes in a british population: EPIC-Norfolk study[J]. Diabetes, 2019,68(12):2315-2326.
[23] Al Muftah WA, Al-Shafai M, Zaghlool SB, et al. Epigenetic associations of type 2 diabetes and BMI in an Arab population[J]. Clin Epigenetics, 2016,8:13.
[24] Kulkarni H, Kos MZ, Neary J, et al. Novel epigenetic determinants of type 2 diabetes in Mexican-American families[J]. Hum Mol Genet, 2015,24(18):5330-5344.
[25] Galmozzi A, Kok BP, Kim AS, et al. PGRMC2 is an intracellular haem chaperone critical for adipocyte function[J]. Nature, 2019,576(7785):138-142.
[26] Alhawiti NM, Al Mahri S, Aziz MA, et al. TXNIP in metabolic regulation: Physiological role and therapeutic outlook[J]. Curr Drug Targets, 2017,18(9):1095-1103.
[27] Waldhart AN, Dykstra H, Peck AS, et al. Phosphorylation of TXNIP by AKT mediates acute influx of glucose in response to insulin[J]. Cell Rep, 2017,19(10):2005-2013.
[28] Gencheva M, Kato M, Newo AN, et al. Contribution of deah-box protein DHX16 in human pre-mRNA splicing[J]. Biochem J, 2010,429(1):25-32.
pmid: 20423332
[1] Xue LOU,Li LIAO,Xing-jun LI,Nan WANG,Shuang LIU,Ruo-mei CUI,Jian XU. Methylation status and expression of TWEAK gene promoter region in peripheral blood of patients with rheumatoid arthritis [J]. Journal of Peking University (Health Sciences), 2021, 53(6): 1020-1025.
[2] SONG Qin-feng,LI Hong-tian,YANG Jing,YUAN Peng-bo,CHENG Zhi-hao,LIU Jian-meng,ZHAO Yang-yu. Thyroid function of first-trimester twin pregnant women and its association with preterm delivery [J]. Journal of Peking University (Health Sciences), 2021, 53(3): 473-478.
[3] WU Zhen-tian, GAO Wen-jing, WANG Bi-qi, CAO Wei-hua, LV Jun, YU Can-qing， PANG Zeng-chang, CONG Li-ming, WANG Hua, WU Xian-ping, LI Li-ming. Correlation between blood pressure and DNA methylation in adult twins [J]. Journal of Peking University(Health Sciences), 2018, 50(3): 387-394.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
[1] Author. English Title Test[J]. Journal of Peking University(Health Sciences), 2010, 42(1): 1 -10 .
[2] . [J]. Journal of Peking University(Health Sciences), 2009, 41(2): 188 -191 .
[3] . [J]. Journal of Peking University(Health Sciences), 2009, 41(3): 376 -379 .
[4] . [J]. Journal of Peking University(Health Sciences), 2009, 41(4): 459 -462 .
[5] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 319 -322 .
[6] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 333 -336 .
[7] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 337 -340 .
[8] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 225 -328 .
[9] . [J]. Journal of Peking University(Health Sciences), 2007, 39(4): 346 -350 .
[10] . [J]. Journal of Peking University(Health Sciences), 2007, 39(4): 361 -364 .
Copyright © Journal of Peking University (Health Sciences), All Rights Reserved.
Powered by Beijing Magtech Co. Ltd