Introduction
Gestational diabetes mellitus (GDM) is described by elevated glucose levels identified for the first time during pregnancy.Reference Casagrande, Linder and Cowie1,Reference Landon and Gabbe2 GDM affects approximately 7.6% of pregnancies in the United States,Reference Casagrande, Linder and Cowie1 and the prevalence has been rising in recent decades.Reference Chen, Wang and Ji3–Reference Getahun, Nath, Ananth, Chavez and Smulian5 GDM is associated with complications at birth including macrosomiaReference He, Qin, Hu, Zhu, Tian and Li6,Reference Rice and Landon7 and long-term consequences in exposed offspring such as type 2 diabetesReference Damm, Houshmand-Oeregaard, Kelstrup, Lauenborg, Mathiesen and Clausen8 and greater adiposity later in life.Reference Wright, Rifas-Shiman, Rich-Edwards, Taveras, Gillman and Oken9,Reference Regnault, Gillman, Rifas-Shiman, Eggleston and Oken10
The developmental origins of health and disease hypothesisReference O'Donnell and Meaney11 suggests that long-term adverse outcomes may result from fetal programming due to exposures during in utero development. Prior studies found that hypothalamic–pituitary–adrenal (HPA) axis activity during childhood may be programmed during fetal development.Reference Reynolds12–Reference Davis, Waffarn and Sandman21 Several animal studies also have implicated maternal hyperglycemia in stunted hippocampal neurodevelopment,Reference Hami, Sadr-Nabavi, Sankian, Balali-Mood and Haghir22,Reference Piazza, Segabinazi and de Meireles23 which could potentially dysregulate the inhibitory role of the hippocampus in HPA axis regulation.Reference Jacobson and Sapolsky24 However, few human studies have assessed offspring cortisol levels in the context of exposure to GDM or maternal hyperglycemia.Reference Jack-Roberts, Maples and Kalkan25–Reference Mina, Lahti and Drake31 Furthermore, results of prior studies have been conflicting and have primarily been conducted in small cohorts.Reference Jack-Roberts, Maples and Kalkan25–Reference Mina, Lahti and Drake31 Additionally, some of them have not accounted for potential confounding, such as maternal body mass index (BMI) and socioeconomic status.Reference Jack-Roberts, Maples and Kalkan25–Reference Ma, Liu and Wu28,Reference Mina, Lahti and Drake31
The aim of this study was to examine the associations of maternal glucose markers in pregnancy with offspring cord blood glucocorticoids at birth and hair cortisol levels in mid-childhood and early adolescence using data from two longitudinal pre-birth cohorts. We hypothesized that high levels of maternal glycemia in pregnancy would be associated with higher levels of cord blood glucocorticoids and offspring hair cortisol levels.
Methods
Project Viva and Gen3G cohorts
Project Viva is an ongoing prospective prebirth cohort study of pregnant women and their children.Reference Oken, Baccarelli and Gold32 We recruited participants from obstetric clinics of Atrius Harvard Vanguard Medical Associates in eastern Massachusetts at their initial prenatal visits between 1999 and 2002. Inclusion criteria included singleton gestation, the ability to complete questionnaires written in English, plans to reside in eastern Massachusetts following delivery, and being at ≤22 weeks of gestation at enrollment. All study participants provided written informed consent at recruitment and postnatal follow-up visits, and children and teens provided verbal assent. The institutional review board at Harvard Pilgrim Health Care approved the study.Reference Oken, Baccarelli and Gold32
Of 2,128 live births included in Project Viva, we excluded 16 women with preexisting type 1 or type 2 diabetes and 45 with no data about pregnancy glucose status. We also excluded 1,139 mother–child pairs with missing values on all the cortisol outcomes. The 1,139 mother–child pairs included 251 nonwhite participants who had hair cortisol data available but no data available for any of the cord blood outcomes. Because of differences in hair growth rates and hair textures between individuals of different races that affect the cortisol assay,Reference Loussouarn33 hair cortisol analyses only included white children. In total, we included 928 pairs in the Project Viva study who were eligible for our current analyses (Supplementary Table S1).
The Genetics of Glucose Regulation in Gestation and Growth (Gen3G) cohortReference Guillemette, Allard and Lacroix34 is a prospective prebirth cohort of pregnant women and their children that began recruitment in 2010. We recruited Gen3G participants in Sherbrooke, Canada during a first trimester prenatal visit. Exclusion criteria included multiple pregnancies, known pregestational diabetes (type 1 or 2), diabetes diagnosed at first trimester, drug and/or alcohol abuse, and any other major medical conditions that could affect glucose regulation. All participants provided written informed consent prior to enrollment, and the institutional review board at Centre Hospitalier Universitaire de Sherbrooke approved the study. The Gen3G cohort is mostly of European descent and did not include enough nonwhite children with hair cortisol measurements available to be included in the present study. Therefore, analyses of hair cortisol in the Gen3G cohort only included white children. Gen3G analyses additionally excluded 39 children with undetectable hair cortisol values and those with missing second trimester oral glucose tolerance test (OGTT) results. In total, we included 313 children from the Gen3G cohort in the current analyses.
Diagnosis of GDM
In Project Viva, participants completed a two-step diagnostic procedure as part of routine clinical practice, in which GDM was diagnosed based on the Carpenter-Coustan criteria.Reference Carpenter and Coustan35 The diagnostic procedure consisted of a 50-g non-fasting glucose challenge test (GCT) followed by a 100-g, 3-h OGTT for those who failed the GCT (1-h glucose > 140 mg/dL). We deemed women who passed the initial GCT to have normal glucose tolerance (NGT). We characterized women who had normal glucose measurements at all four time points during the OGTT after previously failing the GCT as having isolated hyperglycemia (IH). We categorized women with one abnormal glucose measurement during the OGTT as having impaired glucose tolerance (IGT), and women with two or more abnormal measurements with GDM.
In the Gen3G study, participants underwent a fasting, 2-h, 75-g OGTT between the 24th and 28th week of gestation to diagnose GDM based on the diagnostic criteria recommended by the International Association of the Diabetes and Pregnancy Study Groups (IADPSG).Reference Metzger, Gabbe and Persson36 In this study, we categorized women into NGT and GDM groups based on these criteria.
Cord blood glucocorticoids
Cord blood glucocorticoid levels were measured only in the Project Viva study.Reference Huh, Andrew, Rich-Edwards, Kleinman, Seckl and Gillman37 A research staff member collected a sample of umbilical cord blood from each participant at delivery. We used radioimmunoassay techniques to measure cortisol and cortisone in the cord blood in duplicates. We calculated the ratio of cortisol to cortisone in the cord blood using the mean of duplicate cortisol and cortisone assays.
Hair cortisol
We measured child hair cortisol in both studies: in Project Viva at mid-childhood and early adolescenceReference Petimar, Rifas-Shiman, Hivert, Fleisch, Tiemeier and Oken38,Reference Petimar, Rifas-Shiman, Hivert, Fleisch, Tiemeier and Oken39 and in Gen3G at approximately five years of age. In Project Viva, research staff collected a sample of hair measuring 3 cm from the posterior vertex region of the scalp (as close to the scalp as possible). In the laboratory, lab personnel measured cortisol using liquid chromatography tandem mass spectrometry.Reference Gao, Stalder, Foley, Rauh, Deng and Kirschbaum40 We did not assess the intra-assay coefficient of variation because the study did not collect duplicate hair samples, but previous research has shown that it tends to be low when this method is used.Reference Gao, Stalder, Foley, Rauh, Deng and Kirschbaum40 We assayed the proximal 3 cm so that the cortisol measurements reflected cumulative stress over a period of approximately 3 months based on the assumption that hair grows at approximately 1 cm per month.Reference Russell, Koren, Rieder and Van Uum41 Additional details regarding the measurement of hair cortisol in Project Viva are described elsewhere.Reference Petimar, Rifas-Shiman, Hivert, Fleisch, Tiemeier and Oken39
In the Gen3G study, research staff collected hair samples from the posterior vertex region of the scalp as close to the scalp as possible. For the same reason as above, hair samples were cut to be 3 cm long. Details regarding the measurement of hair cortisol in the Gen3G study are provided in the Supplementary Methods section. In the analyses of the Gen3G data, we imputed hair cortisol values below the detectable limit of 1.0 pg/mg to be 0.5 pg/mg. We used this approach to impute these values because it is a commonly used, straightforward approach.
Other variables
In Project Viva, research staff collected demographic data using questionnaires and interviews during the study visits in pregnancy and at delivery, including maternal age, education, smoking status, household income, and race/ethnicity. We categorized child race/ethnicity based on reported data by the mothers in early childhood. We imputed missing child race/ethnicity values with maternal race/ethnicity. We obtained data on child sex from medical records. We categorized maternal education as a dichotomous variable based on whether the woman graduated college. We dichotomized household income using a cut point of $70,000, which was determined based on the distribution. We calculated maternal prepregnancy BMI using values of height and weight that were reported at study enrollment. We categorized smoking status based on whether the woman never smoked, formerly smoked, or smoked during pregnancy.
In Gen3G, research staff collected baseline demographic data during pregnancy including maternal age, type of employment, smoking status, height, and weight. We obtained data on child sex from medical records. At the 5-year follow-up visit, Gen3G research staff collected maternal education and categorized it into three categories: high school or less, some college or university education, and vocational. We categorized maternal employment into four groups: (1) management, business, or finance, (2) sciences or healthcare, (3) education, law, and social, community, and government, and (4) other. We grouped maternal employment into these categories because of their similarities. We calculated maternal BMI based on height and weight measured at the first visit in pregnancy. We dichotomized maternal smoking status based on whether the woman smoked in pregnancy, based on self-report at the first pregnancy visit.
Statistical Analyses
Statistical modeling
In all analyses, we used linear regression to assess the associations of maternal prenatal glycemia with offspring cortisol biomarkers. We performed all analyses involving categorical exposures using women with NGT as the reference group. To address skewed distribution of the hair cortisol outcomes, we used a logarithmic transformation. We selected covariates for inclusion in these models based on prior research.Reference Petimar, Rifas-Shiman, Hivert, Fleisch, Tiemeier and Oken39,Reference Cohen, Schwartz, Epel, Kirschbaum, Sidney and Seeman42–Reference Rosmalen, Oldehinkel, Ormel, de Winter, Buitelaar and Verhulst47 We assessed effect modification by child sex in the models with categorical exposures with the inclusion of an interaction term and via stratified analyses.
For the cord blood outcomes (in Project Viva), minimally adjusted models were adjusted for child race/ethnicity, sex, and gestational age at birth. Fully adjusted analyses were additionally adjusted for maternal prepregnancy BMI, education, and household income. We also assessed effect modification by mode of delivery in the cord blood analyses using tests for interaction and via stratified analyses.
For the hair cortisol outcomes (in both cohorts), minimally adjusted models were corrected for child age at hair samples collection and sex. Fully adjusted models were additionally corrected for maternal prepregnancy BMI, education, and socioeconomic status. For the analyses of the change in hair cortisol from mid-childhood to early adolescence in Project Viva, we used child age at mid-childhood as a covariate. We used household income as a measure of socioeconomic status in Project Viva, and maternal employment type in Gen3G.
Sensitivity and additional analyses
We conducted sensitivity analyses for the cord blood outcomes in Project Viva by including only white children, to include a similar population that was included in the hair cortisol analyses. Using Project Viva data, we performed additional analyses to assess the extent to which the continuous glucose result from the GCT was associated with offspring cortisol biomarkers. In the Gen3G cohort, we performed additional analyses with maternal HbA1c and glucose levels from the OGTT as continuous exposures. Sensitivity analyses including maternal smoking status in the fully adjusted models did not change the strength or the magnitude of the associations in either cohort, so we decided not to include it in our models.
Missing data
We used multiple imputation with SAS PROC MIANALYZEReference Yuan48 to handle missing data. We did not impute missing data in the creation of the figures or the tabulations of summary statistics. We performed all analyses in SAS (version 9.4), and we created all figures in R (version 3.6.1).
Results
Descriptive statistics
In Project Viva, 928 mother–child pairs had at least one child cortisol outcome available for analysis (Table 1). Of these 928 mothers, 83.3% had NGT, 8.8% had IH, 2.7% had IGT, and 5.2% had GDM. Mean (SD) maternal age was 33.1 (4.5) years and mean (SD) prepregnancy BMI was 24.3 (4.8) kg/m2. A total of 449 (48%) children were male. At the time of collection of the hair cortisol measurements at mid-childhood and early adolescence, mean (SD) child age was 7.8 (0.8) and 13.2 (0.9) years, respectively. The distribution of child hair cortisol at each of these time points is shown in Fig. 1.
a NGT = normal glucose tolerance; IH = isolated hyperglycemia; IGT = impaired glucose tolerance; GDM = gestational diabetes mellitus.
b Hair cortisol outcomes are restricted to white participants (N = 632 mid-childhood, N = 563 early adolescence, N = 452 change).
Analyses involving the Gen3G data included 313 women, 9.3% of whom had GDM (Table 2). Mean (SD) maternal age was 28.5 (4.3) years and mean (SD) prepregnancy BMI was 25.0 (5.8) kg/m2. A total of 180 (58%) children were male. Mean (SD) child age at the time of hair collection was 5.4 (0.3) years. The distribution of child hair cortisol in the Gen3G cohort is shown in Fig. 1.
a Data are presented as mean (SD) or median (IQR) for continuous variables and N (%) for categorical variables.
b Prepregnancy BMI was missing for 4 participants.
c N = 97 participants (N = 88 NGT, N = 9 GDM) had hair cortisol measurements below the detectable limit of 1.0 pg/mg. For these participants, we imputed the values to be 0.5 pg/mg.
Associations between maternal glucose status in pregnancy and offspring cortisol markers
We did not observe any associations between maternal pregnancy glucose status and cord blood outcomes in Project Viva (Table 3). For example, the mean cord blood cortisol to cortisone ratio was 0.12 units (95% CI −0.14 to 0.38) higher in women with GDM compared to women with NGT. Our sensitivity analyses for the cord blood outcomes including only white children showed similar results (Supplementary Table S2). We did not observe any evidence of effect modification by mode of delivery for the cord blood outcomes.
a Reporting β (95% CI) comparing to the reference group, which is women with NGT.
b Adjusted for child age and sex. Cord blood outcomes additionally adjusted for child race/ethnicity. In the cord blood analyses, we adjusted for gestational age instead of child age.
c Minimally adjusted model additionally adjusted for maternal prepregnancy BMI, education (≥college graduate vs. <college graduate), and maternal household income at enrollment (>$70,000 per year vs. ≤$70,000 per year).
d We logarithmically transformed all hair cortisol outcomes in the analyses. Hair cortisol analyses included only white children.
Hair cortisol measurements at mid-childhood and early adolescence in Project Viva did not differ between women with NGT and any of the other groups (Table 3). In the models for the hair cortisol outcomes in Project Viva, the effect estimates were close to zero and the confidence intervals all included the null (Table 3). For example, associations between GDM and child hair cortisol were null in fully adjusted models at mid-childhood (β = −0.56 pg/mg, 95% CI −1.16 to 0.04) and early adolescence (β = −0.24 pg/mg, 95% CI −0.79 to 0.31). We conducted additional analyses to investigate the extent to which continuous glucose values measured during the GCT were associated with offspring cortisol biomarkers. However, we only observed null associations (data not shown). Furthermore, we did not observe any evidence of effect modification by child sex for any of the outcomes that we assessed in Project Viva.
In the Gen3G study, we did not observe any associations between GDM and child hair cortisol at age five (β = 0.09 pg/mg, 95% CI −0.38 to 0.57) in fully adjusted analyses. We also did not observe an association of maternal prenatal HbA1c (β = −0.05%, 95% CI −0.54 to 0.44) with child hair cortisol. Results for any of the individual glucose values during OGTT (Table 4) were similarly null. We did not observe any evidence of effect modification by child sex in the Gen3G cohort.
a Reporting β (95% CI).
b Adjusted for child age and sex.
c Minimally adjusted model additionally adjusted for maternal prepregnancy BMI, education, and maternal employment.
Discussion
Using data from two well phenotyped longitudinal prebirth cohorts, we examined associations of maternal glucose markers in pregnancy with cord blood glucocorticoids and offspring hair cortisol levels in mid-childhood and early adolescence. We did not find an association for any of these cortisol biomarkers in both Project Viva and Gen3G cohorts. These results do not support the hypothesis that prenatal exposure to maternal hyperglycemia is associated with dysfunction of the HPA axis in the child.
Few studies have examined the associations between maternal GDM or hyperglycemia and offspring cortisol. Our null findings are in line with some prior reports.Reference Jack-Roberts, Maples and Kalkan25,Reference Berglund, García-Valdés and Torres-Espinola27,Reference Ma, Liu and Wu28,Reference Krishnaveni, Veena and Jones30 For example, a prior study including 54 mother–child pairsReference Mina, Lahti and Drake31 found that maternal fasting glucose levels were not associated with child salivary cortisol reactivity in response to the Marshmallow Test, which is a test of delayed satisfaction that may be used to induce mild stress in young children. Our results are also consistent with a prior study of adverse maternal metabolic conditions which found that GDM was not associated with cord blood cortisol levels.Reference Berglund, García-Valdés and Torres-Espinola27 However, our results are in contrast with prior research that has identified lower cord blood cortisol levelsReference Chen, Guilmette and Luo26 and salivary cortisol levels at ages 11–13Reference Van Dam, Garrett and Schneider29 in children who were exposed to GDM. Notably, these results run counter to our hypothesis that GDM would be associated with higher child cortisol levels, but they are suggestive of potential dysregulation of HPA axis activity. These conflicting results may be due to differences in how cortisol was measured across studies. However, we are not aware of any previous studies that assessed child hair cortisol in relation to GDM. Hair cortisol is different from salivary cortisol because it measures long-term stress over a period of three months,Reference Russell, Koren, Rieder and Van Uum41 whereas salivary cortisol reflects acute stress assessed over a period of 20–30 min.Reference Puhakka and Peltola49 Given the limited amount of research in this area, our study expands upon prior knowledge regarding the relationship between maternal hyperglycemia and child cortisol levels.
This study has strengths and limitations. First, we used data from two prospective cohorts with follow-up into childhood and early adolescence. We also examined child hair cortisol at multiple time points to assess long-term stress axis activation. However, the small number of women with an abnormal pregnancy glucose status (in categorical exposures) is a major limitation that may have limited our statistical power. For example, our sample size in Project Viva would allow us to detect a mean difference of ≥109 nmol/L in cord blood cortisol from the GDM cases in contrast to NGT with 80% power under the assumption of a two-sided statistical test and a type I error rate of 5%. Additionally, we did not examine the association of GDM with diurnal cortisol variation, which could have produced different results because diurnal cortisol variation does not always reflect long-term HPA axis activity. We were unable to exclude participants on the basis of maternal prenatal glucocorticoid or scalp steroid use and were unable to control for shampoo usage in the hair cortisol analyses due to the lack of data on these factors. Moreover, GDM was diagnosed using different criteria in the Project Viva and Gen3G studies. Thus, we cannot conduct a pooled analysis combining data from both studies. However, we believe that a pooled analysis would have resulted in a null association, because we did not observe an association in either study. Finally, the hair cortisol analyses only included white children, so the results of those analyses cannot be generalized to children of other ethnicities.
In conclusion, the present study adds to the limited body of knowledge regarding child cortisol biomarkers in relation to in utero exposure to GDM and glycemic markers. We did not find any association between GDM and child cortisol levels assessed in childhood and early adolescence. Future studies using hair cortisol and other measures of HPA axis function should be conducted in other study populations to assess the generalizability of our results.
Supplementary materials
For supplementary material for this article, please visit https://doi.org/10.1017/S2040174422000381
Acknowledgments
The authors thank the participants and research staff in the Project Viva and Gen3G studies. The work presented in this article was performed as a partial fulfillment of the requirements for the PhD in Public Health with a concentration in Epidemiology for Nathan Cohen.
Financial support
This work was supported by the National Institutes of Health (E.O., grant number R01 HD 034568), (E.O., grant number UH3 OD023286); Fonds de recherche du Québec – Santé (M.F.H., grant number 20697); Canadian Institute of Health Research (M.F.H., grant number MOP-115071), (L.B., grant number PJT-152989); Programme d'aide financière interne (PAFI) from the Centre de recherche du Centre Hospitalier Universitaire de Sherbrooke (P.P., no grant number); a Diabète Québec grant (P.P.); and a Fonds de Recherche du Québec – Société et culture (FRQSC) postdoctoral fellowship award (S.F.). L.B. is a research scholar from the Fonds de recherche du Québec – Santé.
Conflicts of interest
None.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national guidelines on human experimentation (The U.S. Department of Health and Human Services, Office for Human Research Protections and the Tri-Council Policy Statement in Canada) and with the Helsinki Declaration of 1975, as revised in 2008, and has been approved by the institutional committees (Harvard Pilgrim Health Care Institutional Review Board and Centre Hospitalier Universitaire de Sherbrooke Institutional Review Board).