Gene Target Prediction of Environmental Chemicals Using Coupled Matrix-Matrix Completion.

Human exposure to toxic chemicals presents a huge health burden. Key to understanding chemical toxicity is knowledge of the molecular target(s) of the chemicals. Because a comprehensive safety assessment for all chemicals is infeasible due to limited resources, a robust computational method for discovering targets of environmental exposures is a promising direction for public health research. In this study, we implemented a novel matrix completion algorithm named coupled matrix-matrix completion (CMMC) for predicting direct and indirect exposome-target interactions, which exploits the vast amount of accumulated data regarding chemical exposures and their molecular targets. Our approach achieved an AUC of 0.89 on a benchmark data set generated using data from the Comparative Toxicogenomics Database. Our case studies with bisphenol A and its analogues, PFAS, dioxins, PCBs, and VOCs show that CMMC can be used to accurately predict molecular targets of novel chemicals without any prior bioactivity knowledge. Our results demonstrate the feasibility and promise of computationally predicting environmental chemical-target interactions to efficiently prioritize chemicals in hazard identification and risk assessment.

Cross-sectional associations between phthalates, phenols, and parabens with metabolic syndrome risk during early-to-mid adolescence among a cohort of Mexican youth.

BACKGROUND: Epidemiological studies on children and adults have linked toxicants from plastics and personal care products to metabolic disruption. Yet, the impact of endocrine-disrupting chemicals (EDCs) on adolescent metabolic syndrome (MetS) risk during early and mid-adolescence is unclear. METHODS: To examine the links between exposure to EDCs and MetS risk and its components, cross-sectional data from 344 Mexican youth in early-to-mid adolescence (10-17 years) were analyzed. Urinary biomarker concentrations of phthalates, phenol, and paraben analytes were measured from a single spot urine sample collected in 2015; study personnel obtained anthropometric and metabolic measures. We examined associations between summary phthalates and metabolites, phenol, and paraben analytes with MetS risk z-scores using linear regression, adjusted for specific gravity, sex, age, pubertal status, smoking, alcohol intake, physical activity level, and screen time. As a secondary aim, mediation analysis was conducted to evaluate the role of hormones in the association between summary phthalates with lipids and MetS risk z-scores. RESULTS: The mean (SD) age was 13.2 (1.9) years, and 50.9% were female. Sex-stratified analyses revealed associations between summary phthalates and lipids ratio z-scores, including Σ DEHP [ß = 0.21 (95% CI: 0.04, 0.37; p < 0.01)], phthalates from plastic sources (Σ Plastic) [ß = 0.22 (95% CI: 0.05, 0.39; p < 0.01)], anti-androgenic phthalates (Σ AA) [ß = 0.22 (95% CI: 0.05, 0.39; p < 0.01)], and individual phthalate metabolites (MEHHP, MEOHP, and MECPP) among males. Among females, BPA [ß = 0.24 (95% CI: 0.03, 0.44; p < 0.05)] was positively associated with lipids ratio z-score and one phenol (2,5 DCP) [ß = 0.09 (95% CI: 0.01, 0.18); p < 0.05)] was associated with increased waist circumference z-score. Results showed no evidence of mediation by hormone concentrations in the association between summary phthalates with lipids ratio or MetS risk z-scores. CONCLUSION: Higher EDC exposure was positively associated with serum lipids during adolescence, particularly among males.

Toxicoepigenetics and Environmental Health: Challenges and Opportunities.

The rapidly growing field of toxicoepigenetics seeks to understand how toxicant exposures interact with the epigenome to influence disease risk. Toxicoepigenetics is a promising field of environmental health research, as integrating epigenetics into the field of toxicology will enable a more thorough evaluation of toxicant-induced disease mechanisms as well as the elucidation of the role of the epigenome as a biomarker of exposure and disease and possible mediator of exposure effects. Likewise, toxicoepigenetics will enhance our knowledge of how environmental exposures, lifestyle factors, and diet interact to influence health. Ultimately, an understanding of how the environment impacts the epigenome to cause disease may inform risk assessment, permit noninvasive biomonitoring, and provide potential opportunities for therapeutic intervention. However, the translation of research from this exciting field into benefits for human and animal health presents several challenges and opportunities. Here, we describe four significant areas in which we see opportunity to transform the field and improve human health by reducing the disease burden caused by environmental exposures. These include (1) research into the mechanistic role for epigenetic change in environment-induced disease, (2) understanding key factors influencing vulnerability to the adverse effects of environmental exposures, (3) identifying appropriate biomarkers of environmental exposures and their associated diseases, and (4) determining whether the adverse effects of environment on the epigenome and human health are reversible through pharmacologic, dietary, or behavioral interventions. We then highlight several initiatives currently underway to address these challenges.

Early life social and ecological determinants of global DNA methylation in wild spotted hyenas.

Environmental factors early in life can have lasting influence on the development and phenotypes of animals, but the underlying molecular modifications remain poorly understood. We examined cross-sectional associations among early life socioecological factors and global DNA methylation in 293 wild spotted hyenas (Crocuta crocuta) in the Masai Mara National Reserve, Kenya, grouped according to three age classes (cub, subadult and adult). Explanatory variables of interest included annual maternal rank based on outcomes of dyadic agonistic interactions, litter size, wild ungulate prey density and anthropogenic disturbance in the year each hyena was born based on counts of illegal livestock in the Reserve. The dependent variable of interest was global DNA methylation, assessed via the LUminometric Methylation Assay, which provides a percentage methylation value calculated at CCGG sites across the genome. Among cubs, we observed approximately 2.75% higher CCGG methylation in offspring born to high- than low-ranking mothers. Among cubs and subadults, higher anthropogenic disturbance corresponded with greater %CCGG methylation. In both cubs and adults, we found an inverse association between prey density measured before a hyena was 3 months old and %CCGG methylation. Our results suggest that maternal rank, anthropogenic disturbance and prey availability early in life are associated with later life global DNA methylation. Future studies are required to understand the extent to which these DNA methylation patterns relate to adult phenotypes and fitness outcomes.

Neonatal bloodspot DNA methylation patterns are associated with childhood weight status in the Healthy Families Project.

BACKGROUND: This study measured longitudinal DNA methylation dynamics at growth-related genes during childhood, and then tested whether DNA methylation at various stages of childhood was associated with obesity status. METHODS: Using neonatal bloodspot (n = 132) and matched childhood blood samples (n = 65), DNA methylation was quantified at a repetitive element (long interspersed nuclear element-1 (LINE-1)), two imprinted genes (IGF2, H19), and four non-imprinted genes (LEP, PPARA, ESR1, SREBF1) related to growth and adiposity. Logistic regression was used to test whether neonatal bloodspot DNA methylation at target genes was associated with log odds of obesity (Y/N) in children recruited from three age groups-12-24 months old (n = 40), 3-5 years of age (n = 40), and 10-12 years of age (n = 52). RESULTS: In 3-5 year olds, neonatal bloodspot LINE-1 methylation was negatively associated with obesity (log odds = -0.40, p = 0.04). Across childhood age group in matched blood samples, DNA methylation levels in blood decreased (p < 0.05) at LINE-1, PPARA, ESR1, SREBF1, IGF2, and H19, and increased (p < 0.05) at LEP. CONCLUSIONS: Our results suggest that age-related epigenetic changes occur at growth-related genes in the first decade of life, and that gene-specific neonatal bloodspot DNA methylation may be a useful biomarker of obesity likelihood during childhood.

Correlation between Conjugated Bisphenol A Concentrations and Efflux Transporter Expression in Human Fetal Livers.

Because of its widespread use in the manufacturing of consumer products over several decades, human exposure to bisphenol A (BPA) has been pervasive. Fetuses are particularly sensitive to BPA exposure, with a number of negative developmental and reproductive outcomes observed in rodent perinatal models. Xenobiotic transporters are one mechanism to extrude conjugated and unconjugated BPA from the liver. In this study, the mRNA expression of xenobiotic transporters and relationships with total, conjugated, and free BPA levels were explored utilizing human fetal liver samples. The mRNA expression of breast cancer resistance protein (BCRP) and multidrug resistance-associated transporter (MRP)4, as well as BCRP and multidrug resistance transporter 1 exhibited the highest degree of correlation, with r(2) values of 0.941 and 0.816 (P < 0.001 for both), respectively. Increasing concentrations of conjugated BPA significantly correlated with high expression of MRP1 (P < 0.001), MRP2 (P < 0.05), and MRP3 (P < 0.05) transporters, in addition to the NF-E2-related factor 2 transcription factor (P < 0.001) and its prototypical target gene, NAD(P)H quinone oxidoreductase 1 (P < 0.001). These data demonstrate that xenobiotic transporters may be coordinately expressed in the human fetal liver. This is also the first report of a relationship between environmentally relevant fetal BPA levels and differences in the expression of transporters that can excrete the parent compound and its metabolites.

Gene-specific DNA methylation may mediate atypical antipsychotic-induced insulin resistance.

OBJECTIVES: Atypical antipsychotics (AAPs) carry a significant risk of cardiometabolic side effects, including insulin resistance. It is thought that the insulin resistance resulting from the use of AAPs may be associated with changes in DNA methylation. We aimed to identify and validate a candidate gene associated with AAP-induced insulin resistance by using a multi-step approach that included an epigenome-wide association study (EWAS) and validation with site-specific methylation and metabolomics data. METHODS: Subjects with bipolar disorder treated with AAPs or lithium monotherapy were recruited for a cross-sectional visit to analyze peripheral blood DNA methylation and insulin resistance. Epigenome-wide DNA methylation was analyzed in a discovery sample (n = 48) using the Illumina 450K BeadChip. Validation analyses of the epigenome-wide findings occurred in a separate sample (n = 72) using site-specific methylation with pyrosequencing and untargeted metabolomics data. Regression analyses were conducted controlling for known confounders in all analyses and a mediation analysis was performed to investigate if AAP-induced insulin resistance occurs through changes in DNA methylation. RESULTS: A differentially methylated probe associated with insulin resistance was discovered and validated in the fatty acyl CoA reductase 2 (FAR2) gene of chromosome 12. Functional associations of this DNA methylation site with untargeted phospholipid-related metabolites were also detected. Our results identified a mediating effect of this FAR2 methylation site on AAP-induced insulin resistance. CONCLUSIONS: Going forward, prospective, longitudinal studies assessing comprehensive changes in FAR2 DNA methylation, expression, and lipid metabolism before and after AAP treatment are required to assess its potential role in the development of insulin resistance.

Genetic polymorphisms are associated with hair, blood, and urine mercury levels in the American Dental Association (ADA) study participants.

BACKGROUND/AIMS: Mercury (Hg) is a potent toxicant of concern to the general public. Recent studies suggest that several genes that mediate Hg metabolism are polymorphic. We hypothesize that single nucleotide polymorphisms (SNPs) in such genes may underline inter-individual differences in exposure biomarker concentrations. METHODS: Dental professionals were recruited during the American Dental Association (ADA) 2012 Annual Meeting. Samples of hair, blood, and urine were collected for quantifying Hg levels and genotyping (88 SNPs in classes relevant to Hg toxicokinetics including glutathione metabolism, selenoproteins, metallothioneins, and xenobiotic transporters). Questionnaires were administrated to obtain information on demographics and sources of Hg exposure (e.g., fish consumption and use of dental amalgam). Here, we report results for 380 participants with complete genotype and Hg biomarker datasets. ANOVA and linear regressions were used for statistical analysis. RESULTS: Mean (geometric) Hg levels in hair (hHg), blood (bHg), urine (uHg), and the average estimated Hg intake from fish were 0.62µg/g, 3.75µg/L, 1.32µg/L, and 0.12µg/kg body weight/day, respectively. Out of 88 SNPs successfully genotyped, Hg biomarker levels differed by genotype for 25 SNPs, one of which remained significant following Bonferroni correction in ANOVA. When the associations between sources of Hg exposure and SNPs were analyzed with respect to Hg biomarker concentrations, 38 SNPs had significant main effects and/or gene-Hg exposure source interactions. Twenty-five, 23, and four SNPs showed significant main effects and/or interactions for hHg, bHg, and uHg levels, respectively (p<0.05), and six SNPs (in GCLC, MT1M, MT4, ATP7B, and BDNF) remained significant following Bonferroni correction. CONCLUSION: The findings suggest that polymorphisms in environmentally-responsive genes can influence Hg biomarker levels. Hence, consideration of such gene-environment factors may improve the ability to assess the health risks of Hg more precisely.

Perinatal bisphenol A exposure promotes dose-dependent alterations of the mouse methylome.

BACKGROUND: Environmental factors during perinatal development may influence developmental plasticity and disease susceptibility via alterations to the epigenome. Developmental exposure to the endocrine active compound, bisphenol A (BPA), has previously been associated with altered methylation at candidate gene loci. Here, we undertake the first genome-wide characterization of DNA methylation profiles in the liver of murine offspring exposed perinatally to multiple doses of BPA through the maternal diet. RESULTS: Using a tiered focusing approach, our strategy proceeds from unbiased broad DNA methylation analysis using methylation-based next generation sequencing technology to in-depth quantitative site-specific CpG methylation determination using the Sequenom EpiTYPER MassARRAY platform to profile liver DNA methylation patterns in offspring maternally exposed to BPA during gestation and lactation to doses ranging from 0 BPA/kg (Ctr), 50 µg BPA/kg (UG), or 50 mg BPA/kg (MG) diet (N = 4 per group). Genome-wide analyses indicate non-monotonic effects of DNA methylation patterns following perinatal exposure to BPA, corroborating previous studies using multiple doses of BPA with non-monotonic outcomes. We observed enrichment of regions of altered methylation (RAMs) within CpG island (CGI) shores, but little evidence of RAM enrichment in CGIs. An analysis of promoter regions identified several hundred novel BPA-associated methylation events, and methylation alterations in the Myh7b and Slc22a12 gene promoters were validated. Using the Comparative Toxicogenomics Database, a number of candidate genes that have previously been associated with BPA-related gene expression changes were identified, and gene set enrichment testing identified epigenetically dysregulated pathways involved in metabolism and stimulus response. CONCLUSIONS: In this study, non-monotonic dose dependent alterations in DNA methylation among BPA-exposed mouse liver samples and their relevant pathways were identified and validated. The comprehensive methylome map presented here provides candidate loci underlying the role of early BPA exposure and later in life health and disease status.

Epigenetics: relevance and implications for public health.

Improved understanding of the multilayer regulation of the human genome has led to a greater appreciation of environmental, nutritional, and epigenetic risk factors for human disease. Chromatin remodeling, histone tail modifications, and DNA methylation are dynamic epigenetic changes responsive to external stimuli. Careful interpretation can provide insights for actionable public health through collaboration between population and basic scientists and through integration of multiple data sources. We review key findings in environmental epigenetics both in human population studies and in animal models, and discuss the implications of these results for risk assessment and public health protection. To ultimately succeed in identifying epigenetic mechanisms leading to complex phenotypes and disease, researchers must integrate the various animal models, human clinical approaches, and human population approaches while paying attention to life-stage sensitivity, to generate effective prescriptions for human health evaluation and disease prevention.

Adaptive radiation-induced epigenetic alterations mitigated by antioxidants.

Humans are exposed to low-dose ionizing radiation (LDIR) from a number of environmental and medical sources. In addition to inducing genetic mutations, there is concern that LDIR may also alter the epigenome. Such heritable effects early in life can either be positively adaptive or result in the enhanced formation of diseases, including cancer, diabetes, and obesity. Herein, we show that LDIR significantly increased DNA methylation at the viable yellow agouti (A(vy)) locus in a sex-specific manner (P=0.004). Average DNA methylation was significantly increased in male offspring exposed to doses between 0.7 and 7.6 cGy, with maximum effects at 1.4 and 3.0 cGy (P<0.01). Offspring coat color was concomitantly shifted toward pseudoagouti (P<0.01). Maternal dietary antioxidant supplementation mitigated both the DNA methylation changes and coat color shift in the irradiated offspring. Thus, LDIR exposure during gestation elicits epigenetic alterations that lead to positive adaptive phenotypic changes that are negated with antioxidants, indicating they are mediated in part by oxidative stress. These findings provide evidence that in the isogenic A(vy) mouse model, epigenetic alterations resulting from LDIR play a role in radiation hormesis, bringing into question the assumption that every dose of radiation is harmful.

Perinatal bisphenol A exposure promotes hyperactivity, lean body composition, and hormonal responses across the murine life course.

The development of adult-onset diseases is influenced by perinatal exposure to altered environmental conditions. One such exposure, bisphenol A (BPA), has been associated with obesity and diabetes, and consequently labeled an obesogen. Using an isogenic murine model, we examined the effects of perinatal exposure through maternal diet to 50 ng (n=20), 50 µg (n=21), or 50 mg (n=18) BPA/kg diet, as well as controls (n=20) on offspring energy expenditure, spontaneous activity, and body composition at 3, 6, and 9 mo of age, and hormone levels at 9 and 10 mo of age. Overall, exposed females and males exhibited increased energy expenditure (P<0.001 and 0.001, respectively) throughout the life course. In females, horizontal and vertical activity increased (P=0.07 and 0.06, respectively) throughout the life course. Generally, body composition measures were not different throughout the life course in exposed females or males (all P>0.44), although body fat and weight decreased in exposed females at particular ages (all P<0.08). Milligram-exposed females had improved glucose, insulin, adiponectin, and leptin profiles (all P<0.10). Thus, life-course analysis illustrates that BPA is associated with hyperactive and lean phenotypes. Variability across studies may be attributable to differential exposure duration and timing, dietary fat and phytoestrogen content, or lack of sophisticated phenotyping across the life course.

Sex-Specific Deflection of Age-Related DNA Methylation and Gene Expression in Mouse Heart by Perinatal Toxicant Exposures.

Background: Global and site-specific changes in DNA methylation and gene expression are associated with cardiovascular aging and disease, but how toxicant exposures during early development influence the normal trajectory of these age-related molecular changes, and whether there are sex differences, has not yet been investigated. Objectives: We used an established mouse model of developmental exposures to investigate the effects of perinatal exposure to either lead (Pb) or diethylhexyl phthalate (DEHP), two ubiquitous environmental contaminants strongly associated with CVD, on age-related cardiac DNA methylation and gene expression. Methods: Dams were randomly assigned to receive human physiologically relevant levels of Pb (32 ppm in water), DEHP (25 mg/kg chow), or control water and chow. Exposures started two weeks prior to mating and continued until weaning at postnatal day 21 (3 weeks of age). Approximately one male and one female offspring per litter were followed to 3 weeks, 5 months, or 10 months of age, at which time whole hearts were collected (n ≥ 5 per sex per exposure). Enhanced reduced representation bisulfite sequencing (ERRBS) was used to assess the cardiac DNA methylome at 3 weeks and 10 months, and RNA-seq was conducted at all 3 time points. MethylSig and edgeR were used to identify age-related differentially methylated regions (DMRs) and differentially expressed genes (DEGs), respectively, within each sex and exposure group. Cell type deconvolution of bulk RNA-seq data was conducted using the MuSiC algorithm and publicly available single cell RNA-seq data. Results: Thousands of DMRs and hundreds of DEGs were identified in control, DEHP, and Pb-exposed hearts across time between 3 weeks and 10 months of age. A closer look at the genes and pathways showing differential DNA methylation revealed that the majority were unique to each sex and exposure group. Overall, pathways governing development and differentiation were most frequently altered with age in all conditions. A small number of genes in each group showed significant changes in DNA methylation and gene expression with age, including several that were altered by both toxicants but were unchanged in control. We also observed subtle, but significant changes in the proportion of several cell types due to age, sex, and developmental exposure. Discussion: Together these data show that perinatal Pb or DEHP exposures deflect normal age-related gene expression, DNA methylation programs, and cellular composition across the life course, long after cessation of exposure, and highlight potential biomarkers of developmental toxicant exposures. Further studies are needed to investigate how these epigenetic and transcriptional changes impact cardiovascular health across the life course.

piOxi database: a web resource of germline and somatic tissue piRNAs identified by chemical oxidation.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly expressed and extensively studied from the germline. piRNAs associate with PIWI proteins to maintain DNA methylation for transposon silencing and transcriptional gene regulation for genomic stability. Mature germline piRNAs have distinct characteristics including a 24- to 32-nucleotide length and a 2'-O-methylation signature at the 3' end. Although recent studies have identified piRNAs in somatic tissues, they remain poorly characterized. For example, we recently demonstrated notable expression of piRNA in the murine soma, and while overall expression was lower than that of the germline, unique characteristics suggested tissue-specific functions of this class. While currently available databases commonly use length and association with PIWI proteins to identify piRNA, few have included a chemical oxidation method that detects piRNA based on its 3' modification. This method leads to reproducible and rigorous data processing when coupled with next-generation sequencing and bioinformatics analysis. Here, we introduce piOxi DB, a user-friendly web resource that provides a comprehensive analysis of piRNA, generated exclusively through sodium periodate treatment of small RNA. The current version of piOxi DB includes 435 749 germline and 9828 somatic piRNA sequences robustly identified from M. musculus, M. fascicularis and H. sapiens. The database provides species- and tissue-specific data that are further analyzed according to chromosome location and correspondence to gene and repetitive elements. piOxi DB is an informative tool to assist broad research applications in the fields of RNA biology, cancer biology, environmental toxicology and beyond. Database URL:  https://pioxidb.dcmb.med.umich.edu/.

Placental and Immune Cell DNA Methylation Reference Panel for Bulk Tissue Cell Composition Estimation in Epidemiological Studies.

To distinguish DNA methylation (DNAm) from cell proportion changes in whole placental tissue research, we developed a robust cell type-specific DNAm reference to estimate cell composition. We collated newly collected and existing cell type DNAm profiles quantified via Illumina EPIC or 450k microarrays. To estimate cell composition, we deconvoluted whole placental samples (n=36) with robust partial correlation based on the top 50 hyper- and hypomethylated sites per cell type. To test deconvolution performance, we evaluated RMSE in predicting principal component one of DNAm variation in 204 external placental samples. We analyzed DNAm profiles (n=368,435 sites) from 12 cell types: cytotrophoblasts (n=18), endothelial cells (n=19), Hofbauer cells (n=26), stromal cells (n=21), syncytiotrophoblasts (n=4), six lymphocyte types (n=36), and nucleated red blood cells (n=11). Median cell composition was consistent with placental biology: 60.4% syncytiotrophoblast, 17.1% stromal, 8.8% endothelial, 4.5% cytotrophoblast, 3.9% Hofbauer, 1.7% nucleated red blood cells, and 1.2% neutrophils. Our expanded reference outperformed an existing reference in predicting DNAm variation (15.4% variance explained, IQR=21.61) with cell composition estimates (RMSE:10.51 vs. 11.43, p-value<0.001). This cell type reference can robustly estimate cell composition from whole placental DNAm data to detect important cell types, reveal biological mechanisms, and improve casual inference.

Effects of Developmental Lead and Phthalate Exposures on DNA Methylation in Adult Mouse Blood, Brain, and Liver: A Focus on Genomic Imprinting by Tissue and Sex.

BACKGROUND: Maternal exposure to environmental chemicals can cause adverse health effects in offspring. Mounting evidence supports that these effects are influenced, at least in part, by epigenetic modifications. It is unknown whether epigenetic changes in surrogate tissues such as the blood are reflective of similar changes in target tissues such as cortex or liver. OBJECTIVE: We examined tissue- and sex-specific changes in DNA methylation (DNAm) associated with human-relevant lead (Pb) and di(2-ethylhexyl) phthalate (DEHP) exposure during perinatal development in cerebral cortex, blood, and liver. METHODS: Female mice were exposed to human relevant doses of either Pb (32 ppm) via drinking water or DEHP (5mg/kg-day) via chow for 2 weeks prior to mating through offspring weaning. Whole genome bisulfite sequencing (WGBS) was utilized to examine DNAm changes in offspring cortex, blood, and liver at 5 months of age. Metilene and methylSig were used to identify differentially methylated regions (DMRs). Annotatr and ChIP-enrich were used for genomic annotations and gene set enrichment tests of DMRs, respectively. RESULTS: The cortex contained the majority of DMRs associated with Pb (66%) and DEHP (57%) exposure. The cortex also contained the greatest degree of overlap in DMR signatures between sexes (n=13 and 8 DMRs with Pb and DEHP exposure, respectively) and exposure types (n=55 and 39 DMRs in males and females, respectively). In all tissues, detected DMRs were preferentially found at genomic regions associated with gene expression regulation (e.g., CpG islands and shores, 5' UTRs, promoters, and exons). An analysis of GO terms associated with DMR-containing genes identified imprinted genes to be impacted by both Pb and DEHP exposure. Of these, Gnas and Grb10 contained DMRs across tissues, sexes, and exposures, with some signatures replicated between target and surrogate tissues. DMRs were enriched in the imprinting control regions (ICRs) of Gnas and Grb10, and we again observed a replication of DMR signatures between blood and target tissues. Specifically, we observed hypermethylation of the Grb10 ICR in both blood and liver of Pb-exposed male animals. CONCLUSIONS: These data provide preliminary evidence that imprinted genes may be viable candidates in the search for epigenetic biomarkers of toxicant exposure in target tissues. Additional research is needed on allele- and developmental stage-specific effects, as well as whether other imprinted genes provide additional examples of this relationship. https://doi.org/10.1289/EHP14074.

Associations between Sleep and Physical Activity Behavior Clusters and Epigenetic Age Acceleration in Mexican Adolescents.

INTRODUCTION: Epigenetic aging, a marker of biological aging measured by DNA methylation, may be affected by behaviors, including sleep and physical activity. However, investigations of physical activity and sleep with epigenetic aging among pediatric populations are scant and have not accounted for correlated behaviors. METHODS: The study population included 472 Mexico City adolescents (52% female). Blood collection and 7-day wrist actigraphy (Actigraph GTX-BT) occurred during a follow-up visit when participants were 14.5 (2.09) years. Leukocyte DNA methylation was measured with the Infinium MethylationEPIC array after bisulfite conversion, and 9 epigenetic clocks were calculated. Sleep vs wake time was identified through a pruned dynamic programing algorithm, and physical activity was processed with Chandler cut-offs. Kmeans clustering was used to select actigraphy-assessed physical activity and sleep behavior clusters. Linear regression analyses were used to evaluate adjusted associations between the clusters and epigenetic aging. RESULTS: There were 3 unique clusters: "Short sleep/high sedentary behavior", "Adequate sleep duration and late timing/low moderate or vigorous physical activity (MVPA)", and "Adequate sleep duration/high MVPA". Compared to the "Adequate duration/high MVPA", adolescents with "Adequate duration and late sleep timing/low MVPA" had more accelerated aging for the GrimAge clock (ß = 0.63;95% CI 0.07, 1.19). In pubertal-stratified analyses, more mature adolescents in the "Adequate duration and late sleep timing/low MVPA group" had accelerated epigenetic aging. In contrast, females in the "Short sleep/high sedentary" group had decelerated epigenetic aging for the Wu pediatric clock. CONCLUSIONS: Associations between behavior clusters and epigenetic aging varied by pubertal status and sex. Contrary results in the Wu clock suggest the need for future research on pediatric-specific clocks.

Translational toxicoepigenetic Meta-Analyses identify homologous gene DNA methylation reprogramming following developmental phthalate and lead exposure in mouse and human offspring.

Although toxicology uses animal models to represent real-world human health scenarios, a critical translational gap between laboratory-based studies and epidemiology remains. In this study, we aimed to understand the toxicoepigenetic effects on DNA methylation after developmental exposure to two common toxicants, the phthalate di(2-ethylhexyl) phthalate (DEHP) and the metal lead (Pb), using a translational paradigm that selected candidate genes from a mouse study and assessed them in four human birth cohorts. Data from mouse offspring developmentally exposed to DEHP, Pb, or control were used to identify genes with sex-specific sites with differential DNA methylation at postnatal day 21. Associations of human infant DNA methylation in homologous mouse genes with prenatal DEHP or Pb were examined with a meta-analysis. Differential methylation was observed on 6 cytosines (adjusted-p < 0.05) and 90 regions (adjusted-p < 0.001). This translational approach offers a unique method that can detect conserved epigenetic differences that are developmentally susceptible to environmental toxicants.

Phylogenetic and DNA methylation analysis reveal novel regions of variable methylation in the mouse IAP class of transposons.

BACKGROUND: Select retrotransposons in the long terminal repeat (LTR) class exhibit interindividual variation in DNA methylation that is altered by developmental environmental exposures. Yet, neither the full extent of variability at these "metastable epialleles," nor the phylogenetic relationship underlying variable elements is well understood. The murine metastable epialleles, Avy and CabpIAP, result from independent insertions of an intracisternal A particle (IAP) mobile element, and exhibit remarkably similar sequence identity (98.5%). RESULTS: Utilizing the C57BL/6 genome we identified 10802 IAP LTRs overall and a subset of 1388 in a family that includes Avy and CabpIAP. Phylogenetic analysis revealed two duplication and divergence events subdividing this family into three clades. To characterize interindividual variation across clades, liver DNA from 17 isogenic mice was subjected to combined bisulfite and restriction analysis (CoBRA) for 21 separate LTR transposons (7 per clade). The lowest and highest mean methylation values were 59% and 88% respectively, while methylation levels at individual LTRs varied widely, ranging from 9% to 34%. The clade with the most conserved elements had significantly higher mean methylation across LTRs than either of the two diverged clades (p = 0.040 and p = 0.017). Within each mouse, average methylation across all LTRs was not significantly different (71%-74%, p > 0.99). CONCLUSIONS: Combined phylogenetic and DNA methylation analysis allows for the identification of novel regions of variable methylation. This approach increases the number of known metastable epialleles in the mouse, which can serve as biomarkers for environmental modifications to the epigenome.