==== Front
Clin Epigenetics
Clin Epigenetics
Clinical Epigenetics
1868-7075
1868-7083
BioMed Central London

1524
10.1186/s13148-023-01524-7
Research
Maternal glycemia in pregnancy is longitudinally associated with blood DNAm variation at the FSD1L gene from birth to 5 years of age
Taschereau Amélie 1
Thibeault Kathrine 1
Allard Catherine 2
Juvinao-Quintero Diana 8
Perron Patrice 23
Lutz Sharon M. 45
Bouchard Luigi Luigi.Bouchard@USherbrooke.ca

136
Hivert Marie-France 47
1 grid.86715.3d 0000 0000 9064 6198 Department of Biochemistry and Functional Genomics, Faculty of Medicine and Health Sciences (FMHS), Université de Sherbrooke, Sherbrooke, QC Canada
2 grid.411172.0 0000 0001 0081 2808 Centre de Recherche du Centre hospitalier universitaire de Sherbrooke (CR-CHUS), Sherbrooke, QC Canada
3 grid.86715.3d 0000 0000 9064 6198 Department of Medicine, FMHS, Université de Sherbrooke, Sherbrooke, QC Canada
4 grid.38142.3c 000000041936754X Department of Population Medicine, Harvard Pilgrim Health Care Institute, Harvard Medical School, Boston, MA USA
5 grid.38142.3c 000000041936754X Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA USA
6 grid.420762.5 0000 0000 8794 2033 Clinical Department of Laboratory Medicine, Pavillon des Augustines, Centre intégré universitaire de santé et de services sociaux (CIUSSS) du Saguenay-Lac-Saint-Jean – Hôpital de Chicoutimi, 305 rue St-Vallier, Saguenay, QC G7H 5H6 Canada
7 grid.32224.35 0000 0004 0386 9924 Diabetes Unit, Massachusetts General Hospital, Boston, MA USA
8 grid.38142.3c 000000041936754X Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA USA
29 6 2023
29 6 2023
2023
15 10715 3 2023
23 6 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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.
Background

In utero exposure to maternal hyperglycemia has been associated with an increased risk for the development of chronic diseases in later life. These predispositions may be programmed by fetal DNA methylation (DNAm) changes that persist postnatally. However, although some studies have associated fetal exposure to gestational hyperglycemia with DNAm variations at birth, and metabolic phenotypes in childhood, no study has yet examined how maternal hyperglycemia during pregnancy may be associated with offspring DNAm from birth to five years of age.

Hypothesis

Maternal hyperglycemia is associated with variation in offspring DNAm from birth to 5 years of age.

Methods

We estimated maternal hyperglycemia using the area under the curve for glucose (AUCglu) following an oral glucose tolerance test conducted at 24–30 weeks of pregnancy. We quantified DNAm levels in cord blood (n = 440) and peripheral blood at five years of age (n = 293) using the Infinium MethylationEPIC BeadChip (Illumina). Our total sample included 539 unique dyads (mother–child) with 194 dyads having DNAm at both time-points. We first regressed DNAm M-values against the cell types and child age for each time-point separately to account for the difference by time of measurement for these variables. We then used a random intercept model from the linear mixed model (LMM) framework to assess the longitudinal association between maternal AUCglu and the repeated measures of residuals of DNAm. We adjusted for the following covariates as fixed effects in the random intercept model: maternal age, gravidity, smoking status, child sex, maternal body mass index (BMI) (measured at first trimester of pregnancy), and a binary variable for time-point.

Results

In utero exposure to higher maternal AUCglu was associated with lower offspring blood DNAm levels at cg00967989 located in FSD1L gene (β = − 0.0267, P = 2.13 × 10–8) in adjusted linear regression mixed models. Our study also reports other CpG sites for which DNAm levels were suggestively associated (P < 1.0 × 10–5) with in utero exposure to gestational hyperglycemia. Two of these (cg12140144 and cg07946633) were found in the promotor region of PRDM16 gene (β: − 0.0251, P = 4.37 × 10–07 and β: − 0.0206, P = 2.24 × 10–06, respectively).

Conclusion

Maternal hyperglycemia is associated with offspring DNAm longitudinally assessed from birth to 5 years of age.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-023-01524-7.

Keywords

Epigenetics
DNAm
Fetal programming
Gestational hyperglycemia
Pregnancy
Offspring
http://dx.doi.org/10.13039/100008985 Diabète Québec Fonds de la recherche du Québec en santé (FRQS)Canadian Institute of Health Research (CIHR)PJT-152989 MOP 115071 Bouchard Luigi Hivert Marie-France Fonds de recherche du Québec - Santé (FRQS)20697 Hivert Marie-France issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2023
==== Body
pmcBackground

Gestational diabetes mellitus (GDM) is defined by hyperglycemia first occurring during pregnancy in an individual who did not previously have diabetes [1]. GDM affects approximately 14% of pregnancies worldwide [2] with variation in prevalence depending on ethnicity [3] and the diagnostic criteria applied [4]. Regardless of GDM diagnosis, maternal hyperglycemia along the whole spectrum has been linearly associated with adverse perinatal outcomes in both mother and offspring [5]. Exposure to maternal hyperglycemia in utero is also associated with an increased risk for the development of chronic diseases later in childhood, including obesity [6], metabolic [7] and cardiovascular diseases [8], asthma [9], autism [10], as well as other adverse neurodevelopmental outcomes [11]. The Developmental Origin of Health and Disease (DOHaD) hypothesis is central to the concept of fetal programming and based on lasting fetal metabolic adaptions to environmental/maternal stressors during the critical developmental period [12]. Epigenetic marks are highly malleable during in utero development and in early life but can be stable over time [13]. Epigenetics is one of the most likely molecular mechanisms as part of the DOHAD theory and fetal metabolic programming.

DNA methylation (DNAm) is the most studied epigenetic modification. It involves the addition of a methyl group to the fifth carbon of a cytosine when located just upstream of a guanine [14]. Growing evidence supports that DNAm may be the bridge between prenatal exposure to maternal hyperglycemia and the predisposition to chronic disease in later life. The most comprehensive study consists of a meta-analysis of an epigenome-wide association study (EWAS) investigating associations between maternal hyperglycemia exposure and cord blood DNAm variations. This study reported significant associations (FDR < 0.05) at two loci between cord blood DNAm and the area under the curve of glucose (AUCglu) following an oral glucose tolerance test (OGTT) [15]. Candidate gene studies also reported significant correlations between maternal hyperglycemia exposure during pregnancy and offspring DNAm in genes implicated in chronic disease pathways [16, 17].

However, although the association between maternal hyperglycemia exposure and DNAm measured at birth has been explored in a few studies, no study has yet examined how maternal hyperglycemia may be associated with offspring DNAm from birth to five years of age. To address this gap, we conducted an EWAS examining the association between maternal glycemic response to an OGTT performed during pregnancy, and offspring DNAm assessed longitudinally at birth (cord blood) and at age 5 years (peripheral blood), using modeling approaches that integrate DNAm at both time-points.

Methods

Study population

Our study sample consists of 539 mother–child dyads from the Genetics of Glucose Regulation in Gestation and Growth (Gen3G) prospective cohort. We recruited Gen3G participants in the Sherbrooke City area, Canada, to investigate environmental and genetic determinants of glucose regulation in pregnancy, and their impacts on fetal development and offspring health [18]. Between 2010 and 2013, we invited pregnant women to participate in the study when they attended a routine prenatal blood sampling visit during their first trimester of pregnancy at the Centre Hospitalier Universitaire de Sherbrooke (CHUS), if they planned to deliver at CHUS [18]. Eligible women were ≥ 18 years old, with a singleton pregnancy, and not reporting using medication that alters glucose. We excluded women with pre-existing diabetes or diabetes diagnosed in the first trimester based on laboratory screening (A1c > 6.5% or 1 h-glucose ≥ 10.3 mmol/L post-50 g glucose challenge test [GCT]). Details on Gen3G participants have been published previously [18]. All women provided signed consent, and the CHUS IRB approved all Gen3G protocols.

At enrollment, women provided information on age, ethnicity, smoking habits during pregnancy (never/former smoker vs current smoker) and gravidity (primigravid vs multiparous). During the first (6–14 weeks) and second (24–30 weeks) pregnancy visits, research staff measured weight (in kg) and height (in m) using standardized procedures, and we calculated body mass index (BMI) using the standard formula (kg/m2) [18]. At the first visit, most women completed a non-fasting 50 g-GCT. At the second visit, all women conducted a 75 g-OGTT in order to measure glucose tolerance. At delivery, we abstracted from the medical records information about gestational age at birth (in weeks), child sex, and birth weight (in g). Research staff collected umbilical cord blood samples rapidly after delivery (< 30 min).

We conducted in-person follow-up visits with mother–child pairs at the CHUS Research Center approximately 5 years after birth. During this visit, we collected updated medical history and anthropometric measures, in addition to child blood samples.

Maternal glucose measured during a 75 g-OGTT in the second trimester, and estimation of the area under the curve for glucose

We measured glucose levels (mmol/L) using the glucose hexokinase method (Roche Diagnostics, Indianapolis, US) at the CHUS biochemistry clinical laboratory rapidly after blood collection [18, 19]. We calculated the AUCglu based on the trapezoidal formula [20], to capture a global measure of maternal hyperglycemia in pregnancy that encompasses all three time-points of the OGTT since all of them have been associated with adverse pregnancy outcomes [5], and the AUCglu was the most informative phenotype in prior EWAS investigating dysregulation in child cord blood DNAm [15].

Genome-wide DNA methylation level assessment

Genomic DNA was purified from umbilical cord blood cells and from blood cells, respectively, collected at delivery and 5 years of age using the AllPrep DNA/RNA/Protein Mini kit (QIAGEN, Hilden, Germany). DNA was then bisulfite converted with the EZ-Methylation kit (Zymo Research, CA, US), and subsequently analyzed for DNAm at > 850,000 CpG sites at single-nucleotide resolution using the Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA). We first performed a manual inspection of the quality of the methylation assay looking at probes call rate (> 99%), color balance, staining, extension, hybridization, specificity and bisulfite conversion (probe types I and II) using the BeadArray Controls Reporter software. We transferred raw methylation files into R (version 4.1.1) for further preprocessing using the minfi R package [21]. Our final dataset included 440 cord blood samples and 293 peripheral blood samples collected at the ~ 5-year visit after quality controls. In terms of participants, this represents 539 unique dyads (mother–child pairs), from which 194 dyads had both time-points (cord blood and peripheral blood) included in the analyses (see support Fig. 1).Fig. 1 Flowchart illustrating the participant selection from the Gen3G prospective cohort

We implemented functional normalization for background and dye bias correction [22] of the methylation data, and the regression on correlated probes (RCP) method [23] to adjust for probe-type bias, as previously described [24]. We filtered out probes with null methylation variance and detection P-value > 0.05 in at least 5% of the samples. We annotated CpG sites in the array using the R package IlluminaHumanMethylationEPICmanifest (version 1.0 B4), and we excluded probes located on sex chromosomes, non-CpG probes (“rs” or “ch”), single-nucleotide polymorphism (SNP)-associated probes at the CpG site, the single-base extension, or anywhere within the probe, if the SNP had a minor allele frequency (MAF) ≥ 0.05. In addition, we removed cross-reactive probes as previously described [25]. We used ComBat [26] to correct for technical variability and batch effects, while protecting our variables of interest. The model thus included the three glucose values (fasting glucose, 1 h-glucose, and 2 h-glucose post-OGTT), maternal age, BMI in the first trimester of pregnancy, gravidity, smoking status, and child sex. After quality control, we included up to 719,360 high-quality CpG sites in the analysis using DNAm from combined cord blood and childhood peripheral blood. We calculated the predicted cell counts in cord blood based on the method described by Gervin et al.[27], and the predicted cell counts in child blood based on the method from Salas et al.[28], implemented in the FlowSorted.Blood.EPIC R package.

Statistical analyses

We described baseline characteristics of the overall sample, and of samples at each time-point, using standard descriptive statistics. For each CpG site, we used methylation values measured in the β-value scale, as a range estimating between 0% (unmethylated) and 100% (completely methylated), to transform into M-values prior to analyses based on the statistical validity of this method to identify differential DNAm [29]. Due to differences in the cell types predicted in each time-point (i.e., cord blood = CD8 + T cells, CD4 + T cells, monocytes, natural killer, B cells, granulocytes and nucleated red blood cells (nRBC); peripheral blood at ~ 5 years = CD8 + T cells, CD4 + T cells, monocytes, natural killer, B cells and neutrophils) and different age measures (cord blood: gestational age in weeks; peripheral blood: child age in years), we regressed DNAm M-values against the cell types and age for each time-point separately, and retrieved residuals from this association to perform the longitudinal EWAS. For the linear mixed model (LMM), we considered a random intercept model using the nlme R package [30] to assess the longitudinal association between maternal AUCglu (primary exposure) and repeated measures of residuals of DNAm (in M-values adjusted for cell types and age, respectively, in each DNAm time-point) at ~ 719,360 CpG sites as outcomes. We conducted similar analyses for individual glucose components of maternal AUCglu: fasting glucose, 1 h-glucose, and 2 h-glucose post-OGTT as secondary exposures. We included the following covariates as fixed effects: maternal age, gravidity (primigravid vs multigravida), smoking in the first trimester of pregnancy (smoker versus non-smoker), child sex, maternal BMI in the first trimester of pregnancy, and a binary variable for time-point (i.e., birth = 0, early childhood = 1). We used a random intercept for the subject to account for the correlation of multiple measures on the same individual. We also conducted separate linear models for individual CpG at each of the two time-points (cord blood and 5y peripheral blood) using blood DNAm transformed in the same way (residuals on cell counts and age), for all four maternal glycemic exposures (maternal AUCglu, fasting glucose, 1 h post-OGTT and 2 h post-OGTT) using the same covariates (without the time-point variable). We also explored if we could find differentially methylated regions (DMRs) using the DMRff method [31] using the residuals that were used for the EWAS as the individual-level data. We choose this method as it provides decent power for detecting DMRs and without an inflated Type I error rate [32].

We reported coefficient estimates at the individual CpG site as the mean difference in DNAm residuals on cells and age, per unit increase in the AUCglu. We calculated the genomic inflation factor or lambda (λ) for each EWAS (Additional file 1). We adjusted associations for multiple testing using Bonferroni, deeming epigenome-wide significant associations as P < 6.9 × 10–8, while we considered “suggestive” associations identified with an arbitrary P ≤ 1.0 × 10–5. For each top CpG site (P ≤ 1.0 × 10–5 in at least one model (mixed model, or linear models in cord blood and 5y peripheral blood)), we visually inspected outliers and distribution by plotting the residuals (standardized or not) versus the fitted values, and by using quantile plots. Top CpG sites identified in the EWAS were annotated to their gene using the Illumina annotation (IlluminaHumanMethylationEPICanno.ilm10b2.hg19 R package [33]) and if not annotated, we completed using the closest gene using the UCSC Gene Browser. We visualized the results using Manhattan and quantile plots. All analyses were performed in R (version 4.1.1).

Results

Participants’ characteristics

The characteristics of the mothers and children are shown in Table 1. Briefly, mothers were 28.2 ± 4.2 years old (mean ± SD) and had a median [IQR] BMI of 23.9 [21.6; 27.9] kg/m2 at the first trimester of pregnancy (median 9.4 [8.1; 11.7] weeks). Forty-nine (9.1%) women reported smoking and 180 (33.4%) were primigravid. Median and IQR values of maternal glycemia during the 75 g OGTT (median gestational age; 26.3 [25.9; 27.1] weeks) are indicated in Table 1. The median AUCglu was 12.0 [10.5; 13.4] mmol/L*h and fifty (9.3%) women developed GDM. At delivery, the median gestational age was 39.7 [38.9; 40.4] weeks and 252 of these children were female (46.8%). At the 5-year visit, the children had a median age of 5.2 [5.1; 5.4] years old and 130 (44.4%) were female.Table 1 Characteristics of Gen3G mother–child pairs during pregnancy, delivery, and child at ~ 5 years

	N	Mean ± SD or median [IQR] or N (%)	
Maternal measures – 1st trimester	
Age, years	539	28.2 ± 4.2	
Gestational age, weeks	539	9.4 [8.1; 11.7]	
Primigravid, yes	539	180 (33.4%)	
Smoking status, yes	539	49 (9.1%)	
BMI, kg/m2	539	23.9 [21.6; 27.9]	
Maternal measures – 2nd trimester	
Gestational age, weeks	539	26.3 [25.9; 27.1]	
GDM, yes	539	50 (9.3%)	
Glucose fasting, mmol/L	539	4.2 [3.9; 4.4]	
Glucose 1 h post-OGTT, mmol/L	539	7.1 [6.0; 8.2]	
Glucose 2 h post-OGTT, mmol/L	539	5.7 [4.8; 6.6]	
AUCglu, mmol/L * h	539	12.0 [10.5; 13.4]	
Measures at delivery	
Gestational age, weeks	539	39.7 [38.9; 40.4]	
Sex, girls	539	252 (46.8%)	
Child measures – 5 years	
Age, years	293	5.2 [5.1; 5.4]	
Sex, girls	293	130 (44.4%)	
AUCglu, Area Under the Curve of glucose; BMI, Body Mass Index; GDM, Gestational Diabetes Mellitus; SD, Standard Deviation

Maternal AUCglu is associated with DNAm variations in cord blood and in blood at 5 years of age

We first identified a total of 10 CpG sites (9 loci) at which offspring DNAm levels were associated with exposure to maternal AUCglu at a suggestive threshold of P-value ≤ 10–5 (Fig. 2). The strongest association reached the genome-wide significant threshold and was observed at the cg00967989 located within the proximal promotor of the FSD1L gene (β = − 0.0267, P = 2.13 × 10–8) (Fig. 3a). Among the other loci with suggestive associations, cg07946633 (β = − 0.0206, P = 2.24 × 10–6) and cg12140144 (β = − 0.0251, P = 4.37 × 10–7) were both located within the promotor region of the PRDM16 gene (Fig. 3b). The seven other CpG sites were located at the HES1, CCDC28A, INPP5F, CKB, GP1BB/SEPT5, SAR1B, and RGS3 gene loci. We presented in Table 2 the complete list of CpGs at which offspring blood DNAm levels were associated with exposure to high maternal AUCglu at a suggestive threshold (P ≤ 1.0 × 10–5), their associated genes, genomic location, regression β values and P-values from the main longitudinal EWAS (based on LMM). Using DMRff, we identified two DMRs, one containing the two CpG sites (cg07946633 and cg12140144; adj P = 0.03) previously reported at the PRDM16 gene locus and the second containing the one CpG site (cg00967989; adj P = 0.02), previously reported at the FSD1L gene locus (Additional file 2).Fig. 2 Epigenome-wide association plot (− log10 P-values) for maternal AUCglu. Manhattan plot for the EWAS of maternal AUCglu following a 75 g-OGTT with DNAm measured longitudinally at birth and five years of age (red line: Bonferroni threshold = P-value < 6.9 × 10–8; blue line: Suggestive threshold = P < 1.0 × 10–5)

Fig. 3 FSD1L and PRDM16 gene regional plot association (− log10 P-values) for maternal AUCglu. Genomic regional plots for a FSD1L (on Chr 9) and b PRDM16 (on Chr 1) genes of maternal AUCglu following a 75 g-OGTT with DNAm measured longitudinally at birth and five years of age (Solid line: Bonferroni threshold = P-value < 6.9 × 10–8; dashed line: Suggestive threshold = P < 1.0 × 10–5)

Table 2 CpG sites identified (suggestive P < 10–5) in linear mixed models (LMM) testing associations between maternal AUCglu and DNAm measured in cord blood and blood at 5 years of age

CpGs	Chr	Position	Cord blood DNAm levels in β-value Mean ± SD N = 440	Five years blood DNAm levels in β-value Mean ± SD N = 293	Nearest gene	LMM regression estimates	
cg00967989	9	108,210,147	0.013 ± 0.002	0.015 ± 0.003	FSD1L	β: − 0.0267

P: 2.13 × 10–08

	
cg03703356	14	103,989,368	0.073 ± 0.012	0.132 ± 0.016	CKB	β: − 0.0217

P: 4.25 × 10–07

	
cg12140144	1	2,984,275	0.017 ± 0.004	0.045 ± 0.008	PRDM16	β: − 0.0251

P: 4.37 × 10–07

	
cg07946633	1	2,984,245	0.251 ± 0.034	0.339 ± 0.032	PRDM16	β: − 0.0206

P: 2.24 × 10–06

	
cg08694430	9	116,420,234	0.864 ± 0.032	0.838 ± 0.035	RGS3	β: 0.0324

P: 1.33 × 10–06

	
cg20153537	3	193,917,170	0.949 ± 0.017	0.944 ± 0.007	HES1	β: − 0.0223

P: 2.70 × 10–06

	
cg19978242	10	121,578,846	0.856 ± 0.030	0.886 ± 0.026	INPP5F	β: − 0.0322

P: 4.85 × 10–06

	
cg02357751	22	19,710,880	0.014 ± 0.003	0.030 ± 0.015	GP1BB; SEPT5	β: − 0.0368

P: 7.71 × 10–06

	
cg23517035	5	133,971,253	0.954 ± 0.010	0.939 ± 0.011	SAR1B	β: 0.0212

p: 8.68 × 10–06

	
cg17385945	6	139,095,355	0.010 ± 0.002	0.014 ± 0.003	CCDC28A	β: − 0.0257

p: 9.93 × 10–06

	
Model adjusted for maternal age, gravidity, smoking status, child sex, BMI at first trimester of pregnancy and the binary variable for time-point. Significant result after Bonferroni correction with P-value < 6.9 × 10–8 are in bold

AUCglu, Area Under the Curve of glucose; Chr, Chromosome; CpG, Cytosine–phosphate–Guanine; DNAm, DNA methylation; SD, Standard deviation

We also tested associations between maternal AUCglu and DNAm measured at birth (cord blood) and at 5 years (peripheral blood), separately (Additional file 3). When looking at individual time-points for FSD1L, the signal appears to be primarily driven by its association in childhood, with only nominal association observed at birth. A similar pattern applies to the CpGs localized at the GP1BB/SEPT5 gene locus. In contrast, the signals identified for PRDM16 (cg12140144), SAR1B, CCDC28A, and INPP5F appear to be primarily driven by their associations reported at birth. Interestingly, the signal observed for PRDM16 (cg07946633) appears to be equally driven by its associations at birth and at age 5 years with a similar pattern in HES1, CKB, and RGS3 genes.

Most of the CpG sites associated with maternal AUCglu were associated with 1 h and 2 h post-OGTT glucose levels

Using LMM, we tested the association between offspring DNAm levels (cord blood and 5-year blood), and maternal fasting glucose, 1 h, and 2 h post-OGTT glucose levels as separate exposures (Additional files 4, 5, and 6). From the 10 CpG sites initially associated with maternal AUCglu (P-value ≤ 10–5), DNAm levels at cg00967989, located in the FSD1L gene, were associated with maternal 1 h-glucose (β = − 0.0347, P = 6.68 × 10–8) and suggestively with the 2 h-glucose (β = − 0.0357, P = 2.80 × 10–6) post-OGTT. DNAm at CpG sites within the PRDM16 gene associated with maternal AUCglu were suggestively associated with maternal 1 h-glucose (only cg12140144) and 2 h-glucose (both cg12140144 and cg07946633) post-OGTT. In addition, DNAm levels at cg17385945, cg19978242, cg03703356, and cg08694430 located within the CCDC28A, INPP5F, CKB, and RGS3 genes, respectively, were suggestively associated with the 1 h post-OGTT glucose levels (same direction of effect to what we had found in the analyses with maternal AUCglu). None of the CpG sites initially associated with maternal AUCglu or with the 1 h- and 2 h-glucose was found associated with fasting glucose levels. These results, along with all CpG sites we found significantly associated with fasting glucose, 1 h and 2 h post-OGTT, are summarized in a Venn diagram (Fig. 4) and presented in Additional files 7, 8, and 9.Fig. 4 Venn diagram showing common CpGs between maternal AUCglu, fasting, 1 h and 2 h post-OGTT glucose levels. All CpG sites shown in the Venn diagram were at least suggestively associated (P < 1.0 × 10–5) with maternal glucose measurements. Cg00967989 is located at the FSD1L gene locus; cg12140144 and cg07946633 are located at the PRDM16 gene locus

Discussion

In this study, we conducted a longitudinal EWAS based on blood DNAm from birth and 5 years of age in childhood in relation to maternal glycemic response after a 75 g-OGTT conducted during the second trimester of pregnancy. We found CpG sites at which offspring blood DNAm was associated with maternal AUCglu as well as other glycemic indices measured during pregnancy, including 1 h and 2 h post-OGTT glucose levels, but not fasting glucose levels. The different results we obtained with fasting glucose levels could be explained by the different metabolic characteristics observed in the impaired glucose tolerance (IGT) and impaired glucose fasting (IGF) conditions. For example, the insulin resistance state affects different tissues in the two conditions (peripheral in IGT vs. hepatic in IGF) and only IGT subjects show a defect in late-phase insulin secretion in a non-pregnant population [34]. In addition, the two conditions affect different populations that only partially overlap [35], suggesting likely variations in the impact of these complications on potential epigenetic adaptations. Of note, in our previous EWAS of maternal glycemic traits (PACE cohorts meta-analysis), we did not find cord blood DNAm variations that were robustly associated with maternal fasting glucose, while we identified a few CpG sites at which cord blood DNAm was associated with AUCglu and other glycemic indices measured during pregnancy [15]. We identified one CpG site (cg00967989), located in the promotor region of FSD1L, at which higher maternal AUCglu was significantly associated with lower values of DNAm at birth and in childhood (after Bonferroni correction for multiple testing). This association was also confirmed with maternal 1 h-glucose and the 2 h-glucose post-OGTT (with a less stringent P-value threshold). Our study also reported other CpG sites for which offspring blood DNAm levels were suggestively associated (P-value < 10–5) with maternal AUCglu. Interestingly, two of these were found in the promotor region of PRDM16 with the same direction of association. We also observed nominal associations between maternal AUCglu and offspring DNAm at individual time-points (cord blood DNAm and blood DNAm at 5 years of age).

The FSD1L gene is located on chromosome 9 and encodes for type 2 cystatins. Among other functions, cystatins regulate the activity of endogenous cysteine proteinases which are involved in tumor cell invasion and metastasis [36] potentially via angiogenic properties of FSD1L [37]. Genome-wide association studies (GWAS) also reported significant associations between SNPs near its location and age at menarche in females [38], as well as early reproductive traits contributing to fertility in female cattle [39]. However, the exact role of FSD1L in humans has not been clearly established yet. Thus, although our results might support its implication in metabolic health programming in the context of prenatal exposure to gestational hyperglycemia, further studies will be needed to better understand the exact role of FSD1L in this context.

Our results also suggest lower blood DNAm levels in PRDM16 promotor at birth and 5 years in the offspring as a response to elevated gestational glycemia. This result at PRDM16 is consistent in the direction of association with findings by Côté, Gagné-Ouellet et al. who also reported inverse correlations between lower PRDM16 DNAm levels in the placenta of newborns and higher maternal glycemia [40]. The PRDM16 gene is located on chromosome 1 and encodes for one of the 17 PRDM protein family members. These proteins contain a PR domain involved in the regulation of transcriptional activity mediated by chromatin histone modifications [41]. PRDM16 is a transcription factor primarily involved in brown adipose tissue (BAT) differentiation and the transition from white adipose tissue (WAT) to beige adipose tissue (iBAT) [42]. The BAT is enriched in mitochondria whose inner membrane expresses the uncoupling protein 1 (UCP1). UCP1 allows protons to leak across the mitochondrial membrane, which leads to an increase in heat production at the expense of oxidative phosphorylation and adenosine triphosphate (ATP) production [43]. Consistently, BAT and iBAT have been positively associated with leanness and metabolic health in human adults [44]. More specifically, PRDM16 induces the expression of peroxisome proliferator-activated receptor-y (PPAR-y) and PPAR-y coactivator-1α (PGC-1α) with which it forms a complex that regulates the expression of genes associated with lipid catabolism and thermogenesis [45]. Dysregulation of DNAm and gene expression of PGC-1α was observed in the placenta of GDM-exposed participants [46]. Studies involving different polymorphisms of PRDM16 have also shown significant correlations between genetic variants of the gene and the development of obesity in humans [47, 48]. Moreover, Lui et al. [49] demonstrated lower DNAm levels in the promotor region of PRDM16, measured in omental adipose tissue (OAT) in individuals with obesity compared with normal-weight subjects, suggesting a role for epigenetic regulation of PRDM16 in fat mass redistribution of peripheral iBAT to OAT in obesity. However, our findings are based solely on blood DNAm, which may reflect variations in DNAm in other tissues exposed similarly during in utero environmental factors, but we could not biopsy various tissues in healthy newborns for ethical reasons.

Even though computationally intensive, the use of multilevel modeling to analyze trends in DNAm variation is a better strategy to understand the biological relevance of methylation changes over time, versus analyzing their variation at each time-point separately [50]. In our analysis, we identified one signal at FS1DL with epigenome-wide significance in the longitudinal analysis and in the analysis of blood DNAm at 5 years. With few exceptions, most of the top signals detected in the longitudinal analysis showed a higher level of statistical significance in this versus the time-point specific analysis, with consistent association estimates observed across analyses. These findings support the advantage of using multilevel modeling over the single time-point analysis to capture CpGs with a time-varying association with the exposure, in this case, maternal hyperglycemia in mid-pregnancy. To our knowledge, this is the first study that has used a longitudinal approach to investigate changes in offspring DNAm associated with maternal hyperglycemia. Previous studies have looked at associations between cord blood DNAm in relation to elevated maternal AUCglu or individual glucose traits measured after an OGTT in pregnancy and have followed up these findings by investigating the cross-sectional association of important baseline markers with cardiometabolic phenotypes in infants and adults [15]. Similarly, the metabolic risk of infants exposed to GDM has been assessed in early childhood (3–10 years of age), identifying epigenetic age acceleration in those born to moms with GDM, which was further associated with cardiometabolic risk factors [51]. Future prospective studies including larger sample sizes and multiple time-points over the life course are warranted to investigate variation in offspring DNAm in response to maternal hyperglycemia in pregnancy and to evaluate the implication of these markers in the offspring’s future health.

Strengths and limitations

A strength of our study is the prospective cohort design, which allowed collection of samples to evaluate DNAm levels longitudinally from birth to mid-childhood, and their association with maternal glycemia objectively measured during a full 75 g-OGTT in pregnancy. In addition, we measured DNAm using the Illumina MethylationEPIC microarray which provides a coverage of more than 850,000 CpG sites across the human genome. Some limitations are also noteworthy to mention. One CpG site was significantly associated based on a predetermined statistical threshold after correction for multiple testing; however, other findings with suggestive or nominal significance should be interpreted with caution. We investigated DNAm using blood cells; other tissues (liver, beta-cells, adipose, and brain) might be of more relevance for studies of metabolic epigenetic programming related to exposure to maternal hyperglycemia but these tissues are not easily available in healthy children for obvious ethical reasons. We have excluded participants with pre-existing diabetes prior to pregnancy, but 9% of our included population did develop GDM, which may account partially for the observed associations; however, in the past, we have showed that associations between maternal glycemic traits and offspring DNAm markers were stronger in non-GDM participants [15]. Our sample size of participants who had DNAm measurements at both time-points (birth and 5 years) was relatively small (N = 194). For this reason, we chose a statistical model (LMM) that combines the two time of measurements, which improves power to detect associations. However, we cannot exclude that other associations with smaller effect sizes were not detected. Finally, our study only included participants of European descent, which limits the application of our results to other races and ethnicities.

Conclusions

We found associations between maternal hyperglycemia in pregnancy and offspring blood DNAm levels measured at birth and at 5 years. The association was significant after accounting for multiple testing at a CpG site near the FSD1L gene. Notably, our study also reports two CpG sites, located in the promotor region of PRMD16—a gene implicated in adipose tissue regulation, and that we had previously observed associations between maternal hyperglycemia and placenta DNAm. We also confirmed some of our suggestive associations with other maternal glycemic traits, predominantly with 1 h-glucose and 2 h-glucose post-OGTT. Overall, our study provides results supporting that maternal hyperglycemia might be implicated in offspring epigenetic programming which may last until 5 years after birth.

Supplementary Information

Additional file 1: Genomic inflation factor or lambda for each LMM and linear regression models testing associations between maternal hyperglycemia outcomes and DNAm in cord blood and/or in blood at 5 years of age; Table presenting lambda for each LMM and linear regression models.

Additional file 2: Differentially methylated regions based on results from LMM testing associations between maternal AUCglu and DNAm measured in cord blood and blood at 5 years of age, using the DMRff method; Table presenting the differentially methylated regions, including their chromosome number, genomic position, number of CpGin the region, estimate value, standard error, P-value and adjusted P-value, based on results from LMM testing the association between maternal AUCglu and DNAm measured in cord blood at 5 years of age.

Additional file 3: CpG sites identified based on LMM or linear regression models testing associations between maternal AUCglu and DNAm in cord blood and/or in blood at 5 years of age; Table presenting CpG sites, including their chromosome number, genomic position, associated gene and regression beta and p-value, identified at suggestive P < 10−5 based on LMM or linear regression models testing association between maternal AUCglu and DNAm in cord blood and/or in blood at 5 years of age.

Additional file 4: Epigenome-wide association plot for maternal fasting glucose levels; Figure of the epigenome-wide association plot presenting CpG sites from the LMM associated with maternal fasting glucose.

Additional file 5: Epigenome-wide association plot for maternal 1h-post-OGTT glucose levels; Figure of epigenome-wide association plot presenting CpG sites from the LMM associated with maternal 1h-post-OGTT glucose levels.

Additional file 6: Epigenome-wide association plot for maternal 2h-post-OGTT glucose levels; Figure of epigenome-wide association plot presenting CpG sites from the LMM associated with maternal 2h-post-OGTT glucose levels.

Additional file 7: CpG sites identifiedin linear mixed models testing associations between maternal fasting glucose and DNAm measured in cord blood and blood at 5 years of age; Table presenting CpG sites, including their chromosome number, genomic position, and associated gene, identified at suggestive P < 10−5 in linear mixed models testing associations between maternal fasting glucose and DNAm measured in cord blood and blood at 5 years of age.

Additional file 8: CpG sites identified in linear mixed models testing associations between maternal 1h-glucose post-OGTT and DNAm measured in cord blood and blood at 5 years of age; Table presenting CpG sites, including their chromosome number, genomic position, and associated gene, identified at suggestive P < 10−5 in linear mixed models testing associations between maternal 1h post-OGTT glucose levels and DNAm measured in cord blood and blood at 5 years of age.

Additional file 9: CpG sites identified in linear mixed models testing associations between maternal 2h-glucose post-OGTT and DNAm measured in cord blood and blood at 5 years of age; Table presenting CpG sites, including their chromosome number, genomic position, and associated gene, identified at suggestive P < 10−5 in linear mixed models testing associations between maternal 2h post-OGTT glucose levels and DNAm measured in cord blood and blood at 5 years of age.

Abbreviations

ATP Adenosine triphosphate

AUCglu Area under the curve for glucose

BAT Brown adipose tissue

BMI Body mass index

CHUS Centre hospitalier universitaire de Sherbrooke

CpG Cytosine–phosphate–guanine

Chr Chromosome

DNAm DNA methylation

DOHaD Developmental origin of health and disease hypothesis

EWAS Epigenome-wide association study

FG Fasting glucose

Gen3G Genetics of glucose regulation in gestation and growth cohort

GCT 50-G glucose challenge test

GDM Gestational diabetes mellitus

GWAS Genome-wide association study

iBAT Beige adipose tissue

LMM Linear mixed model

MAF Minor allele frequency

nRBC Nucleated red blood cells

OGTT Oral glucose tolerance test

OAT Omental adipose tissue

PGC-1α PPAR-y coactivator-1α

PPAR-y Peroxisome proliferator-activated receptor-y

QC Quality control

RCP Regression on correlated probes

SD Standard deviation

SNP Single-nucleotide polymorphism

UCP1 Uncoupling protein 1

WAT White adipose tissue

Acknowledgements

We first and foremost thank the participants of the Gen3G cohort for their commitment and support over the last 12 years. We also wish to thank Chamaine Washington for editing the English language.

Author contributions

M-FH, SML, DJQ, LB, and PP conceptualized the study; AT, KT, CA, DJQ, LB, and M-FH wrote the main manuscript text; CA prepared Figs. 2 and 3 and supplementary figures; CA performed the statistical analyses; LB and M-FH supervised the work; All authors reviewed the manuscript.

Funding

This work was supported by a Fonds de recherche du Québec - Santé (FRQS) operating grant (to M-FH, grant #20697); a Canadian Institute of Health Research (CIHR) operating grants (to M-FH grant #MOP 115071 and to LB #PJT-152989); AT was supported by Diabète Québec. KT was recipient of a doctoral research award from Fonds de la recherche du Québec en santé (FRQS). LB is a senior research scholar from the FRQS.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of CIUSSS de l’Estrie-CHUS (protocol code #MP-04-2012-189 and September 2011). Informed consent was obtained from all subjects involved in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher's Note

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

Amélie Taschereau and Kathrine Thibeault contributed equally.

Sharon M. Lutz, Luigi Bouchard and Marie-France Hivert have supervised the work equally.
==== Refs
References

1. Quintanilla Rodriguez BS, Mahdy H. Gestational Diabetes. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing; 2023 [cited 2023 May 17]. Available from: http://www.ncbi.nlm.nih.gov/books/NBK545196/
2. Wang H Li N Chivese T Werfalli M Sun H Yuen L IDF diabetes atlas: estimation of global and regional gestational diabetes mellitus prevalence for 2021 by international association of diabetes in pregnancy study group’s criteria Diabetes Res Clin Pract 2022 183 109050 10.1016/j.diabres.2021.109050 34883186
3. Pu J Zhao B Wang EJ Nimbal V Osmundson S Kunz L Racial/ethnic differences in gestational diabetes prevalence and contribution of common risk factors: racial/ethnic differences in GDM risk factors Paediatr Perinat Epidemiol 2015 29 5 436 443 10.1111/ppe.12209 26201385
4. Gupta Y Kalra B Screening and diagnosis of gestational diabetes mellitus JPMA J Pak Med Assoc 2016 66 9 Suppl 1 S19 21 27582144
5. 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.
6. On behalf of the HAPO Follow-up Study Cooperative Research Group, Lowe WL, Lowe LP, Kuang A, Catalano PM, Nodzenski M, et al. Maternal glucose levels during pregnancy and childhood adiposity in the hyperglycemia and adverse pregnancy outcome follow-up study. Diabetologia. 2019;62(4):598–610.
7. Vääräsmäki M Pouta A Elliot P Tapanainen P Sovio U Ruokonen A Adolescent manifestations of metabolic syndrome among children born to women with gestational diabetes in a general-population birth cohort Am J Epidemiol 2009 169 10 1209 1215 10.1093/aje/kwp020 19363101
8. Li Z Wu Y Du B Yu X Wang H Niu Y Associations of maternal gestational diabetes mellitus with alterations in cardiovascular system in early childhood Diabetes Metab Res Rev 2022 10.1002/dmrr.3551 36565185
9. Adgent MA Gebretsadik T Reedus J Graves C Garrison E Bush N Gestational diabetes and childhood asthma in a racially diverse US pregnancy cohort Genuneit J, editor Pediatr Allergy Immunol 2021 32 6 1190 1196 10.1111/pai.13523 33894077
10. Xiang AH Association of maternal diabetes with autism in offspring JAMA 2017 317 5 537 538 10.1001/jama.2016.20122 28170476
11. Wang P Xie J Jiao XC Ma SS Liu Y Yin WJ Maternal glycemia during pregnancy and early offspring development: a prospective birth cohort study J Clin Endocrinol Metab 2021 106 8 2279 2290 10.1210/clinem/dgab331 33982055
12. Heindel JJ Vandenberg LN Developmental origins of health and disease: a paradigm for understanding disease cause and prevention Curr Opin Pediatr 2015 27 2 248 253 10.1097/MOP.0000000000000191 25635586
13. Wolffe AP Matzke MA Epigenetics: regulation through repression Science 1999 286 5439 481 486 10.1126/science.286.5439.481 10521337
14. Moore LD Le T Fan G DNA methylation and its basic function Neuropsychopharmacol Off Publ Am Coll Neuropsychopharmacol 2013 38 1 23 38 10.1038/npp.2012.112
15. Tobi EW Juvinao-Quintero DL Ronkainen J Ott R Alfano R Canouil M Maternal glycemic dysregulation during pregnancy and neonatal blood dna methylation: meta-analyses of epigenome-wide association studies Diabetes Care 2022 45 3 614 623 10.2337/dc21-1701 35104326
16. Gagné-Ouellet V Houde AA Guay SP Perron P Gaudet D Guérin R Placental lipoprotein lipase DNA methylation alterations are associated with gestational diabetes and body composition at 5 years of age Epigenetics 2017 12 8 616 625 10.1080/15592294.2017.1322254 28486003
17. Allard C Desgagné V Patenaude J Lacroix M Guillemette L Battista MC Mendelian randomization supports causality between maternal hyperglycemia and epigenetic regulation of leptin gene in newborns Epigenetics 2015 10 4 342 351 10.1080/15592294.2015.1029700 25800063
18. Guillemette L Allard C Lacroix M Patenaude J Battista MC Doyon M Genetics of glucose regulation in gestation and growth (Gen3G): a prospective prebirth cohort of mother–child pairs in Sherbrooke, Canada BMJ Open 2016 6 2 e010031 10.1136/bmjopen-2015-010031 26842272
19. Hivert MF Cardenas A Allard C Doyon M Powe CE Catalano PM Interplay of placental DNA methylation and maternal insulin sensitivity in pregnancy Diabetes 2020 69 3 484 492 10.2337/db19-0798 31882564
20. Tai MM A mathematical model for the determination of total area under glucose tolerance and other metabolic curves Diabetes Care 1994 17 2 152 154 10.2337/diacare.17.2.152 8137688
21. Aryee MJ Jaffe AE Corrada-Bravo H Ladd-Acosta C Feinberg AP Hansen KD Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays Bioinforma Oxf Engl 2014 30 10 1363 1369 10.1093/bioinformatics/btu049
22. Fortin JP Labbe A Lemire M Zanke BW Hudson TJ Fertig EJ Functional normalization of 450 k methylation array data improves replication in large cancer studies Genome Biol 2014 15 12 503 10.1186/s13059-014-0503-2 25599564
23. Niu L Xu Z Taylor JA RCP: a novel probe design bias correction method for illumina methylation BeadChip Bioinformatics 2016 32 17 2659 2663 10.1093/bioinformatics/btw285 27153672
24. Juvinao-Quintero DL Starling AP Cardenas A Powe CE Perron P Bouchard L Epigenome-wide association study of maternal hemoglobin A1c in pregnancy and cord blood DNA methylation Epigenomics 2021 13 3 203 218 10.2217/epi-2020-0279 33406918
25. Pidsley R Zotenko E Peters TJ Lawrence MG Risbridger GP Molloy P Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling Genome Biol 2016 17 1 208 10.1186/s13059-016-1066-1 27717381
26. Johnson WE Li C Rabinovic A Adjusting batch effects in microarray expression data using empirical Bayes methods Biostat Oxf Engl 2007 8 1 118 127
27. Gervin K Salas LA Bakulski KM van Zelm MC Koestler DC Wiencke JK Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data Clin Epigenet 2019 11 1 125 10.1186/s13148-019-0717-y
28. Salas LA Koestler DC Butler RA Hansen HM Wiencke JK Kelsey KT An optimized library for reference-based deconvolution of whole-blood biospecimens assayed using the Illumina HumanMethylationEPIC BeadArray Genome Biol 2018 19 1 64 10.1186/s13059-018-1448-7 29843789
29. Du P Zhang X Huang CC Jafari N Kibbe WA Hou L Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis BMC Bioinf 2010 11 587 10.1186/1471-2105-11-587
30. J Pinheiro, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. [Internet]. R package version 3.1-155. 2022. Available from: https://CRAN.R-project.org/package=nlme.
31. Suderman M Staley JR French R Arathimos R Simpkin A Tilling K dmrff: identifying differentially methylated regions efficiently with power and control [Internet] Bioinformatics 2018 10.1101/508556 29931280
32. Lent S Cardenas A Rifas-Shiman SL Perron P Bouchard L Liu CT Detecting differentially methylated regions with multiple distinct associations Epigenomics 2021 13 6 451 464 10.2217/epi-2020-0344 33641349
33. Hansen KD. IlluminaHumanMethylationEPICanno.ilm10b2.hg19: Annotation for Illumina’s EPIC methylation arrays. [Internet]. R package version 0.6.0. 2016. Available from: https://bitbucket.com/kasperdanielhansen/Illumina_EPIC.
34. Abdul-Ghani MA Tripathy D DeFronzo RA Contributions of beta-cell dysfunction and insulin resistance to the pathogenesis of impaired glucose tolerance and impaired fasting glucose Diabetes Care 2006 29 5 1130 1139 10.2337/dc05-2179 16644654
35. Qiao Q Lindström J Valle TT Tuomilehto J Progression to clinically diagnosed and treated diabetes from impaired glucose tolerance and impaired fasting glycaemia Diabet Med J Br Diabet Assoc 2003 20 12 1027 1033 10.1111/j.1464-5491.2003.01054.x
36. Talluri R Shete S Evaluating methods for modeling epistasis networks with application to head and neck cancer Cancer Inform 2015 14 Suppl 2 17 23 25733798
37. Gao Y Liu J Zhao D Diao G A novel prognostic model for identifying the risk of hepatocellular carcinoma based on angiogenesis factors Front Genet 2022 13 857215 10.3389/fgene.2022.857215 35368665
38. Perry JRB Stolk L Franceschini N Lunetta KL Zhai G McArdle PF Meta-analysis of genome-wide association data identifies two loci influencing age at menarche Nat Genet 2009 41 6 648 650 10.1038/ng.386 19448620
39. Tahir MS Porto-Neto LR Gondro C Shittu OB Wockner K Tan AWL Meta-analysis of heifer traits identified reproductive pathways in Bos Indicus cattle Genes 2021 12 5 768 10.3390/genes12050768 34069992
40. Côté S Gagné-Ouellet V Guay SP Allard C Houde AA Perron P PPARGC1α gene DNA methylation variations in human placenta mediate the link between maternal hyperglycemia and leptin levels in newborns Clin Epigenet 2016 8 72 10.1186/s13148-016-0239-9
41. Chi J Cohen P The multifaceted roles of PRDM16: adipose biology and beyond Trends Endocrinol Metab TEM 2016 27 1 11 23 10.1016/j.tem.2015.11.005 26688472
42. Becerril S Gómez-Ambrosi J Martín M Moncada R Sesma P Burrell MA Role of PRDM16 in the activation of brown fat programming. Relevance to the development of obesity Histol Histopathol 2013 28 11 1411 1425 23771475
43. Ishibashi J Seale P Functions of Prdm16 in thermogenic fat cells Temp Austin Tex 2015 2 1 65 72
44. Cypess AM Lehman S Williams G Tal I Rodman D Goldfine AB Identification and importance of brown adipose tissue in adult humans N Engl J Med 2009 360 15 1509 1517 10.1056/NEJMoa0810780 19357406
45. Jiang N Yang M Han Y Zhao H Sun L PRDM16 regulating adipocyte transformation and thermogenesis: a promising therapeutic target for obesity and diabetes Front Pharmacol 2022 13 870250 10.3389/fphar.2022.870250 35462933
46. Wang L Fan H Zhou L Wu Y Lu H Luo J Altered expression of PGC-1 α and PDX1 and their methylation status are associated with fetal glucose metabolism in gestational diabetes mellitus Biochem Biophys Res Commun 2018 501 1 300 306 10.1016/j.bbrc.2018.05.010 29730292
47. Yue H He JW Ke YH Zhang H Wang C Hu WW Association of single nucleotide polymorphism Rs2236518 in PRDM16 gene with BMI in Chinese males Acta Pharmacol Sin 2013 34 5 710 716 10.1038/aps.2012.201 23524569
48. AlAmrani A AbdelKarim M AlZoghaibi M PRDM16 gene polymorphism is associated with obesity and blood lipids profiles in Saudi population J Clin Med 2018 7 6 141 10.3390/jcm7060141 29890628
49. Liu L Chen Y Chen J Lu M Guo R Han J The relationship between PRDM16 promoter methylation in abdominal subcutaneous and omental adipose tissue and obesity Clin Nutr 2021 40 4 2278 2284 10.1016/j.clnu.2020.10.016 33127174
50. Staley JR Suderman M Simpkin AJ Gaunt TR Heron J Relton CL Longitudinal analysis strategies for modelling epigenetic trajectories Int J Epidemiol 2018 47 2 516 525 10.1093/ije/dyy012 29462323
51. Shiau S Wang L Liu H Zheng Y Drong A Joyce BT Prenatal gestational diabetes mellitus exposure and accelerated offspring DNA methylation age in early childhood Epigenetics 2021 16 2 186 195 10.1080/15592294.2020.1790924 32614694