In utero arsenic exposure and early childhood motor development in the New Hampshire Birth Cohort Study

Introduction

High-level prenatal and childhood arsenic (As) exposure characteristic of several regions in Asia (e.g., Bangladesh), may impact motor function. However, the relationship between lower-level arsenic exposure (characteristic of other regions) and motor development is largely unstudied, despite the potential for deficient motor skills in childhood to have adverse long-term consequences. Thus, we sought to investigate the association between prenatal As exposure and motor function among 395 children in the New Hampshire Birth Cohort Study, a rural cohort from northern New England.

Methods

Prenatal exposure was estimated by measuring maternal urine speciated As at 24-28 weeks of gestation using high-performance liquid chromatography (HPLC) inductively coupled plasma mass spectrometry (ICP-MS) and summing inorganic As, monomethylarsonic acid, and dimethylarsinic acid to obtain total urinary As (tAs). Motor function was assessed with the Bruininks-Oseretsky Test of Motor Proficiency, 2nd Edition (BOT-2) at a mean (SD) age of 5.5 (0.4) years.

Results

Children who completed this exam were largely reported as white race (97%), born to married mothers (86%) with a college degree or higher (67%). The median (IQR) gestational urine tAs concentration was 4.0 (5.0) µg/L. Mean (SD) BOT-2 scores were 48.6 (8.4) for overall motor proficiency and 48.2 (9.6) for fine manual control [standard score = 50 (10)], and were 16.3 (5.1) for fine motor integration and 12.5 (4.1) for fine motor precision [standard score = 15 (5)]. We found evidence of a non-linear dose response relationship and used a change-point model to assess the association of tAs with overall motor proficiency and indices of fine motor integration, fine motor precision, and their composite, fine manual control, adjusted for age and sex. In models adjusted for potential confounders, each doubling of urine tAs decreased overall motor proficiency by -3.3 points (95% CI: -6.1, -0.4) for tAs concentrations greater than the change point of 9.5 µg/L and decreased fine motor integration by -4.3 points (95% CI: -8.0, -0.6) for tAs concentrations greater than the change point of 17.0 µg/L.

Discussion

In summary, we found that levels of prenatal As exposure above an empirically-derived threshold (i.e., the change point) were associated with decrements in childhood motor development in a US population.

Prenatal Exposure to Per- and Polyfluoroalkyl Substances and Childhood Autism-related Outcomes

BACKGROUND

Epidemiologic evidence linking prenatal exposure to per- and polyfluoroalkyl substances (PFAS) with altered neurodevelopment is inconclusive, and few large studies have focused on autism-related outcomes. We investigated whether blood concentrations of PFAS in pregnancy are associated with child autism-related outcomes.

METHODS

We included 10 cohorts from the National Institutes of Health (NIH)-funded Environmental influences on Child Health Outcomes (ECHO) program (n = 1,429). We measured 14 PFAS analytes in maternal blood collected during pregnancy; eight analytes met detection criteria for analysis. We assessed quantitative autism-related traits in children via parent report on the Social Responsiveness Scale (SRS). In multivariable linear models, we examined relationships of each PFAS (natural log-transformed) with SRS scores. We further modeled PFAS as a complex mixture using Bayesian methods and examined modification of these relationships by child sex.

RESULTS

Most PFAS in maternal blood were not associated with child SRS T-scores. Perfluorononanoic acid (PFNA) showed the strongest and most consistent association: each 1-unit increase in ln-transformed PFNA was associated with greater autism-related traits (adjusted β [95% confidence interval (CI)] = 1.5 [-0.1, 3.0]). The summed mixture, which included six PFAS detected in >70% of participants, was not associated with SRS T-scores (adjusted β [95% highest posterior density interval] = 0.7 [-1.4, 3.0]). We did not observe consistent evidence of sex differences.

CONCLUSIONS

Prenatal blood concentrations of PFNA may be associated with modest increases in child autism-related traits. Future work should continue to examine the relationship between exposures to both legacy and emerging PFAS and additional dimensional, quantitative measures of childhood autism-related outcomes.

Abstract 6675: Integration of associations of immune profiles in peripheral blood and tumor microenvironment with bladder cancer outcomes

Immune cell profiles in peripheral blood have been associated with bladder cancer outcomes, however, their association with response to immunotherapy and the tumor microenvironment is a major unresolved issue. Although tumor growth can be attenuated via the activation of tumor-infiltrating effector T cells, the relationship between tumor infiltration and immune activation remains unclear. This study explored the interaction between bladder cancer outcomes and immune profiles within peripheral blood and a tumor microenvironment (TME) based on DNA methylation profiles. Peripheral blood and the matched tumor FFPE DNA methylation profiles of 60 non-muscle-invasive bladder cancer (NMIBC) and 12 muscle-invasive bladder cancer (MIBC) patients. Cell-type deconvolution approaches were applied to estimate 12 peripheral immune cell-type proportions and 17 cell-type proportions within TME. We found a positive correlation between dendritic cell proportions in the TME with peripheral CD8T memory cell proportions (r = 0.35, P = 0.003) and a negative correlation between dendritic cell proportions in the TME with peripheral regulatory T cell proportions (r = -0.28, P = 0.021). In addition, monocyte cell proportions in TME had a positive correlation with peripheral B memory (r = 0.37, P = 0.002) and CD8T memory cell proportions (r = 0.43, P = 0.0002). To investigate associations of bladder cancer outcomes with immune cell profiles, using Cox proportional hazard models, we observed an association between the fraction of dendritic cells and the hazard of death (HR = 1.27, 95% CI = 1.06-1.53). Further, a high endothelial cell proportion was significantly associated with an increased hazard of death and tumor recurrence (HR = 1.06, 95% CI = 1.01-1.13) in TME. In addition, the peripheral neutrophil-to-lymphocyte ratio (HR = 1.49, 95% CI = 1.01-2.22), monocyte (HR = 1.17, 95% CI = 1.05-1.31), neutrophil (HR = 1.04, 95% CI = 1.01-1.07), and basophil (HR = 1.35, 95% CI = 1.01-1.81) cell proportions were associated with an increased hazard of death and tumor recurrence. Our results integrated the information on bladder cancer outcomes and cell profiles in TME and peripheral blood, providing biomarkers for estimating bladder cancer prognosis using genome-scale DNA methylation measures.

Citation Format: Ji-Qing Chen, Lucas A. Salas, John K. Wiencke, Devin C. Koestler, Annette M. Molinaro, Angeline S. Andrew, John D. Seigne, Margaret R. Karagas, Karl T. Kelsey, Brock C. Christensen. Integration of associations of immune profiles in peripheral blood and tumor microenvironment with bladder cancer outcomes. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6675.

Health Care Utilization During the COVID-19 Pandemic Among Individuals Born Preterm

Importance

Limited data exist on pediatric health care utilization during the COVID-19 pandemic among children and young adults born preterm.

Objective

To investigate differences in health care use related to COVID-19 concerns during the pandemic among children and young adults born preterm vs those born at term.

Design, Setting, and Participants

In this cohort study, questionnaires regarding COVID-19 and health care utilization were completed by 1691 mother-offspring pairs from 42 pediatric cohorts in the National Institutes of Health Environmental Influences on Child Health Outcomes Program. Children and young adults (ages 1-18 years) in these analyses were born between 2003 and 2021. Data were recorded by the August 31, 2021, data-lock date and were analyzed between October 2021 and October 2022.

Exposures

Premature birth (<37 weeks' gestation).

Main Outcomes and Measures

The main outcome was health care utilization related to COVID-19 concerns (hospitalization, in-person clinic or emergency department visit, phone or telehealth evaluations). Individuals born preterm vs term (≥37 weeks' gestation) and differences among preterm subgroups of individuals (<28 weeks', 28-36 weeks' vs ≥37 weeks' gestation) were assessed. Generalized estimating equations assessed population odds for health care used and related symptoms, controlling for maternal age, education, and psychiatric disorder; offspring history of bronchopulmonary dysplasia (BPD) or asthma; and timing and age at COVID-19 questionnaire completion.

Results

Data from 1691 children and young adults were analyzed; among 270 individuals born preterm, the mean (SD) age at survey completion was 8.8 (4.4) years, 151 (55.9%) were male, and 193 (71.5%) had a history of BPD or asthma diagnosis. Among 1421 comparison individuals with term birth, the mean (SD) age at survey completion was 8.4 (2.4) years, 749 (52.7%) were male, and 233 (16.4%) had a history of BPD or asthma. Preterm subgroups included 159 individuals (58.5%) born at less than 28 weeks' gestation. In adjusted analyses, individuals born preterm had a significantly higher odds of health care utilization related to COVID-19 concerns (adjusted odds ratio [aOR], 1.70; 95% CI, 1.21-2.38) compared with term-born individuals; similar differences were also seen for the subgroup of individuals born at less than 28 weeks' gestation (aOR, 2.15; 95% CI, 1.40-3.29). Maternal history of a psychiatric disorder was a significant covariate associated with health care utilization for all individuals (aOR, 1.44; 95% CI, 1.17-1.78).

Conclusions and Relevance

These findings suggest that during the COVID-19 pandemic, children and young adults born preterm were more likely to have used health care related to COVID-19 concerns compared with their term-born peers, independent of a history of BPD or asthma. Further exploration of factors associated with COVID-19-related health care use may facilitate refinement of care models.

Changes in Body Mass Index Among School-Aged Youths Following Implementation of the Healthy, Hunger-Free Kids Act of 2010

Importance

The prevalence of obesity among youths 2 to 19 years of age in the US from 2017 to 2018 was 19.3%; previous studies suggested that school lunch consumption was associated with increased obesity. The Healthy, Hunger-Free Kids Act of 2010 (HHFKA) strengthened nutritional standards of school-based meals.

Objective

To evaluate the association between the HHFKA and youth body mass index (BMI).

Design, Setting, and Participants

This cohort study was conducted using data from the Environmental Influences on Child Health Outcomes program, a nationwide consortium of child cohort studies, between January 2005 and March 2020. Cohorts in the US of youths aged 5 to 18 years with reported height and weight measurements were included.

Exposures

Full implementation of the HHFKA.

Main Outcomes and Measures

The main outcome was annual BMI z-score (BMIz) trends before (January 2005 to August 2016) and after (September 2016 to March 2020) implementation of the HHFKA, adjusted for self-reported race, ethnicity, maternal education, and cohort group. An interrupted time-series analysis design was used to fit generalized estimating equation regression models.

Results

A total of 14 121 school-aged youths (7237 [51.3%] male; mean [SD] age at first measurement, 8.8 [3.6] years) contributing 26 205 BMI measurements were included in the study. Overall, a significant decrease was observed in the annual BMIz in the period following implementation of the HHFKA compared with prior to implementation (-0.041; 95% CI, -0.066 to -0.016). In interaction models to evaluate subgroup associations, similar trends were observed among youths 12 to 18 years of age (-0.045; 95% CI, -0.071 to -0.018) and among youths living in households with a lower annual income (-0.038; 95% CI, -0.063 to -0.013).

Conclusions and Relevance

In this cohort study, HHFKA implementation was associated with a significant decrease in BMIz among school-aged youths in the US. The findings suggest that school meal programs represent a key opportunity for interventions to combat the childhood obesity epidemic given the high rates of program participation and the proportion of total calories consumed through school-based meals.

Human placental microRNAs dysregulated by cadmium exposure predict neurobehavioral outcomes at birth

BACKGROUND

Prenatal cadmium (Cd) exposure has been implicated in both placental toxicity and adverse neurobehavioral outcomes. Placental microRNAs (miRNAs) may function to developmentally program adverse pregnancy and newborn health outcomes in response to gestational Cd exposure.

METHODS

In a subset of the Rhode Island Child Health Study (RICHS, n = 115) and the New Hampshire Birth Cohort Study (NHBCS, = 281), we used small RNA sequencing and trace metal analysis to identify Cd-associated expression of placental miRNAs using negative binomial generalized linear models. We predicted mRNAs targeted by Cd-associated miRNAs and relate them to neurobehavioral outcomes at birth through the integration of transcriptomic data and summary scores from the NICU Network Neurobehavioral Scale (NNNS).

RESULTS

Placental Cd concentrations are significantly associated with the expression level of five placental miRNAs in NHBCS, with similar effect sizes in RICHS. These miRNA target genes overrepresented in nervous system development, and their expression is correlated with NNNS metrics suggestive of atypical neurobehavioral outcomes at birth.

CONCLUSIONS

Gestational Cd exposure is associated with the expression of placental miRNAs. Predicted targets of these miRNAs are involved in nervous system development and may also regulate placental physiology, allowing their dysregulation to modify developmental programming of early life health outcomes.

IMPACT

This research aims to address the poor understanding of the molecular mechanisms governing adverse pregnancy and newborn health outcomes in response to Gestational cadmium (Cd) exposure. Our results outline a robust relationship between Cd-associated placental microRNA expression and NICU Network Neurobehavioral Scales (NNNS) at birth indicative of atypical neurobehavior. This study utilized healthy mother-infant cohorts to describe the role of Cd-associated dysregulation of placental microRNAs as a potential mechanism by which adverse neurobehavioral outcomes are developmentally programmed.

Perinatal Outcomes during versus Prior to the COVID-19 Pandemic and the Role of Maternal Depression and Perceived Stress: A Report from the ECHO Program

OBJECTIVE

 We sought to evaluate the impact of the coronavirus disease 2019 (COVID-19) pandemic on perinatal outcomes while accounting for maternal depression or perceived stress and to describe COVID-specific stressors, including changes in prenatal care, across specific time periods of the pandemic.

STUDY DESIGN

 Data of dyads from 41 cohorts from the National Institutes of Health Environmental influences on Child Health Outcomes Program (N = 2,983) were used to compare birth outcomes before and during the pandemic (n = 2,355), and a partially overlapping sample (n = 1,490) responded to a COVID-19 questionnaire. Psychosocial stress was defined using prenatal screening for depression and perceived stress. Propensity-score matching and general estimating equations with robust variance estimation were used to estimate the pandemic's effect on birth outcomes.

RESULTS

 Symptoms of depression and perceived stress during pregnancy were similar prior to and during the pandemic, with nearly 40% of participants reporting mild to severe stress, and 24% reporting mild depression to severe depression. Gestations were shorter during the pandemic (B =  - 0.33 weeks, p = 0.025), and depression was significantly associated with shortened gestation (B =  - 0.02 weeks, p = 0.015) after adjustment. Birth weights were similar (B =  - 28.14 g, p = 0.568), but infants born during the pandemic had slightly larger birth weights for gestational age at delivery than those born before the pandemic (B = 0.15 z-score units, p = 0.041). More women who gave birth early in the pandemic reported being moderately or extremely distressed about changes to their prenatal care and delivery (45%) compared with those who delivered later in the pandemic. A majority (72%) reported somewhat to extremely negative views of the impact of COVID-19 on their life.

CONCLUSION

 In this national cohort, we detected no effect of COVID-19 on prenatal depression or perceived stress. However, experiencing the COVID-19 pandemic in pregnancy was associated with decreases in gestational age at birth, as well as distress about changes in prenatal care early in the pandemic.

KEY POINTS

· COVID-19 was associated with shortened gestations.. · Depression was associated with shortened gestations.. · However, stress during the pandemic remained unchanged.. · Most women reported negative impacts of the pandemic..

Impact of sedentary behavior and emotional support on prenatal psychological distress and birth outcomes during the COVID-19 pandemic

Abstract.

BACKGROUND

Studies have reported mixed findings regarding the impact of the coronavirus disease 2019 (COVID-19) pandemic on pregnant women and birth outcomes. This study used a quasi-experimental design to account for potential confounding by sociodemographic characteristics.

METHODS

Data were drawn from 16 prenatal cohorts participating in the Environmental influences on Child Health Outcomes (ECHO) program. Women exposed to the pandemic (delivered between 12 March 2020 and 30 May 2021) (n = 501) were propensity-score matched on maternal age, race and ethnicity, and child assigned sex at birth with 501 women who delivered before 11 March 2020. Participants reported on perceived stress, depressive symptoms, sedentary behavior, and emotional support during pregnancy. Infant gestational age (GA) at birth and birthweight were gathered from medical record abstraction or maternal report.

RESULTS

After adjusting for propensity matching and covariates (maternal education, public assistance, employment status, prepregnancy body mass index), results showed a small effect of pandemic exposure on shorter GA at birth, but no effect on birthweight adjusted for GA. Women who were pregnant during the pandemic reported higher levels of prenatal stress and depressive symptoms, but neither mediated the association between pandemic exposure and GA. Sedentary behavior and emotional support were each associated with prenatal stress and depressive symptoms in opposite directions, but no moderation effects were revealed.

CONCLUSIONS

There was no strong evidence for an association between pandemic exposure and adverse birth outcomes. Furthermore, results highlight the importance of reducing maternal sedentary behavior and encouraging emotional support for optimizing maternal health regardless of pandemic conditions.

Birth Outcomes in Relation to Prenatal Exposure to Per- and Polyfluoroalkyl Substances and Stress in the Environmental Influences on Child Health Outcomes (ECHO) Program

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are persistent and ubiquitous chemicals associated with risk of adverse birth outcomes. Results of previous studies have been inconsistent. Associations between PFAS and birth outcomes may be affected by psychosocial stress.

OBJECTIVES

We estimated risk of adverse birth outcomes in relation to prenatal PFAS concentrations and evaluate whether maternal stress modifies those relationships.

METHODS

We included 3,339 participants from 11 prospective prenatal cohorts in the Environmental influences on the Child Health Outcomes (ECHO) program to estimate the associations of five PFAS and birth outcomes. We stratified by perceived stress scale scores to examine effect modification and used Bayesian Weighted Sums to estimate mixtures of PFAS.

RESULTS

We observed reduced birth size with increased concentrations of all PFAS. For a 1-unit higher log-normalized exposure to perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorononanoic acid (PFNA), and perfluorohexane sulfonic acid (PFHxS), we observed lower birthweight-for-gestational-age z-scores of β = - 0.15 [95% confidence interval (CI): - 0.27 , - 0.03 ], β = - 0.14 (95% CI: - 0.28 , - 0.002 ), β = - 0.22 (95% CI: - 0.23 , - 0.10 ), β = - 0.06 (95% CI: - 0.18 , 0.06), and β = - 0.25 (95% CI: - 0.37 , - 0.14 ), respectively. We observed a lower odds ratio (OR) for large-for-gestational-age: OR PFNA = 0.56 (95% CI: 0.38, 0.83), OR PFDA = 0.52 (95% CI: 0.35, 0.77). For a 1-unit increase in log-normalized concentration of summed PFAS, we observed a lower birthweight-for-gestational-age z-score [- 0.28 ; 95% highest posterior density (HPD): - 0.44 , - 0.14 ] and decreased odds of large-for-gestational-age (OR = 0.49 ; 95% HPD: 0.29, 0.82). Perfluorodecanoic acid (PFDA) explained the highest percentage (40%) of the summed effect in both models. Associations were not modified by maternal perceived stress.

DISCUSSION

Our large, multi-cohort study of PFAS and adverse birth outcomes found a negative association between prenatal PFAS and birthweight-for-gestational-age, and the associations were not different in groups with high vs. low perceived stress. This study can help inform policy to reduce exposures in the environment and humans. https://doi.org/10.1289/EHP10723.

Application of novel breast biospecimen cell type adjustment identifies shared DNA methylation alterations in breast tissue and milk with breast cancer risk factors

BACKGROUND

DNA methylation patterning is cell-type-specific and altered DNA methylation is well established to occur early in breast carcinogenesis, affecting non-cancerous, histopathologically normal breast tissue. Previous work assessing risk factor-associated alterations to DNA methylation in breast tissue has been limited, with even less published research in breast milk, a non-invasively obtained biospecimen containing sloughed mammary epithelial cells that may identify early alterations indicative of cancer risk.

METHODS

Here, we present a novel library for the estimation of the cellular composition of breast tissue and milk and subsequent assessment of cell type-independent alterations to DNA methylation associated with established breast cancer risk factors in solid breast tissue (n = 95) and breast milk (n = 48) samples using genome-scale DNA methylation measures from the Illumina HumanMethylation450 array.

RESULTS

We identified 772 hypermethylated CpGs (P < 0.01) associated with age consistent between breast tissue and breast milk samples. Age-associated hypermethylated CpG loci were significantly enriched for CpG island shore regions known to be important for regulating gene expression. Among the overlapping hypermethylated loci mapping to genes, a differentially methylated region was identified in the promoter region of SFRP2, a gene observed to undergo promoter hypermethylation in breast cancer.

CONCLUSIONS

Our findings suggest the potential to identify epigenetic biomarkers of breast cancer risk in noninvasively obtained, tissue-specific breast milk specimens.

IMPACT

This work demonstrates the potential of using breast milk as a non-invasive biomarker of breast cancer risk, improving our ability to detect early stage disease and lowering the overall disease burden.

Analysis of Pregnancy Complications and Epigenetic Gestational Age of Newborns

Importance

Preeclampsia, gestational hypertension, and gestational diabetes, the most common pregnancy complications, are associated with substantial morbidity and mortality in mothers and children. Little is known about the biological processes that link the occurrence of these pregnancy complications with adverse child outcomes; altered biological aging of the growing fetus up to birth is one molecular pathway of increasing interest.

Objective

To evaluate whether exposure to each of these 3 pregnancy complications (gestational diabetes, gestational hypertension, and preeclampsia) is associated with accelerated or decelerated gestational biological age in children at birth.

Design, Setting, and Participants

Children included in these analyses were born between 1998 and 2018 and spanned multiple geographic areas of the US. Pregnancy complication information was obtained from maternal self-report and/or medical record data. DNA methylation measures were obtained from blood biospecimens collected from offspring at birth. The study used data from the national Environmental Influences on Child Health Outcomes (ECHO) multisite cohort study collected and recorded as of the August 31, 2021, data lock date. Data analysis was performed from September 2021 to December 2022.

Exposures

Three pregnancy conditions were examined: gestational hypertension, preeclampsia, and gestational diabetes.

Main Outcomes and Measures

Accelerated or decelerated biological gestational age at birth, estimated using existing epigenetic gestational age clock algorithms.

Results

A total of 1801 child participants (880 male [48.9%]; median [range] chronological gestational age at birth, 39 [30-43] weeks) from 12 ECHO cohorts met the analytic inclusion criteria. Reported races included Asian (49 participants [2.7%]), Black (390 participants [21.7%]), White (1026 participants [57.0%]), and other races (92 participants [5.1%]) (ie, American Indian or Alaska Native, Native Hawaiian or other Pacific Islander, multiple races, and other race not specified). In total, 524 participants (29.0%) reported Hispanic ethnicity. Maternal ages ranged from 16 to 45 years of age with a median of 29 in the analytic sample. A range of maternal education levels, from less than high school (260 participants [14.4%]) to Bachelor's degree and above (629 participants [34.9%]), were reported. In adjusted regression models, prenatal exposure to maternal gestational diabetes (β, -0.423; 95% CI, -0.709 to -0.138) and preeclampsia (β, -0.513; 95% CI, -0.857 to -0.170), but not gestational hypertension (β, 0.003; 95% CI, -0.338 to 0.344), were associated with decelerated epigenetic aging among exposed neonates vs those who were unexposed. Modification of these associations, by sex, was observed with exposure to preeclampsia (β, -0.700; 95% CI, -1.189 to -0.210) and gestational diabetes (β, -0.636; 95% CI, -1.070 to -0.200), with associations observed among female but not male participants.

Conclusions and Relevance

This US cohort study of neonate biological changes related to exposure to maternal pregnancy conditions found evidence that preeclampsia and gestational diabetes delay biological maturity, especially in female offspring.

Early-life exposure to per- and polyfluoroalkyl substances and infant gut microbial composition

Human milk is rich in essential nutrients and immune-activating compounds but is also a source of toxicants including per- and polyfluoroalkyl substances (PFAS). Evidence suggests that immune-related effects of PFAS may, in part, be due to alterations of the microbiome. We aimed to identify the association between milk PFAS exposure and the infant gut microbiome.

Methods

PFAS [perfluorooctane sulfonic acid (PFOS) and perfluorooctanoate (PFOA)] were quantified in milk from ~6 weeks postpartum using high-performance liquid chromatography with tandem mass spectrometry. A molar sum (ΣPFAS) was calculated. Caregivers collected infant stool samples at 6 weeks (n = 116) and/or 1 year postpartum (n = 119). Stool DNA underwent metagenomic sequencing. We estimated the association of PFAS with diversity and relative abundances of species with linear regression. Single- and multi-PFAS models adjusted for potential confounders in complete case analyses and with imputed missing covariate data for 6-week and 1-year microbiomes separately. We assessed sensitive populations with stratification.

Results

PFOS and PFOA were detected in 94% and 83% of milk samples, respectively. PFOS was associated with increased diversity at 6 weeks among infants fed exclusively human milk [β = 0.24 per PFOS doubling, (95% CI = 0.03, 0.45), P = 0.03] and born to primiparous mothers [β = 0.37 (0.06, 0.67), P = 0.02]. Estimates were strongest in multi-PFAS models and among complete cases. ΣPFAS was associated with Bacteroides vulgatus relative abundance at 1 year [(β = -2.34% per doubling (-3.63, -1.05), FDR q = 0.099].

Conclusions

PFAS may increase infant gut microbiome diversity and alter the relative abundance of biologically relevant bacteria. Additional analyses may identify related health outcomes.

Impact of the COVID-19 pandemic on children's sleep habits: an ECHO study

BACKGROUND

Sleep in childhood is affected by behavioral, environmental, and parental factors. We propose that these factors were altered during the COVID-19 pandemic. This study investigates sleep habit changes during the pandemic in 528 children 4-12 years old in the US, leveraging data from the Environmental Influences on Child Health Outcomes (ECHO) Program.

METHODS

Data collection occurred in July 2019-March 2020 (pre-pandemic) and two pandemic periods: December 2020-April 2021 and May-August 2021. Qualitative interviews were performed in 38 participants.

RESULTS

We found no changes in sleep duration, but a shift to later sleep midpoint during the pandemic periods. There was an increase in latency at the first pandemic collection period but no increase in the frequency of bedtime resistance, and a reduced frequency of naps during the pandemic. Qualitative interviews revealed that parents prioritized routines to maintain sleep duration but were more flexible regarding timing. Children from racial/ethnic minoritized communities slept less at night, had later sleep midpoint, and napped more frequently across all collection periods, warranting in-depth investigation to examine and address root causes.

CONCLUSIONS

The COVID-19 pandemic significantly impacted children sleep, but parental knowledge of the importance of sleep might have played a significant protective role.

IMPACT

During the COVID-19 pandemic, US children changed their sleep habits, going to bed and waking up later, but their sleep duration did not change. Sleep latency was longer. Parental knowledge of sleep importance might have played a protective role. Regardless of data collection periods, children from racial/ethnic minoritized communities slept less and went to bed later. This is one of the first study on this topic in the US, including prospective pre-pandemic qualitative and quantitative data on sleep habits. Our findings highlight the pandemic long-term impact on childhood sleep. Results warrants further investigations on implications for overall childhood health.

Trends in Screen Time Use Among Children During the COVID-19 Pandemic, July 2019 Through August 2021

Importance

The COVID-19 pandemic led to widespread lockdowns and school closures that may have affected screen time among children. Although restrictions were strongest early in the pandemic, it is unclear how screen time changed as the pandemic progressed.

Objective

To evaluate change in children's screen time from before the pandemic to during the pandemic, from July 2019 through August 2021.

Design, Setting, and Participants

This is a longitudinal cohort study with repeated measures of screen time collected before the pandemic and during 2 pandemic periods. Children aged 4 to 12 years and their parent were enrolled in 3 pediatric cohorts across 3 states in the US participating in the Environmental Influences of Child Health Outcomes (ECHO) Program. Data analysis was performed from November 2021 to July 2022.

Exposures

COVID-19 pandemic period: prepandemic (July 2019 to March 2020), pandemic period 1 (December 2020 to April 2021), and pandemic period 2 (May 2021 to August 2021).

Main Outcomes and Measures

The primary outcomes were total, educational (not including remote school), and recreational screen time assessed via the ECHO Child Media Use questionnaire. Linear mixed-effects models were used for screen time adjusted for child's age, number of siblings, sex, race, ethnicity, and maternal education.

Results

The cohort included 228 children (prepandemic mean [SD] age, 7.0 [2.7] years; 100 female [43.9%]) with screen time measured during the prepandemic period and at least once during the pandemic period. Prepandemic mean (SD) total screen time was 4.4 (3.9) hours per day and increased 1.75 hours per day (95% CI, 1.18-2.31 hours per day) in the first pandemic period and 1.11 hours per day (95% CI, 0.49-1.72 hours per day) in the second pandemic period, in adjusted models. Prepandemic mean (SD) recreational screen time was 4.0 (3.5) hours per day and increased 0.89 hours per day (95% CI, 0.39-1.39 hours per day) in the first pandemic period and 0.70 hours per day (95% CI, 0.16-1.25 hours per day) in the second pandemic period. Prepandemic mean (SD) educational screen time was 0.5 (1.2) hours per day (median [IQR], 0.0 [0.0-0.4] hours per day) and increased 0.93 hours per day (95% CI, 0.67-1.19 hours per day) in the first pandemic period and 0.46 hours per day (95% CI, 0.18-0.74 hours per day) in the second pandemic period.

Conclusions and Relevance

These findings suggest that screen time among children increased during the COVID-19 pandemic and remained elevated even after many public health precautions were lifted. The long-term association of increased screen time during the COVID-19 pandemic with children's health needs to be determined.

Characterizing changes in behaviors associated with chemical exposures during the COVID-19 pandemic

The COVID-19 pandemic-and its associated restrictions-have changed many behaviors that can influence environmental exposures including chemicals found in commercial products, packaging and those resulting from pollution. The pandemic also constitutes a stressful life event, leading to symptoms of acute traumatic stress. Data indicate that the combination of environmental exposure and psychological stress jointly contribute to adverse child health outcomes. Within the Environmental influences on Child Health Outcomes (ECHO)-wide Cohort, a national consortium initiated to understand the effects of environmental exposures on child health and development, our objective was to assess whether there were pandemic-related changes in behavior that may be associated with environmental exposures. A total of 1535 participants from nine cohorts completed a survey via RedCap from December 2020 through May 2021. The questionnaire identified behavioral changes associated with the COVID-19 pandemic in expected directions, providing evidence of construct validity. Behavior changes reported by at least a quarter of the respondents include eating less fast food and using fewer ultra-processed foods, hair products, and cosmetics. At least a quarter of respondents reported eating more home cooked meals and using more antibacterial soaps, liquid soaps, hand sanitizers, antibacterial and bleach cleaners. Most frequent predictors of behavior change included Hispanic ethnicity and older age (35 years and older). Respondents experiencing greater COVID-related stress altered their behaviors more than those not reporting stress. These findings highlight that behavior change associated with the pandemic, and pandemic-related psychological stress often co-occur. Thus, prevention strategies and campaigns that limit environmental exposures, support stress reduction, and facilitate behavioral change may lead to the largest health benefits in the context of a pandemic. Analyzing biomarker data in these participants will be helpful to determine if behavior changes reported associate with measured changes in exposure.

Concentrations of Per- and Polyfluoroalkyl Substances in Paired Maternal Plasma and Human Milk in the New Hampshire Birth Cohort

Per- and polyfluoroalkyl substances (PFAS) are environmentally persistent endocrine-disrupting chemicals associated with long-term health outcomes. PFAS are transferred from maternal blood to human milk, an important exposure source for infants, and understanding of this transfer is evolving. We characterized concentrations of 10 PFAS in human milk (n = 426) and compared milk-to-plasma concentrations of 9 PFAS among a subset of women with paired samples (n = 294) from the New Hampshire Birth Cohort Study using liquid chromatography-isotope dilution tandem mass spectrometry. We examined the relationship between perfluorooctanoate (PFOA) and perfluorooctane sulfonate (PFOS) in plasma versus milk and fit linear regression models to assess relationships between milk PFOA and PFOS and participant characteristics. The median plasma PFOA concentration was 0.94 ng/mL (interquartile range, IQR, 0.59-1.34) and that of PFOS was 2.60 ng/mL (IQR 1.80-3.90); the median milk PFOA concentration was 0.017 ng/mL (IQR 0.012-0.027) and that of PFOS was 0.024 ng/mL (IQR 0.016-0.036). PFOA and PFOS plasma and milk concentrations showed correlations of ρ = 0.83 and 0.77, respectively (p < 0.001). Parity, previous lactation, week of milk collection, and body mass index were inversely associated with milk PFAS. We estimate that even among our general population cohort, some infants (∼6.5%) are exposed to amounts of PFAS via milk that may have long-term health impacts.

Associations of maternal plasma and umbilical cord plasma metabolomics profiles with birth anthropometric measures

BACKGROUND

The metabolomics profiles of maternal plasma during pregnancy and cord plasma at birth might influence fetal growth and birth anthropometry. The objective was to examine how maternal plasma and umbilical cord plasma metabolites are associated with newborn anthropometric measures, a known predictor of future health outcomes.

METHODS

Pregnant women between 24 and 28 weeks of gestation were recruited as part of a prospective cohort study. Blood samples from 413 women at enrollment and 787 infant cord blood samples were analyzed using the Biocrates AbsoluteIDQ® p180 kit. Multivariable linear regression models were used to examine associations of cord and maternal metabolites with infant anthropometry at birth.

RESULTS

In cord blood samples from this rural cohort from New Hampshire of largely white residents, 13 metabolites showed negative associations, and 10 metabolites showed positive associations with birth weight Z-score. Acylcarnitine C5 showed negative association, and 4 lysophosphatidylcholines showed positive associations with birth length Z-score. Maternal blood metabolites did not significantly correlate with birth weight and length Z-scores.

CONCLUSIONS

Consistent findings were observed for several acylcarnitines that play a role in utilization of energy sources, and a lysophosphatidylcholine that is part of oxidative stress and inflammatory response pathways in cord plasma samples.

IMPACT

The metabolomics profiles of maternal plasma during pregnancy and cord plasma at birth may influence fetal growth and birth anthropometry. This study examines the independent effects of maternal gestational and infant cord blood metabolomes across different classes of metabolites on birth anthropometry. Acylcarnitine species were negatively associated and glycerophospholipids species were positively associated with weight and length Z-scores at birth in the cord plasma samples, but not in the maternal plasma samples. This study identifies lipid metabolites in infants that possibly may affect early growth.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

Human milk-associated bacterial communities associate with the infant gut microbiome over the first year of life

Introduction

Microbial communities inhabiting the human infant gut are important for immune system development and lifelong health. One critical exposure affecting the bacterial colonization of the infant gut is consumption of human milk, which contains diverse microbial communities and prebiotics. We hypothesized that human milk-associated microbial profiles are associated with those of the infant gut.

Methods

Maternal-infant dyads enrolled in the New Hampshire Birth Cohort Study (n = 189 dyads) contributed breast milk and infant stool samples collected approximately at 6 weeks, 4 months, 6 months, 9 months, and 12 months postpartum (n = 572 samples). Microbial DNA was extracted from milk and stool and the V4-V5 region of the bacterial 16S rRNA gene was sequenced.

Results

Clustering analysis identified three breast milk microbiome types (BMTs), characterized by differences in Streptococcus, Staphylococcus, Pseudomonas, Acinetobacter, and microbial diversity. Four 6-week infant gut microbiome types (6wIGMTs) were identified, differing in abundances of Bifidobacterium, Bacteroides, Clostridium, Streptococcus, and Escherichia/Shigella, while two 12-month IGMTs (12mIGMTs) differed primarily by Bacteroides presence. At 6 weeks, BMT was associated with 6wIGMT (Fisher's exact test value of p = 0.039); this association was strongest among infants delivered by Cesarean section (Fisher's exact test value of p = 0.0028). The strongest correlations between overall breast milk and infant stool microbial community structures were observed when comparing breast milk samples to infant stool samples collected at a subsequent time point, e.g., the 6-week breast milk microbiome associated with the 6-month infant gut microbiome (Mantel test Z-statistic = 0.53, value of p = 0.001). Streptoccous and Veillonella species abundance were correlated in 6-week milk and infant stool, and 4- and 6-month milk Pantoea species were associated with infant stool Lachnospiraceae genera at 9 and 12  months.

Discussion

We identified clusters of human milk and infant stool microbial communities that were associated in maternal-infant dyads at 6 weeks of life and found that milk microbial communities were more strongly associated with infant gut microbial communities in infants delivered operatively and after a lag period. These results suggest that milk microbial communities have a long-term effect on the infant gut microbiome both through sharing of microbes and other molecular mechanisms.

Abstract P038: Transcriptome-wide association study identifies novel genes associated with bladder cancer risk

Background: Genome-wide association studies (GWAS) have identified several common germline variants associated with bladder cancer risk. However, the genetic architecture of bladder cancer risk is still largely unknown. Therefore, we conducted a transcriptome-wide association study (TWAS) to identify gene expression patterns associated with bladder cancer risk. Methods: Individual-level genotype data from bladder cancer cases and controls was accessed from dbGaP (phs000346.v2.p2). We included 3,630 muscle invasive and non-muscle invasive bladder cancer cases and 3,356 controls from the Spanish Bladder Cancer Study (1,145 cases and 1,083 controls) and the New England Bladder Cancer Study (Maine and Vermont, 2,485 cases and 2,273 controls) with imputed genotype data for 22 million single-nucleotide polymorphisms (SNPs). We used PrediXcan to predict gene expression in whole blood using reference data from the Depression Genes and Networks (DGN). Covariate-adjusted logistic regression was used to estimate the associations between gene expression and bladder cancer risk. Results were validated using data from 280 non-muscle invasive bladder cancer cases and 353 controls of the New Hampshire Bladder Cancer Study. Results: We identified five genes associated with bladder cancer risk, including a previously reported locus (1p13.3: GSTM1) and four other loci (4p16.3: FAM53A, 7q22.3: RP11-325F22.2, 8p23.1: RP11-297N6.4, 11q24.1: JHY). Among them, higher predicted expression of GSTM1 in whole blood was associated with lower risk of bladder cancer (Discovery dataset OR = 0.78, 95% CI = 0.67, 0.92; Validation dataset OR = 0.58, 95% CI = 0.34, 0.98). Higher predicted expression of FAM53A (Discovery OR = 1.34, 95% CI = 1.08, 1.68; Validation OR = 3.79, 95% CI = 1.76, 8.29), RP11-325F22.2 (Discovery OR = 1.34, 95% CI = 1.09, 1.65; Validation OR = 2.11, 95% CI = 1.05, 4.30), RP11-297N6.4 (Discovery OR = 1.78, 95% CI = 1.22, 2.60; Validation OR = 6.32, 95% CI = 1.74, 23.34) and JHY (Discovery OR = 1.94, 95% CI = 1.30, 2.91; Validation OR = 6.78, 95% CI = 1.85, 25.76) in whole blood was associated with higher risk of bladder cancer. Conclusion: Our TWAS identified five genes associated with bladder cancer risk: GSTM1, a previously reported locus associated with bladder cancer risk, is related to metabolic detoxification. Four novel loci were identified: FAM53A (related to neural tube development and was previously reported to be related to breast cancer progression), JHY (related to axoneme structure), RP11-325F22.2 (lincRNA), and RP11-297N6.4 (lncRNA). These loci require further validation and may help to reveal the underlying mechanism of the relationship between gene expression and bladder cancer risk.

Citation Format: Siting Li, Jiang Gui, Margaret R. Karagas, Michael N. Passarelli. Transcriptome-wide association study identifies novel genes associated with bladder cancer risk. [abstract]. In: Proceedings of the AACR Special Conference: Precision Prevention, Early Detection, and Interception of Cancer; 2022 Nov 17-19; Austin, TX. Philadelphia (PA): AACR; Can Prev Res 2023;16(1 Suppl): Abstract nr P038.

Association of Severe Bronchiolitis during Infancy with Childhood Asthma Development: An Analysis of the ECHO Consortium

Objective: Many studies have shown that severe (hospitalized) bronchiolitis during infancy is a risk factor for developing childhood asthma. However, the population subgroups at the highest risk remain unclear. Using large nationwide pediatric cohort data, namely the NIH Environmental influences on Child Health Outcomes (ECHO) Program, we aimed to quantify the longitudinal relationship of bronchiolitis hospitalization during infancy with asthma in a generalizable dataset and to examine potential heterogeneity in terms of major demographics and clinical factors. Methods: We analyzed data from infants (age <12 months) enrolled in one of the 53 prospective cohort studies in the ECHO Program during 2001−2021. The exposure was bronchiolitis hospitalization during infancy. The outcome was a diagnosis of asthma by a physician by age 12 years. We examined their longitudinal association and determined the potential effect modifications of major demographic factors. Results: The analytic cohort consisted of 11,762 infants, 10% of whom had bronchiolitis hospitalization. Overall, 15% subsequently developed asthma. In the Cox proportional hazards model adjusting for 10 patient-level factors, compared with the no-bronchiolitis hospitalization group, the bronchiolitis hospitalization group had a significantly higher rate of asthma (14% vs. 24%, HR = 2.77, 95%CI = 2.24−3.43, p < 0.001). There was significant heterogeneity by race and ethnicity (Pinteraction = 0.02). The magnitude of the association was greater in non-Hispanic White (HR = 3.77, 95%CI = 2.74−5.18, p < 0.001) and non-Hispanic Black (HR = 2.39, 95%CI = 1.60−3.56; p < 0.001) infants, compared with Hispanic infants (HR = 1.51, 95%CI = 0.77−2.95, p = 0.23). Conclusions: According to the nationwide cohort data, infants hospitalized with bronchiolitis are at a higher risk for asthma, with quantitative heterogeneity in different racial and ethnic groups.

Umbilical cord blood immune cell profiles in relation to the infant gut microbiome

During infancy, the interplay between the developing immune system and the microbiome is critical. We examined whether blood immune cell composition at birth in the umbilical cord (inferred by DNA methylation profiling) related to the early infant gut microbiome (assessed by 16S rRNA gene sequencing) among 73 infants in the New Hampshire Birth Cohort Study. We used generalized estimating equations and controlled for false discovery rate to select microbial taxa associated with immune cells. We found associations between the infant gut microbiome and immune cells, including a positive association between B cells and Enterobacter, a negative association between natural killer cells and Bifidobacterium, and a positive association between granulocytes and Bifidobacterium. Our findings give clues that immune profiles at the time of birth as measured in umbilical cord blood are associated with the development of the gut microbiome in early life.

COVID-19 Pandemic Experiences and Symptoms of Pandemic-Associated Traumatic Stress Among Mothers in the US

Importance

The primary outcomes of the COVID-19 pandemic on the mental health of women with children remain largely unknown.

Objectives

To identify and describe clusters of mothers of children participating in the Environmental influences on Child Health Outcomes (ECHO) Program that characterize pandemic-associated hardships, coping mechanisms, and behaviors, and to evaluate associations between pandemic-associated hardships, coping strategies, and behavior changes with pandemic-associated traumatic stress symptoms.

Design, Setting, and Participants

This multicenter cohort study investigated experiences during the COVID-19 pandemic between April 2020 and August 2021 among maternal caregivers of children participating in the ECHO Program. Data from self-identified mothers of ECHO-enrolled children from 62 US cohorts were included in analyses. Data were analyzed from November 2021 to July 2022.

Exposures

The primary exposures were pandemic-associated changes in mothers' health, health care utilization, work and finances, coping strategies, and health-associated behaviors. Exposures were assessed via a self-reported questionnaire designed by ECHO investigators.

Main Outcomes and Measures

The primary outcome was the total symptoms score of pandemic-associated traumatic stress (PTS), defined as the number of items endorsed at least sometimes or more frequently, from a 10-item self-report measure.

Results

The study surveyed 11 473 mothers (mean [SD] age, 37.8 [7.4] years; 342 American Indian [2.98%], 378 Asian [3.29%], 1701 Black [14.83%], and 7195 White [62.71%]; 2184 with Hispanic/Latina ethnicity [19.04%]) and identified 2 clusters that best characterized their COVID-19 pandemic experiences-one characterized by higher life disruptions (eg, to work and health care), higher social isolation, more coping behaviors to mitigate the outcomes of the pandemic, and more changes to their health behavior routines (high change [1031 mothers]) and the other characterized by lower changes (low change [3061 mothers]). The high change cluster was more socioeconomically advantaged and reported higher PTS (mean [SD] number of symptoms, 3.72 [2.44] vs 2.51 [2.47]). Across both clusters, higher pandemic-associated hardships, coping mechanisms, and behavior changes were associated with higher PTS, and these associations were greater in the low change cluster.

Conclusions and Relevance

In this study of more than 11 000 US mothers, associations between socioeconomic factors, stressful life events, and mental health sequelae were complex. Accordingly, programs, policies, and practices targeting mental health during public health crises such as the COVID-19 pandemic should consider the range and configuration of hardships in designing the most effective interventions to mitigate long-term outcomes.

Alterations in Microbial-Associated Fecal Metabolites in Relation to Arsenic Exposure Among Infants

In utero and early life exposure to inorganic arsenic (iAs) alters immune response in experimental animals and is associated with an increased risk of infant infections. iAs exposure is related to differences in the gut microbiota diversity, community structure, and the relative abundance of individual microbial taxa both in laboratory and human studies. Metabolomics permits a direct measure of molecular products of microbial and host metabolic processes. We conducted NMR metabolomics analysis on infant stool samples and quantified the relative concentrations of 34 known microbial-related metabolites. We examined these metabolites in relation to both in utero and infant log2 urinary total arsenic concentrations (utAs, the sum of iAs and iAs metabolites) collected at approximately 6 weeks of age using linear regression models, adjusted for infant sex, age at sample collection, type of delivery (vaginal vs. cesarean section), feeding mode (breast milk vs. any formula), and specific gravity. Increased fecal butyrate (b = 214.24), propionate (b = 518.33), cholate (b = 8.79), tryptophan (b= 14.23), asparagine (b = 28.80), isoleucine (b = 65.58), leucine (b = 95.91), malonate (b = 50.43), and uracil (b = 36.13), concentrations were associated with a doubling of infant utAs concentrations (p< 0.05). These associations were largely among infants who were formula fed. No clear associations were observed with maternal utAs and infant fecal metabolites. Metabolomic analyses of infant stool samples lend further evidence that the infant gut microbiota is sensitive to As exposure, and these effects may have functional consequences.

Impact of antibiotics on off-target infant gut microbiota and resistance genes in cohort studies

BACKGROUND

Young children are frequently exposed to antibiotics, with the potential for collateral consequences to the gut microbiome. The impact of antibiotic exposures to off-target microbes (i.e., bacteria not targeted by treatment) and antibiotic resistance genes (ARGs) is poorly understood.

METHODS

We used metagenomic sequencing data from paired stool samples collected prior to antibiotic exposure and at 1 year from over 200 infants and a difference-in-differences approach to assess the relationship between subsequent exposures and the abundance or compositional diversity of microbes and ARGs while adjusting for covariates.

RESULTS

By 1 year, the abundance of multiple species and ARGs differed by antibiotic exposure. Compared to infants never exposed to antibiotics, Bacteroides vulgatus relative abundance increased by 1.72% (95% CI: 0.19, 3.24) while Bacteroides fragilis decreased by 1.56% (95% CI: -4.32, 1.21). Bifidobacterium species also exhibited opposing trends. ARGs associated with exposure included class A beta-lactamase gene CfxA6. Among infants attending day care, Escherichia coli and ARG abundance were both positively associated with antibiotic use.

CONCLUSION

Novel findings, including the importance of day care attendance, were identified through considering microbiome data at baseline and post-intervention. Thus, our study design and approach have important implications for future studies evaluating the unintended impacts of antibiotics.

IMPACT

The impact of antibiotic exposure to off-target microbes and antibiotic resistance genes in the gut is poorly defined. We quantified these impacts in two cohort studies using a difference-in-differences approach. Novel to microbiome studies, we used pre/post-antibiotic data to emulate a randomized controlled trial. Compared to infants unexposed to antibiotics between baseline and 1 year, the relative abundance of multiple off-target species and antibiotic resistance genes was altered. Infants who attended day care and were exposed to antibiotics within the first year had a higher abundance of Escherichia coli and antibiotic resistance genes; a novel finding warranting further investigation.

Association between fat-infiltrated axillary lymph nodes on screening mammography and cardiometabolic disease

Objective

Ectopic fat deposition within and around organs is a stronger predictor of cardiometabolic disease status than body mass index (BMI). Fat deposition within the lymphatic system is poorly understood. This study examined the association between the prevalence of cardiometabolic disease and ectopic fat deposition within axillary lymph nodes (LNs) visualized on screening mammograms.

Methods

A cross-sectional study was conducted on 834 women presenting for full-field digital screening mammography. The status of fat-infiltrated LNs was assessed based on the size and morphology of axillary LNs from screening mammograms. The prevalence of cardiometabolic disease was retrieved from the electronic medical records, including type 2 diabetes mellitus (T2DM), hypertension, dyslipidemia, high blood glucose, cardiovascular disease, stroke, and non-alcoholic fatty liver disease.

Results

Fat-infiltrated axillary LNs were associated with a high prevalence of T2DM among all women (adjusted odds ratio: 3.92, 95% CI: [2.40, 6.60], p-value < 0.001) and in subgroups of women with and without obesity. Utilizing the status of fatty LNs improved the classification of T2DM status in addition to age and BMI (1.4% improvement in the area under the receiver operating characteristic curve).

Conclusion

Fat-infiltrated axillary LNs visualized on screening mammograms were associated with the prevalence of T2DM. If further validated, fat-infiltrated axillary LNs may represent a novel imaging biomarker of T2DM in women undergoing screening mammography.

Gestational per- and polyfluoroalkyl substances exposure and infant body mass index trajectory in the New Hampshire Birth Cohort Study

Per- and polyfluoroalkyl substances (PFAS) are environmentally persistent, potential metabolic disruptors of concern for infants. Mothers participating in the New Hampshire Birth Cohort Study (NHBCS) provided a plasma sample during pregnancy to measure concentrations of seven PFAS, and infant weight and length were abstracted from well-child visits between birth and 12 months. Sex-specific growth patterns of child body mass index (BMI) were fit using a growth mixture model (GMM) and the relative risk ratios (RRR) and 95% Confidence Intervals (95% CI) for the association of maternal plasma PFAS with BMI growth patterns during infancy were estimated by using multinomial logistic model for the group probabilities in the GMM. Four growth patterns were identified: Group 1) a steep increase in BMI during the first 6 months, then a leveling off; Group 2) a gradual increase in BMI across the year; Group 3) a steep increase in BMI during months 1-3, then stable BMI; and Group 4) a gradual increase in BMI with plateau around 3 months (reference group). For boys, higher maternal pregnancy perfluorooctanoate concentrations were associated with a 60% decreased chance of being in group 3 as compared to group 4, after adjusting for potential confounding variables (RRR = 0.4; 95% CI: 0.1, 0.9). For girls, higher maternal perfluorooctane sulfonate (PFOS) concentrations during pregnancy were associated with a higher likelihood of following the growth pattern of groups 2 (RRR = 2.5; 95% CI: 1.0, 6.1) and 3 (RRR = 2.8; 95% CI: 1.0, 7.6) as compared to group 4, adjusting for potential confounding variables. In this cohort, sex-specific associations of maternal plasma PFAS concentrations during pregnancy with growth patterns during the first year of life were observed, with greater BMI growth observed among infant girls born to mothers with higher pregnancy concentrations of PFOS.

Associations of Neighborhood Opportunity and Social Vulnerability With Trajectories of Childhood Body Mass Index and Obesity Among US Children

Importance

Physical and social neighborhood attributes may have implications for children's growth and development patterns. The extent to which these attributes are associated with body mass index (BMI) trajectories and obesity risk from childhood to adolescence remains understudied.

Objective

To examine associations of neighborhood-level measures of opportunity and social vulnerability with trajectories of BMI and obesity risk from birth to adolescence.

Design, Setting, and Participants

This cohort study used data from 54 cohorts (20 677 children) participating in the Environmental Influences on Child Health Outcomes (ECHO) program from January 1, 1995, to January 1, 2022. Participant inclusion required at least 1 geocoded residential address and anthropometric measure (taken at the same time or after the address date) from birth through adolescence. Data were analyzed from February 1 to June 30, 2022.

Exposures

Census tract-level Child Opportunity Index (COI) and Social Vulnerability Index (SVI) linked to geocoded residential addresses at birth and in infancy (age range, 0.5-1.5 years), early childhood (age range, 2.0-4.8 years), and mid-childhood (age range, 5.0-9.8 years).

Main Outcomes and Measures

BMI (calculated as weight in kilograms divided by length [if aged <2 years] or height in meters squared) and obesity (age- and sex-specific BMI ≥95th percentile). Based on nationwide distributions of the COI and SVI, Census tract rankings were grouped into 5 categories: very low (<20th percentile), low (20th percentile to <40th percentile), moderate (40th percentile to <60th percentile), high (60th percentile to <80th percentile), or very high (≥80th percentile) opportunity (COI) or vulnerability (SVI).

Results

Among 20 677 children, 10 747 (52.0%) were male; 12 463 of 20 105 (62.0%) were White, and 16 036 of 20 333 (78.9%) were non-Hispanic. (Some data for race and ethnicity were missing.) Overall, 29.9% of children in the ECHO program resided in areas with the most advantageous characteristics. For example, at birth, 26.7% of children lived in areas with very high COI, and 25.3% lived in areas with very low SVI; in mid-childhood, 30.6% lived in areas with very high COI and 28.4% lived in areas with very low SVI. Linear mixed-effects models revealed that at every life stage, children who resided in areas with higher COI (vs very low COI) had lower mean BMI trajectories and lower risk of obesity from childhood to adolescence, independent of family sociodemographic and prenatal characteristics. For example, among children with obesity at age 10 years, the risk ratio was 0.21 (95% CI, 0.12-0.34) for very high COI at birth, 0.31 (95% CI, 0.20-0.51) for high COI at birth, 0.46 (95% CI, 0.28-0.74) for moderate COI at birth, and 0.53 (95% CI, 0.32-0.86) for low COI at birth. Similar patterns of findings were observed for children who resided in areas with lower SVI (vs very high SVI). For example, among children with obesity at age 10 years, the risk ratio was 0.17 (95% CI, 0.10-0.30) for very low SVI at birth, 0.20 (95% CI, 0.11-0.35) for low SVI at birth, 0.42 (95% CI, 0.24-0.75) for moderate SVI at birth, and 0.43 (95% CI, 0.24-0.76) for high SVI at birth. For both indices, effect estimates for mean BMI difference and obesity risk were larger at an older age of outcome measurement. In addition, exposure to COI or SVI at birth was associated with the most substantial difference in subsequent mean BMI and risk of obesity compared with exposure at later life stages.

Conclusions and Relevance

In this cohort study, residing in higher-opportunity and lower-vulnerability neighborhoods in early life, especially at birth, was associated with a lower mean BMI trajectory and a lower risk of obesity from childhood to adolescence. Future research should clarify whether initiatives or policies that alter specific components of neighborhood environment would be beneficial in preventing excess weight in children.

A meta-analysis of pre-pregnancy maternal body mass index and placental DNA methylation identifies 27 CpG sites with implications for mother-child health

Higher maternal pre-pregnancy body mass index (ppBMI) is associated with increased neonatal morbidity, as well as with pregnancy complications and metabolic outcomes in offspring later in life. The placenta is a key organ in fetal development and has been proposed to act as a mediator between the mother and different health outcomes in children. The overall aim of the present work is to investigate the association of ppBMI with epigenome-wide placental DNA methylation (DNAm) in 10 studies from the PACE consortium, amounting to 2631 mother-child pairs. We identify 27 CpG sites at which we observe placental DNAm variations of up to 2.0% per 10 ppBMI-unit. The CpGs that are differentially methylated in placenta do not overlap with CpGs identified in previous studies in cord blood DNAm related to ppBMI. Many of the identified CpGs are located in open sea regions, are often close to obesity-related genes such as GPX1 and LGR4 and altogether, are enriched in cancer and oxidative stress pathways. Our findings suggest that placental DNAm could be one of the mechanisms by which maternal obesity is associated with metabolic health outcomes in newborns and children, although further studies will be needed in order to corroborate these findings.

Inter-method reliability of silicone exposome wristbands and urinary biomarker assays in a pregnancy cohort

Silicone wristbands act as passive environmental samplers capable of detecting and measuring concentrations of a variety of chemicals. They offer a noninvasive method to collect complex exposure data in large-scale epidemiological studies. We evaluated the inter-method reliability of silicone wristbands and urinary biomarkers in the New Hampshire Birth Cohort Study. A subset of study participants (n = 96) provided a urine sample and wore a silicone wristband for 7 days at approximately 12 gestational weeks. Women were instructed to wear the wristbands during all their normal activities. Concentrations of urinary compounds and metabolites in the urine and parent compounds in wristbands were compared. High detection rates were observed for triphenyl phosphate (76.0%) and benzophenone (78.1%) in wristbands, although the distribution of corresponding urinary concentrations of chemicals did not differ according to whether chemicals were detected or not detected in wristbands. While detected among only 8.3% of wristbands, median urinary triclosan concentrations were higher among those with triclosan detected in wristbands (9.04 ng/mL) than without (0.16 ng/mL). For most chemicals slight to fair agreement was observed across exposure assessment methods, potentially due to low rates of detection in the wristbands for chemicals where observed urinary concentrations were relatively low as compared to background concentrations in the general population. Our findings support the growing body of research in support of deploying silicone wristbands as an important exposure assessment tool.

Contribution of gut bacteria to arsenic metabolism in the first year of life in a prospective birth cohort

Gut bacteria are at the interface of environmental exposures and their impact on human systems, and may alter host absorption, metabolism, and excretion of toxic chemicals. We investigated whether arsenic-metabolizing bacterial gene pathways related to urinary arsenic concentrations. In the New Hampshire Birth Cohort Study, urine and stool samples were obtained at six weeks (n = 186) and one year (n = 190) of age. Inorganic arsenic (iAs), monomethylarsonic acid (MMA), dimethylarsinic acid (DMA), and arsenobetaine (AsB) were quantified in infant urine samples using high-performance liquid chromatography with inductively coupled plasma mass spectrometry. Total arsenic exposure (tAs) was summarized as Σ(iAs, MMA, DMA) and log2-transformed. Fecal microbial DNA underwent metagenomic sequencing and the relative abundance of bacterial gene pathways were grouped as KEGG Orthologies (KOs) using BioBakery algorithms. Arsenic metabolism KOs with >80% prevalence were log2-transformed and modeled continuously using linear regression, those with <10% were not evaluated and those with 10-80% prevalence were analyzed dichotomously (detect/non-detect) using logistic regression. In the first set of models, tAs was regressed against KO relative abundance or detection adjusting for age at sample collection and child's sex. Effect modification by delivery mode was assessed in stratified models. In the second set of models, the association between the relative abundance/detection of the KOs and arsenic speciation (%iAs, %MMA, %DMA) was quantified with linear regression. Urinary tAs was associated with the increased relative abundance/detection odds of several arsenic-related KOs, including K16509, an arsenate reductase transcriptional regulator, with stronger associations among six-week-olds than one-year-olds. K16509 was also associated with decreased %MMA and increased %DMA at six weeks and one year. Notably, many associations were stronger among operatively-delivered than vaginally-delivered infants. Our findings suggest associations between arsenic-metabolizing bacteria in the infant gut microbiome and urinary arsenic excretion.

Dietary Exposure to Toxic Elements and the Health of Young Children: Methodological Considerations and Data Needs

Concerns have been raised regarding toxic-element (arsenic, cadmium, lead, and mercury) contamination of commercially available infant foods around the world. Young children are vulnerable to the effects of toxic elements, based on higher absorption levels and potentially poorer detoxification capacities. Toxic-element exposures in early life exact high societal costs, but it is unclear how much dietary exposure to these elements contributes to adverse health outcomes. Well-designed epidemiological studies conducted in different geographical and socioeconomic contexts need to estimate dietary toxicant exposure in young children and to determine whether causal links exist between toxicants in children's diets and health outcomes. This commentary outlines the methodological considerations and data needs to advance such research.

Exposure to melamine and its derivatives and aromatic amines among pregnant women in the United States: The ECHO Program

BACKGROUND

Melamine, melamine derivatives, and aromatic amines are nitrogen-containing compounds with known toxicity and widespread commercial uses. Nevertheless, biomonitoring of these chemicals is lacking, particularly during pregnancy, a period of increased susceptibility to adverse health effects.

OBJECTIVES

We aimed to measure melamine, melamine derivatives, and aromatic amine exposure in pregnant women across the United States (U.S.) and evaluate associations with participant and urine sample collection characteristics.

METHODS

We measured 43 analytes, representing 45 chemicals (i.e., melamine, three melamine derivatives, and 41 aromatic amines), in urine from pregnant women in nine diverse ECHO cohorts during 2008-2020 (N = 171). To assess relations with participant and urine sample collection characteristics, we used generalized estimating equations to estimate prevalence ratios (PRs) for analytes dichotomized at the detection limit, % differences (%Δ) for continuous analytes, and 95% confidence intervals. Multivariable models included age, race/ethnicity, marital status, urinary cotinine, and year of sample collection.

RESULTS

Twelve chemicals were detected in >60% of samples, with near ubiquitous detection of cyanuric acid, melamine, aniline, 4,4'-methylenedianiline, and a composite of o-toluidine and m-toluidine (99-100%). In multivariable adjusted models, most chemicals were associated with higher exposures among Hispanic and non-Hispanic Black participants. For example, concentrations of 3,4-dichloroaniline were higher among Hispanic (%Δ: +149, 95% CI: +17, +431) and non-Hispanic Black (%Δ: +136, 95% CI: +35, +311) women compared with non-Hispanic White women. We observed similar results for ammelide, o-/m-toluidine, 4,4'-methylenedianiline, and 4-chloroaniline. Most chemicals were positively associated with urinary cotinine, with strongest associations observed for o-/m-toluidine (%Δ: +23; 95% CI: +16, +31) and 3,4-dichloroaniline (%Δ: +25; 95% CI: +17, +33). Some chemicals exhibited annual trends (e.g., %Δ in melamine per year: -11; 95% CI: -19, -1) or time of day, seasonal, and geographic variability.

DISCUSSION

Exposure to melamine, cyanuric acid, and some aromatic amines was ubiquitous in this first investigation of these analytes in pregnant women. Future research should expand biomonitoring, identify sources of exposure disparities by race/ethnicity, and evaluate potential adverse health effects.

Adaptive-mixture-categorization (AMC)-based g-computation and its application to trace element mixtures and bladder cancer risk

Several new statistical methods have been developed to identify the overall impact of an exposure mixture on health outcomes. Weighted quantile sum (WQS) regression assigns the joint mixture effect weights to indicate the overall association of multiple exposures, and quantile-based g-computation is a generalized version of WQS without the restriction of directional homogeneity. This paper proposes an adaptive-mixture-categorization (AMC)-based g-computation approach that combines g-computation with an optimal exposure categorization search using the F statistic. AMC-based g-computation reduces variance within each category and retains the variance between categories to build more powerful predictors. In a simulation study, the performance of association analysis was improved using categorizing by AMC compared with quantiles. We applied this method to assess the association between a mixture of 12 trace element concentrations measured from toenails and the risk of non-muscle invasive bladder cancer. Our findings suggested that medium-level (116.7-145.5 μg/g) vs. low-level (39.5-116.2 μg/g) of toenail zinc had a statistically significant positive association with bladder cancer risk.

Estimating cumulative spatial risk over time with low-rank kriging multiple membership models

Many health outcomes result from accumulated exposures to one or more environmental factors. Accordingly, spatial risk studies have begun to consider multiple residential locations of participants, acknowledging that participants move and thus are exposed to environmental factors in several places. However, novel methods are needed to estimate cumulative spatial risk for disease while accounting for other risk factors. To this end, we propose a Bayesian model (LRK-MMM) that embeds a multiple membership model (MMM) into a low-rank kriging (LRK) model in order to estimate cumulative spatial risk at the point level while allowing for multiple residential locations per subject. The LRK approach offers a more computationally efficient means to analyze spatial risk in case-control study data at the point level compared with a Bayesian generalized additive model, and as increased precision in spatial risk estimates by analyzing point locations instead of administrative areas. Through a simulation study, we demonstrate the efficacy of the model and its improvement upon an existing multiple membership model that uses area-level spatial random effects to estimate risk. The results show that our proposed method provides greater spatial sensitivity (improvements ranging from 0.12 to 0.54) and power (improvements ranging from 0.02 to 0.94) to detect regions of elevated risk for disease across a range of exposure scenarios. Finally, we apply our model to case-control data from the New England bladder cancer study to estimate cumulative spatial risk while adjusting for many covariates.

Selenium-associated differentially expressed microRNAs and their targeted mRNAs across the placental genome in two U.S. birth cohorts

Selenium is an important micronutrient for foetal development. MicroRNAs play an important role in the function of the placenta, in communication between the placenta and maternal systems, and their expression can be altered through environmental and nutritional cues. To investigate the associations between placental selenium concentration and microRNA expression in the placenta, our observational study included 393 mother-child pairs from the New Hampshire Birth Cohort Study (NHBCS) and the Rhode Island Child Health Study (RICHS). Placental selenium concentrations were quantified using inductively coupled plasma mass spectrometry, and microRNA transcripts were measured using RNA-seq. We fit negative binomial additive models for assessing the association between selenium and microRNAs. We used the microRNA Data Integration Portal (mirDIP) to predict the target mRNAs of the differentially expressed microRNAs and verified the relationships between miRNA and mRNA targets in a subset of samples using existing whole transcriptome data (N = 199). We identified a non-monotonic association between selenium concentration and the expression of miR-216a-5p/miR-217-5p cluster (effective degrees of freedom, EDF = 2.44 and 2.08; FDR = 3.08 × 10-5) in placenta. Thirty putative target mRNAs of miR-216a-5p and/or miR-217-5p were identified computationally and empirically and were enriched in selenium metabolic pathways (driven by selenoprotein coding genes, TXNRD2 and SELENON). Our findings suggest that selenium influences placental microRNA expression. Further, miR-216a-5p and its putative target mRNAs could be the potential mechanistic targets of the health effect of selenium.

Arsenic exposure and respiratory outcomes during childhood in the INMA study

Ingested inorganic arsenic (iAs) is a human carcinogen that is also linked to other adverse health effects, such as respiratory outcomes. Yet, among populations consuming low-arsenic drinking water, the impact of iAs exposure on childhood respiratory health is still uncertain. For a Spanish child study cohort (INfancia y Medio Ambiente—INMA), low-arsenic drinking water is usually available and ingestion of iAs from food is considered the major source of exposure. Here, we explored the association between iAs exposure and children’s respiratory outcomes assessed at 4 and 7 years of age (n = 400). The summation of 4-year-old children’s urinary iAs, monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) was used as a biomarker of iAs exposure (∑As) (median of 4.92 μg/L). Children’s occurrence of asthma, eczema, sneeze, wheeze, and medication for asthma and wheeze at each assessment time point (i.e., 4- and 7-year) was assessed with maternal interviewer-led questionnaires. Crude and adjusted Poisson regression models using Generalized Estimating Equation (GEE) were performed to account for the association between natural logarithm transformed (ln) urinary ∑As in μg/L at 4 years and repeated assessments of respiratory symptoms at 4 and 7 years of age. The covariates included in the models were child sex, maternal smoking status, maternal level of education, sub-cohort, and children’s consumption of vegetables, fruits, and fish/seafood. The GEE—splines function using Poisson regression showed an increased trend of the overall expected counts of respiratory symptoms with high urinary ∑As. The adjusted expected counts (95% confidence intervals) at ln-transformed urinary ∑As 1.57 (average concentration) and 4.00 (99^th percentile concentration) were 0.63 (0.36, 1.10) and 1.33 (0.61, 2.89), respectively. These exploratory findings suggest that even relatively low-iAs exposure levels, relevant to the Spanish and other populations, may relate to an increased number of respiratory symptoms during childhood.

Racial/ethnic disparities in subjective sleep duration, sleep quality, and sleep disturbances during pregnancy: an ECHO study

In the United States, racial/ethnic minoritized groups experience worse sleep than non-Hispanic Whites (nHW), but less is known about pregnant people. This is a key consideration since poor sleep during pregnancy is common and associated with increased risk of adverse perinatal outcomes. This study reports the prevalence of subjective sleep measures in a multi-racial/ethnic pregnant population from the Environmental influences on Child Health Outcomes (ECHO) program. Participants' self-reported race and ethnicity were grouped into: nHW, non-Hispanic Black/African American (nHB/AA), Hispanic, non-Hispanic Asian (nHA). Analyses examined trimester-specific (first (T1), second (T2), third (T3)) nocturnal sleep duration, quality, and disturbances (Pittsburgh Sleep Quality Index and ECHO maternal sleep health questionnaire). Linear or multinomial regressions estimated the associations between race/ethnicity and each sleep domain by trimester, controlling for body mass index and age, with nHW as reference group. We repeated analyses within maternal education strata. nHB/AA participants reported shorter sleep duration (T2: β = -0.55 [-0.80,-0.31]; T3: β = -0.65 [-0.99,-0.31]) and more sleep disturbances (T2: β = 1.92 [1.09,2.75]; T3: β = 1.41 [0.09,2.74]). Hispanic participants reported longer sleep duration (T1: β = 0.22 [0.00004,0.44]; T2: β = 0.61 [0.47,0.76]; T3: β = 0.46 [0.22,0.70]), better sleep quality (Reference group: Very good. Fairly good T1: OR = 0.48 [0.32,0.73], T2: OR = 0.36 [0.26,0.48], T3: OR = 0.31 [0.18,0.52]. Fairly bad T1: OR = 0.27 [0.16,0.44], T2: OR = 0.46 [0.31, 0.67], T3: OR = 0.31 [0.17,0.55]), and fewer sleep disturbances (T2: β = -0.5 [-1.0,-0.12]; T3: β = -1.21 [-2.07,-0.35]). Differences persisted within the high-SES subsample. Given the stark racial/ethnic disparities in perinatal outcomes and their associations with sleep health, further research is warranted to investigate the determinants of these disparities.

Changes in BMI During the COVID-19 Pandemic

BACKGROUND AND OBJECTIVES

Experts hypothesized increased weight gain in children associated with the coronavirus disease 2019 (COVID-19) pandemic. Our objective was to evaluate whether the rate of change of child body mass index (BMI) increased during the COVID-19 pandemic compared with prepandemic years.

METHODS

The study population of 1996 children ages 2 to 19 years with at least 1 BMI measure before and during the COVID-19 pandemic was drawn from 38 pediatric cohorts across the United States participating in the Environmental Influences on Child Health Outcomes-wide cohort study. We modeled change in BMI using linear mixed models, adjusting for age, sex, race, ethnicity, maternal education, income, baseline BMI category, and type of BMI measure. Data collection and analysis were approved by the local institutional review board of each institution or by the central Environmental Influences on Child Health Outcomes institutional review board.

RESULTS

BMI increased during the COVID-19 pandemic compared with previous years (0.24 higher annual gain in BMI during the pandemic compared with previous years, 95% confidence interval 0.02 to 0.45). Children with BMI in the obese range compared with the healthy weight range were at higher risk for excess BMI gain during the pandemic, whereas children in higher-income households were at decreased risk of BMI gain.

CONCLUSIONS

One effect of the COVID-19 pandemic is an increase in annual BMI gain during the COVID-19 pandemic compared with the 3 previous years among children in our national cohort. This increased risk among US children may worsen a critical threat to public health and health equity.

Associations of maternal urinary arsenic concentrations during pregnancy with childhood cognitive abilities: The HOME study

Arsenic exposure during pregnancy may increase the risk for intellectual deficits in children, but limited data exist from prospective epidemiologic studies, particularly at low arsenic exposure levels. We investigated the association between prenatal maternal urinary arsenic concentrations and childhood cognitive abilities in the Health Outcomes and Measures of the Environment (HOME) Study. We used anion exchange chromatography coupled with inductively coupled plasma mass spectrometry detection to measure arsenic species content in pregnant women's urine. The summation of inorganic arsenic (iAs), monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) refers to ∑As. We assessed children's cognitive function (n = 260) longitudinally at 1-, 2-, and 3-years using Bayley Scales of Infant and Toddler Development, at 5 years using Wechsler Preschool and Primary Scale of Intelligence, and at 8 years using Wechsler Intelligence Scale for Children. We observed a modest decrease in mental development index and full-scale intelligence quotient at ages 3 and 5 years with each doubling of ∑As with estimated score (ß) differences and 95% confidence interval (CI) of -1.8 from -4.1 to 0.5 and -2.5 from -5.1 to 0.0, respectively. This trend was stronger and reached statistical significance among children whose mothers had lower iAs methylation capacity and low urinary arsenobetaine concentrations. Our findings suggest that arsenic exposure levels relevant to the general US population may affect children's cognitive abilities.

Sex-specific relationships of the infant microbiome and early-childhood behavioral outcomes

BACKGROUND

A link between the gut microbiome and behavior is hypothesized, but most previous studies are cross-sectional or in animal models. The modifying role of host sex is poorly characterized. We aimed to identify sex-specific prospective associations between the early-life gut microbiome and preschool-age neurobehavior.

METHODS

In a prospective cohort, gut microbiome diversity and taxa were estimated with 16S rRNA sequencing at 6 weeks, 1 year, and 2 years. Species and gene pathways were inferred from metagenomic sequencing at 6 weeks and 1 year. When subjects were 3 years old, parents completed the Behavioral Assessment System for Children, second edition (BASC-2). A total of 260 children contributed 523 16S rRNA and 234 metagenomics samples to analysis. Models adjusted for sociodemographic characteristics.

RESULTS

Higher diversity at 6 weeks was associated with better internalizing problems among boys, but not girls [β_Boys = -1.86 points/SD Shannon diversity; 95% CI (-3.29, -0.42), p_Boys = 0.01, β_Girls = 0.22 (-1.43, 1.87), p_Girls = 0.8, p_interaction = 0.06]. Among other taxa-specific associations, Bifidobacterium at 6 weeks was associated with Adaptive Skills scores in a sex-specific manner. We observed relationships between functional features and BASC-2 scores, including vitamin B6 biosynthesis pathways and better Depression scores.

CONCLUSIONS

This study advances our understanding of microbe-host interactions with implications for childhood behavioral health.

IMPACT

This is one of the first studies to examine the early-life microbiome and neurobehavior, and the first to examine prospective sex-specific associations. Infant and early-childhood microbiomes relate to neurobehavior including anxiety, depression, hyperactivity, and social behaviors in a time- and sex-specific manner. Our findings suggest future studies should evaluate whether host sex impacts the relationship between the gut microbiome and behavioral health outcomes.

Metagenomic analysis reveals associations between salivary microbiota and body composition in early childhood

Several studies have shown that body mass index is strongly associated with differences in gut microbiota, but the relationship between body weight and oral microbiota is less clear especially in young children. We aimed to evaluate if there is an association between child growth and the saliva microbiome. We hypothesized that associations between growth and the saliva microbiome would be moderate, similarly to the association between growth and the gut microbiome. For 236 toddlers participating in the New Hampshire Birth Cohort Study, we characterized the association between multiple longitudinal anthropometric measures of body height, body weight and body mass. Body Mass Index (BMI) z-scores were calculated, and dual-energy x-ray absorptiometry (DXA) was used to estimate body composition. Shotgun metagenomic sequencing of saliva samples was performed to taxonomically and functionally profile the oral microbiome. We found that within-sample diversity was inversely related to body mass measurements while community composition was not associated. Although the magnitude of associations were small, some taxa were consistently associated with growth and modified by sex. Certain taxa were associated with decreased weight or growth (including Actinomyces odontolyticus and Prevotella melaninogenica) or increased growth (such as Streptococcus mitis and Corynebacterium matruchotii) across anthropometric measures. Further exploration of the functional significance of this relationship will enhance our understanding of the intersection between weight gain, microbiota, and energy metabolism and the potential role of these relationships on the onset of obesity-associated diseases in later life.

Health Behavior Changes during the COVID-19 Pandemic: A Longitudinal Analysis among Children

This longitudinal study compared children’s health behaviors before the COVID-19 pandemic versus during the pandemic. This analysis examined the association between individual-level characteristics and health behavior change. Four prospective cohort studies in the Environmental influences on Child Health Outcomes (ECHO) Program contributed data. Children aged 4−12 years and their caregivers were recruited in California, Colorado, North Dakota, and New Hampshire. Dietary intake, physical activity, screen time, and sleep duration were assessed with questionnaires pre-pandemic and during the pandemic. The final sample included 347 children: 47% female and 62% non-Hispanic White. Compared with pre-pandemic, weekday screen time duration was higher during the pandemic (3.0 vs. 4.5 h, p < 0.001). Unadjusted increases in screen time duration differed by race and ethnicity: 1.3 h/day for non-Hispanic White children, 2.3 h/day for Hispanic children, and 5.3 h/day for non-Hispanic Black children. Overall, no changes occurred in sugar-sweetened beverage (SSB) intake (p = 0.26), discretionary food intake (p = 0.93), and physical activity (p = 0.15). Sleep duration increased by 30 min among children who did not meet sleep recommendations pre-pandemic. Child sex and maternal education level were not associated with health behavior change. The pandemic may have exacerbated disparities in some health behaviors. Families may need support to re-establish healthy routines.

Information enhanced model selection for Gaussian graphical model with application to metabolomic data

In light of the low signal-to-noise nature of many large biological data sets, we propose a novel method to learn the structure of association networks using Gaussian graphical models combined with prior knowledge. Our strategy includes two parts. In the first part, we propose a model selection criterion called structural Bayesian information criterion, in which the prior structure is modeled and incorporated into Bayesian information criterion. It is shown that the popular extended Bayesian information criterion is a special case of structural Bayesian information criterion. In the second part, we propose a two-step algorithm to construct the candidate model pool. The algorithm is data-driven and the prior structure is embedded into the candidate model automatically. Theoretical investigation shows that under some mild conditions structural Bayesian information criterion is a consistent model selection criterion for high-dimensional Gaussian graphical model. Simulation studies validate the superiority of the proposed algorithm over the existing ones and show the robustness to the model misspecification. Application to relative concentration data from infant feces collected from subjects enrolled in a large molecular epidemiological cohort study validates that metabolic pathway involvement is a statistically significant factor for the conditional dependence between metabolites. Furthermore, new relationships among metabolites are discovered which can not be identified by the conventional methods of pathway analysis. Some of them have been widely recognized in biological literature.

The relationship between the gut microbiome and the risk of respiratory infections among newborns

Background

Emerging evidence points to a critical role of the developing gut microbiome in immune maturation and infant health; however, prospective studies are lacking.

Methods

We examined the occurrence of infections and associated symptoms during the first year of life in relation to the infant gut microbiome at six weeks of age using bacterial 16S rRNA V4-V5 gene sequencing (N = 465) and shotgun metagenomics (N = 185). We used generalized estimating equations to assess the associations between longitudinal outcomes and 16S alpha diversity and metagenomics species.

Results

Here we show higher infant gut microbiota alpha diversity was associated with an increased risk of infections or respiratory symptoms treated with a prescription medicine, and specifically upper respiratory tract infections. Among vaginally delivered infants, a higher alpha diversity was associated with an increased risk of all-cause wheezing treated with a prescription medicine and diarrhea involving a visit to a health care provider. Positive associations were specifically observed with Veillonella species among all deliveries and Haemophilus influenzae among cesarean-delivered infants.

Conclusion

Our findings suggest that intestinal microbial diversity and the relative abundance of key taxa in early infancy may influence susceptibility to respiratory infection, wheezing, and diarrhea.

Diet quality, common genetic polymorphisms, and bladder cancer risk in a New England population-based study

PURPOSE

We examined the interaction between common genetic bladder cancer variants, diet quality, and bladder cancer risk in a population-based case-control study conducted in New England.

METHODS

At the time of enrollment, 806 bladder cancer cases and 974 controls provided a DNA sample and completed a diet history questionnaire. Diet quality was assessed using the 2010 Alternate Healthy Eating Index (AHEI-2010) score. Single nucleotide polymorphisms (SNPs) reported in genome-wide association studies to be associated with bladder cancer risk were combined into a polygenic risk score and also examined individually for interaction with the AHEI-2010. Adjusted odds ratios (OR) and 95% confidence intervals (CI) were calculated using logistic regression.

RESULTS

A 1-standard deviation increase in polygenic risk score was associated with higher bladder cancer risk (OR, 1.34; 95% CI 1.21-1.49). Adherence to the AHEI-2010 was not associated with bladder cancer risk (OR, 0.99; 95% CI 0.98-1.00) and the polygenic risk score did not appear to modify the association between the AHEI-2010 and bladder cancer risk. In single-SNP analyses, rs8102137 (bladder cancer risk allele, C) modified the association between the AHEI-2010 total score and bladder cancer risk, with the strongest evidence for the AHEI-2010 long chain fat guideline (OR for TT, 0.92; 95% CI 0.87-0.98; OR for CT, 1.02; 95% CI 0.96-1.08; OR for CC, 1.03; 95% CI 0.93-1.14; p for interaction, 0.02).

CONCLUSIONS

In conclusion, rs8102137 near the cyclin E1 gene ( CCNE1 ) may be involved in gene-diet interactions for bladder cancer risk.

Dietary Exposure to Essential and Non-essential Elements During Infants' First Year of Life in the New Hampshire Birth Cohort Study

Even the low levels of non-essential elements exposure common in the US may have health consequences especially early in life. However, little is known about the infant's dynamic exposure to essential and non-essential elements. This study aims to evaluate exposure to essential and non-essential elements during infants' first year of life and to explore the association between the exposure and rice consumption. Paired urine samples from infants enrolled in the New Hampshire Birth Cohort Study (NHBCS) were collected at approximately 6 weeks (exclusively breastfed) and at 1 year of age after weaning (n = 187). A further independent subgroup of NHBCS infants with details about rice consumption at 1 year of age also was included (n = 147). Urinary concentrations of 8 essential (Co, Cr, Cu, Fe, Mn, Mo, Ni, and Se) and 9 non-essential (Al, As, Cd, Hg, Pb, Sb, Sn, V, and U) elements were determined as a measure of exposure. Several essential (Co, Fe, Mo, Ni, and Se) and non-essential (Al, As, Cd, Hg, Pb, Sb, Sn, and V) elements had higher concentrations at 1 year than at 6 weeks of age. The highest increases were for urinary As and Mo with median concentrations of 0.20 and 1.02 µg/L at 6 weeks and 2.31 and 45.36 µg/L at 1 year of age, respectively. At 1 year of age, As and Mo urine concentrations were related to rice consumption. Further efforts are necessary to minimize exposure to non-essential elements while retaining essential elements to protect and promote children's health.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12403-022-00489-x.

MarZIC: A Marginal Mediation Model for Zero-Inflated Compositional Mediators with Applications to Microbiome Data

BACKGROUND

The human microbiome can contribute to pathogeneses of many complex diseases by mediating disease-leading causal pathways. However, standard mediation analysis methods are not adequate to analyze the microbiome as a mediator due to the excessive number of zero-valued sequencing reads in the data and that the relative abundances have to sum to one. The two main challenges raised by the zero-inflated data structure are: (a) disentangling the mediation effect induced by the point mass at zero; and (b) identifying the observed zero-valued data points that are not zero (i.e., false zeros).

METHODS

We develop a novel marginal mediation analysis method under the potential-outcomes framework to address the issues. We also show that the marginal model can account for the compositional structure of microbiome data.

RESULTS

The mediation effect can be decomposed into two components that are inherent to the two-part nature of zero-inflated distributions. With probabilistic models to account for observing zeros, we also address the challenge with false zeros. A comprehensive simulation study and the application in a real microbiome study showcase our approach in comparison with existing approaches.

CONCLUSIONS

When analyzing the zero-inflated microbiome composition as the mediators, MarZIC approach has better performance than standard causal mediation analysis approaches and existing competing approach.