Extracellular Vesicle-Encapsulated microRNAs as Novel Biomarkers of Lung Health.

Rationale: Early detection of respiratory diseases is critical to facilitate delivery of disease-modifying interventions. Extracellular vesicle-enriched microRNAs (EV-miRNAs) may represent reliable markers of early lung injury. Objectives: Evaluate associations of plasma EV-miRNAs with lung function. Methods: The prospective NAS (Normative Aging Study) collected plasma EV-miRNA measurements from 1996-2015 and spirometry every 3-5 years through 2019. Associations of EV-miRNAs with baseline lung function were modeled using linear regression. To complement the individual miRNA approach, unsupervised machine learning was used to identify clusters of participants with distinct EV-miRNA profiles. Associations of EV-miRNA profiles with multivariate latent longitudinal lung function trajectories were modeled using log binomial regression. Biological functions of significant EV-miRNAs were explored using pathway analyses. Results were replicated in an independent sample of NAS participants and in the HEALS (Health Effects of Arsenic Longitudinal Study). Measurements and Main Results: In the main cohort of 656 participants, 51 plasma EV-miRNAs were associated with baseline lung function (false discovery rate-adjusted P value < 0.05), 28 of which were replicated in the independent NAS sample and/or in the HEALS cohort. A subset of participants with distinct EV-miRNA expression patterns had increased risk of declining lung function over time, which was replicated in the independent NAS sample. Significant EV-miRNAs were shown in pathway analyses to target biological pathways that regulate respiratory cellular immunity, the lung inflammatory response, and airway structural integrity. Conclusions: Plasma EV-miRNAs may represent a robust biomarker of subclinical lung injury and may facilitate early identification and treatment of patients at risk of developing overt lung disease.

Prenatal lead exposure, telomere length in cord blood, and DNA methylation age in the PROGRESS prenatal cohort.

BACKGROUND: Lead is a ubiquitous pollutant with deleterious effects on human health and remains a major current public health concern in developing countries. This heavy metal may interfere with nucleic acids via oxidative stress or epigenetic changes that affect biological markers of aging, e.g., telomere length and DNA methylation (DNAm). Telomere shortening associates with biological age in newborns, and DNA methylation at specific CpG sites can be used to calculate "epigenetic clocks". OBJECTIVE: The aim of this study was to examine the associations of prenatal lead exposures with telomere length and DNA-methylation-based predictors of age in cord blood. DESIGN: The study included 507 mother-child pairs from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) study, a birth cohort in Mexico City. Maternal blood (second trimester, third trimester and at delivery) and bone lead levels (one month postpartum) were measured using inductively coupled plasma-mass spectrometry and X-ray fluorescence, respectively. Cord blood leukocyte telomere length was measured using quantitative PCR and apparent age by DNA methylation biomarkers, i.e., Horvath's DNA methylation age and the Knight's predictor of gestational age. RESULTS: Average maternal age was 28.5 ± 5.5 years, and 51.5% reported low socioeconomic status. Children's mean telomere length was 1.2 ± 1.3 relative units, and mean DNA methylation ages using the Horvath's and Knight's clocks were -2.6 ± 0.1 years and 37.9 ± 1.4 weeks (mean ± SD), respectively. No significant associations were found between maternal blood and bone lead concentrations with telomere length and DNAm age in newborns. CONCLUSION: We found no associations of prenatal lead exposure with telomere length and DNA methylation age biomarkers.

Public drinking water contaminant estimates for birth cohorts in the Environmental Influences on Child Health Outcomes (ECHO) Cohort.

BACKGROUND: The United States Environmental Protection Agency (USEPA) regulates over 80 contaminants in community water systems (CWS), including those relevant to infant health outcomes. Multi-cohort analyses of the association between measured prenatal public water contaminant concentrations and infant health outcomes are sparse in the US. OBJECTIVE: Our objectives were to (1) develop Zip Code Tabulation Area (ZCTA)-level CWS contaminant concentrations for participants in the Environmental Influences on Child Health Outcomes (ECHO) Cohort and (2) evaluate regional, seasonal, and sociodemographic inequities in contaminant concentrations at the ZCTA-level. The ECHO Cohort harmonizes data from over 69 extant pregnancy and pediatric cohorts across the US. METHODS: We used CWS estimates derived from the USEPA's Six-Year Review 3 (2006-2011) to develop population-weighted, average concentrations for 10 contaminants across 7640 ZCTAs relevant to the ECHO Cohort. We evaluated contaminant distributions, exceedances of regulatory thresholds, and geometric mean ratios (with corresponding percent changes) associated with ZCTA sociodemographic characteristics via spatial lag linear regression models. RESULTS: We observed significant regional variability in contaminant concentrations across the US. ZCTAs were most likely to exceed the maximum contaminant level for arsenic (n = 100, 1.4%) and the health-protective threshold for total trihalomethanes (n = 3584, 64.0%). A 10% higher proportion of residents who were American Indian/Alaskan Native and Hispanic/Latino was associated with higher arsenic (11%, 95% CI: 7%, 15%; and 2%, 95% CI: 0%, 3%, respectively) and uranium (15%, 95% CI: 10%, 21%; and 9%, 95% CI: 6%, 12%, respectively) concentrations. IMPACT: Nationwide epidemiologic analyses evaluating the association between US community water system contaminant concentration estimates and associated adverse birth outcomes in cohort studies are sparse because public water contaminant concentration estimates that can be readily linked to participant addresses are not available. We developed Zip Code Tabulation Area (ZCTA)-level CWS contaminant concentrations that can be linked to participants in the Environmental Influences on Child Health Outcomes (ECHO) Cohort and evaluated regional, seasonal, and sociodemographic inequities in contaminant concentrations for these ZCTAs. Future epidemiologic studies can leverage these CWS exposure estimates in the ECHO Cohort to evaluate associations with relevant infant outcomes.

Inverse associations of cord blood mitochondrial DNA copy number with childhood adiposity.

OBJECTIVE: The objective of this study was to examine associations between umbilical cord mitochondrial DNA copy number (mtDNAcn) and adiposity across childhood. METHODS: In a prospective birth cohort of Dominican and African American children from New York City, New York (1998-2006), mtDNAcn was measured in cord blood. Children (N = 336) were evaluated for their height, weight, and bioimpedance at age 5, 7, 9, and 11 years. We used linear mixed-effects models to assess associations of mtDNAcn tertiles in cord blood with child BMI, BMI z scores, fat mass index, and body fat percentage. Latent class growth models and interactions between mtDNAcn and child age or child age2 were used to assess associations between age and adiposity trajectories. RESULTS: BMI was, on average, 1.5 kg/m2 higher (95% CI: 0.58, 2.5) in individuals with mtDNAcn in the low- compared with the middle-mtDNAcn tertile. Results were similar for BMI z score, fat mass index, and body fat percentage. Moreover, children in the low-mtDNAcn group had increased odds of being in an "increasing" or "high-stable" adiposity class. CONCLUSIONS: Lower mtDNAcn at birth may predict greater childhood adiposity, highlighting the potential key role of perinatal mitochondrial function in adiposity during development.

Mitochondrial DNA copy number dynamics and associations with the prenatal environment from birth through adolescence in a population of Dominican and African American children.

Mitochondrial DNA copy number (mtDNAcn) dynamics throughout childhood are poorly understood. We profiled mtDNAcn from birth through adolescence and evaluated how the prenatal environment influences mtDNAcn across childhood. Data were collected from children from New York City followed through 18 years. Using duplexed qRT-PCR, we quantified mtDNAcn relative to nuclear DNA in blood collected from the umbilical cord (n = 450), children aged 5-7 (n = 510), and adolescents aged 15-18 (n = 278). We examined mtDNAcn across childhood with linear mixed-effects models (LMM). Relative mtDNAcn was lowest at birth (mean ± SD: 0.67 ± 0.35) and increased in childhood (1.24 ± 0.50) then slightly declined in adolescence (1.13 ± 0.44). We observed no differences in mtDNAcn by sex or race/ethnicity. mtDNAcn was positively associated with prenatal environmental tobacco smoke exposure (0.077 [ 0.01, 0.14] change in relative mtDNAcn) but negatively associated with maternal completion of high school (-0.066 [-0.13, 0.00]), with the receipt of public assistance at birth (-0.074 [-0.14, -0.01]), and when mother born outside the U.S (-0.061 [-0.13, 0.003]). Infant birth outcomes were not associated with mtDNAcn. MtDNAcn levels were dynamic through childhood and associated with some prenatal factors, underscoring the need for the investigation of longitudinal mtDNAcn for human health research.

Longitudinal transcriptomic analysis of mouse sciatic nerve reveals pathways associated with age-related muscle pathology.

BACKGROUND: Sarcopenia, the age-associated decline in skeletal muscle mass and strength, has long been considered a disease of muscle only, but accumulating evidence suggests that sarcopenia could originate from the neural components controlling muscles. To identify early molecular changes in nerves that may drive sarcopenia initiation, we performed a longitudinal transcriptomic analysis of the sciatic nerve, which governs lower limb muscles, in aging mice. METHODS: Sciatic nerve and gastrocnemius muscle were obtained from female C57BL/6JN mice aged 5, 18, 21 and 24 months old (n = 6 per age group). Sciatic nerve RNA was extracted and underwent RNA sequencing (RNA-seq). Differentially expressed genes (DEGs) were validated using quantitative reverse transcription PCR (qRT-PCR). Functional enrichment analysis of clusters of genes associated with patterns of gene expression across age groups (adjusted P-value < 0.05, likelihood ratio test [LRT]) was performed. Pathological skeletal muscle aging was confirmed between 21 and 24 months by a combination of molecular and pathological biomarkers. Myofiber denervation was confirmed with qRT-PCR of Chrnd, Chrng, Myog, Runx1 and Gadd45ɑ in gastrocnemius muscle. Changes in muscle mass, cross-sectional myofiber size and percentage of fibres with centralized nuclei were analysed in a separate cohort of mice from the same colony (n = 4-6 per age group). RESULTS: We detected 51 significant DEGs in sciatic nerve of 18-month-old mice compared with 5-month-old mice (absolute value of fold change > 2; false discovery rate [FDR] < 0.05). Up-regulated DEGs included Dbp (log2 fold change [LFC] = 2.63, FDR < 0.001) and Lmod2 (LFC = 7.52, FDR = 0.001). Down-regulated DEGs included Cdh6 (LFC = -21.38, FDR < 0.001) and Gbp1 (LFC = -21.78, FDR < 0.001). We validated RNA-seq findings with qRT-PCR of various up- and down-regulated genes including Dbp and Cdh6. Up-regulated genes (FDR < 0.1) were associated with the AMP-activated protein kinase signalling pathway (FDR = 0.02) and circadian rhythm (FDR = 0.02), whereas down-regulated DEGs were associated with biosynthesis and metabolic pathways (FDR < 0.05). We identified seven significant clusters of genes (FDR < 0.05, LRT) with similar expression patterns across groups. Functional enrichment analysis of these clusters revealed biological processes that may be implicated in age-related changes in skeletal muscles and/or sarcopenia initiation including extracellular matrix organization and an immune response (FDR < 0.05). CONCLUSIONS: Gene expression changes in mouse peripheral nerve were detected prior to disturbances in myofiber innervation and sarcopenia onset. These early molecular changes we report shed a new light on biological processes that may be implicated in sarcopenia initiation and pathogenesis. Future studies are warranted to confirm the disease modifying and/or biomarker potential of the key changes we report here.

Breast milk-derived extracellular vesicle miRNAs are associated with maternal asthma and atopy.

Background: Breast milk-derived extracellular vesicle (EV) miRNAs may program child health outcomes associated with maternal asthma and atopy. The authors investigated associations between maternal asthma/atopy and EV miRNAs in the Programming of Intergenerational Stress Mechanisms cohort. Methods: Breast milk-derived EV miRNAs collected 6.1 ± 5.9 weeks postnatally (n = 80 mothers) were profiled using the TaqMan OpenArray Human MicroRNA Panel. The authors assessed associations using adjusted robust regression. Results: Nine EV miRNAs were associated with asthma during pregnancy (a priori criteria: nominal p < 0.05; |Bregression| >0.2). miR-1290 was associated with asthma and atopy during pregnancy (p < 0.05; |Bregression| >0.2). Enriched Kyoto Encyclopedia of Genes and Genomes pathways included TGF-ß signaling and extracellular matrix-receptor interaction (false discovery rate <0.05). Conclusion: In this study, maternal asthma and atopy were associated with breast milk-derived EV miRNAs. Additional studies are needed to understand whether EV miRNAs have direct effects on infant and child health.

Maternal asthma is associated with child health outcomes, although the biological mechanisms involved are not fully understood. miRNAs are small molecules involved in regulating gene expression. miRNAs packaged into membrane-bound particles called extracellular vesicles (EVs) are present in human breast milk and may pass from mother to infant to signal which genes to translate into proteins. This study investigated the extent to which maternal asthma and atopy influenced levels of 130 EV miRNAs measured in breast milk. Nine EV miRNAs were associated with maternal asthma during pregnancy, and one EV miRNA was associated with maternal atopy. miRNAs associated with asthma target genes in pathways related to asthma; however, future research is needed to determine whether changes in breast milk-derived EV miRNAs impact child health.

Associations of Childhood and Perinatal Blood Metals with Children's Gut Microbiomes in a Canadian Gestation Cohort.

BACKGROUND: The gut microbiome is important in modulating health in childhood. Metal exposures affect multiple health outcomes, but their ability to modify bacterial communities in children is poorly understood. OBJECTIVES: We assessed the associations of childhood and perinatal blood metal levels with childhood gut microbiome diversity, structure, species, gene family-inferred species, and potential pathway alterations. METHODS: We assessed the gut microbiome using 16S rRNA gene amplicon sequencing and shotgun metagenomic sequencing in stools collected from 6- to 7-year-old children participating in the GESTation and Environment (GESTE) cohort study. We assessed blood metal concentrations [cadmium (Cd), manganese (Mn), mercury (Hg), lead (Pb), selenium (Se)] at two time points, namely, perinatal exposures at delivery (N=70) and childhood exposures at the 6- to 7-y follow-up (N=68). We used multiple covariate-adjusted statistical models to determine microbiome associations with continuous blood metal levels, including linear regression (Shannon and Pielou alpha diversity indexes), permutational multivariate analysis of variance (adonis; beta diversity distance matrices), and multivariable association model (MaAsLin2; phylum, family, species, gene family-inferred species, and pathways). RESULTS: Children's blood Mn and Se significantly associated with microbiome phylum [e.g., Verrucomicrobiota (coef=-0.305, q=0.031; coef=0.262, q=0.084, respectively)] and children's blood Mn significantly associated with family [e.g., Eggerthellaceae (coef=-0.228, q=0.052)]-level differences. Higher relative abundance of potential pathogens (e.g., Flavonifractor plautii), beneficial species (e.g., Bifidobacterium longum, Faecalibacterium prausnitzii), and both potentially pathogenic and beneficial species (e.g., Bacteriodes vulgatus, Eubacterium rectale) inferred from gene families were associated with higher childhood or perinatal blood Cd, Hg, and Pb (q<0.1). We found significant negative associations between childhood blood Pb and acetylene degradation pathway abundance (q<0.1). Finally, neither perinatal nor childhood metal concentrations were associated with children's gut microbial inter- and intrasubject diversity. DISCUSSION: Our findings suggest both long- and short-term associations between metal exposure and the childhood gut microbiome, with stronger associations observed with more recent exposure. Future epidemiologic analyses may elucidate whether the observed changes in the microbiome relate to children's health. https://doi.org/10.1289/EHP9674.

Association of Prenatal Acetaminophen Exposure Measured in Meconium With Adverse Birth Outcomes in a Canadian Birth Cohort.

Background: The small number of studies examining the association of prenatal acetaminophen with birth outcomes have all relied on maternal self-report. It remains unknown whether prenatal acetaminophen exposure measured in a biological specimen is associated with birth outcomes. Objectives: To investigate the association of acetaminophen measured in meconium with birthweight, gestational age, preterm birth, size for gestational age, gestational diabetes, preeclampsia, and high blood pressure. Methods: This birth cohort from Sherbrooke, QC, Canada, included 773 live births. Mothers with no thyroid disease enrolled at their first prenatal care visit or delivery. Acetaminophen was measured in meconium for 393 children at delivery. We tested associations of prenatal acetaminophen with birthweight, preterm birth, gestational age, small and large for gestational age, gestational diabetes, preeclampsia, and high blood pressure. We imputed missing data via multiple imputation and used inverse probability weighting to account for confounding and selection bias. Results: Acetaminophen was detected in 222 meconium samples (56.5%). Prenatal acetaminophen exposure was associated with decreased birthweight by 136 g (ß = -136; 95% CI [-229, -43]), 20% increased weekly hazard of delivery (hazard ratio = 1.20; 95% CI [1.00, 1.43]), and over 60% decreased odds of being born large for gestational age (odds ratio = 0.38; 95% CI [0.20, 0.75]). Prenatal acetaminophen was not associated with small for gestational age, preterm birth, or any pregnancy complications. Conclusion: Prenatal acetaminophen was associated with adverse birth outcomes. Although unobserved confounding and confounding by indication are possible, these results warrant further investigation into adverse perinatal effects of prenatal acetaminophen exposure.

In Utero Exposure to Caffeine and Acetaminophen, the Gut Microbiome, and Neurodevelopmental Outcomes: A Prospective Birth Cohort Study.

Pregnant individuals are exposed to acetaminophen and caffeine, but it is unknown how these exposures interact with the developing gut microbiome. We aimed to determine whether acetaminophen and/or caffeine relate to the childhood gut microbiome and whether features of the gut microbiome alter the relationship between acetaminophen/caffeine and neurodevelopment. Forty-nine and 85 participants provided meconium and stool samples at 6-7, respectively, for exposure and microbiome assessment. Fecal acetaminophen and caffeine concentrations were quantified, and fecal DNA underwent metagenomic sequencing. Caregivers and study staff assessed the participants' motor and cognitive development using standardized scales. Prenatal exposures had stronger associations with the childhood microbiome than concurrent exposures. Prenatal acetaminophen exposure was associated with a trend of lower gut bacterial diversity in childhood [ß = -0.17 Shannon Index, 95% CI: (-0.31, -0.04)] and was marginally associated with differences in the relative abundances of features of the gut microbiome at the phylum (Firmicutes, Actinobacteria) and gene pathway levels. Among the participants with a higher relative abundance of Proteobacteria, prenatal exposure to acetaminophen and caffeine was associated with lower scores on WISC-IV subscales. Acetaminophen during bacterial colonization of the naïve gut is associated with lasting alterations in childhood microbiome composition. Future studies may inform our understanding of downstream health effects.

Nanoparticle Tracking Analysis for the Quantification and Size Determination of Extracellular Vesicles.

The physiological and pathophysiological roles of extracellular vesicles (EVs) have become increasingly recognized, making the EV field a quickly evolving area of research. There are many different methods for EV isolation, each with distinct advantages and disadvantages that affect the downstream yield and purity of EVs. Thus, characterizing the EV prep isolated from a given source by a chosen method is important for interpretation of downstream results and comparison of results across laboratories. Various methods exist for determining the size and quantity of EVs, which can be altered by disease states or in response to external conditions. Nanoparticle tracking analysis (NTA) is one of the prominent technologies used for high-throughput analysis of individual EVs. Here, we present a detailed protocol for quantification and size determination of EVs isolated from mouse perigonadal adipose tissue and human plasma using a breakthrough technology for NTA representing major advances in the field. The results demonstrate that this method can deliver reproducible and valid total particle concentration and size distribution data for EVs isolated from different sources using different methods, as confirmed by transmission electron microscopy. The adaptation of this instrument for NTA will address the need for standardization in NTA methods to increase rigor and reproducibility in EV research.

Placental mitochondrial DNA mutational load and perinatal outcomes: Findings from a multi-ethnic pregnancy cohort.

Mitochondria fuel placental activity, with mitochondrial dysfunction implicated in several perinatal complications. We investigated placental mtDNA mutational load using NextGen sequencing in relation to birthweight and gestational length among 358 mother-newborn pairs. We found that higher heteroplasmy, especially in the hypervariable displacement loop region, was associated with shorter gestational length. Results were similar among male and female pregnancies, but stronger in magnitude among females. With regard to growth, we observed that higher mutational load was associated with lower birthweight-for-gestational age (BWGA) among females, but higher BWGA among males. These findings support potential sex-differential fetal biological strategies for coping with increased heteroplasmies.

Isolation and characterization of extracellular vesicles in saliva of children with asthma.

AIM: To confirm the presence of extracellular vesicles (EVs) in cell-free saliva (CFS) of children with asthma and describe the isolated EV population. METHODS: A pooled sample of CFS EVs isolated from 180 participants using ExoQuick-TC was examined in downstream analyses. Transmission electron microscopy (TEM) was used to confirm the presence of EVs. Nanoparticle tracking analysis (NTA) and single particle interferometric reflectance imaging sensing (SP-IRIS) with fluorescence were used for sizing, counting, and phenotyping of EVs. Capillary immunoassays were used for protein quantitation. RESULTS: TEM confirmed the presence of EVs of diverse sizes, indicating the prep contained a heterogeneous population of EVs. Capillary immunoassays confirmed the presence of EV-associated proteins (CD9, CD63, CD81, ICAM-1, and ANXA5) and indicated limited cellular contamination. As others have also reported, there were discrepancies in the EV sizing and enumeration across platforms. Fluorescent NTA detected particles with a mode diameter of ~90 nm, whereas SP-IRIS reported sizes of ~55-60 nm that more closely approximated the TEM results. Consistent with protein immunoassay results, SP-IRIS with fluorescence showed that the majority of these EVs were CD9- and CD63-positive, with little expression of CD81. CONCLUSION: EVs from CFS can be isolated using a high-throughput method that can be scaled to large epidemiological studies. To our knowledge, we are the first to characterize CFS EVs from patients with asthma. The use of CFS EVs as potential novel biomarkers in asthma warrants further investigation and opens a new avenue of research for future studies.

Associations between maternal lifetime stressors and negative events in pregnancy and breast milk-derived extracellular vesicle microRNAs in the programming of intergenerational stress mechanisms (PRISM) pregnancy cohort.

Maternal stress is associated with adverse child health. Breast milk microRNAs encapsulated in extracellular vesicles (EVs) are involved in mother-infant biochemical communication during early-life programming. We leverage the PRogramming of Intergenerational Stress Mechanisms (PRISM) pregnancy cohort to investigate associations between maternal stress and breast milk EV-microRNAs. Lifetime stress and negative life events (NLEs) during pregnancy were assessed using the Life Stressor Checklist-Revised (LSCR) and the Crisis in Family Systems-Revised surveys, respectively. RNA was extracted from breast milk EVs (N = 80; collected 6.1 ± 5.9 weeks postnatally), and microRNAs were profiled using the TaqMan OpenArray Human miRNA panel. Associations between stress scores and detection (yes/no) of 173 microRNAs identified in 20-80% of samples were assessed using logistic regression; associations with expression levels of 205 EV-microRNAs identified in >50% of samples were assessed using linear regression. In adjusted models, detection of 60 and 44 EV-microRNAs was associated with higher LSCR and NLE scores, respectively (p < 0.05). Expression level of 8 and 17 EV-microRNAs was associated with LSCR and NLE scores, respectively, at our a priori criteria of p < 0.05 and |Bregression|>0.2. Enriched KEGG pathways for microRNAs associated with stress scores included fatty acid metabolism and the Hippo signaling pathway. Maternal lifetime stress and NLEs during pregnancy were both associated with detection and expression level of breast milk EV-microRNAs, although associations with microRNA profiles differed between stress measures. Further research is needed to identify biological pathways impacted by associated microRNAs and investigate relationships with child health outcomes.Abbreviations: EV: extracellular vesicle; PRISM: PRogramming of Intergenerational Stress Mechanisms pregnancy cohort; LSCR: Life Stressor Checklist-Revised survey; NLE: negative life event; CRISYS-R: Crisis in Family Systems-Revised survey; KEGG: Kyoto Encyclopaedia of Genes and Genomes; NYC: New York City; SD: standard deviation; IQR: interquartile range; Cq: relative cycle threshold values; PCA: principal component analysis.

Maternal anxiety and depression in pregnancy and DNA methylation of the NR3C1 glucocorticoid receptor gene.

Aim: To quantify associations of anxiety and depression during pregnancy with differential cord blood DNA methylation of the glucorticoid receptor (NR3C1). Materials & methods: Pregnancy anxiety, trait anxiety and depressive symptoms were collected using the Pregnancy Related Anxiety Scale, State-Trait Anxiety Index and Edinburgh Postnatal Depression Scale, respectively. NR3C1 methylation was determined at four methylation sites. Results: DNA methylation of CpG1 in the NR3C1 CpG island shore was higher in infants born to women with high pregnancy anxiety (ß 2.54, 95% CI: 0.49-4.58) and trait anxiety (ß 1.68, 95% CI: 0.14-3.22). No significant association was found between depressive symptoms and NR3C1 methylation. Conclusion: We found that maternal anxiety was associated with increased NR3C1 CpG island shore methylation.