Common pregnancy and perinatal complications are associated with offspring cardiometabolic risk factors. These complications may influence multiple metabolic traits in the offspring and these associations might differ with offspring age.
We used data from eight population-based cohort studies to examine and compare associations of pre-eclampsia (PE), gestational hypertension (GH), gestational diabetes (GD), preterm birth (PTB), small (SGA) and large (LGA) for gestational age (vs. appropriate size for gestational age (AGA)) with up to 167 plasma/serum-based nuclear magnetic resonance-derived metabolic traits encompassing lipids, lipoproteins, fatty acids, amino acids, ketones, glycerides/phospholipids, glycolysis, fluid balance, and inflammation. Confounder-adjusted regression models were used to examine associations (adjusted for maternal education, parity age at pregnancy, ethnicity, pre/early pregnancy body mass index and smoking, and offspring sex and age at metabolic trait assessment), and results were combined using meta-analysis by five age categories representing different periods of the offspring life course: neonates (cord blood), infancy (mean ages: 1.1–1.6 years), childhood (4.2–7.5 years); adolescence (12.0–16.0 years), and adulthood (22.0–67.8 years).
Offspring numbers for each age category/analysis varied from 8925 adults (441 PTB) to 1181 infants (135 GD); 48.4% to 60.0% were females. Pregnancy complications (PE, GH, GD) were each associated with up to three metabolic traits in neonates (P≤0.001) with some evidence of persistence to older ages. PTB and SGA were associated with 32 and 12 metabolic traits in neonates respectively, which included an adjusted standardised mean difference of −0.89 standard deviation (SD) units for albumin with PTB (95% CI: −1.10 to −0.69, P=1.3×10−17) and −0.41 SD for total lipids in medium HDL with SGA (95% CI: −0.56 to −0.25, P=2.6×10−7), with some evidence of persistence to older ages. LGA was inversely associated with 19 metabolic traits including lower levels of cholesterol, lipoproteins, fatty acids, and amino acids, with associations emerging in adolescence, (e.g. −0.11 SD total fatty acids, 95% CI: −0.18 to −0.05, P=0.0009), and attenuating with older age across adulthood.
These reassuring findings suggest little evidence of wide-spread and long-term impact of common pregnancy and perinatal complications on offspring metabolic traits, with most associations only observed for newborns rather than older ages, and for perinatal rather than pregnancy complications.
There are widespread changes in maternal circulating metabolites during pregnancy, which return to normal after pregnancy . These alterations are likely to be important for maternal health, normal foetal growth, and development . Maternal metabolic profiles are associated with different common pregnancy and perinatal complications including pre-eclampsia (PE), gestational hypertension (GH), gestational diabetes (GD), preterm birth (PTB) and small (SGA) and large for gestational age (LGA) [3,4,5,6,7,8,9]. PE, GH, PTB, and SGA tend to relate to placental pathologies/foetal growth restriction [10, 11], whereas GD and LGA relate to foetal overgrowth [12, 13]. Both foetal growth restriction and overgrowth might have long-lasting metabolic effects, which may in turn increase cardiovascular disease (CVD) risk [14,15,16].
Studies indicate common pregnancy/perinatal complications associated with cardiovascular disease in the offspring [17,18,19,20,21,22]. However, to the best of our knowledge, effects of common pregnancy/perinatal complications on offspring metabolic traits, and whether these change with age, have not been examined. Identifying whether effects on offspring metabolism are short lived, persist across life, emerge later, or strengthen/weaken with age can improve our understanding of CVD aetiology and may inform the timing of interventions. Therefore, the aim of this study was to examine associations of common pregnancy and perinatal complications related to foetal growth restriction (PE, GH, PTB, and SGA) and foetal overgrowth (GD, LGA) with targeted metabolomic profiles across the offspring life course and investigate whether associations differ by offspring age at assessment of metabolic traits.
This study was carried out by following a pre-specified analysis plan and code developed by AE and DAL (https://osf.io/vfd7g) and is reported in accordance with The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement guidelines for cohort studies .
Participating cohorts were recruited from the EU Child Cohort Network (EUCCN) ; a consortium of European and Australian pregnancy/birth cohorts. Studies were included if they had data on (i) at least one pregnancy/perinatal complication, (ii) offspring metabolic profiles measured at any age in plasma/serum by high-throughput proton nuclear magnetic resonance (NMR)-based targeted metabolomics platform (the most widely used platform across EUCCN cohorts) , and (iii) prespecified confounders.
Eight cohorts were eligible, and all agreed to participate in this analysis. These were the UK-based Avon Longitudinal Study of Parents and Children (ALSPAC) [26,27,28,29,30], Born in Bradford Study (BiB) [31, 32], Finland-based Young Finns Study (YFS) [33,34,35], Northern Finland Birth Cohort 1966 (NFBC1966) [35,36,37], Northern Finland Birth Cohort 1986 (NFBC1986) , and Helsinki Birth Cohort Study (HBCS) , and the Australia-based Barwon Infant Study (BIS) [36, 37] and Longitudinal Study of Australian Children’s Child Health CheckPoint (CheckPoint) [39, 40]. A ninth cohort, the Generation R study from the Netherlands, had metabolic traits measured by mass spectrometry in cord blood and non-fasting blood samples at mean age 9.8 years and was used to replicate results for any overlapping NMR-derived traits [41,42,43]. More detail on included cohorts is in Additional file 1: Supplemental Methods.
Ethics, consent, and permissions
All cohorts had ethical approval from their relevant local or national ethics committees and study participants provided informed consent or assent to participate in the respective cohorts and secondary data analyses. Details on ethics approvals and consent for each cohort can be found in the Additional file 1: Supplemental Methods.
Pregnancy and perinatal complications
Six common pregnancy/perinatal complications related to foetal growth restriction (PE, GH, PTB, and SGA) and foetal overgrowth (GD, LGA) were included. Data on pregnancy/perinatal complications were extracted from medical records or reported in questionnaires. Description of how these were recorded in each cohort is in Additional file 1: Supplemental Methods. Data harmonisation has been previously described . Briefly, PE was defined as elevated blood pressure >20 weeks gestation (≥140 mmHg systolic or ≥90 mmHg diastolic), and proteinuria (>0.3g per 24 h), or by HELLP syndrome (haemolysis, elevated liver enzymes, and low platelets) . GH was defined as new onset hypertension after 20 weeks of gestation, with previously normal blood pressure, without proteinuria or manifestations of PE. GD was defined as glucose intolerance with onset or first diagnosis in pregnancy and continuing past 24–28 weeks of gestation. This was based on a 75-g oral glucose tolerance test comprising fasting and 2h post-load samples at around 26–28 weeks gestation in BiB, and by extraction from health records or questionnaire responses in the other cohorts. PTB was defined as ≤37 completed weeks (or ≤259 days) at birth. SGA and LGA were defined based on the World Health Organisation foetal growth charts, using the 5th and 95th percentiles as cut-offs, respectively .
Offspring exposed to PE, GH, PTB, or GD, were compared to those not exposed to the specific complication. Offspring born SGA and LGA, were each compared to those born appropriate size for gestational age (AGA) offspring, i.e. with SGA/LGA excluded in turn.
Offspring NMR-derived metabolic traits and age categories
A proton NMR-based targeted metabolomics platform  was used to quantify up to 250 offspring metabolic traits (including derived variables) in plasma/serum samples in the eight participating cohorts. The NMR platform uses a single experimental setup to simultaneously quantify metabolic traits from each plasma/serum sample. Metabolic traits were quantified in absolute concentration units or ratios and included circulating lipoprotein lipids and subclasses, fatty acids and their compositions, amino acids and traits related to glycolysis, ketone bodies, fluid balance, and an inflammatory marker. The manufacturer’s standard quality control procedures were performed in all cohorts .
Traits were analysed using non-fasting samples in cord blood and infancy, semi-fasting samples in childhood and adolescence, and fasting samples at older ages. Differences in fasting status were because younger participants were not asked to fast prior to clinic visits due to compliance/ethical issues. Metabolic trait ratios were excluded because of challenges in their interpretation, leaving up to 167 metabolic traits in the analysis. Description of the methods and ages at the assessment of metabolic traits in each cohort is in Additional file 1: Supplemental Methods. The metabolic traits available in each cohort are listed in Additional file 2: Data Set 1.
Metabolic traits measured at all available ages from each cohort were included. Cohort results were combined into five pre-specified age categories for meta-analysis, chosen to reflect key life course periods and to maximise the number of participants, with results from at least two cohorts available for each age group. The age categories were neonates (cord blood), infancy (mean age 1.1 to 1.6 years), childhood (mean age 4.2 to 7.5 years), adolescence (mean age 12.0 to 16.0 years), and adulthood (mean age 22.0 to 67.8 years).
To estimate unconfounded associations of pregnancy/perinatal complications with offspring metabolic traits, we identified and adjusted for potential confounders, i.e. factors that could plausibly cause pregnancy/perinatal complications and influence offspring metabolism, and avoided adjustment for mediators on causal path of any effect (e.g. offspring adiposity) and other sources of collider bias . The identified confounders were maternal education (the most available and consistent indicator of early life socioeconomic position across cohorts), ethnicity, age at pregnancy/birth of offspring, parity, pre/early pregnancy body mass index (BMI) and smoking in pregnancy. Offspring sex and age at metabolic trait assessment were included as adjustments to improve modelling precision. Details on how confounders were measured in each cohort are provided in Additional file 1: Supplemental Methods. Harmonised variables were derived as described previously .
ALSPAC, BiB, BIS, and YFS were able to adjust for all confounders. NFBC1966, NFBC1986, and HBCS did not adjust for ethnicity, but they were predominantly white ethnicity, HBCS was unable to adjust for smoking, and CheckPoint was unable to adjust for BMI or parity.
Associations between each of the six pregnancy/perinatal complications and each offspring NMR-derived metabolic trait from all available timepoints were examined in each cohort by fitting adjusted (for confounders plus offspring age and sex) linear regression models (with robust standard errors). Analyses were restricted to those with complete data on the relevant pregnancy/perinatal complication, metabolic trait, and confounders. The impact of missing data was explored by comparing the characteristics of included offspring with those that were excluded due to missing data (Additional file 3: Table S1). To allow comparison of results across different pregnancy/perinatal complications, traits, and ages, metabolic traits were analysed using cohort-specific standard deviation (SD) units (mean=0, SD=1).
Cohort-specific results were then combined using meta-analysis in five age categories for neonates, infancy, childhood, adolescence, and adulthood, using a random effects model to allow for between-cohort heterogeneity. Variability in the meta-analysis results that were due to between-cohort heterogeneity was measured by computing the I2 statistic [48, 49]. Where evidence of substantial between-cohort heterogeneity was found, we inspected each cohort’s results to identify the reason for heterogeneity. A P-value threshold of P≤0.001 was selected to identify statistically robust associations between the pregnancy/perinatal complications and metabolic traits. This was chosen instead of a Bonferroni-corrected P-value because many metabolic traits are highly correlated  and so independent tests were not performed, yet our threshold is still more stringent than a conventional P-value threshold. For association that reached this threshold in one age category, we highlighted the equivalent association in all other age categories, to explore changes with age in the context of different numbers of participants for each age category (Table 1).
We further investigated the change in associations with older age by fitting confounder-adjusted natural cubic spline mixed effects trajectory models  in 4980 ALSPAC offspring with up to 4 repeated NMR-based assessments from 7 to 26 years. This analysis included all offspring with complete data on the relevant pregnancy/perinatal complication and confounders, and at least one of the 4 repeated measures of metabolic traits (i.e. those with incomplete outcome measurements were included). Because ALSPAC data spans ages 7–26 years, trajectory analysis was only done for the metabolic traits showing a meta-analysis association in childhood, adolescence, or adulthood (and not in neonates or infancy). An interaction between the pregnancy/perinatal complication and age was included to allow different trajectories for exposed/nonexposed offspring. The predicted mean trajectories and differences in metabolic traits were obtained .
Lastly, we sought support for associations identified with NMR-derived traits by performing a replication analysis in the Generation R study. Replication was done for all NMR-derived traits that were available from the mass spectrometry platform, using regression models with similar adjustment for confounders.
A total of 24,864 offspring from 8 cohorts were included in this study for the age sub-group meta-analysis (Table 1). Offspring numbers in each age category analysis varied from 8925 adults (441 PTB) to 1181 infants (135 GD) (Table 2). Offspring birth years were from 1934 to 2013, and 48.4% to 60.0% were female (Table 1). The proportion of offspring exposed to pregnancy/perinatal complications in each cohort ranged from 1.9% to 3.7% PE, 1.5% to 14.2% GH, 0.5% to 15.3% GD, 4.6% to 10.4% PTB, 3.5% to 9.2% SGA, and 6.2% to 13.9% LGA (Table 1). Those excluded due to missing data had lower maternal education, younger maternal age, and higher prevalence of maternal pregnancy smoking and non-white ethnicity than those included in the analysis (Additional file 3: Table S1).
Pre-eclampsia, gestational hypertension, preterm birth and small for gestational age associations with offspring metabolic traits
PE, GH, PTB, and SGA were associated with 1, 1, 32, and 12 metabolic traits, respectively, with all but 3 of these associations observed for neonates only. PE was associated with lower levels of aromatic amino acid phenylalanine in infancy (mean difference: −0.44 SD, 95% CI: −0.67 to −0.22, P=0.0002), with differences close to zero in neonates, children, and adults, but there was some evidence in favour of lower levels in adolescents (mean difference: −0.17 SD, 95% CI: −0.34 to −0.01, P=0.03) (Fig. 1, Additional file 4: Table S2). GH was inversely associated with the ketone acetate in infants (mean difference: −0.43 SD, 95% CI: −0.59 to −0.28, P=3.7×10−8), although equivalent differences for other age groups were either close to zero or were imprecisely estimated (Fig. 1, Additional file 4: Table S2).
PTB was inversely associated with total lipids in small HDL, concentration of small HDL particles, degree of unsaturation, glucose, fluid balance markers creatinine and albumin, and inflammatory marker glycoprotein acetyls (GlycA), and positively associated with cholesterol measures, lipoprotein subclasses, glycerides, phospholipids, Apolipoprotein B, saturated fatty acids, the α-amino acid glutamine, and the aromatic amino acid tyrosine, with all associations observed in neonates (Fig. 2). Associations with neonate metabolic traits for PTB (vs. not PTB) ranged in magnitude from 0.22 SD (95% CI: 0.09 to 0.35, P=0. 0007) for total cholines to −0.89 SD (95% CI: −1.10 to −0.69, P=1.3×10−17) for albumin. For most results, equivalent differences across older age categories were close to zero though there was some evidence of persistence to older ages for some metabolic traits, e.g. the mean difference in total lipids in small LDL in neonates and adults was 0.33 SD (95% CI: 0.16 to 0.50, P=0.0001), and 0.14 SD (95% CI: 0.04 to 0.24, P=0.004), respectively (Fig. 2, Additional file 4: Table S2).
SGA (vs. AGA) in neonates was inversely associated with total cholesterol in HDL, total cholesterol in HDL2, total cholesterol in HDL3, total lipids in medium HDL, concentration of medium HDL particles, apolipoprotein A-I, and histidine, and positively associated with total cholesterol in VLDL, total lipids in very small VLDL, concentration of very small VLDL particles, and omega-3 fatty acids (Fig. 2). Differences ranged in magnitude from −0.21 SD (95% CI: −0.33 to −0.08, P=0.001) for histidine to −0.41 SD (95% CI: −0.56 to −0.25, P=2.6×10−7) for total lipids in medium HDL (Additional file 4: Table S2). Most were reduced at older ages but there was some evidence for higher levels of total cholesterol in VLDL, total lipids in very small VLDL, concentration of very small VLDL particles, and omega-3 fatty acids during adolescence, e.g. mean differences in total lipids in very small VLDL in neonates and adolescence were 0.34 SD (95% CI: 0.20 to 0.47, P=1.4×10−6), and 0.17 SD (95% CI: 0.06 to 0.29, P=0.003), respectively. SGA was also inversely associated with amino acid alanine in childhood (−0.25 SD, 95% CI: −0.38 to −0.11, P=0.0003), with no clear differences in alanine at other age groups (Fig. 2).
Of the NMR-derived metabolic traits that PE, GH, PTB, and SGA were associated with, four were found among mass spectroscopy measures in the Generation R Study and included for replication (PTB: tyrosine, sphingomyelins, glutamine; SGA: histidine). Consistent with the pooled difference for PTB neonates in NMR-derived tyrosine (0.34 SD, 95% CI: 0.19 to 0.49, P=7.0×10−6) and sphingomyelins (0.39 SD, 95% CI: 0.23 to 0.54, P=6.8×10−7), PTB was also associated with higher tyrosine (0.82 SD, 95% CI: 0.40 to 1.24, P=0.0001, n=725 (29 PTB)) and sphingomyelins (0.49 SD, 95% CI: 0.07 to 0.90, P=0.02) in neonates in Generation R. In contrast, the meta-analysis association of PTB with glutamine (0.56 SD, 95% CI: 0.36 to 0.76, P=2.4×10−8) did not replicate in Generation R (−0.09, 95% CI: −0.51 to 0.34, P=0.7). Lastly, the mean difference in histidine for SGA (vs. AGA) neonates was similar in the meta-analysis (−0.21 SD, 95% CI: −0.33 to −0.08, P=0.001) and Generation R but this result was imprecisely estimated (−0.25 SD, 95% CI: −0.61 to 0.11, P=0.2, n=651 (35 SGA)).
Gestational diabetes and large for gestational age associations with offspring metabolic traits
GD and LGA were associated with 3 and 19 metabolic traits, respectively. GD was associated with smaller LDL particle size (mean difference: −0.25 SD, 95% CI: −0.39 to −0.10, P=0.0007) and with lower isoleucine (mean difference: −0.27 SD, 95% CI: −0.41 to −0.14, P=0.00008) in neonates, with differences in both metabolic traits close to zero for older ages (Fig. 3). GD was positively associated with glucose in infants (mean difference: 0.35 SD, 95% CI: 0.18 to 0.52, P=0.00005), with no difference in glucose found for other age categories (Fig. 3).
Of the nineteen associations of LGA with offspring metabolic traits, none were observed in neonates or infants, two were observed in children, sixteen in adolescents, and one in adults (Fig. 3). All were inverse associations and represented lower levels of cholesterol, fatty acids, lipoprotein subclasses, three branched-chain amino acids, and the inflammatory marker glycoprotein acetyls for LGA (vs. AGA). Associations ranged in magnitude from −0.10 SD (95% CI: −0.16 to −0.04, P=0.001) for leucine in adolescents to −0.19 SD (95% CI: −0.29 to −0.09, P=0.0003) for valine in adulthood (Additional file 4: Table S2). For associations seen in adolescence, equivalent associations in adults were slightly attenuated and had wider confidence intervals, and for some traits, there was also evidence for lower levels in childhood, e.g. difference in total lipids in small HDL in childhood, adolescence and adulthood was −0.14 SD (95% CI: −0.24 to −0.04, P=0.004), −0.13 SD (95% CI: −0.19 to −0.07, P=.00006), and −0.08 SD (95% CI: −0.16 to 0.01, P=0.1), respectively.
Of the 22 associations identified for GD and LGA, only one metabolic trait overlapped with mass spectroscopy measures in the Generation R Study (LGA and child total triglycerides). The inverse association in our meta-analysis (difference in child total triglycerides for LGA vs. AGA (−0.15 SD, 95% CI: −0.24 to −0.06, P=0.001) was weaker and imprecisely estimated in Generation R (−0.04 SD, 95% CI: −0.40 to 0.32, P=0.8, n=339 total with n=75 LGA),
Overlap in metabolic trait associations
Comparing identified pregnancy/perinatal complication–metabolic trait associations showed that PE/GH were each associated with a different metabolic trait, PTB/SGA were associated with total HDL cholesterol, total HDL2 cholesterol, and total HDL3 cholesterol in neonates (higher with PTB and lower with SGA), whereas GD (neonate) and LGA (adolescent) were inversely associated with isoleucine. Additionally, SGA (neonates) and LGA (adolescents) were both inversely associated with concentration of medium HDL particles, PTB (neonate) and LGA (adolescence) were both inversely associated with GlycA and total lipids in small HDL, and both PTB (inverse association: neonates) and GD (positive association: infants) were associated with glucose.
For most (72%) of the total 3787 results there was little to no between-cohort heterogeneity (I2≤25%, with I2=0% for 86% of these), and 8% of the results showed evidence of substantial or high heterogeneity between cohorts (I2≥75%) (Additional file 5: Data Set 2). Four of these were results that met the P≤0.001 threshold and all were for PTB neonates (total lipids in small and very large HDL, and concentrations of small and very large HDL particles). Inspecting results from the cohorts contributing to this age (BiB and BIS) revealed a consistent direction of association in both but mean differences in BiB double those in BIS (e.g. mean difference in total lipids in small HDL was −0.86 SD in BiB and −0.48 SD in BIS. Among other results with substantial or high heterogeneity, 97 were for LGA–adult metabolic traits. Further investigation revealed that this was because LGA was inversely associated with adult biomarkers in ALSPAC (age 24.5 years), with smaller positive associations seen in NFBC1966 (age 46.6 years). Notably, results for NFBC1966 at age 31.2 years were between the ALSPAC age 24.5 years results and the NFBC1966 age 46.6 years results, suggesting a possible age effect (Fig. 4).
ALSPAC age-change trajectory analysis
Of the twenty associations identified with metabolic traits beyond infancy, one was for SGA, nineteen were for LGA, and all except one for LGA were included in trajectory analysis from childhood to adulthood in the ALSPAC cohort. Consistent with the meta-analysis, SGA (vs. AGA) was associated with lower alanine at age 7 years, and this difference was reduced with increasing child age, with alanine slightly higher in SGA from mid adolescence (Additional file 6: Figure S1). Similarly, most associations between LGA and metabolic traits appeared to change with age from childhood to adulthood e.g. the inverse association in the meta-analysis between LGA (vs. AGA) and triglycerides in VLDL in adolescence was attenuated with older age (Fig. 5, Additional file 7: Figure S2).
We examined the association between common pregnancy/perinatal complications and NMR-derived metabolic profile across the offspring life course in eight population-based cohort studies. Pregnancy complications (PE, GH, GD) were related to only a few metabolic traits mostly in neonates, with little evidence of persistence. More extensive disruptions in neonate metabolic traits were seen for perinatal complications, mostly for PTB, with little evidence of persistence to older ages, except for differences in metabolic traits with LGA, which mostly emerged during adolescence and appeared to weaken by older adulthood.
To the best of our knowledge, ours is the first study to simultaneously investigate short- and long-term effects of common pregnancy/perinatal complications on offspring metabolomics. We found that most pregnancy/perinatal complications were associated with metabolic traits in neonates only with little evidence of persistence of associations beyond early life. This is somewhat in contradiction to commonly expressed views on the importance of the pregnancy period for offspring long-term health. However, our results are consistent with Mendelian randomisation studies that show that intrauterine exposures related to variation in birth weight are unlikely to causally influence offspring metabolic health [53, 54]. Most of the associations with metabolic traits in our study were observed in preterm compared with term/post-term born neonates, which appears to support earlier findings of widespread gestational age-dependent effects on metabolites from an untargeted metabolomics analysis in 298 neonates .
Our finding of inverse associations of LGA with adolescent metabolic traits is consistent with a study of 18,288 adolescents and adults with the same metabolomics platform as our study  which showed a higher mean birth weight across the distribution was associated with healthier lipid profile, including lower triglycerides, and other lipids that we see inverse associations of LGA within our study. These findings may reflect an interplay with offspring adiposity and puberty whereby LGA offspring have higher prepubertal body fat compared to AGA  with this leading to earlier puberty for LGA [57, 58]. Earlier puberty in turn might influence metabolic traits  and explain these findings. For example, a decreasing body fat with advanced puberty  in the earlier maturing LGA adolescents might explain their healthier metabolic profile at this age . This could also explain why the association of LGA with healthier metabolic profile in adolescence was weakened or even reversed by early midlife (because all other offspring would have already completed puberty). Shared lifestyles between the mothers and offspring might also contribute to these findings [62, 63].
Strengths of this work include the larger sample size compared with previous studies and the examination of metabolic traits during important life periods including whether these change over the life course. Limitations include the small numbers of exposed offspring despite the relatively larger sample size, (e.g. 18.4% exposed to GH — the most common complication). This also meant that we had little power (and so were unable) to examine sex interactions or explore associations separately by sex. Residual confounding due to use of crude harmonised variables across cohorts [64, 65] and unmeasured confounders (e.g. maternal health) could influence our findings. Complete case analysis may have reduced precision of estimates and introduced bias due to missing data and limited the generalisability of our findings due to the difference in characteristics between included and excluded participants. Bias due to missing data on metabolic traits is reduced in the ALSPAC trajectory analysis because all offspring with incomplete outcome measures were included, under the missing at random assumption (i.e. the probability of a missing outcome measurement depends on the observed values of the outcome, conditional on the covariates in the model).
Attrition could introduce a selection bias which may contribute to null/weaker associations in adults and reduced sensitivity to detect associations with small effect sizes, if older offspring were healthier (with a lower prevalence of complications) than those lost to follow-up. The P-value threshold used to identify robust associations was somewhat arbitrary, so results require replication in other studies. Approaches to select an effective number of tests have been developed  but to the best of our knowledge have not been extended to multi-cohort analyses, like our study. We have provided all coefficients and exact P-values, which readers can use to apply their own P-value thresholds (Additional file 5: Data Set 2). Only one of five cohorts included in the GD analysis had universal diagnostic testing at the time of pregnancy therefore, some misclassification of GD is possible. Direct measurements of neonatal body fat were not examined and may have provided greater insight than birthweight-derived exposures (SGA and LGA). Finally, there was a wide age gap for adults, especially for preterm birth where 67-year-olds were included, which might limit the interpretation of results in adults.
Our results offer reassurance that PE, GH, and GD do not result in widespread metabolic disruption in offspring and that more widespread disruptions for PTB and SGA are mostly confined to neonates. Differences in offspring metabolic traits for LGA require further exploration to establish why they primarily arose during adolescence.
Availability of data and materials
The datasets used during the current study are available to researchers by request from each cohort.
Liang L, Rasmussen M-LH, Piening B, et al. Metabolic Dynamics and Prediction of Gestational Age and Time to Delivery in Pregnant Women. Cell. 2020;181(7):1680–92.e15. https://doi.org/10.1016/j.cell.2020.05.002.
Kelly RS, Giorgio RT, Chawes BL, et al. Applications of Metabolomics in the Study and Management of Preeclampsia; A Review of the Literature. Metabolomics. 2017;13(7):86. https://doi.org/10.1007/s11306-017-1225-8.
McBride N, Yousefi P, Sovio U, et al. Do Mass Spectrometry-Derived Metabolomics Improve the Prediction of Pregnancy-Related Disorders? Findings from a UK Birth Cohort with Independent Validation. Metabolites. 2021;11(8). https://doi.org/10.3390/metabo11080530.
McBride N, Yousefi P, White SL, et al. Do nuclear magnetic resonance (NMR)-based metabolomics improve the prediction of pregnancy-related disorders? Findings from a UK birth cohort with independent validation. BMC Med. 2020;18(1):366. https://doi.org/10.1186/s12916-020-01819-z.
Taylor K, Ferreira DLS, West J, Yang T, Caputo M, Lawlor DA. Differences in Pregnancy Metabolic Profiles and Their Determinants between White European and South Asian Women: Findings from the Born in Bradford Cohort. Metabolites. 2019;9(9):190. https://doi.org/10.3390/metabo9090190.
Parikh NI, Gonzalez JM, Anderson CAM, et al. Adverse Pregnancy Outcomes and Cardiovascular Disease Risk: Unique Opportunities for Cardiovascular Disease Prevention in Women: A Scientific Statement From the American Heart Association. Circulation. 2021;143(18):e902–e16. https://doi.org/10.1161/CIR.0000000000000961.
Brosens I, Pijnenborg R, Vercruysse L, Romero R. The “Great Obstetrical Syndromes” are associated with disorders of deep placentation. Am J Obstet Gynecol. 2011;204(3):193–201. https://doi.org/10.1016/j.ajog.2010.08.009.
Kaitu'u-Lino TJ, MacDonald TM, Cannon P, et al. Circulating SPINT1 is a biomarker of pregnancies with poor placental function and fetal growth restriction. Nat Commun. 2020;11(1):2411. https://doi.org/10.1038/s41467-020-16346-x.
Brand JS, West J, Tuffnell D, et al. Gestational diabetes and ultrasound-assessed fetal growth in South Asian and White European women: findings from a prospective pregnancy cohort. BMC Med. 2018;16(1):203. https://doi.org/10.1186/s12916-018-1191-7.
Li M, Hinkle SN, Grantz KL, et al. Glycaemic status during pregnancy and longitudinal measures of fetal growth in a multi-racial US population: a prospective cohort study. Lancet Diabetes Endocrinol. 2020;8(4):292–300. https://doi.org/10.1016/s2213-8587(20)30024-3.
Bavineni M, Wassenaar TM, Agnihotri K, Ussery DW, Lüscher TF, Mehta JL. Mechanisms linking preterm birth to onset of cardiovascular disease later in adulthood. Eur Heart J. 2019;40(14):1107–12. https://doi.org/10.1093/eurheartj/ehz025.
Alsnes IV, Vatten LJ, Fraser A, et al. Hypertension in Pregnancy and Offspring Cardiovascular Risk in Young Adulthood: Prospective and Sibling Studies in the HUNT Study (Nord-Trøndelag Health Study) in Norway. Hypertension. 2017;69(4):591–8. https://doi.org/10.1161/hypertensionaha.116.08414.
Yu Y, Arah OA, Liew Z, et al. Maternal diabetes during pregnancy and early onset of cardiovascular disease in offspring: population based cohort study with 40 years of follow-up. BMJ. 2019;367:l6398. https://doi.org/10.1136/bmj.l6398.
Jaddoe VWV, de Jonge LL, Hofman A, Franco OH, Steegers EAP, Gaillard R. First trimester fetal growth restriction and cardiovascular risk factors in school age children: population based cohort study. BMJ. 2014;348:g14. https://doi.org/10.1136/bmj.g14.
Kurbasic A, Fraser A, Mogren I, et al. Maternal Hypertensive Disorders of Pregnancy and Offspring Risk of Hypertension: A Population-Based Cohort and Sibling Study. Am J Hypertens. 2019;32(4):331–4. https://doi.org/10.1093/ajh/hpy176.
Geelhoed JJ, Fraser A, Tilling K, et al. Preeclampsia and gestational hypertension are associated with childhood blood pressure independently of family adiposity measures: the Avon Longitudinal Study of Parents and Children. Circulation. 2010;122(12):1192–9. https://doi.org/10.1161/circulationaha.110.936674.
von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet. 2007;370(9596):1453–7. https://doi.org/10.1016/s0140-6736(07)61602-x.
Jaddoe VWV, Felix JF, Andersen AN, et al. The LifeCycle Project-EU Child Cohort Network: a federated analysis infrastructure and harmonized data of more than 250,000 children and parents. Eur J Epidemiol. 2020;35(7):709–24. https://doi.org/10.1007/s10654-020-00662-z.
Boyd A, Golding J, Macleod J, et al. Cohort Profile: the ‘children of the 90s’--the index offspring of the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2013;42(1):111–27. https://doi.org/10.1093/ije/dys064.
Fraser A, Macdonald-Wallis C, Tilling K, et al. Cohort Profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int J Epidemiol. 2013;42(1):97–110. https://doi.org/10.1093/ije/dys066.
Northstone K, Lewcock M, Groom A, et al. The Avon Longitudinal Study of Parents and Children (ALSPAC): an update on the enrolled sample of index children in 2019. Wellcome Open Res. 2019;4:51. https://doi.org/10.12688/wellcomeopenres.15132.1.
Bell JA, Bull CJ, Gunter MJ, et al. Early Metabolic Features of Genetic Liability to Type 2 Diabetes: Cohort Study With Repeated Metabolomics Across Early Life. Diabetes Care. 2020:dc192348. https://doi.org/10.2337/dc19-2348.
Wright J, Small N, Raynor P, et al. Cohort Profile: the Born in Bradford multi-ethnic family cohort study. Int J Epidemiol. 2013;42(4):978–91. https://doi.org/10.1093/ije/dys112.
Kaikkonen JE, Würtz P, Suomela E, et al. Metabolic profiling of fatty liver in young and middle-aged adults: Cross-sectional and prospective analyses of the Young Finns Study. Hepatology. 2017;65(2):491–500. https://doi.org/10.1002/hep.28899.
Nordström T, Miettunen J, Auvinen J, et al. Cohort Profile: 46 years of follow-up of the Northern Finland Birth Cohort 1966 (NFBC1966). Int J Epidemiol. 2021. https://doi.org/10.1093/ije/dyab109.
Santos Ferreira DL, Williams DM, Kangas AJ, et al. Association of pre-pregnancy body mass index with offspring metabolic profile: Analyses of 3 European prospective birth cohorts. PLoS Med. 2017;14(8):e1002376. https://doi.org/10.1371/journal.pmed.1002376.
Eriksson JG, Sandboge S, Salonen MK, Kajantie E, Osmond C. Long-term consequences of maternal overweight in pregnancy on offspring later health: findings from the Helsinki Birth Cohort Study. Ann Med. 2014;46(6):434–8. https://doi.org/10.3109/07853890.2014.919728.
Clifford SA, Davies S, Wake M. Child Health CheckPoint: cohort summary and methodology of a physical health and biospecimen module for the Longitudinal Study of Australian Children. BMJ Open. 2019;9(Suppl 3):3–22. https://doi.org/10.1136/bmjopen-2017-020261.
Ellul S, Wake M, Clifford SA, et al. Metabolomics: population epidemiology and concordance in Australian children aged 11-12 years and their parents. BMJ Open. 2019;9(Suppl 3):106–17. https://doi.org/10.1136/bmjopen-2017-020900.
Voerman E, Jaddoe VWV, Uhl O, et al. A population-based resource for intergenerational metabolomics analyses in pregnant women and their children: the Generation R Study. Metabolomics. 2020;16(4):43. https://doi.org/10.1007/s11306-020-01667-1.
Pinot de Moira A, Haakma S, Strandberg-Larsen K, et al. The EU Child Cohort Network’s core data: establishing a set of findable, accessible, interoperable and re-usable (FAIR) variables. Eur J Epidemiol. 2021. https://doi.org/10.1007/s10654-021-00733-9.
Kiserud T, Piaggio G, Carroli G, et al. The World Health Organization Fetal Growth Charts: A Multinational Longitudinal Study of Ultrasound Biometric Measurements and Estimated Fetal Weight. PLoS Med. 2017;14(1):e1002220. https://doi.org/10.1371/journal.pmed.1002220.
Elhakeem A, Hughes RA, Tilling KM, et al. Using linear and natural cubic splines, SITAR, and latent trajectory models to characterise nonlinear longitudinal growth trajectories in cohort studies. BMC Med Res Methodol. 2022;22(1):68. https://doi.org/10.1186/s12874-022-01542-8
Moen GH, Brumpton B, Willer C, et al. Mendelian randomization study of maternal influences on birthweight and future cardiometabolic risk in the HUNT cohort. Nat Commun. 2020;11(1):5404. https://doi.org/10.1038/s41467-020-19257-z.
Bond TA, Richmond RC, Karhunen V, et al. Exploring the causal effect of maternal pregnancy adiposity on offspring adiposity: Mendelian randomization using polygenic risk scores. BMC Med. 2022;20(1):34.
Taal HR, Vd Heijden AJ, Steegers EA, Hofman A, Jaddoe VW. Small and large size for gestational age at birth, infant growth, and childhood overweight. Obesity (Silver Spring). 2013;21(6):1261–8. https://doi.org/10.1002/oby.20116.
Di Giovanni I, Marcovecchio ML, Chiavaroli V, de Giorgis T, Chiarelli F, Mohn A. Being born large for gestational age is associated with earlier pubertal take-off and longer growth duration: a longitudinal study. Acta Paediatr. 2017;106(1):61–6. https://doi.org/10.1111/apa.13633.
Perak AM, Lancki N, Kuang A, et al. Associations of Maternal Cardiovascular Health in Pregnancy With Offspring Cardiovascular Health in Early Adolescence. JAMA. 2021;325(7):658–68. https://doi.org/10.1001/jama.2021.0247.
Dhana K, Haines J, Liu G, et al. Association between maternal adherence to healthy lifestyle practices and risk of obesity in offspring: results from two prospective cohort studies of mother-child pairs in the United States. BMJ. 2018;362:k2486. https://doi.org/10.1136/bmj.k2486.
Brand JS, Gaillard R, West J, et al. Associations of maternal quitting, reducing, and continuing smoking during pregnancy with longitudinal fetal growth: Findings from Mendelian randomization and parental negative control studies. PLoS Med. 2019;16(11):e1002972. https://doi.org/10.1371/journal.pmed.1002972.
Macdonald-Wallis C, Tilling K, Fraser A, Nelson SM, Lawlor DA. Established preeclampsia risk factors are related to patterns of blood pressure change in normal term pregnancy: findings from the Avon Longitudinal Study of Parents and Children. J Hypertens. 2011;29(9):1703–11. https://doi.org/10.1097/HJH.0b013e328349eec6.
We thank all cohort members and researchers who participated in the study. Cohort-specific acknowledgements can be found in Table S3.
This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreements No. 733206 (LifeCycle), No. 874739 (LongITools) and No. 101021566 (ART-HEALTH), and the British Heart Foundation (CH/F/20/90003 and AA/18/1/34219). AE and DAL work in a unit that is supported by the University of Bristol and UK Medical Research Council (MC_UU_00011/6). Cohort funders are listed in Additional file 8: Table S3. The funders had no role in the design and conduct of the study; management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. This research reflects only the authors’ view, and the European Commission is not responsible for any use that may be made of the information it contains.
Authors and Affiliations
MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, UK
Ahmed Elhakeem, Tim Cadman & Deborah A. Lawlor
Population Health Science, Bristol Medical School, University of Bristol, Bristol, UK
Ahmed Elhakeem & Deborah A. Lawlor
Research Unit of Population Health, Faculty of Medicine, University of Oulu, Oulu, Finland
AE and DAL designed the study and developed the analysis plan with input from all authors. AE undertook the analysis in ALSPAC, BiB and undertook the meta-analysis. JR undertook analysis in the NFBC cohorts. TM1 undertook analysis in the BIS cohort. KL undertook analysis in the CheckPoint cohort. TM2 undertook analysis in NFBC cohort. BHM undertook analysis in YFS. RW undertook analysis in Generation R. AE wrote the first draft of the manuscript. JR TM1, KL, TM2, BHM, RW, TC, TY, DB, JGE, MRJ, RG, VWVJ, TL, OTR, RS, MW, JW, SS, and DAL provided feedback on the draft. All authors read and approved the final manuscript.
All cohorts had ethical approval from their relevant local or national ethics committees [ALSPAC: ALSPAC Law and Ethics committee and local research ethics committees (NHS Haydock REC: 10/H1010/70); BiB: Bradford Research Ethics Committee (07/H1302/112), YFS: Hospital District of Southwest Finland (ETMK:68/1801/2017), NFBC66 and NFBC86: Ethics Committee of Northern Ostrobothnia Hospital District, Finland (EETTMK: 108/2017); HBCS: Ethics Committee of the Hospital District of Helsinki and Uusimaa (023/98, 344/E3/2000, and HUS/2020/2016) and that of the National Public Health Institute, Helsinki (May 19, 1995); BIS: Barwon Health Human Research Ethics Committee (10/24); CheckPoint: Royal Children's Hospital Melbourne Human Research Ethics Committee (33225D) and Australian Institute of Family Studies Ethics Committee (14-26); Generation R: Medical Ethical Committee of the Erasmus MC, University Medical Center in Rotterdam (MEC-2012-165-NL40020.078.12)]. All study participants provided informed consent or assent to participate in the respective cohorts and secondary data analyses. Cohort descriptions and details on ethics approvals and consent for each cohort can be found in the Additional file 1: Supplemental Methods.
Consent for publication
DAL reported grants from national and international government and charity funders, Roche Diagnostics, and Medtronic Ltd for work unrelated to this publication. DAL also declares that she is an editor for BMC Medicine. The other authors report no conflicts.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Cohort-specific acknowledgements and funding statements.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Elhakeem, A., Ronkainen, J., Mansell, T. et al. Effect of common pregnancy and perinatal complications on offspring metabolic traits across the life course: a multi-cohort study.
BMC Med21, 23 (2023). https://doi.org/10.1186/s12916-022-02711-8