Signaling Events Downstream of AHR Activation That Contribute to Toxic Responses: The Functional Role of an AHR-Dependent Long Noncoding RNA (slincR) Using the Zebrafish Model

Leveraging integrative toxicogenomic approach towards development of stressor-centric adverse outcome pathway networks for plastic additives

Plastics are widespread pollutants found in atmospheric, terrestrial and aquatic ecosystems due to their extensive usage and environmental persistence. Plastic additives, that are intentionally added to achieve specific functionality in plastics, leach into the environment upon plastic degradation and pose considerable risk to ecological and human health. Limited knowledge concerning the presence of plastic additives throughout plastic life cycle has hindered their effective regulation, thereby posing risks to product safety. In this study, we leveraged the adverse outcome pathway (AOP) framework to understand the mechanisms underlying plastic additives-induced toxicities. We first identified an exhaustive list of 6470 plastic additives from chemicals documented in plastics. Next, we leveraged heterogenous toxicogenomics and biological endpoints data from five exposome-relevant resources, and identified associations between 1287 plastic additives and 322 complete and high quality AOPs within AOP-Wiki. Based on these plastic additive-AOP associations, we constructed a stressor-centric AOP network, wherein the stressors are categorized into ten priority use sectors and AOPs are linked to 27 disease categories. We visualized the plastic additives-AOP network for each of the 1287 plastic additives and made them available in a dedicated website: https://cb.imsc.res.in/saopadditives/ . Finally, we showed the utility of the constructed plastic additives-AOP network by identifying highly relevant AOPs associated with benzo[a]pyrene (B[a]P), bisphenol A (BPA), and bis(2-ethylhexyl) phthalate (DEHP) and thereafter, explored the associated toxicity pathways in humans and aquatic species. Overall, the constructed plastic additives-AOP network will assist regulatory risk assessment of plastic additives, thereby contributing towards a toxic-free circular economy for plastics.

Molecular mechanisms of environmental pollutant-induced cartilage damage: from developmental disorders to osteoarthritis

The objective of the present study was to review the molecular mechanisms of the adverse effects of environmental pollutants on chondrocytes and extracellular matrix (ECM). Existing data demonstrate that both heavy metals, including cadmium (Cd), lead (Pb), and arsenic (As), as well as organic pollutants, including polychlorinated dioxins and furans (PCDD/Fs) and polychlorinated biphenyls (PCB), bisphenol A, phthalates, polycyclic aromatic hydrocarbons (PAH), pesticides, and certain other organic pollutants that target cartilage ontogeny and functioning. Overall, environmental pollutants reduce chondrocyte viability through the induction apoptosis, senescence, and inflammatory response, resulting in cell death and impaired ECM production. The effects of organic pollutants on chondrocyte development and viability were shown to be mediated by binding to the aryl hydrocarbon receptor (AhR) signaling and modulation of non-coding RNA expression. Adverse effects of pollutant exposures were observed in articular and growth plate chondrocytes. These mechanisms also damage chondrocyte precursors and subsequently hinder cartilage development. In addition, pollutant exposure was shown to impair chondrogenesis by inhibiting the expression of Sox9 and other regulators. Along with altered Runx2 signaling, these effects also contribute to impaired chondrocyte hypertrophy and chondrocyte-to-osteoblast trans-differentiation, resulting in altered endochondral ossification. Several organic pollutants including PCDD/Fs, PCBs and PAHs, were shown to induce transgenerational adverse effects on cartilage development and the resulting skeletal deformities. Despite of epidemiological evidence linking human environmental pollutant exposure to osteoarthritis or other cartilage pathologies, the data on the molecular mechanisms of adverse effects of environmental pollutant exposure on cartilage tissue were obtained from studies in laboratory rodents, fish, or cell cultures and should be carefully extrapolated to humans, although they clearly demonstrate that cartilage should be considered a putative target for environmental pollutant toxicity.

Exposure to the aryl hydrocarbon receptor agonist dioxin disrupts formation of the muscle, nerves, and vasculature in the developing jaw

Human exposure to environmental pollutants can disrupt embryonic development and impact juvenile and adult health outcomes by adversely affecting cell and organ function. Notwithstanding, environmental contamination continues to increase due to industrial development, insufficient regulations, and the mobilization of pollutants as a result of extreme weather events. Dioxins are a class of structurally related persistent organic pollutants that are highly toxic, carcinogenic, and teratogenic. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent dioxin compound and has been shown to induce toxic effects in developing organisms by activating the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor targeted by multiple persistent organic pollutants. Contaminant-induced AHR activation results in malformations of the craniofacial cartilages and neurocranium; however, the mechanisms mediating these phenotypes are not well understood. In this study, we utilized the optically transparent zebrafish model to elucidate novel cellular targets and potential transcriptional targets underlying TCDD-induced craniofacial malformations. To this end, we exposed zebrafish embryos at 4 h post fertilization to TCDD and employed a mixed-methods approach utilizing immunohistochemistry staining, transgenic reporter lines, fixed and in vivo confocal imaging, and timelapse microscopy to determine the targets mediating TCDD-induced craniofacial phenotypes. Our data indicate that embryonic TCDD exposure reduced jaw and pharyngeal arch Sox10+ chondrocytes and Tcf21+ pharyngeal mesoderm progenitors. Exposure to TCDD correspondingly led to a reduction in collagen type II deposition in Sox10+ domains. Embryonic TCDD exposure impaired development of tissues derived from or guided by Tcf21+ progenitors, namely: nerves, muscle, and vasculature. Specifically, TCDD exposure disrupted development of the hyoid and mandibular arch muscles, decreased neural innervation of the jaw, resulted in compression of cranial nerves V and VII, and led to jaw vasculature malformations. Collectively, these findings reveal novel structural targets and potential transcriptional targets of TCDD-induced toxicity, showcasing how contaminant exposures lead to congenital craniofacial malformations.

AOP Report: Aryl Hydrocarbon Receptor Activation Leads to Early-Life Stage Mortality via Sox9 Repression-Induced Craniofacial and Cardiac Malformations

The aryl hydrocarbon receptors (Ahrs) are evolutionarily conserved ligand-dependent transcription factors that are activated by structurally diverse endogenous compounds as well as environmental chemicals such as polycyclic aromatic hydrocarbons and halogenated aromatic hydrocarbons. Activation of the Ahr leads to several transcriptional changes that can cause developmental toxicity resulting in mortality. Evidence was assembled and evaluated for two novel adverse outcome pathways (AOPs) which describe how Ahr activation (molecular initiating event) can lead to early-life stage mortality (adverse outcome), via either SOX9-mediated craniofacial malformations (AOP 455) or cardiovascular toxicity (AOP 456). Using a key event relationship (KER)-by-KER approach, we collected evidence using both a narrative search and a systematic review based on detailed search terms. Weight of evidence for each KER was assessed to inform overall confidence of the AOPs. The AOPs link to previous descriptions of Ahr activation and connect them to two novel key events (KEs), increase in slincR expression, a newly characterized long noncoding RNA with regulatory functions, and suppression of SOX9, a critical transcription factor implicated in chondrogenesis and cardiac development. In general, confidence levels for KERs ranged between medium and strong, with few inconsistencies, as well as several opportunities for future research identified. While the majority of KEs have only been demonstrated in zebrafish with 2,3,7,8-tetrachlorodibenzo-p-dioxin as an Ahr activator, evidence suggests that the two AOPs likely apply to most vertebrates and many Ahr-activating chemicals. Addition of the AOPs into the AOP-Wiki (https://aopwiki.org/) helps expand the growing Ahr-related AOP network to 19 individual AOPs, of which six are endorsed or in progress and the remaining 13 relatively underdeveloped. Environ Toxicol Chem 2023;42:2063-2077. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

A CRISPR-Cas9 mutation in sox9b long intergenic noncoding RNA (slincR) affects zebrafish development, behavior, and regeneration

The role of long noncoding RNAs (lncRNAs) regulators of toxicological responses to environmental chemicals is gaining prominence. Previously, our laboratory discovered an lncRNA, sox9b long intergenic noncoding RNA (slincR), that is activated by multiple ligands of aryl hydrocarbon receptor (AHR). Within this study, we designed a CRISPR-Cas9-mediated slincR zebrafish mutant line to better understand its biological function in presence or absence of a model AHR ligand, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). The slincRosu3 line contains an 18 bp insertion within the slincR sequence that changes its predicted mRNA secondary structure. Toxicological profiling showed that slincRosu3 is equally or more sensitive to TCDD for morphological and behavioral phenotypes. Embryonic mRNA-sequencing showed differential responses of 499 or 908 genes in slincRosu3 in absence or presence of TCDD Specifically, unexposed slincRosu3 embryos showed disruptions in metabolic pathways, suggesting an endogenous role for slincR. slincRosu3 embryos also had repressed mRNA levels of sox9b-a transcription factor that slincR is known to negatively regulate. Hence, we studied cartilage development and regenerative capacity-both processes partially regulated by sox9b. Cartilage development was disrupted in slincRosu3 embryos both in presence and absence of TCDD. slincRosu3 embryos also displayed a lack of regenerative capacity of amputated tail fins, accompanied by a lack of cell proliferation. In summary, using a novel slincR mutant line, we show that a mutation in slincR can have widespread impacts on gene expression and structural development endogenously and limited, but significant impacts in presence of AHR induction that further highlights its importance in the developmental process.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and subchronic (4 wk) TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of lncRNAs. TCDD dysregulated >4000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 lncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cell-cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolism. Networks were validated by the striking enrichments that predicted regulatory lncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

Polycyclic aromatic hydrocarbons exposure and plasma lncRNA signature: A profile and functional analysis

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental pollutants that pose detrimental effects on human health, and the exploration of the associations of PAHs exposure with long non-coding RNA (lncRNA) may provide novel clues to the underlying mechanisms. In the present study, we detected 10 urinary PAHs metabolites by GC-MS and plasma lncRNAs levels by Human LncRNA Array v4 among 230 participants from two panels (160 in the Shiyan panel and 70 in the Wuhan-Zhuhai panel). We applied linear regression models to assess the associations between PAHs metabolites and lncRNAs separately in each panel and combined the results using fixed-effect meta-analysis. To explore the potential origin of PAHs-related lncRNAs in plasma, we estimated their tissue-specificity and associations between lncRNAs levels in plasma and leukocytes. Leukocytes mRNA sequencing data and RNA binding proteins were utilized to explore implicated pathways of identified lncRNAs. We found that urinary 1-hydroxyphenanthrene (1-OH-Phe) was inversely associated with 8 lncRNAs and positively associated with 1 lncRNA, as well as 9-hydroxyphenanthrene (9-OH-Phe) was inversely associated with 11 lncRNAs (FDR < 0.1). Tissue specificity analysis using Genome Tissue Expression database suggested that several identified lncRNAs might specifically express in organs targeted by PAHs exposure (lung, liver, heart, kidney, and brain). Besides, plasma levels of 1-OH-Phe related ENSG00000260616 and 9-OH-Phe related STARD4-AS1 were inversely associated with their intra-leukocytes levels (P value < 0.05). Notably, STARD4-AS1 was positively associated with the expression levels of its neighboring protein-coding gene (CAMK4 and STARD4) in leukocytes and were involved in pathways related to cellular response to DNA damage, which we further confirmed using DNA damage biomarker, 8-hydroxydeoxyguanosine. Functional analysis also revealed vital pathways related to cytokine-mediated signaling and glucose homeostasis. Our findings provided novel insights into plausible biological mechanisms underlying the adverse effects of PAHs exposure.

Exposure to the persistent organic pollutant 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD, dioxin) disrupts development of the zebrafish inner ear

Dioxins are a class of highly toxic and persistent environmental pollutants that have been shown through epidemiological and laboratory-based studies to act as developmental teratogens. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the most potent dioxin congener, has a high affinity for the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. TCDD-induced AHR activation during development impairs nervous system, cardiac, and craniofacial development. Despite the robust phenotypes previously reported, the characterization of developmental malformations and our understanding of the molecular targets mediating TCDD-induced developmental toxicity remains limited. In zebrafish, TCDD-induced craniofacial malformations are produced, in part, by the downregulation of SRY-box transcription factor 9b ( sox9b ), a member of the SoxE gene family. sox9b , along with fellow SoxE gene family members sox9a and sox10 , have important functions in the development of the otic placode, the otic vesicle, and, ultimately, the inner ear. Given that sox9b in a known target of TCDD and that transcriptional interactions exist among SoxE genes, we asked whether TCDD exposure impaired the development of the zebrafish auditory system, specifically the otic vesicle, which gives rise to the sensory components of the inner ear. Using immunohistochemistry, in vivo confocal imaging, and time-lapse microscopy, we assessed the impact of TCDD exposure on zebrafish otic vesicle development. We found exposure resulted in structural deficits, including incomplete pillar fusion and altered pillar topography, leading to defective semicircular canal development. The observed structural deficits were accompanied by reduced collagen type II expression in the ear. Together, our findings reveal the otic vesicle as a novel target of TCDD-induced toxicity, suggest that the function of multiple SoxE genes may be affected by TCDD exposure, and provide insight into how environmental contaminants contribute to congenital malformations.

Highlights

The zebrafish ear is necessary to detect changes in motion, sound, and gravity.Embryos exposed to TCDD lack structural components of the developing ear.TCDD exposure impairs formation of the fusion plate and alters pillar topography.The semicircular canals of the ear are required to detect changes in movement.Following TCDD exposure embryos fail to establish semicircular canals.

SLAMR, a synaptically targeted lncRNA, facilitates the consolidation of contextual fear memory

LncRNAs are involved in critical processes for cell homeostasis and function. However, it remains largely unknown whether and how the transcriptional regulation of long noncoding RNAs results in activity-dependent changes at the synapse and facilitate formation of long-term memories. Here, we report the identification of a novel lncRNA, SLAMR, that becomes enriched in CA1- but not in CA3-hippocampal neurons upon contextual fear conditioning. SLAMR is transported to dendrites via the molecular motor KIF5C and recruited to the synapse in response to stimulation. Loss of function of SLAMR reduced dendritic complexity and impaired activity dependent changes in spine structural plasticity. Interestingly, gain of function of SLAMR enhanced dendritic complexity, and spine density through enhanced translation. Analyses of the SLAMR interactome revealed its association with CaMKIIα protein through a 220-nucleotide element and its modulation of CaMKIIα activity. Furthermore, loss-of-function of SLAMR in CA1 selectively impairs consolidation but neither acquisition, recall, nor extinction of fear memory and spatial memory. Together, these results establish a new mechanism for activity dependent changes at the synapse and consolidation of contextual fear.

Environmentally relevant uptake, elimination, and metabolic changes following early embryonic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin in zebrafish

Dioxin and dioxin-like compounds are ubiquitous environmental contaminants that induce toxicity by binding to the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. The zebrafish model has been used to define the developmental toxicity observed following exposure to exogenous AHR ligands such as the potent agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (dioxin, TCDD). While the model has successfully identified cellular targets of TCDD and molecular mechanisms mediating TCDD-induced phenotypes, fundamental information such as the body burden produced by standard exposure models is still unknown. We performed targeted gas chromatography (GC) high-resolution mass spectrometry (HRMS) in tandem with non-targeted liquid chromatography (LC) HRMS to quantify TCDD uptake, model the elimination dynamics of TCDD, and determine how TCDD exposure affects the zebrafish metabolome. We found that 50 ppt, 10 ppb, and 1 ppb waterborne exposures to TCDD during early embryogenesis produced environmentally relevant body burdens: 38 ± 4.34, 26.6 ± 1.2, and 8.53 ± 0.341 pg/embryo, respectively, at 24 hours post fertilization. TCDD exposure was associated with the dysregulation of metabolic pathways that are associated with the AHR signaling pathway as well as pathways shown to be affected in mammals following TCDD exposure. In addition, we discovered that TCDD exposure affected several metabolic pathways that are critical for brain development and function including glutamate metabolism, chondroitin sulfate biosynthesis, and tyrosine metabolism. Together, these data demonstrate that existing exposure methods produce environmentally relevant body burdens of TCDD in zebrafish and provide insight into the biochemical pathways impacted by toxicant-induced AHR activation.

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

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.

Polycyclic aromatic hydrocarbons, long non-coding RNA expression, and DNA damage in coke oven workers

Exposure to polycyclic aromatic hydrocarbons (PAHs) was associated with DNA damage, while the roles of long non-coding RNAs (lncRNAs) in the associations were unclear. We aimed to assess the association of lncRNA NR_024564 with urinary monohydroxy PAHs (OH-PAHs) and DNA damage among 332 coke oven workers. We determined 12 OH-PAHs by gas chromatography-mass spectrometry, and the expression level of NR_024564 by droplet digital RT-PCR and DNA damage by the comet assay. In total participants, we found that NR_024564 was not significantly associated with OH-PAHs or comet parameters. However, among workers with ≥ 20 working years, multiple OH-PAHs including urinary 1-hydroxyphenanthrene (1-OHPh), 2-OHPh, 3-OHPh, 9-OHPh, 1‑hydroxypyrene, and total PAH metabolites were related to increased comet parameters. Moreover, NR_024564 was significantly associated with 2-OHPh and four comet parameters. Each 1% increase in 2-OHPh was associated with 0.35% reduction (95% CI: 0.16%, 0.55%) in NR_024564 (P-_FDR = 0.005), and 2-OHPh was marginally interacted with working years in relation to NR_024564 decrease. Also, each 1% increment of NR_024564 was related to 0.04-0.13% decrease of Olive tail moment, percent DNA in the comet tail, tail length, and tail moment (all P-_FDR < 0.05). Furthermore, low NR_024564 level combined with high levels of 1-OHPh and 2-OHPh or ≥ 20 working years was positively associated with the comet parameters among the total participants. Our results indicated that NR_024564 might be linked to the adverse associations of PAHs with the DNA damage of coke oven workers who worked for ≥ 20 years.

The Ahr2-Dependent wfikkn1 Gene Influences Zebrafish Transcriptome, Proteome, and Behavior

The aryl hydrocarbon receptor (AHR) is required for vertebrate development and is also activated by exogenous chemicals, including polycyclic aromatic hydrocarbons (PAHs) and TCDD. AHR activation is well-understood, but roles of downstream molecular signaling events are largely unknown. From previous transcriptomics in 48-hours post fertilization (hpf) zebrafish exposed to several PAHs and TCDD, we found wfikkn1 was highly co-expressed with cyp1a (marker for AHR activation). Thus, we hypothesized wfikkn1's role in AHR signaling, and showed that wfikkn1 expression was Ahr2 (zebrafish ortholog of human AHR)-dependent in developing zebrafish exposed to TCDD. To functionally characterize wfikkn1, we made a CRISPR-Cas9 mutant line with a 16-bp deletion in wfikkn1's exon, and exposed wildtype and mutants to DMSO or TCDD. 48-hpf mRNA sequencing revealed over 700 genes that were differentially expressed (p < 0.05, log2FC > 1) between each pair of treatment combinations, suggesting an important role for wfikkn1 in altering both the 48-hpf transcriptome and TCDD-induced expression changes. Mass spectrometry-based proteomics of 48-hpf wildtype and mutants revealed 325 significant differentially expressed proteins. Functional enrichment demonstrated wfikkn1 was involved in skeletal muscle development and played a role in neurological pathways after TCDD exposure. Mutant zebrafish appeared morphologically normal but had significant behavior deficiencies at all life stages, and absence of Wfikkn1 did not significantly alter TCDD-induced behavior effects at all life stages. In conclusion, wfikkn1 did not appear to be significantly involved in TCDD's overt toxicity but is likely a necessary functional member of the AHR signaling cascade.

Analysis of Gene-Environment Interactions Related to Developmental Disorders

Various genetic and environmental factors are associated with developmental disorders (DDs). It has been suggested that interaction between genetic and environmental factors (G × E) is involved in the etiology of DDs. There are two major approaches to analyze the interaction: genome-wide and candidate gene-based approaches. In this mini-review, we demonstrate how these approaches can be applied to reveal the G × E related to DDs focusing on zebrafish and mouse models. We also discuss novel approaches to analyze the G × E associated with DDs.

Gene co-expression network analysis in zebrafish reveals chemical class specific modules

BACKGROUND

Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response.

RESULTS

Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure.

CONCLUSIONS

Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

Polycyclic aromatic hydrocarbons associated long non-coding RNAs and heart rate variability in coke oven workers

Epidemiological studies have showed that polycyclic aromatic hydrocarbons (PAHs) were associated with heart rate variability (HRV), but the role of long non-coding RNAs (lncRNAs) in the association is unknown. We aimed to identify PAHs-related lncRNAs and assess their associations with HRV among coke oven workers. Differential lncRNAs expression between 12 exposed workers and 12 controls was tested by Human 8X60k LncRNA Arrays in discovery stage, then selected NR_024564 was validated in 353 workers using droplet digital RT-PCR. Microarray results showed that 1234 lncRNAs were downregulated with 805 lncRNAs upregulated in exposed group (≥ 2-fold change). In validation stage, no significant association was observed between NR_024564 and PAH exposure or HRV in total subjects, while urinary 2-hydroxyfluorene (2-OHFlu) was inversely related to root mean square successive difference (RMSSD). However, in current smokers, NR_024564 was inversely related to urinary 2-OHFlu, 2-hydroxyphenanthrene, 1-hydroxypyrene (1-OHP), and total PAHs metabolites (ΣOH-PAHs), of which 1-OHP accounted for the strongest estimation for interaction with smoking status (P_interaction = 0.011). Also, the positive associations of NR_024564 with RMSSD and high frequency power showed an interaction with smoking status (P_interaction = 0.034 and 0.023, respectively). Also, urinary 2-OHFlu and ΣOH-PAHs were inversely associated with RMSSD in current smokers. In addition, elevated NR_024564 was dose-responsive related to increased RMSSD in above high-PAHs groups among smokers (all P_trend < 0.05). Our results revealed that NR_024564 and its interactions with smoking status might act as novel mechanisms regulating the adverse effects of PAHs on HRV.

A mechanism linking perinatal 2,3,7,8 tetrachlorodibenzo-p-dioxin exposure to lower urinary tract dysfunction in adulthood

Benign prostatic hyperplasia/lower urinary tract dysfunction (LUTD) affects nearly all men. Symptoms typically present in the fifth or sixth decade and progressively worsen over the remainder of life. Here, we identify a surprising origin of this disease that traces back to the intrauterine environment of the developing male, challenging paradigms about when this disease process begins. We delivered a single dose of a widespread environmental contaminant present in the serum of most Americans [2,3,7,8 tetrachlorodibenzo-p-dioxin (TCDD), 1 µg/kg], and representative of a broader class of environmental contaminants, to pregnant mice and observed an increase in the abundance of a neurotrophic factor, artemin, in the developing mouse prostate. Artemin is required for noradrenergic axon recruitment across multiple tissues, and TCDD rapidly increases prostatic noradrenergic axon density in the male fetus. The hyperinnervation persists into adulthood, when it is coupled to autonomic hyperactivity of prostatic smooth muscle and abnormal urinary function, including increased urinary frequency. We offer new evidence that prostate neuroanatomical development is malleable and that intrauterine chemical exposures can permanently reprogram prostate neuromuscular function to cause male LUTD in adulthood.

Summary: We describe a new mechanism of benign prostate disease, initiated by fetal chemical exposure, which durably increases prostatic noradrenergic axon density and causes smooth muscle hyperactivity and urinary voiding dysfunction.

Phenotypically Anchored mRNA and miRNA Expression Profiling in Zebrafish Reveals Flame Retardant Chemical Toxicity Networks

The ubiquitous use of flame retardant chemicals (FRCs) in the manufacture of many consumer products leads to inevitable environmental releases and human exposures. Studying toxic effects of FRCs as a group is challenging since they widely differ in physicochemical properties. We previously used zebrafish as a model to screen 61 representative FRCs and showed that many induced behavioral and teratogenic effects, with aryl phosphates identified as the most active. In this study, we selected 10 FRCs belonging to diverse physicochemical classes and zebrafish toxicity profiles to identify the gene expression responses following exposures. For each FRC, we executed paired mRNA-micro-RNA (miR) sequencing, which enabled us to study mRNA expression patterns and investigate the role of miRs as posttranscriptional regulators of gene expression. We found widespread disruption of mRNA and miR expression across several FRCs. Neurodevelopment was a key disrupted biological process across multiple FRCs and was corroborated by behavioral deficits. Several mRNAs (e.g., osbpl2a) and miRs (e.g., mir-125b-5p), showed differential expression common to multiple FRCs (10 and 7 respectively). These common miRs were also predicted to regulate a network of differentially expressed genes with diverse functions, including apoptosis, neurodevelopment, lipid regulation and inflammation. Commonly disrupted transcription factors (TFs) such as retinoic acid receptor, retinoid X receptor, and vitamin D regulator were predicted to regulate a wide network of differentially expressed mRNAs across a majority of the FRCs. Many of the differential mRNA-TF and mRNA-miR pairs were predicted to play important roles in development as well as cancer signaling. Specific comparisons between TBBPA and its derivative TBBPA-DBPE showed contrasting gene expression patterns that corroborated with their phenotypic profiles. The newer generation FRCs such as IPP and TCEP produced distinct gene expression changes compared to the legacy FRC BDE-47. Our study is the first to establish a mRNA-miR-TF regulatory network across a large group of structurally diverse FRCs and diverse phenotypic responses. The purpose was to discover common and unique biological targets that will help us understand mechanisms of action for these important chemicals and establish this approach as an important tool for better understanding toxic effects of environmental contaminants.

A Review of the Functional Roles of the Zebrafish Aryl Hydrocarbon Receptors

Over the last 2 decades, the zebrafish (Danio rerio) has emerged as a stellar model for unraveling molecular signaling events mediated by the aryl hydrocarbon receptor (AHR), an important ligand-activated receptor found in all eumetazoan animals. Zebrafish have 3 AHRs-AHR1a, AHR1b, and AHR2, and studies have demonstrated the diversity of both the endogenous and toxicological functions of the zebrafish AHRs. In this contemporary review, we first highlight the evolution of the zebrafish ahr genes, and the characteristics of the receptors including developmental and adult expression, their endogenous and inducible roles, and the predicted ligands from homology modeling studies. We then review the toxicity of a broad spectrum of AHR ligands across multiple life stages (early stage, and adult), discuss their transcriptomic and epigenetic mechanisms of action, and report on any known interactions between the AHRs and other signaling pathways. Through this article, we summarize the promising research that furthers our understanding of the complex AHR pathway through the extensive use of zebrafish as a model, coupled with a large array of molecular techniques. As much of the research has focused on the functions of AHR2 during development and the mechanism of TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) toxicity, we illustrate the need to address the considerable knowledge gap in our understanding of both the mechanistic roles of AHR1a and AHR1b, and the diverse modes of toxicity of the various AHR ligands.

Aryl Hydrocarbon Receptor Mediates Larval Zebrafish Fin Duplication Following Exposure to Benzofluoranthenes

The aryl hydrocarbon receptor (AHR) mediates developmental toxicity of several xenobiotic classes including polycyclic aromatic hydrocarbons. Using embryonic zebrafish, we previously identified 4 polycyclic aromatic hydrocarbons that caused a novel phenotype among AHR ligands-growth of a lateral, duplicate caudal fin fold. The window of sensitivity to the most potent inducer of this phenotype, benzo[k]fluoranthene (BkF), was prior to 36 h postfertilization (hpf), although the phenotype was not manifest until 60 hpf. AHR dependency via Ahr2 was demonstrated using morpholino knockdown. Hepatocyte ablation demonstrated that hepatic metabolism of BkF was not required for the phenotype, nor was it responsible for the window of sensitivity. RNA sequencing performed on caudal trunk tissue from BkF-exposed animals collected at 48, 60, 72, and 96 hpf showed upregulation of genes associated with AHR activation, appendage development, and tissue patterning. Genes encoding fibroblast growth factor and bone morphogenic protein ligands, along with retinaldehyde dehydrogenase, were prominently upregulated. Gene Ontology term analysis revealed that upregulated genes were enriched for mesoderm development and fin regeneration, whereas downregulated genes were enriched for Wnt signaling and neuronal development. MetaCore (Clarivate Analytics) systems analysis of orthologous human genes predicted that R-SMADs, AP-1, and LEF1 regulated the expression of an enriched number of gene targets across all time points. Our results demonstrate a novel aspect of AHR activity with implications for developmental processes conserved across vertebrate species.

Using Zebrafish to Analyze the Genetic and Environmental Etiologies of Congenital Heart Defects

Congenital heart defects (CHDs) are among the most common human birth defects. However, the etiology of a large proportion of CHDs remains undefined. Studies identifying the molecular and cellular mechanisms that underlie cardiac development have been critical to elucidating the origin of CHDs. Building upon this knowledge to understand the pathogenesis of CHDs requires examining how genetic or environmental stress changes normal cardiac development. Due to strong molecular conservation to humans and unique technical advantages, studies using zebrafish have elucidated both fundamental principles of cardiac development and have been used to create cardiac disease models. In this chapter we examine the unique toolset available to zebrafish researchers and how those tools are used to interrogate the genetic and environmental contributions to CHDs.

Casting a wide net: use of diverse model organisms to advance toxicology

Toxicology - the study of how chemicals interact with biological systems - has clear relevance to human health and disease. Persistent exposure to natural and synthetic chemicals is an unavoidable part of living on our planet; yet, we understand very little about the effects of exposure to the vast majority of chemicals. While epidemiological studies can provide strong statistical inference linking chemical exposure to disease, research in model systems is essential to elucidate the mechanisms of action and to predict outcomes. Most research in toxicology utilizes a handful of mammalian models that represent a few distinct branches of the evolutionary tree. This narrow focus constrains the understanding of chemical-induced disease processes and systems that have evolved in response to exposures. We advocate for casting a wider net in environmental toxicology research to utilize diverse model systems, including zebrafish, and perform more mechanistic studies of cellular responses to chemical exposures to shift the perception of toxicology as an applied science to that of a basic science. This more-inclusive perspective will enrich the field and should remain central to research on chemical-induced disease.

Concentration-dependent transcriptome of zebrafish embryo for environmental chemical assessment

Mechanistic information is essential to screen and predict the adverse effects of a large number of chemicals during early-life exposure. Concentration-dependent omics can capture the extent of perturbations of biological pathways or processes and provide information on the mechanism of toxicity. However, the application of concentration-dependent transcriptome to assess the developmental toxicity of environmental chemicals is still limited. Here, twelve chemicals representing five different modes of action (MOAs) were tested by the concentration-dependent reduced zebrafish transcriptome approach (CRZT) in combination with a phenotype-based high content screen (PHCS). The responsiveness, sensitivity and mechanistic differentiation of CRZT were validated in comparison with PHCS. First, PHCS identified 10 chemicals with obvious embryotoxicity (LD₅₀ range: 2.11-70.68 μM), while the potencies of the biological pathways perturbed by 12 chemicals (POD_path20 range: 0.002-2.1 μM) were demonstrated by CRZT. Second, although the potency of the transcriptome perturbations was positively correlated with lethality (LD₅₀) (R² = 0.64, P-value < 0.05) for most tested chemicals, BbF was non-embryotoxic but was the most potent on the perturbance of biological pathways. Finally, the profiles of the perturbed biological processes and the transcriptome potency (POD_path20) captured by CRZT could effectively classify most chemicals corresponding to their known MOAs. In summary, CRZT could significantly improve testing the developmental toxicity of environmental chemicals.

Translational Toxicology in Zebrafish

A major goal of translational toxicology is to identify adverse chemical effects and determine whether they are conserved or divergent across experimental systems. Translational toxicology encompasses assessment of chemical toxicity across multiple life stages, determination of toxic mode-of-action, computational prediction modeling, and identification of interventions that protect or restore health following toxic chemical exposures. The zebrafish is increasingly used in translational toxicology because it combines the genetic and physiological advantages of mammalian models with the higher-throughput capabilities and genetic manipulability of invertebrate models. Here, we review recent literature demonstrating the power of the zebrafish as a model for addressing all four activities of translational toxicology. Important data gaps and challenges associated with using zebrafish for translational toxicology are also discussed.

Use of Zebrafish in Drug Discovery Toxicology

Unpredicted human safety events in clinical trials for new drugs are costly in terms of human health and money. The drug discovery industry attempts to minimize those events with diligent preclinical safety testing. Current standard practices are good at preventing toxic compounds from being tested in the clinic; however, false negative preclinical toxicity results are still a reality. Continual improvement must be pursued in the preclinical realm. Higher-quality therapies can be brought forward with more information about potential toxicities and associated mechanisms. The zebrafish model is a bridge between in vitro assays and mammalian in vivo studies. This model is powerful in its breadth of application and tractability for research. In the past two decades, our understanding of disease biology and drug toxicity has grown significantly owing to thousands of studies on this tiny vertebrate. This Review summarizes challenges and strengths of the model, discusses the 3Rs value that it can deliver, highlights translatable and untranslatable biology, and brings together reports from recent studies with zebrafish focusing on new drug discovery toxicology.

Mechanisms of Developmental Toxicity of Dioxins and Related Compounds

Dioxins and related compounds induce morphological abnormalities in developing animals in an aryl hydrocarbon receptor (AhR)-dependent manner. Here we review the studies in which 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is used as a prototypical compound to elucidate the pathogenesis of morphological abnormalities. TCDD-induced cleft palate in fetal mice involves a delay in palatogenesis and dissociation of fused palate shelves. TCDD-induced hydronephrosis, once considered to be caused by the anatomical obstruction of the ureter, is now separated into TCDD-induced obstructive and non-obstructive hydronephrosis, which develops during fetal and neonatal periods, respectively. In the latter, a prostaglandin E₂ synthesis pathway and urine concentration system are involved. TCDD-induced abnormal development of prostate involves agenesis of the ventral lobe. A suggested mechanism is that AhR activation in the urogenital sinus mesenchyme by TCDD modulates the wingless-type MMTV integration site family (WNT)/β-catenin signaling cascade to interfere with budding from urogenital sinus epithelium. TCDD exposure to zebrafish embryos induces loss of epicardium progenitor cells and heart malformation. AHR2-dependent downregulation of Sox9b expression in cardiomyocytes is a suggested underlying mechanism. TCDD-induced craniofacial malformation in zebrafish is considered to result from the AHR2-dependent reduction in SRY-box 9b (SOX9b), probably partly via the noncoding RNA slincR, resulting in the underdevelopment of chondrocytes and cartilage.