Associations of maternal early-pregnancy blood glucose and insulin concentrations with DNA methylation in newborns
Clinical Epigenetics volume 12, Article number: 134 (2020)
Intrauterine exposure to a disturbed maternal glucose metabolism is associated with adverse offspring outcomes. DNA methylation is a potential mechanism underlying these associations. We examined whether maternal early-pregnancy glucose and insulin concentrations are associated with newborn DNA methylation. In a population-based prospective cohort study among 935 pregnant women, maternal plasma concentrations of non-fasting glucose and insulin were measured at a median of 13.1 weeks of gestation (95% range 9.4–17.4). DNA methylation was measured using the Infinium HumanMethylation450 BeadChip (Ilumina). We analyzed associations of maternal early-pregnancy glucose and insulin concentrations with single-CpG DNA methylation using robust linear regression models. Differentially methylated regions were analyzed using the dmrff package in R. We stratified the analyses on normal weight versus overweight or obese women. We also performed a look-up of CpGs and differently methylated regions from previous studies to be associated with maternal gestational diabetes, hyperglycemia or hyperinsulinemia, or with type 2 diabetes in adults.
Maternal early-pregnancy glucose and insulin concentrations were not associated with DNA methylation at single CpGs nor with differentially methylated regions in the total group. In analyses stratified on maternal BMI, maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG (cg03617420, XKR6) among normal weight women and at another (cg12081946, IL17D) among overweight or obese women. No stratum-specific associations were found for maternal early-pregnancy insulin concentrations. The two CpGs were not associated with birth weight or childhood glycemic measures (p values > 0.1). Maternal early-pregnancy insulin concentrations were associated with one CpG known to be related to adult type 2 diabetes. Enrichment among nominally significant findings in our maternal early-pregnancy glucose concentrations was found for CpGs identified in a previous study on adult type 2 diabetes.
Maternal early-pregnancy glucose concentrations, but not insulin concentrations, were associated with DNA methylation at one CpG each in the subgroups of normal weight and of overweight or obese women. No associations were present in the full group. The role of these CpGs in mechanisms underlying offspring health outcomes needs further study. Future studies should replicate our results in larger samples with early-pregnancy information on maternal fasting glucose metabolism.
The prevalence of gestational diabetes is rising worldwide and has been reported to complicate up to 25% of all pregnancies [1, 2]. This rise is partly due to the increasing prevalence of obesity among women of reproductive age and depends on screening tools and diagnostic criteria [2,3,4]. Intrauterine exposure to maternal gestational diabetes or impaired glucose tolerance measured in mid-pregnancy and late pregnancy is associated with increased risks of adverse maternal and fetal perinatal outcomes and of diabetes and obesity in the offspring [5,6,7,8]. These associations of increased risks on perinatal outcomes are already present for higher maternal glucose concentrations below the threshold of gestational diabetes [9, 10]. Additionally, the associations are stronger among women who are overweight or obese at the start of their pregnancy [2, 3]. Women who develop hyperglycemia and gestational diabetes may already have suboptimal glucose metabolism earlier in pregnancy. The first trimester of pregnancy is a critical period for embryonic and placental growth and development . As such, impaired maternal glucose metabolism may already exert negative effects in that early stage. Early-pregnancy glucose metabolism has been shown to be associated with altered fetal growth, adverse birth outcomes, and childhood glucose metabolism, but not with other childhood cardiometabolic outcomes after adjusting for maternal pre-pregnancy BMI [9, 12, 13]. Thus, early pregnancy may be an important time window for the effects of suboptimal maternal glycemic measures and as such an influential period for future interventions.
The mechanisms underlying these associations are unknown. DNA methylation has been suggested as a potential mechanism linking adverse exposures during pregnancy and impaired offspring health [14, 15]. Previous studies using candidate-gene approaches suggested that maternal gestational diabetes is associated with epigenetic modifications in placenta and cord blood at loci relevant to growth, energy homeostasis, and diabetes mellitus [14, 16,17,18]. Epigenome-wide association studies (EWAS) of gestational diabetes or maternal glucose concentrations showed varying results, with no clear pattern of associations [1, 15, 19,20,21,22,23,24,25]. The inconsistent results of candidate-gene studies and EWAS may be due to differences in study design. The studies varied in their exposure definition: gestational diabetes as binary exposure or glucose concentrations after an oral challenge test, in the tissues in which DNA methylation was measured: placenta or blood, and in the extent of adjustment for covariates, with most not adjusting for cell-type heterogeneity. Also, the majority had limited sample sizes [15, 19, 21,22,23,24,25]. It is not known whether maternal glucose and insulin concentrations across the full range in early pregnancy are associated with cord blood DNA methylation and whether these associations differ between normal weight versus overweight or obese women. Insight into these associations and their underlying mechanisms is important, as maternal blood glucose metabolism can be a target for preventive interventions to improve child health outcomes.
We hypothesized that maternal early-pregnancy glucose and insulin concentrations are associated with offspring DNA methylation at birth. Therefore, we conducted an epigenome-wide association analysis in a population-based prospective cohort study, with maternal glucose and insulin concentrations measured at a median of 13.1 weeks of gestation (95% range 9.4–17.4). As a secondary analysis, we stratified on maternal body mass index (BMI) categories to observe if maternal BMI modifies the studied associations. We also examined whether differential DNA methylation at any CpGs found to be associated with maternal glucose or insulin concentrations in cord blood persisted in peripheral blood of 10-year-old children. To obtain further insight into the potential significance of the observed DNA methylation changes, we conducted exploratory analyses on the associations of identified CpGs with offspring health outcomes. We also performed a look-up in our results of CpGs identified to be associated with maternal glucose metabolism during pregnancy or with type 2 diabetes in adults in previous literature.
The population characteristics for the total group and stratified on maternal BMI in women with normal weight versus women with overweight or obesity are shown in Table 1. The mean maternal early-pregnancy glucose concentration was 4.4 mmol/l (standard deviation 0.8) and the median maternal early-pregnancy insulin concentration was 126.3 pmol/l (95% range 19.9–774.6). Gestational diabetes was diagnosed in 12 (1.3%) women and 59 (6.3%) women were obese. Non-response analyses showed that mothers without data on early-pregnancy glucose and insulin measurements delivered more often female children (Additional file 1: Table S1).
Associations of maternal early-pregnancy glucose and insulin concentrations with DNA methylation at birth
After Bonferroni (p value cutoff < 1.0 × 10−7) or false-discovery rate (FDR) correction, we did not observe associations of maternal early-pregnancy glucose or insulin concentrations with offspring DNA methylation in cord blood. These models were adjusted for gestational age at glucose/insulin measurement, maternal age, educational level, parity, smoking, pre-pregnancy BMI, child sex, cell-type proportions, and batch. The results of both analyses are presented in Additional file 2: Figure S1a and Figure S1b. The CpGs with p values < 1.0 × 10−4 for both analyses are shown in Additional file 3: Table S2 and Table S3. A model without adjustment for maternal pre-pregnancy BMI showed largely similar results (Additional file 4: Table S4 and Table S5).
The analyses stratified on maternal BMI showed that among normal weight women, maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG (cg03617420 in XKR6; effect estimate = 7.3 × 10−3 (standard error (SE) 1.3 × 10−3), p value = 7.4 × 10−9) (Fig. 1a). DNA methylation at this CpG was not significantly associated with glucose concentrations in overweight or obese women (effect estimate = − 2.4 × 10−3 (SE 3.2 × 10−3), p value = 0.46). Among overweight or obese women maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG (cg12081946 in IL17D; effect estimate = − 3.4 × 10−2 (SE 5.6 × 10−3), p value = 8.9 × 10−10) (Fig. 1b). DNA methylation at this CpG was not significantly associated with glucose concentrations in normal weight women (effect estimate = 3.1 × 10−3 (SE 4.1 × 10−3), p value = 0.45). Maternal early-pregnancy insulin concentrations were not associated with DNA methylation in cord blood in normal weight and in overweight or obese women (Additional file 5: Figure S2a and Figure S2b).
Neither maternal early-pregnancy glucose nor insulin concentrations were associated with differentially methylated regions in cord blood, analyzed using the dmrff package in R . The differentially methylated regions with p values < 1.0 × 10−4 are shown in Additional file 6: Table S6 and Table S7.
Maternal glucose concentrations were not associated with DNA methylation levels at the two identified CpGs in peripheral blood of 10-year-old children (p values > 0.1) (Additional file 7: Table S8). In exploratory analyses, the two identified CpGs were not associated with birth weight or childhood glucose concentrations, which both were previously found to be associated with maternal early-pregnancy glucose concentrations in our data (p values > 0.1) [9, 13].
Look-up of maternal glucose metabolism and adult type 2 diabetes-associated CpGs
In a look-up in our results of CpGs and DMRs identified in previous studies to be associated with maternal glucose metabolism or with adult type 2 diabetes, we found that one CpG, cg1680945 at MDN1, known to be related to adult type 2 diabetes was also significantly associated with maternal early-pregnancy insulin concentrations (effect estimate = − 3.3 × 10−3 (SE 1.1 × 10−3), p value = 2.2 × 10−3) . The look-up of other previously described CpGs and DMRs in the maternal early-pregnancy glucose and insulin EWAS results showed no associations (Additional file 8: Table S9; Additional file 9: Table S10; Additional File 10: Table S11; Additional file 11: Table S12) [1, 14, 17, 18, 22, 23, 25, 27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42]. We found enrichment for findings from one previous study on adult type 2 diabetes among the 24,935 nominally significant CpGs from the maternal early-pregnancy glucose EWAS results (Fisher combined probability p value = 0.04) . No evidence for enrichment of the CpGs from other previous studies among the 24,935 nominally significant CpGs from the maternal early-pregnancy glucose cord blood analysis, nor among the 19,418 nominally significant CpGs from the maternal early-pregnancy insulin cord blood analysis was found (lowest Fisher combined probability p value = 0.15 in maternal early-pregnancy glucose EWAS results and p value = 0.12 in insulin EWAS results) [1, 14, 17, 18, 22, 23, 25, 27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42].
In this population-based EWAS we did not observe associations of maternal early-pregnancy glucose and insulin concentrations across the full spectrum with offspring cord blood DNA methylation in the full group. However, after stratification on maternal BMI, maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG each among normal weight and among overweight or obese women. Associations of DNA methylation at these CpG sites did not persist in 10-year-old children. Also, we did not find associations with offspring health outcomes. Maternal early-pregnancy insulin was associated with one CpG known from a previous adult type 2 diabetes-associated study. Also, we found enrichment of CpGs identified in a previous EWAS on adult type 2 diabetes among our maternal early-pregnancy glucose EWAS results. Overall, our results constitute a first step toward a better understanding of a potential role of DNA methylation underlying the associations of maternal glycemic traits in early pregnancy with offspring health outcomes.
Interpretation of main findings
Gestational diabetes or impaired glucose tolerance diagnosed in the second half of pregnancy increases the risks of adverse birth outcomes, of obesity and diabetes in the offspring [5, 6, 8]. It has been suggested that women who develop gestational diabetes or hyperglycemia later in pregnancy already have suboptimal glucose metabolism before or in early pregnancy [43, 44]. Maternal glycemic measures in early pregnancy have been described to be associated with altered fetal growth and glucose metabolism in childhood, but not with child adiposity, lipid levels, and blood pressure after adjustment for maternal pre-pregnancy BMI [9, 12, 13]. Thus, early pregnancy may already be a critical period for the effects of maternal glucose concentrations on offspring birth outcomes and glycemic health in childhood. The associations of maternal glucose metabolism with offspring outcomes may be explained by differential DNA methylation. Therefore, we hypothesized that maternal early-pregnancy glucose and insulin concentrations are associated with offspring DNA methylation levels at birth and that these associations may be different for normal weight and overweight women.
Results from a recent meta-analysis from seven pregnancy cohorts among 3677 mother-newborn pairs showed that gestational diabetes was not associated with differential methylation at a single CpG level, but it was associated with lower cord blood methylation levels within two specific regions . In the current population-based EWASs, we did not find any associations in the full group, but maternal early-pregnancy glucose concentrations were associated with DNA methylation at cg03617420 (XKR6) among normal weight women, and at cg12081946 (IL17D) among overweight or obese women. The effect estimates of both CpGs were in opposite directions for normal weight and overweight or obese women, which could imply a modifying effect of maternal BMI, as we hypothesized.
XKR6, XK-related 6 gene, is located on chromosome 8 and is classified as a member of the Kell blood group complex subunit-related family. Genetic variants in this gene have previously been associated with type 2 diabetes, lipid concentrations, systolic blood pressure, and kidney function, among others [45,46,47,48]. XKR6 is expressed in many tissues, most strongly in testis and lymphocytes, but also in the cerebellum and pancreas, among others. IL17D, interleukin 17D, is part of the cytokine family and located on chromosome 13 and has been previously associated with autoimmune and inflammatory diseases . Autoimmune processes are part of the pathogenesis of type 1 diabetes mellitus . Genetic variants close to IL17D have been associated with PR segment on electrocardiogram . IL17D is most strongly expressed in the brain and skeletal muscle. Based on Roadmap Epigenomics Data Complete Collection extracted from the UCSC Genome Browser, both cg03617420 and cg12081946 coincide with DNAseI hypersensitivity clusters and transcription factor binding regions, indicating a location in potential regulatory elements.
DNA methylation levels at these CpGs have not been previously described in relation to maternal early-pregnancy glucose concentrations, and our results need replication in larger groups. We found that DNA methylation at one CpG, cg1680945 (MDN1), which was associated with adult type 2 diabetes in a previous study, was also associated with maternal early-pregnancy insulin concentrations . We also found enrichment for CpGs identified in a previous EWAS on adult type 2 diabetes among the nominally significant CpGs from the maternal early-pregnancy glucose EWAS . Maternal early-pregnancy glucose and insulin concentrations were not associated with any of the other previously reported maternal gestational diabetes, hyperglycemia or hyperinsulinemia, or adult type 2 diabetes-associated CpGs or differently methylated regions [1, 14, 17, 18, 22, 23, 25, 28,29,30,31,32,33,34,35,36,37,38,39, 41, 42].
The lack of identified associations in our total study group may have multiple reasons. Our study population is relatively healthy with on average lean women, with limited variability in maternal early-pregnancy glucose and insulin concentrations and with a low percentage of women who developed gestational diabetes. Associations of maternal early-pregnancy glucose and insulin concentrations with DNA methylation may be more apparent when using fasting glucose and insulin concentrations or among higher risk populations, as observed in studies in women with gestational diabetes [22, 31, 34]. Besides this, the moderate sample size of this study also limits the power to detect smaller differences. Another possibility is that associations of maternal early-pregnancy glucose and insulin concentrations with DNA methylation in offspring are more apparent in other tissues than cord blood, such as the placental tissue, body fat, skeletal muscle, liver, or pancreas. As shown in the analysis stratified on BMI, for some CpGs, associations may not be apparent in the full group as the directions of effect may be opposite between certain subgroups in the study population.
Maternal glucose concentrations were not associated with DNA methylation in child peripheral blood at age 10 years at the two CpGs identified in cord blood. This could imply a temporary effect of maternal glucose metabolism on offspring DNA methylation or it may be due to the relatively small sample size at age 10 years. Further studies should repeat these analyses among different age-groups in children to replicate our findings at birth and explore persistence of differential DNA methylation at the two identified CpGs. Also, the two identified CpGs were not associated with birth weight or childhood glucose concentrations in exploratory analyses. This may be due to the relatively small sample size or it may indicate that DNA methylation at these sites does not represent a biological pathway linking maternal glucose levels to birth weight or childhood glucose concentrations. Further studies, including studies among high-risk populations, are needed to examine the exact pathways involved in the associations of maternal glycemic measures during pregnancy with adverse birth outcomes and with offspring health outcomes, such as diabetes at later ages, in more detail.
This study suggests that maternal glycemic traits are associated with DNA methylation and that these associations may differ between mothers with overweight/obesity and those without. The role of the identified differential DNA methylation in pathways to offspring health needs further study. This is a first step toward discovering the underlying biological pathways and, if confirmed, it emphasizes the first trimester being a potentially important window of pregnancy for intervention studies to improve child health outcomes. Further, larger studies, with maternal early-pregnancy fasting blood samples, are needed to replicate our results.
Major strengths of this study are the population-based prospective design and the fact that we have information on maternal plasma glucose and insulin concentrations in early pregnancy in combination with cord blood DNA methylation. In addition to single-CpG analyses, differential methylated regions were also evaluated. We were able to adjust for a large number of potential confounders and for estimated cell-type proportions. The relatively small number of mothers with gestational diabetes (1.3% versus 2–5% in the general Dutch population ) may be due to the fact that information on gestational diabetes was taken from medical records and there was no structural testing of all pregnant women. However, the low number of mothers with gestational diabetes and obesity included in the sample may also indicate a selection toward a healthy, non-diabetic and lean population that might influence the generalizability of our findings and may have limited our statistical power to detect significant associations. Glucose and insulin concentrations were measured once during early pregnancy. Future studies are needed to measure maternal glucose and insulin concentrations at multiple time points during pregnancy to observe whether normal glucose and insulin concentrations in early pregnancy will worsen or maintain normal during pregnancy and whether patterns of glucose and insulin concentrations during pregnancy may be more informative than single measurements. The blood samples in the study are non-fasting. They were collected after a fasting time of at least 30 min. Since glucose and insulin concentrations vary during the day and are sensitive to carbohydrate intake, this may have led to non-differential misclassification of glucose and insulin concentrations. However, it has been suggested that maternal non-fasting glucose concentrations may better reflect the normal physiological state in pregnancy [10, 53]. Ideally, data on oral glucose tolerance tests would have been included, but these are not available in the Generation R Study. Blood samples were collected and processed in a standardized way, but time from sampling to freezing could be up to four hours. This may have affected the measured glucose concentrations. DNA methylation was measured in blood, which may differ from methylation patterns in other tissues. The study participants are of European ancestry and therefore, the findings might not be generalizable to other populations.
Maternal early-pregnancy blood glucose and insulin concentrations were not associated with differential DNA methylation at birth in the full group. However, maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG in XKR6 among normal weight women and at another CpG in IL17D among overweight or obese women. Their role in mechanisms underlying offspring health outcomes needs further study. These results await confirmation by future studies in larger samples with early-pregnancy information on maternal fasting glucose metabolism and exploring potential tissue-specific methylation effects in the offspring.
This study was embedded in the Generation R Study, a population-based prospective cohort from early fetal life onwards, based in Rotterdam, the Netherlands . The study has been approved by the Medical Ethical Committee of the Erasmus MC, University Medical Center Rotterdam (MEC 198.782/2001/31). Written informed consent was obtained for all participants . In total, 8879 women were enrolled during pregnancy (response rate at baseline: 61%), of whom, 6186 had measurements of glucose and insulin concentrations available. DNA methylation was measured in cord blood of a randomly selected European-ancestry subset of n = 1396 mothers. Out of these mothers, n = 945 had measurements on early-pregnancy glucose metabolism available. We excluded women with pre-existing diabetes (n = 3); twin pregnancies and in case of multiple (non-twin) children per mother, we excluded one of each sibling pair, based on data completeness or, if equal, randomly (n = 7). The population for analysis of this study comprised 935 mother-newborn pairs (Fig. 2).
Maternal glucose and insulin concentrations assessment
Blood samples were collected once in early pregnancy at 13.1 median weeks’ gestation (95% range 9.4–17.4), as described previously. After at least 30 min of fasting, venous blood samples were collected from pregnant women by research nurses and temporally stored at room temperature. Samples were minimally 30 min fasting. As we did not have information on the exact fasting duration, we consider all samples random. The time-interval of 30 min was chosen because of the design of the study, in which it was not possible to obtain fasting samples from all pregnant women. At least every 3 h, blood samples were transported to a dedicated laboratory facility of the regional laboratory in Rotterdam, the Netherlands (Star-MDC), for further processing and storage . Glucose (mmol/l) was measured with the c702 module on a Cobas 8000 analyzer. Insulin (pmol/l) was measured with electrochemiluminescence immunoassay on a Cobas e411 analyzer. Quality control samples demonstrated intra- and inter-assay CVs of 0.9% and 1.2% for glucose concentrations and of 1.3% and 2.5% for insulin concentrations, respectively.
Information on pre-existing diabetes was obtained from self-reported questionnaires and on gestational diabetes from medical records after delivery. Gestational diabetes was diagnosed by a community midwife or an obstetrician according to Dutch midwifery and obstetric guidelines at the time of inclusion into the study, using the following criteria: either a random glucose concentration > 11.0 mmol/l, a fasting glucose ≥ 7.0 mmol/l, or a fasting glucose between 6.1 and 6.9 mmol/l with a subsequent abnormal glucose tolerance test .
DNA was extracted from cord blood using the salting-out method. Five hundred nanograms of DNA per sample underwent bisulfite conversion using the EZ-96 DNA Methylation kit (Shallow) (Zymo Research Corporation, Irvine, CA, USA). Samples were plated randomly onto 96-well plates. Samples were processed with the Illumina Infinium HumanMethylation450 (450 k) BeadChip (Illumina Inc., San Diego, CA, USA). Quality control of analyzed samples was performed using standardized criteria. Quality control and normalization of the array data was performed according to the Control Probe Adjustment and reduction of global CORrelation (CPACOR) workflow using R [56, 57]. Probes that had a detection p value ≥ 1E−16 were set to missing per array. Next, the intensity values were quantile normalized for each of the six probe-type categories separately: type II red/green, type I methylated red/green, and type I unmethylated red/green. Beta values were calculated as proportion of methylated intensity value to the sum of (methylated and unmethylated intensities plus 100). Arrays with observed technical problems such as failed bisulfite conversion, hybridization or extension, and arrays with a sex mismatch were removed from subsequent analyses. Additionally, only arrays with a call rate > 95% per sample were processed further. Probes on the X and Y chromosomes were excluded from the analyses. Outlying methylation beta values were excluded using the following method: values < (25th percentile – 3*interquartile range (3IQR)) and values > (75th percentile +3IQR) were removed . For each analysis, we excluded plates with fewer than 3 samples because of convergence issues. This did not lead to exclusions in the analyses of the full group, but led to exclusion of 4 and 14 participants, in the normal weight and overweight/obese stratum, respectively. For all CpGs and differentially methylated regions, the official gene name of the nearest gene was noted using Illumina’s annotation information and we enhanced the annotation provided by Illumina with the UCSC Genome Browser build hg19 using the CpG location [59, 60].
Information on maternal age, pre-pregnancy weight, educational level, and parity was obtained from questionnaires at enrolment . Maternal smoking during pregnancy was assessed by questionnaires in pregnancy. We measured maternal height at enrolment without shoes and heavy clothing. Pre-pregnancy BMI was calculated (self-reported pre-pregnancy weight in kilograms divided by height measured at enrolment in meters, squared). Information on gestational age at birth, child sex, and birth weight was obtained from medical records. To adjust for batch effects, plate number was included as a covariate in the analyses. We estimated leukocyte subtypes using a cord blood-specific reference . This method estimates the relative proportions of six white blood cell subtypes (CD4+ T lymphocytes, CD8+ T lymphocytes, natural killer cells, B lymphocytes, monocytes, and granulocytes) and nucleated red blood cells.
First, non-response analysis was conducted among participants with singleton children and information available on cord blood DNA methylation, comparing participants with to those without data on maternal early-pregnancy glucose metabolism available, using Student’s t tests, Mann-Whitney tests, and chi-square tests. Second, we used robust linear regression models in an EWAS framework to assess the associations of maternal early-pregnancy glucose and insulin concentrations with single-CpG DNA methylation in cord blood . Maternal early-pregnancy insulin had a skewed distribution and was natural log-transformed for the analyses. The analyses were performed in two models: first model—adjusted for gestational age at glucose/insulin measurement, maternal age at intake, educational level, parity, smoking, child sex, cell-type proportions, and batch; and a secondary (main) model additionally adjusted for pre-pregnancy BMI. Since maternal obesity enhances the effect of higher glucose and insulin concentrations on adverse offspring outcomes, the effect on DNA methylation could be modified by maternal BMI. Therefore, as a secondary analysis, we stratified women into two strata of normal weight and overweight or obese women and repeated the main linear regression models in these strata. Included covariates were based on previous studies and factors known to be strongly associated with DNA methylation [22, 63,64,65]. Multiple testing was accounted for using Bonferroni correction, with CpGs with a p value < 1.0 × 10−7 considered significant. Additionally, we planned a priori to also report results using FDR correction for multiple testing, using the method by Benjamini and Hochberg . Third, we identified differentially methylated regions using the dmrff package (https://github.com/perishky/dmrff), which identifies differentially methylated regions by combining EWAS summary statistics from nearby CpGs . Significant differentially methylated regions were defined as regions spanning a set of CpG sites with at most 500 bp between consecutive sites with nominal EWAS p values < 0.05 and effect estimates with the same direction. Fourth, we examined whether associations of any CpGs identified in cord blood persisted in peripheral blood of 10-year-old children, using the main model additionally adjusted for child age at measurement. Missing covariate data were multiple-imputed using the Markov chain Monte Carlo method. All analyses were performed using R version 3.4.3 . All authors had access to the study data and reviewed and approved the final manuscript.
Associations of identified CpGs with offspring health outcomes
Exploratory analyses were performed in the relevant strata of maternal BMI to examine associations of identified CpGs with offspring birth weight and childhood glucose concentrations measured at age 10 years. We chose these outcomes because we have previously found them to be associated with maternal early-pregnancy glucose concentrations [9, 13]. We ran linear regression models using gestational age and sex-adjusted birth weight standard deviation scores (SDS) and childhood glucose concentrations as outcomes. Birth weight SDS were calculated based on the Niklasson reference charts, using Growth Analyzer (version 3.5; Dutch Growth Research Foundation, Rotterdam, the Netherlands) . Models were adjusted for maternal age, educational level, parity, and smoking, as well as for plate number and the seven cell types from the cord blood reference . Childhood glucose models were additionally adjusted for child sex and age at glucose measurement.
Look-up of previously identified CpGs
We performed a look-up in our maternal early-pregnancy glucose and insulin results of previously described maternal glucose metabolism-associated CpGs and DMRs and of previously described adult type 2 diabetes-associated CpGs in adults. The PubMed search terms are described in Additional file 12: Note 1. We took those studies into account that (1) included more than 50 participants in total, (2) measured DNA methylation in cord blood or peripheral blood, (3) were epigenome-wide studies or candidate-gene studies, and (4) reported p values for single CpGs or differential methylated regions. Significance was determined based on a Bonferroni corrected p value < 0.05/number of tested CpGs per reference study. We also evaluated enrichment of these CpGs among CpGs with a p value < 0.05 in our EWAS results using a hypergeometric test.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Body mass index
Differentially methylated region
Epigenome-wide association study
Howe CG, Cox B, Fore R, Jungius J, Kvist T, Lent S, et al. Maternal gestational diabetes mellitus and newborn DNA methylation: findings from the pregnancy and childhood epigenetics consortium. Diabetes Care. 2020;43(1):98–105.
Poston L, Caleyachetty R, Cnattingius S, Corvalan C, Uauy R, Herring S, et al. Preconceptional and maternal obesity: epidemiology and health consequences. Lancet Diabetes Endocrinol. 2016;4(12):1025–36.
Catalano PM, McIntyre HD, Cruickshank JK, McCance DR, Dyer AR, Metzger BE, et al. The hyperglycemia and adverse pregnancy outcome study: associations of GDM and obesity with pregnancy outcomes. Diabetes Care. 2012;35(4):780–6.
Fernandez-Twinn DS, Hjort L, Novakovic B, Ozanne SE, Saffery R. Intrauterine programming of obesity and type 2 diabetes. Diabetologia. 2019;62(10):1789–801.
Kawasaki M, Arata N, Miyazaki C, Mori R, Kikuchi T, Ogawa Y, et al. Obesity and abnormal glucose tolerance in offspring of diabetic mothers: a systematic review and meta-analysis. PLoS One. 2018;13(1):e0190676.
Lowe WL Jr, Scholtens DM, Kuang A, Linder B, Lawrence JM, Lebenthal Y, et al. Hyperglycemia and Adverse Pregnancy Outcome Follow-up Study (HAPO FUS): maternal gestational diabetes mellitus and childhood glucose metabolism. Diabetes Care. 2019;42(3):372–80.
Riskin-Mashiah S, Younes G, Damti A, Auslender R. First-trimester fasting hyperglycemia and adverse pregnancy outcomes. Diabetes Care. 2009;32(9):1639–43.
Yu Y, Arah OA, Liew Z, Cnattingius S, Olsen J, Sorensen HT, et al. Maternal diabetes during pregnancy and early onset of cardiovascular disease in offspring: population based cohort study with 40 years of follow-up. Bmj. 2019;367:l6398.
Geurtsen ML, van Soest EEL, Voerman E, Steegers EAP, Jaddoe VWV, Gaillard R. High maternal early-pregnancy blood glucose levels are associated with altered fetal growth and increased risk of adverse birth outcomes. Diabetologia. 2019;62(10):1880–90.
Hapo Study Cooperative Research Group, Metzger BE, Lowe LP, Dyer AR, Trimble ER, Chaovarindr U, et al. Hyperglycemia and adverse pregnancy outcomes. N Engl J Med. 2008;358(19):1991–2002.
Jaddoe VW, de Jonge LL, Hofman A, Franco OH, Steegers EA, Gaillard R. First trimester fetal growth restriction and cardiovascular risk factors in school age children: population based cohort study. Bmj. 2014;348:g14.
Daraki V, Georgiou V, Papavasiliou S, Chalkiadaki G, Karahaliou M, Koinaki S, et al. Metabolic profile in early pregnancy is associated with offspring adiposity at 4 years of age: the Rhea pregnancy cohort Crete, Greece. PLoS One. 2015;10(5):e0126327.
Wahab RJ, Voerman E, Jansen PW, Oei EHG, Steegers EAP, Jaddoe VWV, et al. Maternal glucose concentrations in early pregnancy and cardiometabolic risk factors in childhood. Obesity (Silver Spring). 2020;28(5):985–93.
El Hajj N, Pliushch G, Schneider E, Dittrich M, Muller T, Korenkov M, et al. Metabolic programming of MEST DNA methylation by intrauterine exposure to gestational diabetes mellitus. Diabetes. 2013;62(4):1320–8.
Ruchat SM, Houde AA, Voisin G, St-Pierre J, Perron P, Baillargeon JP, et al. Gestational diabetes mellitus epigenetically affects genes predominantly involved in metabolic diseases. Epigenetics. 2013;8(9):935–43.
Gagne-Ouellet V, Houde AA, Guay SP, Perron P, Gaudet D, Guerin R, et al. Placental lipoprotein lipase DNA methylation alterations are associated with gestational diabetes and body composition at 5 years of age. Epigenetics. 2017;12(8):616–25.
Mansell T, Ponsonby AL, Collier F, Burgner D, Vuillermin P, Lange K, et al. Genetic variation, intrauterine growth, and adverse pregnancy conditions predict leptin gene DNA methylation in blood at birth and 12 months of age. Int J Obes. 2019.
Mansell T, Ponsonby AL, Januar V, Novakovic B, Collier F, Burgner D, et al. Early-life determinants of hypoxia-inducible factor 3A gene (HIF3A) methylation: a birth cohort study. Clin Epigenetics. 2019;11(1):96.
Cardenas A, Gagne-Ouellet V, Allard C, Brisson D, Perron P, Bouchard L, et al. Placental DNA methylation adaptation to maternal glycemic response in pregnancy. Diabetes. 2018;67(8):1673–83.
Elliott HR, Sharp GC, Relton CL, Lawlor DA. Epigenetics and gestational diabetes: a review of epigenetic epidemiology studies and their use to explore epigenetic mediation and improve prediction. Diabetologia. 2019.
Finer S, Mathews C, Lowe R, Smart M, Hillman S, Foo L, et al. Maternal gestational diabetes is associated with genome-wide DNA methylation variation in placenta and cord blood of exposed offspring. Hum Mol Genet. 2015;24(11):3021–9.
Haertle L, El Hajj N, Dittrich M, Muller T, Nanda I, Lehnen H, et al. Epigenetic signatures of gestational diabetes mellitus on cord blood methylation. Clin Epigenetics. 2017;9:28.
Quilter CR, Cooper WN, Cliffe KM, Skinner BM, Prentice PM, Nelson L, et al. Impact on offspring methylation patterns of maternal gestational diabetes mellitus and intrauterine growth restraint suggest common genes and pathways linked to subsequent type 2 diabetes risk. FASEB J. 2014;28(11):4868–79.
Reichetzeder C, Dwi Putra SE, Pfab T, Slowinski T, Neuber C, Kleuser B, et al. Increased global placental DNA methylation levels are associated with gestational diabetes. Clin Epigenetics. 2016;8:82.
Yang IV, Zhang W, Davidson EJ, Fingerlin TE, Kechris K, Dabelea D. Epigenetic marks of in utero exposure to gestational diabetes and childhood adiposity outcomes: the EPOCH study. Diabet Med. 2018;35(5):612–20.
Suderman M, Staley JR, French R, Arathimos R, Simpkin A, Tilling K. dmrff: identifying differentially methylated regions efficiently with power and control. bioRxiv. 2018; https://www.biorxiv.org/content/biorxiv/early/2018/12/31/508556.full.pdf:508556.
Kulkarni H, Kos MZ, Neary J, Dyer TD, Kent JW Jr, Göring HHH, et al. Novel epigenetic determinants of type 2 diabetes in Mexican-American families. Hum Mol Genet. 2015;24(18):5330–44.
Al Muftah WA, Al-Shafai M, Zaghlool SB, Visconti A, Tsai P-C, Kumar P, et al. Epigenetic associations of type 2 diabetes and BMI in an Arab population. Clin Epigenetics. 2016;8:13.
Arpón A, Santos JL, Milagro FI, Cataldo LR, Bravo C, Riezu-Boj J-I, et al. Insulin sensitivity is associated with lipoprotein lipase (LPL) and catenin delta 2 (CTNND2) DNA methylation in peripheral white blood cells in non-diabetic young women. Int J Mol Sci. 2019;20(12):2928.
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 Diabet Endocrinol. 2015;3(7):526–34.
Chen P, Piaggi P, Traurig M, Bogardus C, Knowler WC, Baier LJ, et al. Differential methylation of genes in individuals exposed to maternal diabetes in utero. Diabetologia. 2017;60(4):645–55.
Florath I, Butterbach K, Heiss J, Bewerunge-Hudler M, Zhang Y, Schöttker B, et al. Type 2 diabetes and leucocyte DNA methylation: an epigenome-wide association study in over 1,500 older adults. Diabetologia. 2016;59(1):130–8.
Hidalgo B, Irvin MR, Sha J, Zhi D, Aslibekyan S, Absher D, et al. Epigenome-wide association study of fasting measures of glucose, insulin, and HOMA-IR in the Genetics of Lipid Lowering Drugs and Diet Network study. Diabetes. 2014;63(2):801–7.
Hjort L, Novakovic B, Grunnet LG, Maple-Brown L, Damm P, Desoye G, et al. Diabetes in pregnancy and epigenetic mechanisms-how the first 9 months from conception might affect the child's epigenome and later risk of disease. Lancet Diabetes Endocrinol. 2019;7(10):796–806.
Kriebel J, Herder C, Rathmann W, Wahl S, Kunze S, Molnos S, et al. Association between DNA methylation in whole blood and measures of glucose metabolism: KORA F4 study. PLoS One. 2016;11(3):e0152314-e.
Lin X, Lim IY, Wu Y, Teh AL, Chen L, Aris IM, et al. Developmental pathways to adiposity begin before birth and are influenced by genotype, prenatal environment and epigenome. BMC Med. 2017;15(1):50.
Liu J, Carnero-Montoro E, van Dongen J, Lent S, Nedeljkovic I, Ligthart S, et al. An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis. Nat Commun. 2019;10(1):2581.
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;48(1):58–70.
Soriano-Tárraga C, Jiménez-Conde J, Giralt-Steinhauer E, Mola-Caminal M, Vivanco-Hidalgo RM, Ois A, et al. Epigenome-wide association study identifies TXNIP gene associated with type 2 diabetes mellitus and sustained hyperglycemia. Hum Mol Genet. 2016;25(3):609–19.
Toperoff G, Aran D, Kark JD, Rosenberg M, Dubnikov T, Nissan B, et al. Genome-wide survey reveals predisposing diabetes type 2-related DNA methylation variations in human peripheral blood. Hum Mol Genet. 2012;21(2):371–83.
Walaszczyk E, Luijten M, Spijkerman AMW, Bonder MJ, Lutgers HL, Snieder H, et al. DNA methylation markers associated with type 2 diabetes, fasting glucose and HbA(1c) levels: a systematic review and replication in a case-control sample of the Lifelines study. Diabetologia. 2018;61(2):354–68.
Yuan W, Xia Y, Bell CG, Yet I, Ferreira T, Ward KJ, et al. An integrated epigenomic analysis for type 2 diabetes susceptibility loci in monozygotic twins. Nat Commun. 2014;5:5719.
Sletner L, Jenum AK, Yajnik CS, Morkrid K, Nakstad B, Rognerud-Jensen OH, et al. Fetal growth trajectories in pregnancies of European and South Asian mothers with and without gestational diabetes, a population-based cohort study. PLoS One. 2017;12(3):e0172946.
Sovio U, Murphy HR, Smith GC. Accelerated fetal growth prior to diagnosis of gestational diabetes mellitus: a prospective cohort study of nulliparous women. Diabetes Care. 2016;39(6):982–7.
Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, Schadt EE, et al. Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet. 2009;41(1):56–65.
Mahajan A, Taliun D, Thurner M, Robertson NR, Torres JM, Rayner NW, et al. Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat Genet. 2018;50(11):1505–13.
Wain LV, Vaez A, Jansen R, Joehanes R, van der Most PJ, Erzurumluoglu AM, et al. Novel blood pressure locus and gene discovery using genome-wide association study and expression data sets from blood and the kidney. Hypertension. 2017; HYPERTENSIONAHA.117.09438.
Wuttke M, Li Y, Li M, Sieber KB, Feitosa MF, Gorski M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet. 2019;51(6):957–72.
Matsuzaki G, Umemura M. Interleukin-17 family cytokines in protective immunity against infections: role of hematopoietic cell-derived and non-hematopoietic cell-derived interleukin-17 s. Microbiol Immunol. 2018;62(1):1–13.
Giwa AM, Ahmed R, Omidian Z, Majety N, Karakus KE, Omer SM, et al. Current understandings of the pathogenesis of type 1 diabetes: genetics to environment. World J Diabetes. 2020;11(1):13–25.
Verweij N, Mateo Leach I, van den Boogaard M, van Veldhuisen DJ, Christoffels VM, LifeLines Cohort S, et al. Genetic determinants of P wave duration and PR segment. Circ Cardiovasc Genet. 2014;7(4):475–81.
van Leeuwen M, Prins SM, de Valk HW, Evers IM, Visser GH, Mol BW. [Gestational diabetes mellitus: treatment reduces the risk of complications] Diabetes gravidarum. Behandeling vermindert kans op complicaties. Ned Tijdschr Geneeskd. 2011;155:A2291.
Clausen T, Burski TK, Oyen N, Godang K, Bollerslev J, Henriksen T. Maternal anthropometric and metabolic factors in the first half of pregnancy and risk of neonatal macrosomia in term pregnancies. A prospective study. Eur J Endocrinol. 2005;153(6):887–94.
Kooijman MN, Kruithof CJ, van Duijn CM, Duijts L, Franco OH, Van IMH, et al. The Generation R Study: design and cohort update 2017. Eur J Epidemiol. 2016;31(12):1243–64.
Kruithof CJ, Kooijman MN, van Duijn CM, Franco OH, de Jongste JC, Klaver CC, et al. The Generation R Study: Biobank update 2015. Eur J Epidemiol. 2014;29(12):911–27.
Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan ST, et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015;16:37.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.
Tukey JW. Exploratory data analysis. Reading: Addison-Wesley; 1977.
Triche T. FDb.InfiniumMethylation.hg19: annotation package for Illumina Infinium DNA methylation probes. R package version 220. 2014.
Infinium HumanMethylation450K v1.2 Product Files. 2014. http://support.illumina.com/downloads/infinium_humanmethylation450_product_files.html.
Jaddoe VW, van Duijn CM, Franco OH, van der Heijden AJ, van Iizendoorn MH, de Jongste JC, et al. The Generation R Study: design and cohort update 2012. Eur J Epidemiol. 2012;27(9):739–56.
Gervin K, Salas LA, Bakulski KM, van Zelm MC, Koestler DC, Wiencke JK, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenetics. 2019;11(1):125.
Hunt KJ, Schuller KL. The increasing prevalence of diabetes in pregnancy. Obstet Gynecol Clin N Am. 2007;34(2):173–99 vii.
Joubert BR, Felix JF, Yousefi P, Bakulski KM, Just AC, Breton C, et al. DNA Methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96.
Sharp GC, Salas LA, Monnereau C, Allard C, Yousefi P, Everson TM, et al. Maternal BMI at the start of pregnancy and offspring epigenome-wide DNA methylation: findings from the pregnancy and childhood epigenetics (PACE) consortium. Hum Mol Genet. 2017;26(20):4067–85.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Niklasson A, Ericson A, Fryer JG, Karlberg J, Lawrence C, Karlberg P. An update of the Swedish reference standards for weight, length and head circumference at birth for given gestational age (1977-1981). Acta Paediatr Scand. 1991;80(8-9):756–62.
The Generation R Study is conducted by the Erasmus Medical Center in close collaboration with the School of Law and Faculty of Social Sciences of the Erasmus University Rotterdam; the Municipal Health Service Rotterdam area, Rotterdam; the Rotterdam Homecare Foundation, Rotterdam; and the Stichting Trombosedienst & Artsenlaboratorium Rijnmond (STAR-MDC), Rotterdam. We gratefully acknowledge the contribution of children and parents, general practitioners, hospitals, midwives, and pharmacies in Rotterdam. The study protocol was approved by the Medical Ethical Committee of the Erasmus Medical Centre, Rotterdam. Written informed consent was obtained for all participants. The generation and management of the Illumina 450K methylation array data (EWAS data) for the Generation R Study was executed by the Human Genotyping Facility of the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, the Netherlands. We thank Mr. Michael Verbiest, Ms. Mila Jhamai, Ms. Sarah Higgins, Mr. Marijn Verkerk, and Dr. Lisette Stolk for their help in creating the EWAS database. We thank Dr. A. Teumer for his work on the quality control and normalization scripts.
The general design of the Generation R Study is made possible by financial support from the Erasmus Medical Center, Rotterdam; the Erasmus University Rotterdam; the Netherlands Organization for Health Research and Development; and the Ministry of Health, Welfare, and Sport. The EWAS data was funded by a grant to VWVJ from the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO), Netherlands Consortium for Healthy Aging (NCHA; project nr. 050-060-810), by funds from the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, and by a grant from the National Institute of Child and Human Development (R01HD068437). VWVJ received a Consolidator Grant from the European Research Council (ERC-2014-CoG-648916). This project received funding from the European Union’s Horizon 2020 research and innovation programme LifeCycle (733206) and from the European Joint Programming Initiative “A Healthy Diet for a Healthy Life” (JPI HDHL, NutriPROGRAM project, ZonMw the Netherlands no.529051022 and Precise project no. P75416). RG received funding from the Dutch Heart Foundation (grant number 2017 T013), the Dutch Diabetes Foundation (grant number 2017.81.002), and the Netherlands Organization for Health Research and Development (NWO, ZonMW, grant number 543003109).
Ethics approval and consent to participate
The study has been approved by the Medical Ethical Committee of the Erasmus MC, University Medical Center Rotterdam (MEC 198.782/2001/31). Written informed consent was obtained for all participants (49).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of mother and child characteristics between participants included and non-participants.
Epigenome-wide association study results of maternal early-pregnancy glucose concentrations and DNA methylation in cord blood. Figure S1b. Epigenome-wide association study results of maternal early-pregnancy insulin concentrations and DNA methylation in cord blood.
CpGs with p-values <1.0 x 10-4 from epigenome-wide association study of maternal early-pregnancy glucose concentrations and DNA methylation. Table S3. CpGs with p-values <1.0 x 10-4 from epigenome-wide association study of maternal early-pregnancy insulin concentrations and DNA methylation.
CpGs with p-values <1.0 x 10-4 from epigenome-wide association study of maternal early-pregnancy glucose concentrations and DNA methylation - model without pre-pregnancy BMI adjustment. Table S5. CpGs with p-values <1.0 x 10-4 from epigenome-wide association study of maternal early-pregnancy insulin concentrations and DNA methylation - model without pre-pregnancy BMI adjustment.
Epigenome-wide association study results of maternal early-pregnancy insulin concentrations and DNA methylation in cord blood in normal weight women. Figure S2b. Epigenome-wide association study results of maternal early-pregnancy insulin concentrations and DNA methylation in cord blood in overweight or obese women.
Differentially Methylated Regions with p-values <1.0 x 10-4 associated with maternal early-pregnancy glucose concentrations. Table S7. Differentially Methylated Regions with p-values <1.0 x 10-4 associated maternal early-pregnancy insulin concentrations.
Look-up of two CpGs identified in cord blood in corresponding maternal BMI strata in epigenome-wide association study results of maternal early-pregnancy glucose concentrations and DNA methylation in 10-year-old children.
Look-up results of maternal glucose metabolism associated CpGs in maternal early-pregnancy glucose EWAS results.
Look-up results of maternal glucose metabolism associated CpGs in maternal early-pregnancy insulin EWAS results.
Look-up results of adult type 2 diabetes associated CpGs in maternal early-pregnancy glucose EWAS results.
Look-up results of adult type 2 diabetes associated CpGs in maternal early-pregnancy insulin EWAS results.
About this article
Cite this article
Geurtsen, M.L., Jaddoe, V.W.V., Gaillard, R. et al. Associations of maternal early-pregnancy blood glucose and insulin concentrations with DNA methylation in newborns. Clin Epigenet 12, 134 (2020). https://doi.org/10.1186/s13148-020-00924-3
- Maternal glucose
- Maternal hyperglycemia
- Gestational diabetes
- Maternal insulin
- Diabetes mellitus
- DNA methylation
- Differentially methylated regions