Evaluating changes in firefighter urinary metabolomes after structural fires: an untargeted, high resolution approach.

Firefighters have elevated rates of urinary tract cancers and other adverse health outcomes, which may be attributable to environmental occupational exposures. Untargeted metabolomics was applied to characterize this suite of environmental exposures and biological changes in response to occupational firefighting. 200 urine samples from 100 firefighters collected at baseline and two to four hours post-fire were analyzed using untargeted liquid-chromatography and high-resolution mass spectrometry. Changes in metabolite abundance after a fire were estimated with fixed effects linear regression, with false discovery rate (FDR) adjustment. Partial least squares discriminant analysis (PLS-DA) was also used, and variable important projection (VIP) scores were extracted. Systemic changes were evaluated using pathway enrichment for highly discriminating metabolites. Metabolome-wide-association-study (MWAS) identified 268 metabolites associated with firefighting activity at FDR q < 0.05. Of these, 20 were annotated with high confidence, including the amino acids taurine, proline, and betaine; the indoles kynurenic acid and indole-3-acetic acid; the known uremic toxins trimethylamine n-oxide and hippuric acid; and the hormone 7a-hydroxytestosterone. Partial least squares discriminant analysis (PLS-DA) additionally implicated choline, cortisol, and other hormones. Significant pathways included metabolism of urea cycle/amino group, alanine and aspartate, aspartate and asparagine, vitamin b3 (nicotinate and nicotinamide), and arginine and proline. Firefighters show a broad metabolic response to fires, including altered excretion of indole compounds and uremic toxins. Implicated pathways and features, particularly uremic toxins, may be important regulators of firefighter's increased risk for urinary tract cancers.

Associations Between Epigenetic Age Acceleration and microRNA Expression Among U.S. Firefighters.

Epigenetic changes may be biomarkers of health. Epigenetic age acceleration (EAA), the discrepancy between epigenetic age measured via epigenetic clocks and chronological age, is associated with morbidity and mortality. However, the intersection of epigenetic clocks with microRNAs (miRNAs) and corresponding miRNA-based health implications have not been evaluated. We analyzed DNA methylation and miRNA profiles from blood sampled among 332 individuals enrolled across 2 U.S.-based firefighter occupational studies (2015-2018 and 2018-2020). We considered 7 measures of EAA in leukocytes (PhenoAge, GrimAge, Horvath, skin-blood, and Hannum epigenetic clocks, and extrinsic and intrinsic epigenetic age acceleration). We identified miRNAs associated with EAA using individual linear regression models, adjusted for sex, race/ethnicity, chronological age, and cell type estimates, and investigated downstream effects of associated miRNAs with miRNA enrichment analyses and genomic annotations. On average, participants were 38 years old, 88% male, and 75% non-Hispanic white. We identified 183 of 798 miRNAs associated with EAA (FDR q < 0.05); 126 with PhenoAge, 59 with GrimAge, 1 with Horvath, and 1 with the skin-blood clock. Among miRNAs associated with Horvath and GrimAge, there were 61 significantly enriched disease annotations including age-related metabolic and cardiovascular conditions and several cancers. Enriched pathways included those related to proteins and protein modification. We identified miRNAs associated with EAA of multiple epigenetic clocks. PhenoAge had more associations with individual miRNAs, but GrimAge and Horvath had greater implications for miRNA-associated pathways. Understanding the relationship between these epigenetic markers could contribute to our understanding of the molecular underpinnings of aging and aging-related diseases.

Repeat measures of DNA methylation in an inception cohort of firefighters.

OBJECTIVES: Firefighters face exposures associated with adverse health outcomes including risk for multiple cancers. DNA methylation, one type of epigenetic regulation, provides a potential mechanism linking occupational hazards to adverse health outcomes. We hypothesised that DNA methylation profiles would change in firefighters after starting their service and that these patterns would be associated with occupational exposures (cumulative fire-hours and fire-runs). METHODS: We profiled DNA methylation with the Infinium MethylationEPIC in blood leucocytes at two time points in non-smoking new recruits: prior to live fire training and 20-37 months later. Linear mixed effects models adjusted for potential confounders were used to identify differentially methylated CpG sites over time using data from 50 individuals passing all quality control. RESULTS: We report 680 CpG sites with altered methylation (q value <0.05) including 60 with at least a 5% methylation difference at follow-up. Genes with differentially methylated CpG sites were enriched in biological pathways related to cancers, neurological function, cell signalling and transcription regulation. Next, linear mixed effects models were used to determine associations between occupational exposures with methylation at the 680 loci. Of these, more CpG sites were associated with fire-runs (108 for all and 78 for structure-fires only, q<0.05) than with fire-hours (27 for all fires and 1 for structure fires). These associations were independent of time since most recent fire, suggesting an impact of cumulative exposures. CONCLUSIONS: Overall, this study provides evidence that DNA methylation may be altered by fireground exposures, and the impact of this change on disease development should be evaluated.

Differential DNA Methylation by Hispanic Ethnicity Among Firefighters in the United States.

Firefighters are exposed to a variety of environmental hazards and are at increased risk for multiple cancers. There is evidence that risks differ by ethnicity, yet the biological or environmental differences underlying these differences are not known. DNA methylation is one type of epigenetic regulation that is altered in cancers. In this pilot study, we profiled DNA methylation with the Infinium MethylationEPIC in blood leukocytes from 31 Hispanic white and 163 non-Hispanic white firefighters. We compared DNA methylation (1) at 12 xenobiotic metabolizing genes and (2) at all loci on the array (>740 000), adjusting for confounders. Five of the xenobiotic metabolizing genes were differentially methylated at a raw P-value <.05 when comparing the 2 ethnic groups, yet were not statistically significant at a 5% false discovery rate (q-value <.05). In the epigenome-wide analysis, 76 loci exhibited DNA methylation differences at q < .05. Among these, 3 CpG sites in the promoter region of the biotransformation gene SULT1C2 had lower methylation in Hispanic compared to non-Hispanic firefighters. Other differentially methylated loci included genes that have been implicated in carcinogenesis in published studies (FOXK2, GYLTL1B, ZBTB16, ARHGEF10, and more). In this pilot study, we report differential DNA methylation between Hispanic and non-Hispanic firefighters in xenobiotic metabolism genes and other genes with functions related to cancer. Epigenetic susceptibility by ethnicity merits further study as this may alter risk for cancers linked to toxic exposures.

An epigenome-wide association study of ambient pyrethroid pesticide exposures in California's central valley.

BACKGROUND: Pyrethroid pesticide use is increasing worldwide, although the full extent of associated health effects is unknown. An epigenome-wide association study (EWAS) with exploratory pathway analysis may help identify potential pyrethroid-related health effects. METHODS: We performed an exploratory EWAS of chronic ambient pyrethroid exposure using control participants' blood in the Parkinson's Environment and Genes Study in the Central Valley of California (N = 237). We estimated associations of living and working near agricultural pyrethroid pesticide applications in the past 5 years (binary) with site-specific differential methylation, and used a false discovery rate (FDR) cut off of 0.05 for significance. We controlled for age, sex, education, cell count, and an ancestral marker for Hispanic ethnicity. We normalized methylation values for Type I/II probe bias using Beta-Mixture Quantile (BMIQ) normalization, filtered out cross-reactive probes, and evaluated for remaining bias with Surrogate Variable Analysis (SVA). We also evaluated the effects of controlling for cell count and normalizing for Type I/II probe bias by comparing changes in effect estimates and p-values for the top hits across BMIQ and GenomeStudio normalization methods, and controlling for cell count. To facilitate broader interpretation, we annotated genes to the CpG sites and performed gene set overrepresentation analysis, using genes annotated to CpG sites that were associated with pyrethroids at a raw p < 0.05, and controlling for background representation of CpG sites on the chip. We did this for both a biological process context (Gene Ontology terms) using missMethyl, and a disease set context using WebGestalt. For these gene set overrepresentation analyses we also used an FDR cut off of 0.05 for significance of gene sets. RESULTS: After controlling for cell count and applying BMIQ normalization, 4 CpG sites were differentially methylated in relation to pyrethroid exposures. When using GenomeStudio's Illumina normalization, 415 CpG sites were differentially methylated, including all four identified with the BMIQ method. In the gene set overrepresentation analyses, we identified 6 GO terms using BMIQ normalization, and 76 using Illumina normalization, including the 6 identified by BMIQ. For disease sets, we identified signals for Alzheimer's disease, leukemia and several other cancers, diabetes, birth defects, and other diseases, for both normalization methods. We identified minimal changes in effect estimates after controlling for cell count, and controlling for cell count generally weakened p-values. BMIQ normalization, however, resulted in different beta coefficients and weakened p-values. CONCLUSIONS: Chronic ambient pyrethroid exposure is associated with differential methylation at CpG sites that annotate to a wide variety of disease states and biological mechanisms that align with prior research. However, this EWAS also implicates several novel diseases for future investigation, and highlights the relative importance of different background normalization methods in identifying associations.

Prenatal exposure to pyrethroid pesticides and childhood behavior and executive functioning.

Several previous studies of pyrethroid biomarkers and behavior have reported associations between concurrent pyrethroid levels and adverse behavioral problems in children. One geospatial study reported associations between prenatal exposure to pyrethroids and autism. However, the association between prenatal pyrethroid biomarkers and childhood behavior is unknown. The Mount Sinai Children's Environmental Health Center is a prospective birth cohort with urinary pyrethroid biomarkers during pregnancy and behavioral measurements at 4, 6, and 7-9 years of age. Primiparous women were enrolled between 1998 and 2002. 162 mother/child pairs with complete exposure and behavioral outcomes data were used to investigate associations between detectable levels of prenatal pyrethroid metabolites and scores on the Behavioral Assessment System for Children and the Behavior Rating Inventory of Executive Function. Overall, detection frequencies of pyrethroid metabolites were low (<30%). In longitudinal mixed models, detectable levels of 3-PBA during pregnancy were associated with worse Internalizing (ß -4.50, 95% CI -8.05, -0.95), Depression (ß -3.21, 95% CI -6.38, -0.05), Somatization (ß -3.22, 95% CI -6.38, -0.06), Behavioral Regulation (ß -3.59, 95% CI -6.97, -0.21), Emotional Control (ß -3.35, 95% CI -6.58, -0.12), Shifting (ß -3.42, 95% CI -6.73, -0.11), and Monitoring (ß -4.08, 95% CI -7.07, -1.08) scales. Detectable levels of cis-DCCA were associated with worse Externalizing (ß -4.74, 95% CI -9.37, -0.10), Conduct Problems (ß -5.35, 95% CI -9.90, -0.81), Behavioral Regulation (ß -6.42, 95% CI -11.39, -1.45), and Inhibitory Control (ß -7.20, 95% CI -12.00, -2.39). Although detection frequencies of pyrethroid metabolites were low, we found suggestive evidence that prenatal exposure to 3-PBA and cis-DCCA may be associated with a variety of behavioral and executive functioning deficits.

Prenatal exposure to organophosphorus pesticides and childhood neurodevelopmental phenotypes.

Prenatal exposure to organophosphorus pesticides (OPs) has been associated with different neurodevelopmental outcomes across different cohorts. A phenotypic approach may address some of these differences by incorporating information across scales and accounting for the complex correlational structure of neurodevelopmental outcomes. Additionally, Bayesian hierarchical modeling can account for confounding by collinear co-exposures. We use this framework to examine associations between prenatal exposure to OPs and behavior, executive functioning, and IQ assessed at age 6-9 years in a cohort of 404 mother/infant pairs recruited during pregnancy. We derived phenotypes of neurodevelopment with a factor analysis, and estimated associations between OP metabolites and these phenotypes in Bayesian hierarchical models for exposure mixtures. We report seven factors: 1) Impulsivity and Externalizing, 2) Executive Functioning, 3) Internalizing, 4) Perceptual Reasoning, 5) Adaptability, 6) Processing Speed, and 7) Verbal Intelligence. These, along with the Working Memory Index, were standardized and scaled so that positive values reflected positive attributes and negative values represented adverse outcomes. Standardized dimethylphosphate metabolites were negatively associated with Internalizing factor scores (ß^ - 0.13, 95% CI - 0.26, 0.00) but positively associated with Executive Functioning factor scores (ß^ 0.18, 95% CI 0.04, 0.31). Standardized diethylphosphate metabolites were negatively associated with the Working Memory Index (ß^ - 0.17, 95% CI - 0.33, - 0.03). Associations with factor scores were generally stronger and more precise than associations with individual instrument-specific items. Factor analysis of outcomes may provide some advantages in etiological studies of childhood neurodevelopment by incorporating information across scales to reduce dimensionality and improve precision.

Prenatal exposure to organophosphate pesticides and reciprocal social behavior in childhood.

Prenatal exposure to organophosphate pesticides (OPs) has been associated with adverse neurodevelopmental outcomes in childhood, including low IQ, pervasive developmental disorder (PDD), attention problems and ADHD. Many of these disorders involve impairments in social functioning. Thus, we investigated the relationship between biomarkers of prenatal OP exposure and impaired reciprocal social behavior in childhood, as measured by the Social Responsiveness Scale (SRS). Using a multi-ethnic urban prospective cohort of mother-infant pairs in New York City recruited between 1998 and 2002 (n=404) we examined the relation between third trimester maternal urinary levels of dialkylphosphate (ΣDAP) OP metabolites and SRS scores among 136 children who returned for the 7-9year visit. Overall, there was no association between OPs and SRS scores, although in multivariate adjusted models, associations were heterogeneous by race and by sex. Among blacks, each 10-fold increase in total diethylphosphates (ΣDEP) was associated with poorer social responsiveness (ß=5.1 points, 95% confidence interval (CI) 0.8, 9.4). There was no association among whites or Hispanics, or for total ΣDAP or total dimethylphosphate (ΣDMP) biomarker levels. Additionally, stratum-specific models supported a stronger negative association among boys for ΣDEPs (ß=3.5 points, 95% CI 0.2, 6.8), with no notable association among girls. Our results support an association of prenatal OP exposure with deficits in social functioning among blacks and among boys, although this may be in part reflective of differences in exposure patterns.