Maternal pre-pregnancy BMI, offspring epigenome-wide DNA methylation, and childhood obesity: findings from the Boston Birth Cohort
BMC Medicine volume 21, Article number: 317 (2023)
Maternal pre-pregnancy obesity is an established risk factor for childhood obesity. Investigating epigenetic alterations induced by maternal obesity during fetal development could gain mechanistic insight into the developmental origins of childhood obesity. While obesity disproportionately affects underrepresented racial and ethnic mothers and children in the USA, few studies investigated the role of prenatal epigenetic programming in intergenerational obesity of these high-risk populations.
This study included 903 mother–child pairs from the Boston Birth Cohort, a predominantly urban, low-income minority birth cohort. Mother-infant dyads were enrolled at birth and the children were followed prospectively to age 18 years. Infinium Methylation EPIC BeadChip was used to measure epigenome-wide methylation level of cord blood. We performed an epigenome-wide association study of maternal pre-pregnancy body mass index (BMI) and cord blood DNA methylation (DNAm). To quantify the degree to which cord blood DNAm mediates the maternal BMI-childhood obesity, we further investigated whether maternal BMI-associated DNAm sites impact birthweight or childhood overweight or obesity (OWO) from age 1 to age 18 and performed corresponding mediation analyses.
The study sample contained 52.8% maternal pre-pregnancy OWO and 63.2% offspring OWO at age 1–18 years. Maternal BMI was associated with cord blood DNAm at 8 CpG sites (genome-wide false discovery rate [FDR] < 0.05). After accounting for the possible interplay of maternal BMI and smoking, 481 CpG sites were discovered for association with maternal BMI. Among them 123 CpGs were associated with childhood OWO, ranging from 42% decrease to 87% increase in OWO risk for each SD increase in DNAm. A total of 14 identified CpG sites showed a significant mediation effect on the maternal BMI-child OWO association (FDR < 0.05), with mediating proportion ranging from 3.99% to 25.21%. Several of these 14 CpGs were mapped to genes in association with energy balance and metabolism (AKAP7) and adulthood metabolic syndrome (CAMK2B).
This prospective birth cohort study in a high-risk yet understudied US population found that maternal pre-pregnancy OWO significantly altered DNAm in newborn cord blood and provided suggestive evidence of epigenetic involvement in the intergenerational risk of obesity.
The prevalence of childhood obesity was 19.7% and affected about 14.7 million children and adolescents in the USA in 2017–2020 , increasing from 18.5% in 2015–2016 . Childhood obesity is more common in racial minority children and has increased dramatically . Obesity remains a serious public health threat that disproportionately affects racial minority children .
Maternal pre-pregnancy obesity is an established risk factor for offspring cardiometabolic dysfunction, including obesity . Given that 25% of conceptions occur among obese women in the USA , maternal obesity before pregnancy is one of the most common intrauterine risk factors for childhood obesity. However, there is still limited mechanistic insight, especially in racial minority children. One plausible hypothesis points to epigenetic modification, as mouse models demonstrated the induction of epigenetic alterations by maternal obesity can affect the later risk of obesity in offspring [7,8,9].
However, evidence from prospective birth cohort studies connecting maternal BMI, fetal DNAm, and childhood BMI remains limited . Only five studies explored the association of maternal BMI with genome-wide DNAm in newborn cord blood [11,12,13,14,15]. The majority of previous studies were based on Illumina HumanMethylation27 BeadChip  or 450 BeadChip [12,13,14]. Notably, none of the previous findings overlapped with one another in terms of gene level. Few studies conducted on African Americans were limited by a small sample size [11, 14]. In addition, previous studies did not consider other maternal lifestyle factors like smoking, which may exacerbate the adverse DNAm alterations induced by maternal BMI, and such evidence is lacking. Most notably, there is a lack of longitudinal follow-up in children from birth to adolescence to estimate the impact of maternal pre-pregnancy BMI-related CpG methylation on offspring obesity.
To fill this gap, we examined the association between maternal pre-pregnancy BMI and epigenome-wide methylation of cord blood-derived DNA (using the latest Illumina Infinium MethylationEPIC BeadChip with more than 850,000 methylation sites) among mother–child pairs enrolled in the Boston Birth Cohort (BBC), a predominantly urban, low-income ethnic minority (Black and Hispanic) sample in Boston, MA. We then investigated the relationships between the identified maternal BMI-associated methylation sites and birthweight/child overweight or obesity (OWO) from 1 to 18 years old and further identified potential pathways by causal mediation analysis.
The present study included 903 mother–child pairs from the BBC. The BBC, a predominantly urban, low-income minority birth cohort, was initiated in 1998 with a rolling enrolment in Boston, MA. Briefly, mothers who delivered singleton live births at the Boston Medical Center (BMC) were invited to participate shortly after giving birth. Twins or triplets and newborns with major birth defects were excluded. The BBC study was designed to over-sample low birthweight (< 2.5 kg) and preterm (< 37 weeks of gestation) births at enrollment. All enrolled mothers completed a standard questionnaire interview to assess maternal socio-demographic characteristics, lifestyle (including smoking, alcohol consumption, and dietary habits), and reproductive and medical history.
Cord blood samples were obtained at delivery by trained nursing staff of the labor and delivery (L&D) service. Blood samples were then immediately stored in a designated − 20° freezers in L&D for study collection. Cord blood was centrifuged (0° Celsius for 13 min) and fractionated into plasma, white blood cells, and red blood cells. All biospecimens are stored in − 80° freezers pending analysis.
Clinical information on mothers and newborns, including birth outcomes, was obtained from their medical records. Beginning in 2004 under a separate but linked IRB protocol, maternal-child pairs were invited to continue their participation in a postnatal follow-up study if the child continued primary or specialty care at BMC. A standardized questionnaire was used to assess postnatal demographic, infant feeding, and environmental information. The study design, data collection methods, and long-term follow-up have been described and published previously .
The study protocol has received initial and annual continuation approval by the Institute Review Boards of Boston Medical Center and the Johns Hopkins Bloomberg School of Public Health. All study mothers provided written informed consent. Study children are not able to give true informed consent legally until they turn 18. Thus, we obtained assent from study children at BMC IRB determined ages. Written informed consent was obtained when the child reached age 18.
DNAm profiling in cord blood
DNA samples were isolated from EDTA-treated peripheral white blood cells and shipped to the University of Minnesota Genomics Center for methylation profiling. The quality of the DNA samples has been demonstrated in our previous study using the Illumina BeadChip . Detailed quality control process of the DNA methylation data has been previously published . Briefly, DNA methylation was measured for 963 cord blood samples (plus 21 replicates) using the Illumina Infinium MethylationEPIC BeadChip (850K). This platform interrogates a total of 865,859 CpG sites. Both cord blood samples and maternal samples (n ~ 420 samples) were randomly placed in each 96-well DNA plate during the DNA methylation profiling. The laboratory staff were blinded to the sample placement. Only methylation data measured in cord blood was used in the present study.
β-value for each CpG site was reported, ranging from 0 to 1, to signify the percentage of DNAm at each CpG site. We used minfi package  to perform existing analytic pipelines. Probes were excluded if they had bead count < 3 in 5% of samples or had > 5% of samples with a detection P > 0.01 (n = 4193) and had an annotated SNP at the measured or extension site or that overlapped with the probe or that potentially cross-hybridized to other genomic locations  (n = 140,271). Samples were excluded due to the following reasons: if they (1) were outliers detected by multidimensional scaling (MDS) analysis or a median log2 intensity value < 10 (n = 13); (2) were sex mixed-up samples (n = 7); (3) had missing rate > 0.02 across probes (n = 1).
Furthermore, MDS plots confirmed male cord blood samples, female cord blood samples, and maternal blood samples, which clustered separately, as expected. Correlations for the duplicates were computed among the 21 pairs of duplicates. The correlation for duplicate measurements on the same sample ranged from 0.990 to 0.996. We then performed the single-sample Noob (ssNoob) methods for background and dye bias correction  and performed quantile normalization to normalize type 1 and type II probes. We further excluded 39 mother–child pairs with missing maternal BMI data. These filtering process resulted in 903 samples with 721,395 CpG sites for the downstream analyses. Additional file 1: Table S1 presents the characteristics of the 903 pairs enrolled in this study compared to the remaining recruited mother–child pairs (n = 7720) in the BBC. The key variables were comparable between the 903 pairs and the remaining ones.
Assessment of perinatal and postnatal variables
Maternal pre-pregnancy weight and height assessments were based on maternal questionnaire interviews at enrollment. Maternal BMI was calculated as weight in kilograms divided by the square of height in meters and then categorized into two groups: non-OWO (< 25 kg/m2) and OWO (≥ 25 kg/m2). Maternal smoking during pregnancy was defined according to the following questions:(1) in the 6 months before you found out you were pregnant, did you smoke/use tobacco? (2) did you smoke/use tobacco in the first three months of pregnancy? (3) did you smoke/use tobacco in the middle 3 months of pregnancy? and (4) did you smoke/use tobacco in the last 3 months of pregnancy. We defined ever-smoking during pregnancy mothers as those who answered “yes” to any of the above questions and defined others as never-smoking during pregnancy mothers.
Birthweight was abstracted from the electronic medical record (EMR). Other important covariates included: child’s sex, maternal age at delivery, maternal race (Black, White, or Hispanic), maternal education level (high school or less, versus others), maternal alcohol consumption (ever drinking alcohol during pregnancy versus others), and parity (not including the index pregnancy, zero versus one or more). Gestational age was determined using an established algorithm based on both the last menstrual period and the result of early ultrasound (< 20 weeks of gestation) to maximize accuracy and minimize missing data. Details of the covariates can be found elsewhere .
Assessment of overweight or obesity (OWO) in childhood
Child weight and height were measured by medical staff during pediatric well-child visits and retrieved from the EMR. Children’s repeated measurements of BMI from birth to age 18 were collected. The first measure was used if repeated measures were taken within one year of age. BMI z-scores and percentiles for each age-window (1 year of age) were calculated based on US national reference data by age and gender . Missing percentiles were imputed using the average of the last and following observed values. OWO was defined as BMI ≥ 85th percentile of age and gender.
Figure 1 presented the overall analysis flowchart. We compared the demographic and clinical characteristics of newborns exposed to maternal pre-pregnancy OWO versus those who were not, with chi-square test for categorical variables and ANOVA for continuous variables. We used false discovery rate (FDR)  < 0.05 to determine epigenome-wide significant CpGs in relation to maternal BMI.
Differentially methylated sites associated with maternal BMI
Linear regression was applied to investigate differentially methylated sites in cord blood that were associated with maternal BMI, with methylation values as dependent variables, maternal BMI as the independent variable, and child’s sex, maternal age at delivery, maternal race, maternal education level, maternal alcohol consumption, parity, and gestational age as covariates. We also adjusted for estimated cord blood cell composition (CD4 + , CD8 + T cells, B cells, monocytes, granulocytes, natural killer cells, and nucleated red blood cells). The cell type proportion was inferred for each cord blood sample, based on external cord blood reference DNA methylation signatures of the constituent cell type from the Illumina Infinium Methylation450 BeadChip by the minfi package . To quantify latent factors, including the effects of unobserved batch effects and other unmeasured confounding factors, we used smart surrogate variable (SV) analysis by the smartSVA package . A total of 54 SVs were generated and included as covariates in the above model.
To account for the possible interplay of maternal BMI and smoking in impacting the methylation level, we performed likelihood ratio test comparing model 1 (with maternal BMI and the cross-product term of maternal BMI x smoking) to model 2 (without these two terms).
The CpG sites associated with pre-pregnancy maternal BMI were annotated to the corresponding genes based on an official EPIC array annotation file from Illumina . Further enrichment analyses were conducted using the Genetic Association Database (GAD) enclosed within the Database for Annotation, Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/).
As a first step, we performed multiple linear regression (for birthweight) or logistic regression (for child OWO) analyses and identified whether the maternal BMI-associated CpG sites (identified in EWAS) were associated with newborn birthweight and child OWO. Covariates adjusted included maternal BMI, smoking, maternal age at delivery, race, education level, alcohol consumption, child’s sex, parity, gestational age, cord blood cell compositions, and all SVs. In the analyses of child OWO, we additionally adjusted for newborn birthweight.
The subsequent analyses focused on CpG sites associated with both maternal BMI and newborn birthweight or child OWO (P < 0.05). Using parametric regression models, we performed causal mediation analysis, achieved by the R “mediation” package . Two models were estimated for each CpG: (1) a model for the mediator (β-value for each CpG site as a continuous variable) conditional on exposure (maternal BMI) and covariates (maternal smoking, age at delivery, race, education level, and alcohol consumption, child’s sex, parity, gestational age, cord blood cell compositions and all SVs); (2) a model for the outcome (newborn birthweight or child OWO) conditional on exposure, the mediator, and covariates. We allowed for the presence of exposure-mediator interactions in the outcome regression model. In our non-smoking restricted analysis, we removed the smoking covariate from the model.
We aimed to calculate how much maternal BMI association with child OWO risk (average total effect) was attributable to the mediating effect of methylation level at a specific locus (average causal mediation effect, ACME). The proportion attributable to the ACME was calculated as ACME divided by the total effect on the log odds scale, with 0 indicating no mediation effect. Stratified mediation analyses by maternal smoking status during pregnancy were also performed. We adjusted for multiple testing by the FDR method. All analyses were conducted using R software, version 3.4.1 (R Project for Statistical Computing).
A total of 903 mother–child pairs were included after data quality control steps; 641 mothers (71.0%) were Black; 478 (52.9%) children were boys. The median value of pre-pregnancy BMI was 26.62 kg/m2, with the interquartile range being 7.86 kg/m2. The study sample contained 52.8% maternal pre-pregnancy OWO and 63.2% offspring OWO at age 1–18 years, respectively. Compared to non-OWO mothers, OWO mothers were more likely to be older, Black, and had a higher proportion of smokers. The children born to OWO mothers had a higher prevalence of OWO from ages 1 to 18 years (Table 1).
Epigenome-wide association analysis of maternal BMI in cord blood
We identified 8 CpG sites, corresponding to 6 genes, associated with maternal BMI at genome-wide significance (FDR < 0.05) (Table 2 and Fig. 2). Figure 2a shows the Manhattan plot and the Quantile–Quantile plot of the association between maternal BMI and cord blood DNAm. The genomic inflation factor was 1.01 indicating little residual confounding. The maximum adjusted hypermethylated difference (standard error, SE) per one unit increase in maternal pre-pregnancy BMI was 0.00026 (4.89E − 05) for cg21701395, which was annotated to TP53INP1 (P = 2.31E − 07). The corresponding maximum hypomethylated adjusted difference (SE) was − 0.00067 (1.30E − 04) for cg10272744, which was annotated to LL22NC01-81G9.3 (P = 3.58E − 07). Other significant hypermethylated CpG sites were identified in the PTEN and KILLIN genes as indicated by the positive coefficients in relation to maternal BMI. The CpG sites in the AAGAB, ALPK1, and ERCC8 genes were significantly hypomethylated in newborns born to mothers with higher BMI as indicated by the negative coefficients. When further dividing mothers into three groups (under/normal weight, overweight, and obesity), we observed obesity mothers having the highest DNA methylation level at cg13694461 and cg06466203, lowest at cg02266725, cg10272744, and cg23151800 (Additional file 2: Fig. S1). Accounting for the possible interplay of maternal BMI and smoking in impacting the methylation level, we performed a likelihood ratio test and identified 481 CpG sites with altered DNAm in newborns that were associated with maternal BMI (FDR < 0.05). Additional genomic information on these CpG sites was provided in Additional file 3: Table S2. Manhattan plot indicated that significant associations were distributed across the genome (Fig. 2b). The elevated tail of the QQ-plot indicated enrichment of robust association and the genomic inflation factor 1.19 indicated small residual confounding or polygenic effect. Enrichment analysis of the 481 CpG sites revealed that their annotated genes were enriched in genes associated with triglycerides, hematocrit, and tobacco use disorder (FDR < 0.05; Additional file 1: Table S3). We identified ten genes with at least two significantly altered DNAm sites: DIP2C, ADARB2, LAMA3, MAP4K4, MCEE, PLEKHG1, PTPRN2, RIN3, SMTNL2, and THSD4.
Maternal BMI-associated CpGs and birthweight
Of the total 481 maternal BMI-associated CpG sites, 26 were associated with newborn birthweight in this study (P value < 0.05). The adjusted differences in newborn birthweight per 1 SD increment of cord blood DNAm ranged from -92.5 g for cg04107082 (P = 0.031) to 64.9 g for cg26779585 (P = 0.003). Then, out of the total 26 CpGs, we estimated the mediation effect of each single CpG on the maternal BMI-birthweight association. After correction for multiple testing, we did not detect any significant mediation effect (Additional file 1: Table S4).
Maternal BMI-associated CpGs and child OWO
Additional file 1: Table S5 presents the associations between the identified 481 maternal BMI-associated CpG sites and child OWO, with adjustment for maternal BMI, birthweight, and other covariates. A total of 123 CpGs showed a significant association with child OWO in at least one age-window (P value < 0.05). The adjusted odds ratio (ORs) per 1 SD increment of cord blood DNAm ranged from 0.58 for cg07686439 (P = 0.005, OWO at 3 years old) to 1.87 for cg08469255 (P = 0.015, OWO at 4 years old).
In the following mediation analysis, we focused on the CpGs in association with child OWO, and found a total of 14 CpG sites showed a significant mediation effect on the maternal BMI-child OWO association in at least one age group (FDR < 0.05) (Additional file 2: Fig. S2). As presented in Table 3, 21.69%, 17.41%, and 14.88% of the maternal BMI-associated child OWO risk at age 3 years were mediated through methylation levels at cg02059896 (P < 0.001), cg14257335 (P = 0.004), and cg15059222 (P = 0.003), respectively. A total of 8 CpG sites mediated the association of maternal BMI with child OWO at age 7–9 years, with the mediating proportions ranging from 4.38% (cg06292624, P = 0.007) to 25.21% (cg14015044, P = 0.009). At the time window of age 13 years, we also observed 5 CpG sites mediating the association of maternal BMI with child OWO. The mediating proportions ranged from 3.99% at cg06292624 (P = 0.012) to 10.57% at cg08391482 (P = 0.011).
To test whether the mediation effect differs among smoking and non-smoking participants, we restricted the mediation analyses to never-smoking during pregnancy mother–child pairs (n = 739). Of note, we did not observe a significant mediation effect for cg04347874, cg09545197, cg08391482, cg02059896, and cg15059222 in the non-smoking group. The main mediation results for the overall sample appeared to be driven by smoking in these mother–child pairs for these 5 CpG sites (Table 3). The rest of CpG sites tended to have comparable mediating effects on the maternal BMI-child OWO association in both groups.
In a population of predominantly Black mother–newborn pairs who are low-income and living in disinvested communities in Boston, Massachusetts, we identified 481 differentially methylated CpG sites in the newborn cord blood associated with maternal pre-pregnancy BMI. This study extends previous work based on the Illumina Infinium Methylation27 or 450 BeadChip, which were biased toward promoter regions and missed more dynamic methylation sites, including the enhancers. To the best of our knowledge, this is the first EWAS using the latest Illumina Infinium MethylationEPIC BeadChip and the largerst EWAS in a racial minority population and the first prospective study to link maternal BMI-associated cord blood DNA methylation with children’s BMI from birth to adolescence. Further mediation analyses revealed multiple CpG sites mediated the impact of maternal BMI on offspring obesity, making contributions to understanding intergenerational transmission of overweight and obesity.
Childhood obesity prevalence in the US was 12.7, 20.7, and 22.2% among 2–5, 6–11, and 12–19 years old in 2017–2020, respectively . The percentage of obese offspring in the present study was 31.4, 34.7, and 32.0% among the age group of 2–5, 6–11, and 12–18 years old. The BBC cohort mainly recruited low-income Hispanic and non-Hispanic Black mother–child pairs, which have the highest rates of obesity in the US. Compared to the cross-sectional study design, we can capture more obese children through repeated measurements across the pediatric age range.
In our EWAS of maternal pre-pregnancy BMI and cord blood DNA methylation, we replicated several associations at the gene level reported in previous studies: PRR16 , PTPRN2, SEPT9, XXYLT1 , and GFI1 . Some of these findings were supported by existing literature and biologically plausible. For example, PTPRN2 encodes a protein involved in regulating insulin secretion . The specific substrates and binding partners suggest its function for metabolism, including obesity and type 2 diabetes [29, 30]. We identified two CpGs of PTPRN2 gene with significant association with maternal BMI. A previous case–control study (252 obese children and 386 controls) showed an epigenetic association of PTPRN2 gene with childhood obesity under certain genetic backgrounds . PTPRN2 has been reported in relation to BMI in GWAS [32,33,34]. Together with our result, these literatures suggest that the methylation level of PTPRN2 might be on the pathway from maternal BMI to offspring obesity.
For the first time, we quantitatively estimated how much of the effects of maternal BMI on offspring long-term obesity risk are mediated through DNAm alterations in newborn cord blood. The mediated proportions were up to > 25% (cg14015044, TNFRSF10C). We summarized the annotated or nearest annotated gene of the identified CpGs with significant mediation effects in our study and the previous GWAS findings (Fig. 3 and Additional file 1: Table S6) [33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58]. Three identified CpG sites were linked to physical measurements in previous GWAS studies. CpG cg06292624 was located at the gene AKAP7, a member of scaffolding proteins that bind to protein kinase A (PKA), which plays a vital role in intracellular PKA-dependent signaling pathways . Both mouse and human models demonstrated a central role for the PKA signaling pathway in regulating energy balance and metabolism across multiple systems . Previous GWAS identified that the AKAP7 SNPs were significantly associated with body weight , height , and hip circumference adjusted for BMI . The methylation level at CpG cg06292624 was found to be associated with the risk of child OWO in our study. Further mediation analysis showed that ~ 5% of the offspring obesity risk related to maternal BMI was mediated through methylation at cg06292624, suggesting that such epigenetic regulation related to maternal BMI might influence the subsequent risk of obesity in offspring.
Four of the total 14 identified CpGs mediating the association between maternal BMI and child’s OWO were mapped to genes that have been reported in association with adulthood metabolic syndrome in previous GWAS studies. CpG cg15059222 maps to the CAMK2B (Calcium/Calmodulin-Dependent Protein Kinase Type II Beta Chain) gene. Among its related pathways are cytokine signaling in the immune system and the immune response IFN gamma signaling pathway . The perturbation of this process may lead to impaired homeostatic responses and possibly disease states, including obesity and diabetes [63, 64]. The methylation level at CpG cg15059222 was found to be associated with the risk of obesity at 3 years old in our study. Further mediation analysis noted that more than 14% of the increased risk of early-onset (at age 3 years) OWO related to maternal pre-pregnancy BMI was mediated through methylation level at cg15059222. SNPs in CAMK2B have been consistently linked to birthweight and risk of type 2 diabetes in different populations [34, 50,51,52,53,54,55,56,57]. Taken together, maternal BMI-induced epigenetic modification of DNA may play an important part in the underlying pathway in the intergenerational transmission of obesity and in adult metabolic syndrome.
There are considerable variations in the degree of mediation across the subgroups of smoking status for several CpG sites. The mediation effect appears to be stronger in smoking mother–child pairs for the following CpGs: cg04347874 (annotated gene: NKX2-1), cg09545197 (KCNA3), cg08391482 (SERBP1), cg02059896 (nearest gene: PIKFYVE), and cg15059222 (CAMK2B) than in the non-smoking mother–child pairs. The combination of smoking and OWO aggravates the adverse effect on plasma cholesterol levels . As animal studies have indicated, maternally driven cholesterol is an important source of cholesterol for the fetus in the first trimester of pregnancy . Fetus cholesterol disturbance-induced DNAm alterations are possibly augmented. Therefore, our study suggested that maternal obesity and smoking may both be involved in the same pathway that causes offspring obesity through maternal cholesterol disturbances. Similarly, maternal overweight and smoking have been found to have a synergistic adverse effect on the development of congenital heart anomalies in offspring . We furthermore suggest that future studies should include interaction calculations when exploring maternal-BMI altered DNAm.
Our study had several strengths. The longitudinal follow-up in children from birth to adolescence and the usage of the latest DNA methylation array that provides comprehensive genome-wide coverage and captures additional enhancer and intergenic regions of the genome allowed us to make insight into possible unrevealed biological pathways of transgenerational amplification of obesity. The BBC, a prospective study with long-term postnatal follow-up (from birth to adolescence), is mainly an under-represented, urban, low-income Black population. Our EWAS findings were based on the largest sample size among a racial minority population. Findings from this study may directly relate to understanding and preventing obesity and reducing health disparities in people of underrepresented race in research and have the highest rates of obesity. In addition, the interplay of maternal smoking and pre-pregnancy BMI in impacting fetal DNAm could be leveraged to offer new insight to identify maternal BMI-associated methylation sites in future studies. Such approaches are infrequently studied beyond locus-specific interactions due to computational burden .
Our study also had some limitations. First, residual confounding due to mother’s nutritional information, substance abuse, exposure to air pollution, children’s nutrition, and physical activity during follow-up is possible, although we have made a comprehensive adjustment for preselected potential confounders, including cell type proportion (primary confounders in the previous study ) and also used the recommended SVA method to remove the unknown confounding effect. The genomic control factor was only 1.01 for analysis without considering smoking interaction and 1.19 for integrating potential interaction effect. It indicated small residual confoundings or polygenic effect. The elevated tail of the QQ-plot also suggested enrichment of true signals among the top CpGs. Second, maternal pre-pregnancy BMI was calculated based on self-reported data by peripartum/postpartum interviews and may be less reliable than such information obtained pre-pregnancy or in early pregnancy. However, a previous study based on the BBC mothers observed self-reporting pre-pregnancy BMI at peripartum/postpartum interviews was highly correlated with the first trimester BMI obtained from the obstetric clinic (Pearson correlation coefficient = 0.87) , indicating the accuracy of reported BMI. Third, both a strength and a limitation of the BBC is that it is conducted in a population of underrepresented race; while findings from the BBC may not be generalizable to other populations, the findings apply to those most impacted by obesity and its long-term consequences. Fourthly, the observational studies are not strong enough to yield standalone causal inference compared to randomized controlled trials. Future studies are warranted to validate our findings and unravel the underlying mechanisms of identified CpGs. Fifthly, DNAm data was available only at birth. Future studies with multiple time points from birth to adolescence are expected to provide insights into the role of dynamic changes in methylation and expression levels in the progress of obesity. Future studies are also expected to expand our research and integrate critical contributors to childhood obesity, including genome, epigenome, and metabolome profile. Lastly, emerging evidence suggested the same topologically associating domains (TADs) demonstrated similar gene expression and epigenetic modification profiles . The DNAm itself is also suggested to have impact on the different layers of chromatin organization . Integrating the effect of TADs and DNAm patterns in the same population is promising to unravel the complex associations.
We presented novel findings on associations of maternal pre-pregnancy BMI with altered DNAm at 481 CpGs in newborn cord blood and the link of altered DNAm with future development of OWO among this predominantly urban, low-income racialized minority population. Our findings highlight the importance of promoting healthy weight and smoking cessation among women who wish to become pregnant. Our findings lent further support for the potential role of epigenetic modification in the intergenerational OWO. Our findings, if futher confirmed, may also help identify theraputic targets for the early-onset obesity. Future studies are warranted to validate and elucidate the functional mechanisms of our findings.
Availability of data and materials
The data underlying this article will be available on reasonable request to the corresponding author and conditional on Institutional Review Board approval.
Average causal mediation effect
Boston Birth Cohort
Boston Medical Center
Body mass index
- CpG sites:
Cytosine-phosphodiester bond-guanine (CpG) dinucleotides
Electronic medical record
Epigenome-wide association study
False discovery rate
Labor and delivery
Overweight or obesity
Single nucleotide polymorphism
Centers for Disease Control and Prevention. Childhood Obesity Facts | Overweight & Obesity | CDC. 2022. https://www.cdc.gov/obesity/data/childhood.html. Accessed 05 May 2023.
Hales CM, Fryar CD, Carroll MD, Freedman DS, Ogden CL. Trends in obesity and severe obesity prevalence in US youth and adults by sex and age, 2007–2008 to 2015–2016. JAMA. 2018;319(16):1723–5.
Madsen KA, Weedn AE, Crawford PB. Disparities in peaks, plateaus, and declines in prevalence of high BMI among adolescents. Pediatrics. 2010;126(3):434–42.
Subica AM. Public health interventions in the era of childhood obesity and other racial health disparities. Am J Public Health. 2018;108(9):1125–7.
Yu Z, Han S, Zhu J, Sun X, Ji C, Guo X. Pre-pregnancy body mass index in relation to infant birth weight and offspring overweight/obesity: a systematic review and meta-analysis. PLoS ONE. 2013;8(4): e61627.
Branum AM, Kirmeyer SE, Gregory ECW. Prepregnancy body mass index by maternal characteristics and state: data from the birth certificate, 2014. Natl Vital Stat Rep. 2016;65(6):1–11.
Jansson N, Rosario FJ, Gaccioli F, et al. Activation of placental mTOR signaling and amino acid transporters in obese women giving birth to large babies. J Clin Endocrinol Metab. 2013;98(1):105–13.
Li CCY, Young PE, Maloney CA, et al. Maternal obesity and diabetes induces latent metabolic defects and widespread epigenetic changes in isogenic mice. Epigenetics. 2013;8(6):602–11.
Masuyama H, Hiramatsu Y. Effects of a high-fat diet exposure in utero on the metabolic syndrome-like phenomenon in mouse offspring through epigenetic changes in adipocytokine gene expression. Endocrinology. 2012;153(6):2823–30.
Reichetzeder C. Overweight and obesity in pregnancy: their impact on epigenetics. Eur J Clin Nutr. 2021;75(12):1710–22.
Liu X, Chen Q, Tsai H-J, et al. Maternal preconception body mass index and offspring cord blood DNA methylation: exploration of early life origins of disease. Environ Mol Mutagen. 2014;55(3):223–30.
Sharp GC, Lawlor DA, Richmond RC, et al. Maternal pre-pregnancy BMI and gestational weight gain, offspring DNA methylation and later offspring adiposity: findings from the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2015;44(4):1288–304.
Sharp GC, Salas LA, Monnereau C, 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.
Martin CL, Jima D, Sharp GC, et al. Maternal pre-pregnancy obesity, offspring cord blood DNA methylation, and offspring cardiometabolic health in early childhood: an epigenome-wide association study. Epigenetics. 2019;14(4):325–40.
Ghildayal N, Fore R, Lutz SM, et al. Early-pregnancy maternal body mass index is associated with common DNA methylation markers in cord blood and placenta: a paired-tissue epigenome-wide association study. Epigenetics. 2021;0(0):1–11.
Pearson C, Bartell T, Wang G, et al. Boston Birth Cohort profile: rationale and study design. Precis Nutr. 2022;1(2): e00011.
Hong X, Hao K, Ji H, et al. Genome-wide approach identifies a novel gene-maternal pre-pregnancy BMI interaction on preterm birth. Nat Commun. 2017;8(1):15608.
Xu R, Hong X, Zhang B, et al. DNA methylation mediates the effect of maternal smoking on offspring birthweight: a birth cohort study of multi-ethnic US mother-newborn pairs. Clin Epigenetics. 2021;13(1):47.
Aryee MJ, Jaffe AE, Corrada-Bravo H, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinforma Oxf Engl. 2014;30(10):1363–9.
Pidsley R, Zotenko E, Peters TJ, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208.
Fortin J-P, Triche TJ Jr, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33(4):558–60.
National Center for Health Statistics. CDC growth charts, United States. 2022. http://www.cdc.gov/growthcharts/. Accessed 05 May 2023.
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.
Bakulski KM, Feinberg JI, Andrews SV, et al. DNA methylation of cord blood cell types: applications for mixed cell birth studies. Epigenetics. 2016;11(5):354–62.
Chen J, Behnam E, Huang J, et al. Fast and robust adjustment of cell mixtures in epigenome-wide association studies with SmartSVA. BMC Genomics. 2017;18(1):413.
Illumina. Infinium MethylationEPIC Product files. 2017. https://support.illumina.com/downloads/infinium-methylationepic-v1-0-product-files.html. Accessed 01 Jan 2022.
Tingley D, Yamamoto T, Hirose K, et al. Mediation: R Package for Causal Mediation Analysis. J Stat Softw. 2014;59(5):1–38.
National Center for Biotechnology Information. PTPRN2 protein tyrosine phosphatase receptor type N2 [Homo sapiens (human)]. 2023. https://www.ncbi.nlm.nih.gov/gene/5799. Accessed 05 May 2023.
Cui L, Yu W-P, DeAizpurua HJ, Schmidli RS, Pallen CJ. Cloning and characterization of islet cell antigen-related protein-tyrosine phosphatase (PTP), a novel receptor-like PTP and autoantigen in insulin-dependent diabetes. J Biol Chem. 1996;271(40):24817–23.
Bence KK, editor. Protein tyrosine phosphatase control of metabolism. New York: Springer; 2013.
Lee S. The association of genetically controlled CpG methylation (cg158269415) of protein tyrosine phosphatase, receptor type N2 (PTPRN2) with childhood obesity. Sci Rep. 2019;9(1):4855.
Zhu Z, Guo Y, Shi H, et al. Shared genetic and experimental links between obesity-related traits and asthma subtypes in U.K. Biobank. J Allergy Clin Immunol. 2020;145(2):537–49.
Kichaev G, Bhatia G, Loh P-R, et al. Leveraging polygenic functional enrichment to improve GWAS power. Am J Hum Genet. 2019;104(1):65–75.
Sakaue S, Kanai M, Tanigawa Y, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53(10):1415–24.
Christakoudi S, Evangelou E, Riboli E, Tsilidis KK. GWAS of allometric body-shape indices in UK Biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci Rep. 2021;11(1):10688.
Liu J, Zhou Y, Liu S, et al. The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease. Hum Genet. 2018;137(6–7):553–67.
Chen M-H, Raffield LM, Mousas A, et al. Trans-ethnic and ancestry-specific blood-cell genetics in 746,667 individuals from 5 global populations. Cell. 2020;182(5):1198-1213.e14.
Thompson A, King K, Morris AP, Pirmohamed M. Assessing the impact of alcohol consumption on the genetic contribution to mean corpuscular volume. Hum Mol Genet. 2021;30(21):2040–51.
Vuckovic D, Bao EL, Akbari P, et al. The polygenic and monogenic basis of blood traits and diseases. Cell. 2020;182(5):1214-1231.e11.
Rietveld CA, Esko T, Davies G, et al. Common genetic variants associated with cognitive performance identified using the proxy-phenotype method. Proc Natl Acad Sci U S A. 2014;111(38):13790–4.
Treutlein J, Cichon S, Ridinger M, et al. Genome-wide association study of alcohol dependence. Arch Gen Psychiatry. 2009;66(7):773–84.
Frank J, Cichon S, Treutlein J, et al. Genome-wide significant association between alcohol dependence and a variant in the ADH gene cluster. Addict Biol. 2012;17(1):171–80.
Ellinghaus E, Stanulla M, Richter G, et al. Identification of germline susceptibility loci in ETV6-RUNX1-rearranged childhood acute lymphoblastic leukemia. Leukemia. 2012;26(5):902–9.
Charles BA, Hsieh MM, Adeyemo AA, et al. Analyses of genome wide association data, cytokines, and gene expression in African-Americans with benign ethnic neutropenia. PLoS ONE. 2018;13(3): e0194400.
Astle WJ, Elding H, Jiang T, et al. The allelic landscape of human blood cell trait variation and links to common complex disease. Cell. 2016;167(5):1415-1429.e19.
Sinnott-Armstrong N, Tanigawa Y, Amar D, et al. Genetics of 35 blood and urine biomarkers in the UK Biobank. Nat Genet. 2021;53(2):185–94.
Donati G, Dumontheil I, Pain O, Asbury K, Meaburn EL. Evidence for specificity of polygenic contributions to attainment in English, maths and science during adolescence. Sci Rep. 2021;11(1):3851.
Richardson TG, Sanderson E, Palmer TM, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable Mendelian randomisation analysis. PLoS Med. 2020;17(3): e1003062.
Vujkovic M, Keaton JM, Lynch JA, et al. Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis. Nat Genet. 2020;52(7):680–91.
Mahajan A, Taliun D, Thurner M, 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.
Chen J, Spracklen CN, Marenne G, et al. The trans-ancestral genomic architecture of glycemic traits. Nat Genet. 2021;53(6):840–60.
Almoguera B, Vazquez L, Mentch F, et al. Identification of Four Novel Loci in Asthma in European American and African American Populations. Am J Respir Crit Care Med. 2017;195(4):456–63.
Plotnikov D, Williams C, Guggenheim JA. Association between birth weight and refractive error in adulthood: a Mendelian randomisation study. Br J Ophthalmol. 2020;104(2):214–9.
Evangelou E, Kerkhof HJ, Styrkarsdottir U, et al. A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip. Ann Rheum Dis. 2014;73(12):2130–6.
Kristiansson K, Perola M, Tikkanen E, et al. Genome-wide screen for metabolic syndrome susceptibility Loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits. Circ Cardiovasc Genet. 2012;5(2):242–9.
Giri A, Hellwege JN, Keaton JM, et al. Trans-ethnic association study of blood pressure determinants in over 750,000 individuals. Nat Genet. 2019;51(1):51–62.
Wojcik GL, Graff M, Nishimura KK, et al. Genetic analyses of diverse populations improves discovery for complex traits. Nature. 2019;570(7762):514–8.
Lee JJ, Wedow R, Okbay A, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–21.
Suryavanshi SV, Jadhav SM, McConnell BK. Polymorphisms/mutations in A-kinase anchoring proteins (AKAPs): role in the cardiovascular system. J Cardiovasc Dev Dis. 2018;5(1):7.
London E, Bloyd M, Stratakis CA. PKA functions in metabolism and resistance to obesity: lessons from mouse and human studies. J Endocrinol. 2020;246(3):R51-64.
Christakoudi S, Evangelou E, Riboli E, Tsilidis KK. GWAS of allometric body-shape indices in U.K. Biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci Rep. 2021;11(1):10688.
Nair JS, DaFonseca CJ, Tjernberg A, et al. Requirement of Ca2+ and CaMKII for Stat1 Ser-727 phosphorylation in response to IFN-gamma. Proc Natl Acad Sci U S A. 2002;99(9):5971–6.
Icer MA, Gezmen-Karadag M. The multiple functions and mechanisms of osteopontin. Clin Biochem. 2018;59:17–24.
Guzik TJ, Skiba DS, Touyz RM, Harrison DG. The role of infiltrating immune cells in dysfunctional adipose tissue. Cardiovasc Res. 2017;113(9):1009–23.
Akbartabartoori M, Lean MEJ, Hankey CR. Smoking combined with overweight or obesity markedly elevates cardiovascular risk factors. Eur J Cardiovasc Prev Rehabil. 2006;13(6):938–46.
Burke KT, Colvin PL, Myatt L, Graf GA, Schroeder F, Woollett LA. Transport of maternal cholesterol to the fetus is affected by maternal plasma cholesterol concentrations in the golden Syrian hamster. J Lipid Res. 2009;50(6):1146–55.
Baardman ME, Kerstjens-Frederikse WS, Corpeleijn E, et al. Combined adverse effects of maternal smoking and high body mass index on heart development in offspring: evidence for interaction? Heart Br Card Soc. 2012;98(6):474–9.
Salas LA, Peres LC, Thayer ZM, et al. A transdisciplinary approach to understand the epigenetic basis of race/ethnicity health disparities. Epigenomics. 2021;13(21):1761–70.
Kumar R, Story RE, Pongracic JA, et al. Maternal pre-pregnancy obesity and recurrent wheezing in early childhood. Pediatr Allergy Immunol Pulmonol. 2010;23(3):183–90.
Akdemir KC, Le VT, Chandran S, et al. Disruption of chromatin folding domains by somatic genomic rearrangements in human cancer. Nat Genet. 2020;52(3):294–305.
Sun Z, Wang Y, Song Z, et al. DNA methylation in transposable elements buffers the connection between three-dimensional chromatin organization and gene transcription upon rice genome duplication. J Adv Res. 2022;42:41–53.
The authors would like to thank all study participants in the Boston Birth Cohort. The authors also acknowledge the nursing staff at Labor and Delivery of the Boston Medical Center and the field team for their contributions to the Boston Birth Cohort.
The Boston Birth Cohort (the parent study) was supported in part by the National Institutes of Health (NIH) grants (R01HD098232-01, R21ES011666, 2R01HD041702, R21HD066471, R01HD086013, R01HD098232, R01ES031272, R21AI154233, R01ES031521, and U01 ES034983) and the Health Resources and Services Administration (HRSA) of the U.S. Department of Health and Human Services (HHS) (UT7MC45949). This information or content and conclusions are those of the authors and should not be construed as the official position or policy of nor should any endorsements be inferred by any funding agencies.
Ethics approval and consent to participate
The work described has been carried out in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) for studies involving humans. The study sample consists of understudied and underrepresented US Black, indigenous, and people of color (BIPOC) mothers and newborns. Written informed consent was obtained from all the participating mothers before their enrollment to the Boston Birth Cohort and assent was obtained from the participating children at institutional review boards (IRB)-determined ages. Written informed consent was obtained from the child at age 18. The informed consent, protocol, and data collection instruments were reviewed and approved by the IRB of the Boston Medical Center (H-23525, H-23237) and Johns Hopkins Bloomberg School of Public Health (3966/CR1489). The privacy rights of human participants have been observed, including obtaining an NIH certificate of confidentiality.
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.
Dr. Jiahui Si is currently affiliated with National Institute of Health Data Science at Peking University, Peking University, Beijing, China and is supported by the China Postdoctoral Science Foundation (88014Y0215).
Table S1. [Characteristics of mother-newborn pairs included versus excluded from this current study]. Table S3. [Gene enrichment analysis of 481 probes]. Table S4. [Association of maternal BMI-associated CpG sites with newborn birthweight]. Table S5. [Association of maternal BMI-associated CpG sites with child OWO at age 1-18 years]. Table S6. [The annotated gene of the identified CpGs with significant mediation effect].
Fig. S1. [Distributions of DNA methylation at 8 maternal-BMI associated CpG sites in cord blood]. Fig. S2. [Number of CpG sites and single age-window child OWO with significant mediation effect].
Table S2. 481 CpG sites significantly associated with maternal pre-pregnancy body mass index or the interplay of maternal pre-pregnancy body mass index and smoking.
About this article
Cite this article
Si, J., Meir, A.Y., Hong, X. et al. Maternal pre-pregnancy BMI, offspring epigenome-wide DNA methylation, and childhood obesity: findings from the Boston Birth Cohort. BMC Med 21, 317 (2023). https://doi.org/10.1186/s12916-023-03003-5
- DNA methylation
- Childhood Obesity
- Maternal pre-pregnancy obesity
- Early life origins of chronic diseases
- Boston Birth Cohort