Search

Per- and polyfluoroalkyl substances (PFAS) are a widespread and persistent class of contaminants posing significant environmental and human health concerns. Comprehensive understanding of the modes of action underlying toxicity among structurally diverse PFAS is mostly lacking. To address this need, we recently reported on our application of developing zebrafish to evaluate a large library of PFAS for developmental toxicity. In the present study, we prioritized 15 bioactive PFAS that induced significant morphological effects and performed RNA-sequencing to characterize early transcriptional responses at a single timepoint (48 h post fertilization) after early developmental exposures (8 h post fertilization). Internal concentrations of 5 of the 15 PFAS were measured from pooled whole fish samples across multiple timepoints between 24-120 h post fertilization, and additional temporal transcriptomics at several timepoints (48-96 h post fertilization) were conducted for Nafion byproduct 2. A broad range of differentially expressed gene counts were identified across the PFAS exposures. Most PFAS that elicited robust transcriptomic changes affected biological processes of the brain and nervous system development. While PFAS disrupted unique processes, we also found that similarities in some functional head groups of PFAS were associated with the disruption in expression of similar gene sets. Body burdens after early developmental exposures to select sulfonic acid PFAS, including Nafion byproduct 2, increased from the 24-96 h post fertilization sampling timepoints and were greater than those of sulfonamide PFAS of similar chain lengths. In parallel, the Nafion byproduct 2-induced transcriptional responses increased between 48 and 96 h post fertilization. PFAS characteristics based on toxicity, transcriptomic effects, and modes of action will contribute to further prioritization of PFAS structures for testing and informed hazard assessment., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Rericha, St. Mary, Truong, McClure, Martin, Leonard, Thunga, Simonich, Waters, Field and Tanguay.)

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 h post-fertilization after exposure to eight retene concentrations. We aimed to identify the most sensitive transcriptomic responses and to determine whether this approach could uncover gene sets uniquely differentially expressed at concentrations which induce a phenotype. We identified a concentration-response relationship for differential gene expression in both number of differentially expressed genes (DEGs) and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-β (TGF-β) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module seven genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module seven genes are largely responsible for retene toxicity., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Wilson, McClure, Waters, Simonich and Tanguay.)