Placental microRNAs relate to early childhood growth trajectories.

BACKGROUND: Poor placental function is a common cause of intrauterine growth restriction, which in turn is associated with increased risks of adverse health outcomes. Our prior work suggests that birthweight and childhood obesity-associated genetic variants functionally impact placental function and that placental microRNA are associated with birthweight. To address the influence of the placenta beyond birth, we assessed the relationship between placental microRNAs and early childhood growth. METHODS: Using the SITAR package, we generated two parameters that describe individual weight trajectories of children (0-5 years) in the New Hampshire Birth Cohort Study (NHBCS, n = 238). Using negative binomial generalized linear models, we identified placental microRNAs that relate to growth parameters (FDR < 0.1), while accounting for sex, gestational age at birth, and maternal parity. RESULTS: Genes targeted by the six growth trajectory-associated microRNAs are enriched (FDR < 0.05) in growth factor signaling (TGF/beta: miR-876; EGF/R: miR-155, Let-7c; FGF/R: miR-155; IGF/R: Let-7c, miR-155), calmodulin signaling (miR-216a), and NOTCH signaling (miR-629). CONCLUSIONS: Growth-trajectory microRNAs target pathways affecting placental proliferation, differentiation and function. Our results suggest a role for microRNAs in regulating placental cellular dynamics and supports the Developmental Origins of Health and Disease hypothesis that fetal environment can have impacts beyond birth. IMPACT: We found that growth trajectory associated placenta microRNAs target genes involved in signaling pathways central to the formation, maintenance and function of placenta; suggesting that placental cellular dynamics remain critical to infant growth to term and are under the control of microRNAs. Our results contribute to the existing body of research suggesting that the placenta plays a key role in programming health in the offspring. This is the first study to relate molecular patterns in placenta, specifically microRNAs, to early childhood growth trajectory.

An integrated -omics analysis of the epigenetic landscape of gene expression in human blood cells.

BACKGROUND: Gene expression can be influenced by DNA methylation 1) distally, at regulatory elements such as enhancers, as well as 2) proximally, at promoters. Our current understanding of the influence of distal DNA methylation changes on gene expression patterns is incomplete. Here, we characterize genome-wide methylation and expression patterns for ~ 13 k genes to explore how DNA methylation interacts with gene expression, throughout the genome. RESULTS: We used a linear mixed model framework to assess the correlation of DNA methylation at ~ 400 k CpGs with gene expression changes at ~ 13 k transcripts in two independent datasets from human blood cells. Among CpGs at which methylation significantly associates with transcription (eCpGs), > 50% are distal (> 50 kb) or trans (different chromosome) to the correlated gene. Many eCpG-transcript pairs are consistent between studies and ~ 90% of neighboring eCpGs associate with the same gene, within studies. We find that enhancers (P < 5e-18) and microRNA genes (P = 9e-3) are overrepresented among trans eCpGs, and insulators and long intergenic non-coding RNAs are enriched among cis and distal eCpGs. Intragenic-eCpG-transcript correlations are negative in 60-70% of occurrences and are enriched for annotated gene promoters and enhancers (P < 0.002), highlighting the importance of intragenic regulation. Gene Ontology analysis indicates that trans eCpGs are enriched for transcription factor genes and chromatin modifiers, suggesting that some trans eCpGs represent the influence of gene networks and higher-order transcriptional control. CONCLUSIONS: This work sheds new light on the interplay between epigenetic changes and gene expression, and provides useful data for mining biologically-relevant results from epigenome-wide association studies.

Placental transcriptome variation associated with season, location, and urinary prenatal pyrethroid metabolites of Thai farm-working women.

Prenatal exposure to pyrethroids is linked to adverse health effects in early life and proper placental function is critical to fetal development. This study explores the impact of prenatal pyrethroid exposure, as well as factors impacting exposure and effect, on the placental transcriptome, to understand pyrethroid exposures' relationship to placental function. The study of Asian Women and their Offspring's Development and Environmental Exposures (SAWASDEE) recruited pregnant farm-working women from two agricultural districts in the Chiang Mai province of Thailand between 2017 and 2019. This cohort was predominantly exposed to cypermethrin (type II), alongside pyrethroids such as cyfluthrin (type II) and permethrin (type I). In 253 participants, maternal urinary pyrethroid metabolites, 3-phenoxybenzoic acid (PBA), cis-3-(2,2-Dichlorovinyl)-2,2-dimethylcyclopropane carboxylic acid (CDCCA), and trans-3-(2,2-Dichlorovinyl)-2,2-dimethylcyclopropane carboxylic acid (TDCCA) were measured in early, middle, and late pregnancy and adjusted for urinary creatinine. The placental transcriptome was analyzed using RNA-Seq. Using generalized linear regression, we identified differentially expressed genes (DEGs) associated with the sum of each metabolite across pregnancy, as well as those associated with location of residence and season of birth. Pathway and upstream transcription factor analyses were performed to examine potential mechanisms associated with DEGs. Notably, TDCCA and CDCCA levels peaked in late pregnancy, with significant regional differences, particularly higher levels in the Fang region. Placental gene expression analysis showed no DEGs associated with individual metabolites at FDR<0.05. However, 251 DEGs by location, implicating immune response and oxidative phosphorylation pathways, were identified, while season of birth was associated with 2585 DEGs, over-represented in fibrosis signaling and metabolism pathways. Finally, transcription factor analysis identified 226 and 282 transcription factors associated with location and season, respectively, related to cell proliferation, differentiation, and the immune system. These alterations may have significant implications for fetal development and other pathologic processes, highlighting the importance of monitoring environmental exposures during pregnancy.

Variation in placental microRNA expression associates with maternal family history of cardiovascular disease.

In the United States, cardiovascular disease is the leading cause of death and the rate of maternal mortality remains among the highest of any industrialized nation. Maternal cardiometabolic health throughout gestation and postpartum is representative of placental health and physiology. Both proper placental functionality and placental microRNA expression are essential to successful pregnancy outcomes, and both are highly sensitive to genetic and environmental sources of variation. Placental pathologies, such as preeclampsia, are associated with maternal cardiovascular health but may also contribute to the developmental programming of chronic disease in offspring. However, the role of more subtle alterations to placental function and microRNA expression in this developmental programming remains poorly understood. We performed small RNA sequencing to investigate microRNA in placentae from the Rhode Island Child Health Study (n = 230). MicroRNA counts were modeled on maternal family history of cardiovascular disease using negative binomial generalized linear models. MicroRNAs were considered to be differentially expressed at a false discovery rate (FDR) less than 0.10. Parallel mRNA sequencing data and bioinformatic target prediction software were then used to identify potential mRNA targets of differentially expressed microRNAs. Nine differentially expressed microRNAs were identified (FDR < 0.1). Bioinformatic target prediction revealed 66 potential mRNA targets of these microRNAs, many of which are implicated in TGFß signaling pathway but also in pathways involving cellular metabolism and immunomodulation. A robust association exists between familial cardiovascular disease and placental microRNA expression which may be implicated in both placental insufficiencies and the developmental programming of chronic disease.

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.

Maternal depression and adverse neighbourhood conditions during pregnancy are associated with gestational epigenetic age deceleration.

Gestational epigenetic age (GEA) acceleration and deceleration can indicate developmental risk and may help elucidate how prenatal exposures lead to offspring outcomes. Depression and neighbourhood conditions during pregnancy are well-established determinants of birth and child outcomes. Emerging research suggests that maternal depression may contribute to GEA deceleration. It is unknown whether prenatal neighbourhood adversity would likewise influence GEA deceleration. This study examined whether maternal depression and neighbourhood conditions independently or jointly contributed to GEA deceleration, and which social and environmental neighbourhood conditions were associated with GEA. Participants were from the Albany Infant and Mother Study (n = 204), a prospective non-probability sampled cohort of higher risk racial/ethnic diverse mother/infant dyads. GEA was estimated from cord blood. Depressive symptoms and census-tract level neighbourhood conditions were assessed during pregnancy. Maternal depression (ß = -0.03, SE = 0.01, p = 0.008) and neighbourhood adversity (ß = -0.32, SE = 0.14, p = 0.02) were independently associated with GEA deceleration, controlling for all covariates including antidepressant use and cell type proportions. Neighbourhood adversity did not modify the association of maternal depression and GEA (ß = 0.003, SE = 0.03, p = 0.92). igher levels of neighbourhood poverty, public assistance, and lack of healthy food access were each associated with GEA deceleration; higher elementary school test scores (an indicator of community tax base) were associated with GEA acceleration (all p < 0.001). The results of this study indicated that maternal depression and neighbourhood conditions were independently and cumulatively associated GEA in this diverse population.

Prenatal lead (Pb) exposure is associated with differential placental DNA methylation and hydroxymethylation in a human population.

Prenatal lead (Pb) exposure is associated with adverse developmental outcomes and to epigenetic alterations such as DNA methylation and hydroxymethylation in animal models and in newborn blood. Given the importance of the placenta in foetal development, we sought to examine how prenatal Pb exposure was associated with differential placental DNA methylation and hydroxymethylation and to identify affected biological pathways linked to developmental outcomes. Maternal (n = 167) and infant (n = 172) toenail and placenta (n = 115) samples for prenatal Pb exposure were obtained from participants in a US birth cohort, and methylation and hydroxymethylation data were quantified using the Illumina Infinium MethylationEPIC BeadChip. An epigenome-wide association study was applied to identify differential methylation and hydroxymethylation associated with Pb exposure. Biological functions of the Pb-associated genes were determined by overrepresentation analysis through ConsensusPathDB. Prenatal Pb quantified from maternal toenail, infant toenail, and placenta was associated with 480, 27, and 2 differentially methylated sites (q < 0.05), respectively, with both increases and decreases associated with exposure. Alternatively, we identified 2, 1, and 14 differentially hydroxymethylated site(s) associated with maternal toenail, infant toenail, and placental Pb, respectively, with most showing increases in hydroxymethylation with exposure. Significantly overrepresented pathways amongst genes associated with differential methylation and hydroxymethylation (q < 0.10) included mechanisms pertaining to nervous system and organ development, calcium transport and regulation, and signalling activities. Our results suggest that both methylation and hydroxymethylation in the placenta can be variable based on Pb exposure and that the pathways impacted could affect placental function.

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.

Extracellular vesicle microRNA in early versus late pregnancy with birth outcomes in the MADRES study.

Circulating miRNA may contribute to the development of adverse birth outcomes. However, few studies have investigated extracellular vesicle (EV) miRNA, which play important roles in intercellular communication, or compared miRNA at multiple time points in pregnancy. In the current study, 800 miRNA were profiled for EVs from maternal plasma collected in early (median: 12.5 weeks) and late (median: 31.8 weeks) pregnancy from 156 participants in the MADRES Study, a health disparity pregnancy cohort. Associations between miRNA and birth weight, birth weight for gestational age (GA), and GA at birth were examined using covariate-adjusted robust linear regression. Differences by infant sex and maternal BMI were also investigated. Late pregnancy measures of 13 miRNA were associated with GA at birth (PFDR<0.050). Negative associations were observed for eight miRNA (miR-4454+ miR-7975, miR-4516, let-7b-5p, miR-126-3p, miR-29b-3p, miR-15a-5p, miR-15b-5p, miR-19b-3p) and positive associations for five miRNA (miR-212-3p, miR-584-5p, miR-608, miR-210-3p, miR-188-5p). Predicted target genes were enriched (PFDR<0.050) in pathways involved in organogenesis and placental development. An additional miRNA (miR-107), measured in late pregnancy, was positively associated with GA at birth in infants born to obese women (PFDR for BMI interaction = 0.011). In primary analyses, the associations between early pregnancy miRNA and birth outcomes were not statistically significant (PFDR≥0.05). However, sex-specific associations were observed for early pregnancy measures of 37 miRNA and GA at birth (PFDR for interactions<0.050). None of the miRNA were associated with fetal growth measures (PFDR≥0.050). Our findings suggest that EV miRNA in both early and late pregnancy may influence gestational duration.

Placental microRNA expression associates with birthweight through control of adipokines: results from two independent cohorts.

MicroRNAs are non-coding RNAs that regulate gene expression post-transcriptionally. In the placenta, the master regulator of foetal growth and development, microRNAs shape the basic processes of trophoblast biology and specific microRNA have been associated with foetal growth. To comprehensively assess the role of microRNAs in placental function and foetal development, we have performed small RNA sequencing to profile placental microRNAs from two independent mother-infant cohorts: the Rhode Island Child Health Study (n = 225) and the New Hampshire Birth Cohort Study (n = 317). We modelled microRNA counts on infant birthweight percentile (BWP) in each cohort, while accounting for race, sex, parity, and technical factors, using negative binomial generalized linear models. We identified microRNAs that were differentially expressed (DEmiRs) with BWP at false discovery rate (FDR) less than 0.05 in both cohorts. hsa-miR-532-5p (miR-532) was positively associated with BWP in both cohorts. By integrating parallel whole transcriptome and small RNA sequencing in the RICHS cohort, we identified putative targets of miR-532. These targets are enriched for pathways involved in adipogenesis, adipocytokine signalling, energy metabolism, and hypoxia response, and included Leptin, which we further demonstrated to have a decreasing expression with increasing BWP, particularly in male infants. Overall, we have shown a robust and reproducible association of miR-532 with BWP, which could influence BWP through regulation of adipocytokines Leptin and Adiponectin.