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.

Coupling Genome-wide Transcriptomics and Developmental Toxicity Profiles in Zebrafish to Characterize Polycyclic Aromatic Hydrocarbon (PAH) Hazard.

Polycyclic Aromatic Hydrocarbons (PAHs) are diverse environmental pollutants associated with adverse human health effects. Many studies focus on the carcinogenic effects of a limited number of PAHs and there is an increasing need to understand mechanisms of developmental toxicity of more varied yet environmentally relevant PAHs. A previous study characterized the developmental toxicity of 123 PAHs in zebrafish. Based on phenotypic responses ranging from complete inactivity to acute mortality, we classified these PAHs into eight bins, selected 16 representative PAHs, and exposed developing zebrafish to the concentration of each PAH that induced 80% phenotypic effect. We conducted RNA sequencing at 48 h post fertilization to identify gene expression changes as a result of PAH exposure. Using the Context Likelihood of Relatedness algorithm, we inferred a network that links the PAHs based on coordinated gene responses to PAH exposure. The 16 PAHs formed two broad clusters: Cluster A was transcriptionally more similar to the controls, while Cluster B consisted of PAHs that were generally more developmentally toxic, significantly elevated cyp1a transcript levels, and induced Ahr2-dependent Cyp1a protein expression in the skin confirmed by gene-silencing studies. We found that cyp1a transcript levels were associated with transcriptomic response, but not with PAH developmental toxicity. While all cluster B PAHs predominantly activated Ahr2, they also each enriched unique pathways like ion transport signaling, which likely points to differing molecular events between the PAHs downstream of Ahr2. Thus, using a systems biology approach, we have begun to evaluate, classify, and define mechanisms of PAH toxicity.

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.

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 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). AHR activation is well-understood, but roles of downstream molecular signaling events are largely unknown. From previous transcriptomics in 48 h postfertilization (hpf) zebrafish exposed to several PAHs and TCDD, we found wfikkn1 was highly coexpressed 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 dimethyl sulfoxide or TCDD. 48-hpf mRNA sequencing revealed over 700 genes that were differentially expressed (p < .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.

Developmental toxicity in zebrafish (Danio rerio) exposed to uranium: A comparison with lead, cadmium, and iron.

Populations of plants and animals, including humans, living in close proximity to abandoned uranium mine sites are vulnerable to uranium exposure through drainage into nearby waterways, soil accumulation, and blowing dust from surface soils. Little is known about how the environmental impact of uranium exposure alters the health of human populations in proximity to mine sites, so we used developmental zebrafish (Danio rerio) to investigate uranium toxicity. Fish are a sensitive target for modeling uranium toxicity, and previous studies report altered reproductive capacity, enhanced DNA damage, and gene expression changes in fish exposed to uranium. In our study, dechorionated zebrafish embryos were exposed to a concentration range of uranyl acetate (UA) from 0 to 3000 µg/L for body burden measurements and developmental toxicity assessments. Uranium was taken up in a concentration-dependent manner by 48 and 120 h post fertilization (hpf)-zebrafish without evidence of bioaccumulation. Exposure to UA was not associated with teratogenic outcomes or 24 hpf behavioral effects, but larvae at 120 hpf exhibited a significant hypoactive photomotor response associated with exposure to 3 µg/L UA which suggested potential neurotoxicity. To our knowledge, this is the first time that uranium has been associated with behavioral effects in an aquatic organism. These results were compared to potential metal co-contaminants using the same exposure paradigm. Similar to uranium exposure, lead, cadmium, and iron significantly altered neurobehavioral outcomes in 120-hpf zebrafish without inducing significant teratogenicity. Our study informs concerns about the potential impacts of developmental exposure to uranium on childhood neurobehavioral outcomes. This work also sets the stage for future, environmentally relevant metal mixture studies. Summary Uranium exposure to developing zebrafish causes hypoactive larval swimming behavior similar to the effect of other commonly occurring metals in uranium mine sites. This is the first time that uranium exposure has been associated with altered neurobehavioral effects in any aquatic organism.

The chemistry and toxicology of vaping.

Vaping is the process of inhaling and exhaling an aerosol produced by an e-cigarette, vape pen, or personal aerosolizer. When the device contains nicotine, the Food and Drug Administration (FDA) lists the product as an electronic nicotine delivery system or ENDS device. Similar electronic devices can be used to vape cannabis extracts. Over the past decade, the vaping market has increased exponentially, raising health concerns over the number of people exposed and a nationwide outbreak of cases of severe, sometimes fatal, lung dysfunction that arose suddenly in otherwise healthy individuals. In this review, we discuss the various vaping technologies, which are remarkably diverse, and summarize the use prevalence in the U.S. over time by youths and adults. We examine the complex chemistry of vape carrier solvents, flavoring chemicals, and transformation products. We review the health effects from epidemiological and laboratory studies and, finally, discuss the proposed mechanisms underlying some of these health effects. We conclude that since much of the research in this area is recent and vaping technologies are dynamic, our understanding of the health effects is insufficient. With the rapid growth of ENDS use, consumers and regulatory bodies need a better understanding of constituent-dependent toxicity to guide product use and regulatory decisions.

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.

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.

BACKGROUND: A structurally diverse group of chemicals, including dioxins [e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)] and polycyclic aromatic hydrocarbons (PAHs), can xenobiotically activate the aryl hydrocarbon receptor (AHR) and contribute to adverse health effects in humans and wildlife. In the zebrafish model, repression of sox9b has a causal role in several AHR-mediated toxic responses, including craniofacial cartilage malformations; however, the mechanism of sox9b repression remains unknown. We previously identified a long noncoding RNA, sox9b long intergenic noncoding RNA (slincR), which is increased (in an AHR-dependent manner) by multiple AHR ligands and is required for the AHR-activated repression of sox9b. OBJECTIVE: Using the zebrafish model, we aimed to enhance our understanding of the signaling events downstream of AHR activation that contribute to toxic responses by identifying: a) whether slincR is enriched on the sox9b locus, b) slincR's functional contributions to TCDD-induced toxicity, c) PAHs that increase slincR expression, and d) mammalian orthologs of slincR. METHODS: We used capture hybridization analysis of RNA targets (CHART), qRT-PCR, RNA sequencing, morphometric analysis of cartilage structures, and hemorrhaging screens. RESULTS: The slincR transcript was enriched at the 5' untranslated region (UTR) of the sox9b locus. Transcriptome profiling and human ortholog analyses identified processes related to skeletal and cartilage development unique to TCDD-exposed controls, and angiogenesis and vasculature development unique to TCDD-exposed zebrafish that were injected with a splice-blocking morpholino targeting slincR. In comparison to TCDD exposed control morphants, slincR morphants exposed to TCDD resulted in abnormal cartilage structures and a smaller percentage of animals displaying the hemorrhaging phenotype. In addition, slincR expression was significantly increased in six out of the sixteen PAHs we screened. CONCLUSION: Our study establishes that in zebrafish, slincR is recruited to the sox9b 5' UTR to repress transcription, can regulate cartilage development, has a causal role in the TCDD-induced hemorrhaging phenotype, and is up-regulated by multiple environmentally relevant PAHs. These findings have important implications for understanding the ligand-specific mechanisms of AHR-mediated toxicity. https://doi.org/10.1289/EHP3281.