Skip to main content

DNA methylation risk score for type 2 diabetes is associated with gestational diabetes

Abstract

Background

Gestational diabetes mellitus (GDM) and type 2 diabetes mellitus (T2DM) share many pathophysiological factors including genetics, but whether epigenetic marks are shared is unknown. We aimed to test whether a DNA methylation risk score (MRS) for T2DM was associated with GDM across ancestry and GDM criteria.

Methods

In two independent pregnancy cohorts, EPIPREG (n = 480) and EPIDG (n = 32), DNA methylation in peripheral blood leukocytes was measured at a gestational age of 28 ± 2. We constructed an MRS in EPIPREG and EPIDG based on CpG hits from a published epigenome-wide association study (EWAS) of T2DM.

Results

With mixed models logistic regression of EPIPREG and EPIDG, MRS for T2DM was associated with GDM: odd ratio (OR)[95% CI]: 1.3 [1.1–1.8], P = 0.002 for the unadjusted model, and 1.4 [1.1–1.7], P = 0.00014 for a model adjusted by age, pre-pregnant BMI, family history of diabetes and smoking status. Also, we found 6 CpGs through a meta-analysis (cg14020176, cg22650271, cg14870271, cg27243685, cg06378491, cg25130381) associated with GDM, and some of their methylation quantitative loci (mQTLs) were related to T2DM and GDM.

Conclusion

For the first time, we show that DNA methylation marks for T2DM are also associated with GDM, suggesting shared epigenetic mechanisms between GDM and T2DM.

Introduction

Type 2 Diabetes Mellitus (T2DM) affects 10.5% of the global population in the age group 20 to 79 years. T2DM threatens health of the individuals and healthcare systems due to its numerous complications and high healthcare cost and is among the top ten causes of death in the world [1].

During pregnancy, insulin resistance increases to maintain adequate glucose flow to the offspring. However, if the pancreatic beta cells cannot compensate with sufficiently high insulin secretion, this can result in hyperglycaemia and Gestational Diabetes Mellitus (GDM), which is defined by hyperglycaemia with first onset during pregnancy [2]. GDM increases the risk of pregnancy complications such as pre-eclampsia, caesarean section, neonatal hypoglycaemia, preterm birth, and foetal macrosomia. Moreover, GDM severely increases the future T2DM risk in the mother [3]. Although the prevalence of GDM is increasing globally, it varies depending on population characteristics and the diagnostic criteria used [4].

Epigenetics is the study of changes in gene expression that are not caused by changes in the DNA sequence. Some epigenome-wide association studies (EWAS) suggest that epigenetic mechanisms contribute to the pathogenesis of T2DM [5,6,7]. Genetic risk scores (GRS) have been increasingly used to assess disease risk, such as for T2DM [8, 9]. Similarly, methylation risk scores (MRS) are increasingly studied, to assess associations with outcomes of interest such as some cancers and kidney disease [10], and show promising potential as a tool to aid the prediction of T2DM and to understand gene-environment interactions. Schrader et al. showed that MRS separated T2DM subjects into different groups and were associated with diabetic complications like cardiovascular disease, chronic kidney disease and retinopathy [11].

Despite similar genetics and pathophysiology between GDM and T2DM [12], genetics only explain a small proportion of overall T2DM risk, and environment is known to play an important role for epigenetics. However, whether epigenetic marks are common between GDM and T2D has not been reported previously. Epigenetics marks common for GDM and T2DM may help us understand these pathological mechanisms better in order to prevent, diagnose or improve treatment of GDM and T2DM. We aimed to test whether an MRS for T2DM is associated with a GDM across ancestry and GDM criteria.

Methods

Study population

The EPIPREG cohort

A subgroup of 480 women was selected from the STORK Groruddalen (STORK G) pregnancy cohort. STORK G included 823 healthy women from different ethnic origins (European, South Asian, African, Middle Eastern and South American) who attended three different Child Health Clinics in the area of Groruddalen, Oslo, Norway, during the 2008–2010 period [13]. Ethnic origin was defined by either the individual’s country of birth or their mother’s country of birth if the latter was born outside of Europe [17]. The EPIPREG subgroup has a total of 312 European subjects (EPIPREG_EU), whereof 73 were diagnosed with GDM, and 168 South Asians (EPIPREG_SA), whereof 68 were diagnosed with GDM. European and South Asian ancestry was determined by genetic principal components [14]. Fasting blood samples were collected, and a 75 g oral glucose tolerance test (OGTT) was offered to all women (universal testing) at week 28 ± 2 of pregnancy. For the present study GDM was classified according to the slightly modified International Association of the Diabetes and Pregnancy (IADPSG) criteria (fasting glucose ≥ 5.1 mmol/l and/or 2-hour glucose ≥ 8.5 mmol/l, as 1-h glucose values were not available).

The Norwegian Regional Committee for Medical Health Research Ethics South East approved the STORK-G study including genetic and epigenetic data (ref. number 2015/1035). Written informed consent from all participants was obtained before any study-related procedure.

EPIDG cohort

A total of 32 (16 GDM, 16 non-GDM) pregnant women were selected from EPIDG cohort, which started in 2019 and is still recruiting participants. The EPIDG cohort has a total of 230 Mediterranean (South European) pregnant women who attended the Unit of Diabetes and Pregnancy at University Hospital Virgen de la Victoria, Málaga, Spain. All participants gave their consent to participate in the study. GDM criteria followed a two-step strategy from the National Diabetes Data Group (NDDG) guideline [15]. The first step was a screening using the O’Sullivan test (50 g glucose overload) between 24 and 28 weeks. Then an oral glucose tolerance test (OGTT)-100gr was performed on those women with a positive result in the O’Sullivan test (≥ 0.7.8 mmol/L). GDM was diagnosed if glucose values were for the OGTT-100 higher than the threshold, at least at two points: fasting > 5.8 mmol/l, after 1 h > 10.6 mmol/l, after 2 h > 9.2 mmol/L, after 3 h > 8.0mmol/l. Pregnant women with normal OGTT-100gr were considered as controls (non-GDM). The 32 samples were selected based on the availability of blood samples and matched by age, gestational age and pre-pregnant BMI. The study was approved by the Institutional review board at the Hospital Universitario Virgen de la Victoria, Málaga.

Samples extraction, DNA isolation and bisulfite conversion

In both cohorts, samples were collected in gestational week 28 ± 2. The samples were either aliquoted and biobanked or subject to routine laboratory analyses. In the EPIPREG sample, DNA from peripheral blood leukocytes was extracted subsequently throughout the data collection, at the Hormone Laboratory, Oslo University Hospital, using a salting out procedure [16], and stored at -80ºC. EZ DNA methylation TM Kit (ZYMO Research, Tustin, CA, USA) was used for the bisulphite conversion of DNA. In EPIDG cohort, DNA from peripheral blood leukocytes was extracted using DNA Blood Mini Kit (Quiagen, Hiden, Germany). Epitect Bisulfite Kit (Qiagen, Germany) was used for bisulphite conversion. DNA methylation in both cohorts was quantified with Infinium MethylationEPIC BeadChip (Illumina, San Diego, USA).

Methylation values extraction

Raw data from both studies were analysed in R. In EPIPREG, the Meffil R package (https://cran.r-project.org/) was used for quality control (QC), normalization and reporting of B-values [14]. During the QC procedure implemented in Meffil R package, we removed six samples who were considered outliers from the methylated to unmethylated ratio comparison (> 3SD), one sample that was an outlier in control probes bisulfite 1 and 2 (> 5SD), and one sample due to sex discrepancy (predicted sex outlier > 5 SD). Probes with a bead count fewer than three and dectection p-value < 0.01. In total 472 of the 480 samples and 864,560 probes passed the QC. Functional normalization was implemented in Meffil R, which takes into account potential batch effects such as slide, row and columns. Proportions of blood cells (namely, CD8T, CD4T, NK, β-cells, monocytes, and neutrophils) were calculated using the Houseman method [27] during the QC [17].

After processing, 472 of the 480 available individuals remained and 810.266 in the EPIPREG sample passed QC. In EPIDG, the full sample and 741.479 probes passed QC (Supplementary Fig. 1).

Methylation risk score

To construct the weighted DNA MRS for T2DM in our sample we used summary data for the CpGs sites discovered in an EWAS of T2DM in five prospective European cohorts (n = 5.859) [7, 18,19,20,21]. The regression coefficients for the EWAS of T2DM was based in β-values. We calculated MRS for each of the discovered CpG sites available in our two cohorts. Due the QC filters of CpGs in each cohort only 42 of 72 CpG could be included in the MRS (Supplementary Table 1). MRS was constructed by multiplying the regression coefficient from the EWAS of T2DM with individual β-values in EPIPREG and EPIDG (EPIPREG-EU, EPIPREG-SA and EPIDG). Thereafter we summarized the individual score for each of the 42 CpGs sites to obtain the MRS. Figure 1 illustrates the workflow, including development of the MRS and the statistical analysis.

Fig. 1
figure 1

Workflow of the methodology followed to construct the MRS and the rest of the statistical analysis. preBMI: previous pregnancy BMI, Fam his diabetes: family history of diabetes. Blue Line indicate the in-house data vs. the public data used in our study

Statistical analysis

All the statistical analyses were performed using Rstudio (3.4.4). Clinical variables were compared between GDM and non-GDM in each group separately (EPIPRE_EU, EPIPREG_SA and EPIDG), depending on if the variable was normally distributed or not, t-test or Whitney U test were used for the continuous variable. A logistic mixed model regression was used to elucidate if MRS was associated with GDM in our three groups. Linear mixed model regressions were performed to test the association between MRS and continuous variables such as HOMA-IR, C-peptide, fasting glucose and 2 h glucose level, using a Bonferroni corrected threshold p < 0.0125 (0.05/4). We used the R packages lme4 and lmerTest [22] to perform mixed models regression. In these models, ancestry (European, South Asian, Mediterranean (EPIDG)) was treated as a random intercept to overcome potential ancestry-related differences in DNA methylation. We performed an unadjusted and an adjusted model by age, pre-pregnant BMI, family history of diabetes and smoking status, since these are covariates that are associated with both epigenetic marks and gestational diabetes.

To identify individual CpG sites associated with GDM in the studied populations, we performed logistic regressions for each CpGs used in the MRS across the three groups, followed by a fixed effects meta-analysis of the three groups for each CpG site using METAL [23]. Since random-effect models may bias the results if smaller studies have a large effect compared to well-performed, larger studies [24], we decided to use a fixed-effect model for a conservative approach, where the case-control design in EPIDG would not count more than the population-based design in EPIPREG.

mQTLs and phenotypes associated

Due to low statistical power for mQTL analysis in EPIPREG, we performed look-ups in GoDMC of the significant CpGs obtained in the meta-analysis to identify methylation quantitative trait loci (mQTL) to see if these methylation sites were influenced by genetic variants. GoDMC include Cis and trans meta-analysis results from genome-wide scans of 420.509 DNA methylation sites. This information come from several projects, with the aim to share data on genetic basis of DNA methylation variation. GoDMC provides a list of SNPs that have been associated with the CpG of interest. We used the LDlink-SNPclip tool [25] to identify SNPs in linkage disequilibrium (r2 > 0.8, MAF = 0.01) and the selected populations were European and South Asians. To assess relevance of the mQTLs to diabetes related phenotypes, we performed explorative look-ups in Phenoscanner [26] of phenotypes nominally associated with the mQTLs. Phenoscanner is a database holding publicly available results from large-scale GWAS. Although the default threshold used by Phenoscanner is p < 1 × 10− 5, we also looked up phenotypes with p-value < 0.05 due to the explorative nature of this look-up [27].

Results

Characteristics of the samples

The clinical characteristics of EPIPREG and EPIDG are represented in Table 1. In the EPIPREG cohort, non-GDM and GDM women differed significantly in pre-pregnant BMI, fasting glucose (Gluc0), 2 h glucose after OGTT (Gluc2), systolic blood pressure (SBP), glycosylated hemoglobin A1C (HbA1C), Homeostatic Model Assessment for insulin resistance (HOMA-IR), triglycerides (TAG), cholesterol and HDL across ethnicities. Only fasting glucose levels differed significantly between GDM cases and controls in the EPIDG cohort (Table 1).

Table 1 Characteristics of the study subjects in each cohort

Methylation risk score

In the mixed models regression, MRS for T2D levels were higher in women with GDM compared to non-GDM (adjusted model O.R: 1.4, 95%C.I: 1.10–1.74, P.val = 0.00014). For continuous traits, the MRS for T2D was significantly associated only with fasting glucose in the adjusted linear mixed models (P.val = 0.0086) (Table 2). When assessing EPIDG, EPIPREG Europeans and EPIPREG South Asians separately, MRS for T2D was associated with higher OR for GDM in EPIPREG South Asians, and followed the same direction of effect in EPIPREG Europeans and EPIDG although not statistically significant (Fig. 2; Supplementary Table 2).

Fig. 2
figure 2

Boxplot MRS for each group: red (SA_nonGDM)- EPIPREG south Asians non-GDM, Yellow (SA_GDM)- EPIPREG South Asians GDM, Green (EU_nonGDM)- EPIPREG Europeans non-GDM, Light Blue (EU_GDM): EPIPREG Europeans GDM, Dark Blue (M_nonGDM)- EPIDG (Mediterraneas) non-GDM, Pink (M_GDM)- EPIDG GDM

mQTLs for CpG sites common between T2DM and GDM

In the meta-analysis of the 42 CpG sites included in the MRS for T2DM, six CpGs were significantly associated with GDM (FDR < 0.05) (Table 3). All the identified CpG sites, except cg25130381, had significant mQTLs from lookups in GoDMC. A total of 23 mQTLs were found (supplementary Table 3). According to lookups in Phenoscanner, only two CpGs had mQTLs associated at p-value < 1 × 10− 5 with phenotypes related to immunoglobulin G, Rheumatoid arthritis, and body composition (Supplementary Table 4). However, when using a liberal threshold (p-value < 0.05), we observed some mQTLs possibly related to diabetes related phenotypes (Table 4).

Table 2 Association between the MRS for T2D and phenotypes in EPIPREG and EPIDG
Table 3 List of the 6 CpGs significantly associated with GDM across the three samples (EPIPREG_EU, EPIPREG_SA, EPIDG)
Table 4 mQTLs and disease associated, according to phenoscanner

Discussion

This is the first study to show that an MRS for T2DM is associated with GDM across our two cohorts, diagnostic criteria and ancestry. In addition, the MRS for T2DM was associated with higher fasting glucose levels. Six CpGs of the 42 included in the MRS, were significantly associated with GDM across EPIPREG Europeans, EPIPREG South Asians and EPIDG Mediterranean. We identified 23 mQTLs linked to the 6 CpG sites. Some of them were associated with T2DM and GDM with nominal significance in lookups of GWAS summary data.

GDM and T2DM share many pathophysiological factors, but the exact underlying mechanisms are largely unknown. There is a need to know how T2DM and GDM are linked to improve the prevention and thus avoid metabolic complications in the future. Genome-wide association studies (GWAS) have been used to investigate the potential link between GDM and T2DM [12], but none have so far tested whether GDM and T2DM share epigenetic marks. An increasing number of studies suggest that environmental factors are associated with epigenetic marks and that those marks may be important to understand interactions between genetics and environmental factors. A few EWAS of GDM in maternal peripheral blood are published so far, but the data is not sufficient for an MRS for GDM due to heterogeneity in diagnosis criteria used as well as their small sample sizes [28,29,30,31].

Our MRS was constructed from a meta-analysis of five European cohorts that identified CpGs associated with incident T2DM [18]. According to the bibliography of the genes included in the MRS, some of them were associated with glucose metabolism [32, 33], T2DM, type 1 diabetes mellitus and GDM [34]. ATP binding cassette subfamily G member 1 (ABCG1) is a gene involved in macrophage cholesterol and phospholipid transport and may regulate cellular lipid homeostasis. ABCG1 has been associated with diabetes and metabolic syndrome [35]. Furthermore, some studies show that GDM affects cholesterol homeostasis through this gene, therefore it could have a role in GDM and cardiovascular events [36].

The mQTLs were nominally associated with T2DM and GDM in lookups, but were interestingly also associated with phenotypes related to body composition, immunoglobulin G and autoimmune diseases associated with inflammation, such as rheumatoid arthritis and Crohn’s disease. These observations may suggest a pleiotropic relationship between these mQTLs for diabetes related CpG sites, but further exploration is necessary to assess the relevance and understand potential mechanisms.

The strengths of this study include the combination of two independent cohorts with case-control design in EPIDG which maximizes differences in GDM and non-GDM versus the population-based design in EPIPREG which includes the full range of values including less severe GDM. This combined with the two ancestries in EPIPREG, strengthen the evidence that these CpG sites may be important to understand the common epigenetic grounds of GDM and T2D. Limitations include limited statistical power in analysis of each sample separately and in associations between separate CpG sites and GDM. Still, the MRS for T2D was associated with GDM in EPIPREG South Asians, with the same direction of effect in EPIPREG Europeans and EPIDG. Further, an MRS for GDM based on an EWAS of GDM instead could have given a more precise prediction of GDM. Unfortunately published EWAS of GDM are still few and have small sample sizes which would result in a less robust MRS [40,41,42]. Also, our cross-sectional study design cannot entangle whether GDM status influence DNA methylation or vice versa. Hence this MRS should be tested in early pregnancy to be considered as potential predictor of GDM.

Genetic risk score reflects the inherited risk but have shown poor performance in complex disease such as diabetes. MRS is thought to reflect environmental triggers of the disease or phenotype and help to understand disease mechanism. Recent studies have shown the utility of both GRS and MRS for clinical prediction of different phenotypes [43]. Although GRS is being used in multiple diseases, their use present some limitations [44]. For example, GRS doesn’t reflect the effect of the environment on the phenotype whereas MRS incorporate both genetic factors and environmental exposures and their variation over time. Undoubtedly, combining both scores would enhance their clinical application for predicting or stratification of risk subjects. The MRS level differed between EPIPREG and EPIDG cohort, suggesting that the MRS may be population specific. Although we cannot rule out that part of these differences may be due to the batch or handling effect, they could also be a result of the EPIDG control group not having negative O’sullivan test, while non-GDM women in EPIPREG had normal fasting and 2 h glucose values (1 h glucose OGTT missing in EPIPREG). Further, the populations in EPIPREG and EPIDG differ largely in many aspects, including the severity of GDM.

To conclude, we are the first to show that MRS for T2DM is significantly associated with GDM, suggesting shared epigenetic mechanisms between GDM and T2DM. This may help explain some of the molecular mechanisms mediating the increased risk of developing T2DM after GDM. Future research should compare whether GRS and MRS or a combination of the two provide better predictions of GDM and T2DM and explore how transcription of the identified genes impact methylation on nearby genes. Furthermore, understanding the role of the genetic variants in disease development may help to improve prevention and management of both GDM and T2DM in the future.

Data availability

The data of this study are protected by Norwegian data protection law. Data, analytical methods and research materials supporting the results of this study may be reasonably requested from CS and SM.

Abbreviations

GDM:

Gestational Diabetes Mellitus

T2DM:

Type 2 Diabetes Mellitus

GRS:

Genetic Risk Score

MRS:

Methylation Risk Score

OGTT:

Oral Glucose Tolerance Test

QC:

Quality Control

mQTL:

Methylation Quantitative Trait Loci

References

  1. Khan MAB, Hashim MJ, King JK, Govender RD, Mustafa H, Kaabi JA. Epidemiology of Type 2 Diabetes – Global Burden of Disease and Forecasted Trends. J Epidemiol Glob Health. 2020 [cited 2022 Feb 2];10(1):107. Available from: https://www.pmc/articles/PMC7310804/.

  2. Diagnostic criteria and classification of hyperglycaemia first detected in pregnancy: a World Health Organization Guideline. Diabetes Res Clin Pract. 2014 [cited 2023 Feb 27];103(3):341–63. Available from: https://www.pubmed.ncbi.nlm.nih.gov/24847517/.

  3. Zhu Y, Zhang C. Prevalence of Gestational Diabetes and Risk of Progression to Type 2 Diabetes: a Global Perspective. Curr Diab Rep. 2016 [cited 2022 Nov 30];16(1):7. Available from: http://www.pmc/articles/PMC6675405/.

  4. Farrar D. Hyperglycemia in pregnancy: prevalence, impact, and management challenges. Int J Womens Health. 2016 [cited 2022 Nov 23];8:519–27. Available from: https://www.pubmed.ncbi.nlm.nih.gov/27703397/.

  5. Wang Z, Peng H, Gao W, Cao W, Lv J, Yu C, et al. Blood DNA methylation markers associated with type 2 diabetes, fasting glucose, and HbA1c levels: an epigenome-wide association study in 316 adult twin pairs. Genomics. 2021;113(6):4206–13.

    Article  CAS  PubMed  Google Scholar 

  6. Chambers JC, Loh M, Lehne B, Drong A, Kriebel J, Motta V et al. Epigenome-wide association of DNA methylation markers in peripheral blood from Indian Asians and Europeans with incident type 2 diabetes: a nested case-control study. lancet Diabetes Endocrinol. 2015 [cited 2022 Dec 22];3(7):526–34. Available from: https://pubmed.ncbi.nlm.nih.gov/26095709/.

  7. Cardona A, Day FR, Perry JRB, Loh M, Chu AY, Lehne B et al. Epigenome-wide association study of incident type 2 diabetes in a British population: EPIC-Norfolk study. Diabetes. 2019 [cited 2022 Dec 22];68(12):2315–26. Available from: https://pubmed.ncbi.nlm.nih.gov/31506343/.

  8. Iwata M, Maeda S, Kamura Y, Takano A, Kato H, Murakami S et al. Genetic Risk Score Constructed Using 14 Susceptibility Alleles for Type 2 Diabetes Is Associated With the Early Onset of Diabetes and May Predict the Future Requirement of Insulin Injections Among Japanese Individuals. Diabetes Care. 2012 [cited 2023 Mar 8];35(8):1763–70. Available from: https://diabetesjournals.org/care/article/35/8/1763/29901/Genetic-Risk-Score-Constructed-Using-14.

  9. Stanislawski MA, Litkowski E, Raghavan S, Harrall KK, Shaw J, Glueck DH et al. Genetic Risk Score for Type 2 Diabetes and Traits Related to Glucose-Insulin Homeostasis in Youth: The Exploring Perinatal Outcomes Among Children (EPOCH) Study. Diabetes Care. 2021 [cited 2023 Mar 8];44(9):2018–24. Available from: https://pubmed.ncbi.nlm.nih.gov/34257098/.

  10. Thompson M, Hill BL, Rakocz N, Chiang JN, Geschwind D, Sankararaman S et al. Methylation risk scores are associated with a collection of phenotypes within electronic health record systems. NPJ genomic Med. 2022 [cited 2022 Dec 2];7(1). Available from: https://pubmed.ncbi.nlm.nih.gov/36008412/.

  11. Schrader S, Perfilyev A, Ahlqvist E, Groop L, Vaag A, Martinell M et al. Novel Subgroups of Type 2 Diabetes Display Different Epigenetic Patterns That Associate With Future Diabetic Complications. Diabetes Care. 2022 [cited 2023 Feb 3];45(7):1621–30. Available from: https://pubmed.ncbi.nlm.nih.gov/35607770/.

  12. Pervjakova N, Moen GH, Borges MC, Ferreira T, Cook JP, Allard C et al. Multi-ancestry genome-wide association study of gestational diabetes mellitus highlights genetic links with type 2 diabetes. Hum Mol Genet. 2022 [cited 2023 Apr 26];31(19):3377–91. Available from: https://academic.oup.com/hmg/article/31/19/3377/6537590.

  13. Jenum AK, Mrøkrid K, Sletner L, Vange S, Torper JL, Nakstad B et al. Impact of ethnicity on gestational diabetes identified with the WHO and the modified International Association of Diabetes and Pregnancy Study Groups criteria: a population-based cohort study. Eur J Endocrinol. 2012 [cited 2023 Jan 16];166(2):317–24. Available from: https://eje.bioscientifica.com/view/journals/eje/166/2/317.xml.

  14. Fragoso-Bargas N, Opsahl JO, Kiryushchenko N, Böttcher Y, Lee-Ødegård S, Qvigstad E et al. Cohort profile: epigenetics in pregnancy (EPIPREG) – Population-based sample of European and south Asian pregnant women with epigenome-wide DNA methylation (850k) in peripheral blood leukocytes. PLoS ONE. 2021;16(8 August).

  15. Classification and diagnosis of diabetes mellitus and other categories of glucose intolerance. National Diabetes Data Group. Diabetes. 1979 [cited 2022 Jun 10];28(12):1039–57. Available from: https://pubmed.ncbi.nlm.nih.gov/510803/.

  16. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988 [cited 2023 Feb 6];16(3):1215. Available from: https://pubmed.ncbi.nlm.nih.gov/3344216/.

  17. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012 [cited 2022 Jun 13];13(1):1–16. Available from: https://bmcbioinformatics.biomedcentral.com/articles/https://doi.org/10.1186/1471-2105-13-86.

  18. Fraszczyk E, Spijkerman AMW, Zhang Y, Brandmaier S, Day FR, Zhou L et al. Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts. Diabetologia. 2022 [cited 2022 Dec 2];65(5):763–76. Available from: https://pubmed.ncbi.nlm.nih.gov/35169870/.

  19. Verschuren WMM, Blokstra A, Picavet HSJ, Smit HA. Cohort profile: the Doetinchem Cohort Study. Int J Epidemiol. 2008 [cited 2023 Sep 7];37(6):1236–41. Available from: https://pubmed.ncbi.nlm.nih.gov/18238821/.

  20. Raum E, Rothenbacher D, Löw M, Stegmaier C, Ziegler H, Brenner H. Changes of cardiovascular risk factors and their implications in subsequent birth cohorts of older adults in Germany: a life course approach. Eur J Cardiovasc Prev Rehabil. 2007 [cited 2023 Sep 7];14(6):809–14. Available from: https://pubmed.ncbi.nlm.nih.gov/18043304/.

  21. Holle R, Happich M, Löwel H, Wichmann HE. KORA–a research platform for population based health research. Gesundheitswesen. 2005 [cited 2023 Sep 7];67 Suppl 1(SUPPL. 1). Available from: https://pubmed.ncbi.nlm.nih.gov/16032513/.

  22. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest Package: Tests in Linear Mixed Effects Models. J Stat Softw. 2017 [cited 2022 Dec 2];82(13):1–26. Available from: https://www.jstatsoft.org/index.php/jss/article/view/v082i13.

  23. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010 [cited 2024 Jan 16];26(17):2190–1. Available from: https://pubmed.ncbi.nlm.nih.gov/20616382/.

  24. Chapter 10. : Analysing data and undertaking meta-analyses | Cochrane Training. [cited 2024 Jan 16]. Available from: https://training.cochrane.org/handbook/current/chapter-10.

  25. LDlink | An Interactive Web Tool for Exploring Linkage. Disequilibrium in Population Groups. [cited 2022 Dec 2]. Available from: https://ldlink.nci.nih.gov/?tab=snpclip.

  26. PhenoScanner. [cited 2022 Dec 2]. Available from: http://www.phenoscanner.medschl.cam.ac.uk/.

  27. Althouse AD. Adjust for Multiple Comparisons? It’s Not That Simple. Ann Thorac Surg. 2016 May 1 [cited 2023 Sep 7];101(5):1644–5. Available from: http://www.annalsthoracicsurgery.org/article/S0003497515018731/fulltext.

  28. Opsahl JO, Moen GH, Qvigstad E, Böttcher Y, Birkeland KI, Sommer C. Epigenetic signatures associated with maternal body mass index or gestational weight gain: a systematic review. J Dev Orig Health Dis. 2021 Jun 1 [cited 2023 May 15];12(3):373–83. Available from: https://pubmed.ncbi.nlm.nih.gov/32873364/.

  29. Wu P, Farrell WE, Haworth KE, Emes RD, Kitchen MO, Glossop JR et al. Maternal genome-wide DNA methylation profiling in gestational diabetes shows distinctive disease-associated changes relative to matched healthy pregnancies. Epigenetics. 2016;0.

  30. Dias S, Adam S, Rheeder P, Louw J, Pheiffer C. Altered Genome-Wide DNA Methylation in Peripheral Blood of South African Women with Gestational Diabetes Mellitus. Int J Mol Sci. 2019 [cited 2020 Apr 20];20(23):5828. Available from: http://www.ncbi.nlm.nih.gov/pubmed/31757015.

  31. Ballesteros M, Gil-Lluís P, Ejarque M, Diaz-Perdigones C, Martinez-Guasch L, Fernández-Veledo S et al. DNA methylation in gestational diabetes and its predictive value for postpartum glucose disturbances. J Clin Endocrinol Metab. 2022 [cited 2022 Aug 9]; Available from: https://pubmed.ncbi.nlm.nih.gov/35914803/.

  32. Cardona A, Day FR, Perry JRB, Loh M, Chu AY, Lehne B et al. Epigenome-wide association study of incident type 2 diabetes in a British population: EPIC-Norfolk study. Diabetes. 2019 [cited 2023 Feb 3];68(12):2315–26. Available from: https://pubmed.ncbi.nlm.nih.gov/31506343/.

  33. Meeks KAC, Henneman P, Venema A, Addo J, Bahendeka S, Burr T et al. Epigenome-wide association study in whole blood on type 2 diabetes among sub-Saharan African individuals: findings from the RODAM study. Int J Epidemiol. 2019 [cited 2023 Feb 3];48(1):58–70. Available from: https://pubmed.ncbi.nlm.nih.gov/30107520/.

  34. Tobi EW, Juvinao-Quintero DL, Ronkainen J, Ott R, Alfano R, Canouil M et al. Maternal Glycemic Dysregulation During Pregnancy and Neonatal Blood DNA Methylation: Meta-analyses of Epigenome-Wide Association Studies. Diabetes Care. 2022 [cited 2022 Nov 30];45(3):614–23. Available from: https://pubmed.ncbi.nlm.nih.gov/35104326/.

  35. Qie R, Chen Q, Wang T, Chen X, Wang J, Cheng R et al. Association of ABCG1 gene methylation and its dynamic change status with incident type 2 diabetes mellitus: the Rural Chinese Cohort Study. J Hum Genet. 2021 [cited 2023 Feb 9];66(4):347–57. Available from: https://pubmed.ncbi.nlm.nih.gov/32968204/.

  36. Tian M, Zhang C, Han H, Liu R, Gao Z, Chen P et al. Effects of the preassembly of benzohydroxamic acid with Fe (III) ions on its adsorption on cassiterite surface. Miner Eng. 2018 [cited 2018 Nov 13];127:32–41. Available from: https://www.sciencedirect.com/science/article/pii/S0892687518303352.

  37. MAP4K2 mitogen-activated protein. kinase kinase kinase kinase 2 [Homo sapiens (human)] - Gene - NCBI. [cited 2023 Feb 9]. Available from: https://www.ncbi.nlm.nih.gov/gene?Db=gene&Cmd=DetailsSearch&Term=5871.

  38. PDGFB platelet derived growth. factor subunit B [Homo sapiens (human)] - Gene - NCBI. [cited 2023 Feb 13]. Available from: https://www.ncbi.nlm.nih.gov/gene/5155#gene-ontology.

  39. Shan Z, Xu C, Wang W, Li W. Enhanced PDGF signaling in gestational diabetes mellitus is involved in pancreatic β-cell dysfunction. Biochem Biophys Res Commun. 2019;516(2):402–7.

    Article  CAS  PubMed  Google Scholar 

  40. Linares-Pineda T, Peña-Montero N, Fragoso-Bargas N, Gutiérrez-Repiso C, Lima-Rubio F, Suarez-Arana M et al. Epigenetic marks associated with gestational diabetes mellitus across two time points during pregnancy. Clin Epigenetics. 2023 [cited 2023 Jul 7];15(1):110. Available from: https://pubmed.ncbi.nlm.nih.gov/37415231/.

  41. Linares-Pineda TM, Peña-Montero N, Gutiérrez-Repiso C, Lima-Rubio F, Sánchez-Pozo A, Tinahones FJ et al. Epigenome wide association study in peripheral blood of pregnant women identifies potential metabolic pathways related to gestational diabetes. Epigenetics. 2023 [cited 2023 May 30];18(1):2211369. Available from: https://pubmed.ncbi.nlm.nih.gov/37192269/.

  42. Moen GH, Sommer C, Prasad RB, Sletner L, Groop L, Qvigstad E et al. MECHANISMS IN ENDOCRINOLOGY: Epigenetic modifications and gestational diabetes: a systematic review of published literature. Eur J Endocrinol. 2017 [cited 2024 Jan 17];176(5):R247–67. Available from: https://pubmed.ncbi.nlm.nih.gov/28232369/.

  43. Nabais MF, Gadd DA, Hannon E, Mill J, McRae AF, Wray NR. An overview of DNA methylation-derived trait score methods and applications. Genome Biol 2023 241. 2023 [cited 2024 Jan 12];24(1):1–23. Available from: https://genomebiology.biomedcentral.com/articles/https://doi.org/10.1186/s13059-023-02855-7.

  44. Abu-El-Haija A, Reddi HV, Wand H, Rose NC, Mori M, Qian E et al. The clinical application of polygenic risk scores: A points to consider statement of the American College of Medical Genetics and Genomics (ACMG). Genet Med. 2023 [cited 2024 Jan 12];25(5). Available from: https://pubmed.ncbi.nlm.nih.gov/36920474/.

Download references

Acknowledgements

The authors thank the women who participated in the STORK Groruddalen study and EPIDG study.

Funding

EPIPREG is supported by the South Eastern Norway Regional Health Authority (grant number: 2019092), and the Norwegian Diabetes Association (grant number: N/A). GH-M. is the recipient of an Australian Research Council Discovery Early Career Award (Project number: DE220101226) funded by the Australian Government and supported by the Research Council of Norway (Project grant: 325640). JO-O is supported by the Norwegian Research Council (Medical Student Research Program, grant 271555/F20).

EPIDG was funded by the Juan Rodés program from “Instituto de Salud Carlos III” (JR20-00040 to MM-V), PFIS program (FI19/00178 to TML-P), and the Nicolas Monardes Program from the “Servicio Andaluz de Salud, Junta de Andalucía”, Spain (RC-0008-2021 to SM). In addition, this study was supported by the research grants from the ISCIII (PI18/01175, PI21/01864) and from “Servicio Andaluz de Salud”, Junta de Andalucía (PI-0283-2018, PI-0419-2019). This study has been co-funded by FEDER funds (“A way to make Europe”).

Author information

Authors and Affiliations

Authors

Contributions

TML-P: methodology, analysis, interpretation of the data and wrote original draft; NB-F: methodology and analysis; MJP: EPIDG funding acquisition and acquired data from the patients and collected samples; MM-V: EPIDG acquired data from the patients and collected samples; AKJ: designed the STORK project; LS: contributed with data acquisition in STORK-G ; SL-Ø: curated EPIPREG data ; JO-O: contributed to the manuscript ;GH-M: contributed with the imputation in EPIPREG data ; EQ: acquired genotype and DNA methylation data in STORK G. ; RB-P: realized the wer laboratory experiments regarding the methylation in EPIPREG; KI-B: participated in the design of the STORK-G, acquired genotype and DNA methylation data in STORK G. ; SM: EPIDG funding acquisition, interpretation of the data and has substantively revised the work, CS: acquired genotype and DNA methylation data in STORK G and she is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Corresponding author

Correspondence to Christine Sommer.

Ethics declarations

Ethics approval and consent to participate

The Norwegian Regional Committee for Medical Health Research Ethics South East approved the STORK-G study including genetic and epigenetic data (ref. number 2015/1035). And the EPIDG study was approved by the Institutional review board at the Hospital Universitario Virgen de la Victoria, Málaga. Written informed consent from all participants was obtained before any study-related procedure.

Consent for publication

All authors revised the manuscript critically and approver the final version.

Conflict of interest

The authors declare there is no conflict of interest.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Supplementary figure 1.

Work flow for sample selection A) in EPIPREG cohort, and B) EPIDG cohort

Supplementary Material 2: Supplementary table 1.

List of the 42 CpGs used in MRS from meta-analysis of 5 Europeans cohorts

Supplementary Material 3: Supplementary table 2.

Logistic model to associate MRS with GDM in each group separately

Supplementary Material 4: Supplementary table 3.

List of the 5 CpGs with their mQTLS

Supplementary Material 5: Supplementary table 4.

Phenotypes associated with 2 CpGs with mQTLs higher p-value

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Linares-Pineda, T.M., Fragoso-Bargas, N., Picón, M.J. et al. DNA methylation risk score for type 2 diabetes is associated with gestational diabetes. Cardiovasc Diabetol 23, 68 (2024). https://doi.org/10.1186/s12933-024-02151-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-024-02151-z

Keywords