Changes in fetal DNA methylation (DNAm) of the leptin (LEP) gene have been associated with exposure to maternal hyperglycemia, but their links with childhood obesity risk are still unclear. We investigated the association between maternal hyperglycemia, placental LEP DNAm (25 5'-C-phosphate-G-3' (CpG) sites), neonatal leptinemia, and adiposity (i.e., BMI and skinfold thickness (ST) (subscapular (SS) + triceps (TR) skinfold measures, and the ratio of SS:TR) at 3-years-old, in 259 mother-child dyads, from Gen3G birth cohort. We conducted multivariate linear analyses adjusted for gestational age at birth, sex of the child, age at follow-up, and cellular heterogeneity. We assessed the causal role of DNAm in the association between maternal glycemia and childhood outcomes, using mediation analysis. We found three CpGs associated with neonatal leptinemia (p ≤ 0.002). Of these, cg05136031 and cg15758240 were also associated with BMI (β = -2.69, p = 0.05) and fat distribution (β = -0.581, p = 0.05) at 3-years-old, respectively. Maternal glycemia was associated with DNAm at cg15758240 (β = -0.01, p = 0.04) and neonatal leptinemia (β = 0.19, p = 0.004). DNAm levels at cg15758240 mediates 0.8% of the association between maternal glycemia and neonatal leptinemia (p < 0.001). Our results support that DNAm regulation of the leptin pathway in response to maternal glycemia might be involved in programming adiposity in childhood.

The aim of this study was to identify placental DNA methylation (DNAm) variations associated with adiposity at 3 years of age. We quantified placental DNAm using the Infinium MethylationEPIC BeadChips. We assessed associations between DNAm at single-CpGs and skinfold thickness using robust linear regression models adjusted for gestational age, child's sex, age at follow-up and cellular heterogeneity. We sought replication of DNAm association with child adiposity in an independent cohort. We quantified placental mRNA levels for annotated gene using qRT-PCR and tested for correlation with DNAm. Lower DNAm at cg22593959 and cg22436429 was associated with higher adiposity (β = -1.18, q = 0.002 and β = -0.82, q = 0.04). The cg22593959 is located in an intergenic region (chr7q31.3), whereas cg22436429 is within the TFAP2E gene (1p34.3). DNAm at cg22593959 and cg22436429 was correlated with mRNA levels at FAM3C (rs = -0.279, p = 0.005) and TFAP2E (rs = 0.216, p = 0.03). In an independent cohort, the association between placental DNAm at cg22593959 and childhood adiposity was of similar strength and direction (β = -3.8 ± 4.1, p = 0.36), yet non-significant. Four genomic regions were also associated with skinfold thickness within FMN1, MAGI2, SKAP2 and BMPR1B genes. We identified placental epigenetic variations associated with adiposity at 3 years of age suggesting that childhood fat accretion patterns might be established during fetal life.