The addition of mammalian cell culture medium impacts nanoparticle toxicity in zebrafish

Engineered nanomaterials (ENMs) are ubiquitous in contemporary applications, yet their environmental and human health impacts remain inadequately understood. This study addresses the challenge of identifying potential risks associated with ENM exposure by highlighting the significant variability in existing research methodologies. Without a systematic collection of toxicological data that encompasses standardized materials, relevant platforms, and assays, the task of identifying potential risks linked to ENM exposure becomes an intricate challenge. In vitro assessments often use media rich in ionic species, such as RPMI and fetal bovine serum (FBS). Zebrafish embryos, known to develop normally in low-ionic environments, were exposed to Cerium Oxide, Zinc Oxide, and Graphene Oxides in different media at varying concentrations. Here, we discovered that zebrafish embryos tolerated a mix of 80 % RPMI, 2 % FBS, and 1 % antibiotic cocktail. The results revealed that adverse effects observed in zebrafish with certain nanomaterials in Ultra-Pure (UP) water were mitigated in cell culture medium, emphasizing the importance of revisiting previously considered non-toxic materials in vitro. The zebrafish results underscore the importance of utilizing a multidimensional in vivo platform to gauge the biological activity of nanomaterials accurately.

PAH bioremediation with Rhodococcus rhodochrous ATCC 21198: Impact of cell immobilization and surfactant use on PAH treatment and post-remediation toxicity

Polycyclic aromatic hydrocarbons (PAHs) are prevalent environmental contaminants that are harmful to ecological and human health. Bioremediation is a promising technique for remediating PAHs in the environment, however bioremediation often results in the accumulation of toxic PAH metabolites. The objectives of this research were to demonstrate the cometabolic treatment of a mixture of PAHs by a pure bacterial culture, Rhodococcus rhodochrous ATCC 21198, and investigate PAH metabolites and toxicity. Additionally, the surfactant Tween ® 80 and cell immobilization techniques were used to enhance bioremediation. Total PAH removal ranged from 70-95% for fluorene, 44-89% for phenanthrene, 86-97% for anthracene, and 6.5-78% for pyrene. Maximum removal was achieved with immobilized cells in the presence of Tween ® 80. Investigation of PAH metabolites produced by 21198 revealed a complex mixture of hydroxylated compounds, quinones, and ring-fission products. Toxicity appeared to increase after bioremediation, manifesting as mortality and developmental effects in embryonic zebrafish. 21198's ability to rapidly transform PAHs of a variety of molecular structures and sizes suggests that 21198 can be a valuable microorganism for catalyzing PAH remediation. However, implementing further treatment processes to address toxic PAH metabolites should be pursued to help lower post-remediation toxicity in future studies.

Interlaboratory Study on Zebrafish in Toxicology: Systematic Evaluation of the Application of Zebrafish in Toxicology's (SEAZIT's) Evaluation of Developmental Toxicity

Embryonic zebrafish represent a useful test system to screen substances for their ability to perturb development. The exposure scenarios, endpoints captured, and data analysis vary among the laboratories who conduct screening. A lack of harmonization impedes the comparison of the substance potency and toxicity outcomes across laboratories and may hinder the broader adoption of this model for regulatory use. The Systematic Evaluation of the Application of Zebrafish in Toxicology (SEAZIT) initiative was developed to investigate the sources of variability in toxicity testing. This initiative involved an interlaboratory study to determine whether experimental parameters altered the developmental toxicity of a set of 42 substances (3 tested in duplicate) in three diverse laboratories. An initial dose-range-finding study using in-house protocols was followed by a definitive study using four experimental conditions: chorion-on and chorion-off using both static and static renewal exposures. We observed reasonable agreement across the three laboratories as 33 of 42 test substances (78.6%) had the same activity call. However, the differences in potency seen using variable in-house protocols emphasizes the importance of harmonization of the exposure variables under evaluation in the second phase of this study. The outcome of the Def will facilitate future practical discussions on harmonization within the zebrafish research community.

Developmental Toxicity Assessment Using Zebrafish-Based High-Throughput Screening

Zebrafish-based high-throughput screening has been extensively used to study toxicological profiles of individual chemicals and mixtures, identify novel toxicants, and study modes of action to prioritize chemicals for further testing and policy decisions. Within this chapter, we describe a protocol for automated zebrafish developmental high-throughput screening in our laboratory, with emphasis on exposure setups, morphological and behavioral readouts, and quality control.

Defining the environmental determinants of dysbiosis at scale with zebrafish

The gut microbiome, critical to maintaining vertebrate homeostasis, is susceptible to a various exposures. In some cases, these exposures induce dysbiosis, wherein the microbiome changes into a state conducive to disease progression. To better prevent, manage, and treat health disorders, we need to define which exposures induce dysbiosis. Contemporary methods face challenges due to the immense diversity of the exposome and the restricted throughput of conventional experimental tools used for dysbiosis evaluation. We propose integrating high-throughput model systems as an augment to traditional techniques for rapid identification of dysbiosis-inducing agents. Although high-throughput screening tools revolutionized areas such as pharmacology and toxicology, their incorporation in gut microbiome research remains limited. One particularly powerful high-throughput model system is the zebrafish, which affords access to scalable in vivo experimentation involving a complex gut microbiome. Numerous studies have employed this model to identify potential dysbiosis triggers. However, its potential could be further harnessed via innovative study designs, such as evaluation of synergistic effects from combined exposures, expansions to the methodological toolkit to discern causal effects of microbiota, and efforts to assess and improve the translational relevance of the model. Ultimately, this burgeoning experimental resource can accelerate the discovery of agents that underlie dysbiotic disorders.

Comparative analysis between zebrafish and an automated live-cell assay to classify developmental neurotoxicant chemicals

Modern toxicology's throughput has dramatically increased due to alternative models, laboratory automation, and machine learning. This has enabled comparative studies across species and assays to prioritize chemical hazard potential and to understand how different model systems might complement one another. However, such comparative studies of high-throughput data are still in their infancy, with more groundwork needed to firmly establish the approach. Therefore, this study aimed to compare the bioactivity of the NIEHS Division of Translational Toxicology's (DTT) 87-compound developmental neurotoxicant (DNT) library in zebrafish and an in vitro high-throughput cell culture system. The early life-stage zebrafish provided a whole animal approach to developmental toxicity assessment. Chemical hits for abnormalities in embryonic zebrafish morphology, mortality, and behavior (ZBEscreen™) were compared with chemicals classified as high-risk by the Cell Health Index (CHI™), which is an outcome class probability from a machine learning classifier using 12 parameters from the SYSTEMETRIC® Cell Health Screen (CHS). The CHS was developed to assess human toxicity risk using supervised machine learning to classify acute cell stress phenotypes in a human leukemia cell line (HL60 cells) following a 4-h exposure to a chemical of interest. Due to the design of the screen, the zebrafish assays were more exhaustive, yielding 86 total bioactive hits, whereas the SYSTEMETRIC® CHS focusing on acute toxicity identified 20 chemicals as potentially toxic. The zebrafish embryonic and larval photomotor response assays (EPR and LPR, respectively) detected 40 of the 47 chemicals not found by the zebrafish morphological screen and CHS. Collectively, these results illustrate the advantages of using two alternative models in tandem for rapid hazard assessment and chemical prioritization and the effectiveness of CHI™ in identifying toxicity within a single multiparametric assay.

Deep autoencoder-based behavioral pattern recognition outperforms standard statistical methods in high-dimensional zebrafish studies

Zebrafish have become an essential tool in screening for developmental neurotoxic chemicals and their molecular targets. The success of zebrafish as a screening model is partially due to their physical characteristics including their relatively simple nervous system, rapid development, experimental tractability, and genetic diversity combined with technical advantages that allow for the generation of large amounts of high-dimensional behavioral data. These data are complex and require advanced machine learning and statistical techniques to comprehensively analyze and capture spatiotemporal responses. To accomplish this goal, we have trained semi-supervised deep autoencoders using behavior data from unexposed larval zebrafish to extract quintessential "normal" behavior. Following training, our network was evaluated using data from larvae shown to have significant changes in behavior (using a traditional statistical framework) following exposure to toxicants that include nanomaterials, aromatics, per- and polyfluoroalkyl substances (PFAS), and other environmental contaminants. Further, our model identified new chemicals (Perfluoro-n-octadecanoic acid, 8-Chloroperfluorooctylphosphonic acid, and Nonafluoropentanamide) as capable of inducing abnormal behavior at multiple chemical-concentrations pairs not captured using distance moved alone. Leveraging this deep learning model will allow for better characterization of the different exposure-induced behavioral phenotypes, facilitate improved genetic and neurobehavioral analysis in mechanistic determination studies and provide a robust framework for analyzing complex behaviors found in higher-order model systems.

3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

Chronic Vitamin E Deficiency Dysregulates Purine, Phospholipid, and Amino Acid Metabolism in Aging Zebrafish Skeletal Muscle

Muscle wasting occurs with aging and may be a result of oxidative stress damage and potentially inadequate protection by lipophilic antioxidants, such as vitamin E. Previous studies have shown muscular abnormalities and behavioral defects in vitamin E-deficient adult zebrafish. To test the hypothesis that there is an interaction between muscle degeneration caused by aging and oxidative damage caused by vitamin E deficiency, we evaluated long-term vitamin E deficiency in the skeletal muscle of aging zebrafish using metabolomics. Zebrafish (55 days old) were fed E+ and E- diets for 12 or 18 months. Then, skeletal muscle samples were analyzed using UPLC-MS/MS. Data were analyzed to highlight metabolite and pathway changes seen with either aging or vitamin E status or both. We found that aging altered purines, various amino acids, and DHA-containing phospholipids. Vitamin E deficiency at 18 months was associated with changes in amino acid metabolism, specifically tryptophan pathways, systemic changes in the regulation of purine metabolism, and DHA-containing phospholipids. In sum, while both aging and induced vitamin E deficiency did have some overlap in altered and potentially dysregulated metabolic pathways, each factor also presented unique alterations, which require further study with more confirmatory approaches.

Review of the Zebrafish as a Model to Investigate per- and Polyfluoroalkyl Substance Toxicity

The existence of thousands of per- and polyfluoroalkyl substances (PFAS) and evidence that some cause adverse health effects has created immense need to better understand PFAS toxicity and to move beyond one-chemical-at-a-time approaches to hazard assessment for this chemical class. The zebrafish model enables rapid assessment of large libraries of PFAS, powerful comparison of compounds in a single in vivo system, and evaluation across life stages and generations, and has led to significant advances in PFAS research in recent years. The focus of this review is to assess contemporary findings regarding PFAS toxicokinetics, toxicity and apical adverse health outcomes, and potential modes of action using the zebrafish model. Much of the peer-reviewed literature has focused on a small subset of PFAS structural subclasses, such as the perfluoroalkyl sulfonic acids and perfluoroalkyl carboxylic acids. However, recent data on more diverse PFAS structures are enabling prioritization of compounds of concern. Structure-activity comparisons and the utilization of modeling and 'omics technologies in zebrafish have greatly contributed to our understanding of the hazard potential for a growing number of PFAS and will surely inform our understanding and predictive capabilities for many more PFAS in the future.

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

The role of long non-coding RNAs (lncRNAs) regulators of toxicological responses to environmental chemicals is gaining prominence. Previously, our laboratory discovered a 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.

11-Cl-BBQ, a select modulator of AhR-regulated transcription, suppresses lung cancer cell growth via activation of p53 and p27Kip1

Aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor and functions as a tumour suppressor in different cancer models. In the present study, we report detailed characterization of 11-chloro-7H-benzimidazo[2,1-a]benzo[de]iso-quinolin-7-one (11-Cl-BBQ) as a select modulator of AhR-regulated transcription (SMAhRT) with anti-cancer actions. Treatment of lung cancer cells with 11-Cl-BBQ induced potent and sustained AhR-dependent anti-proliferative effects by promoting G1 phase cell cycle arrest. Investigation of 11-Cl-BBQ-induced transcription in H460 cells with or without the AhR expression by RNA-sequencing revealed activation of p53 signalling. In addition, 11-Cl-BBQ suppressed multiple pathways involved in DNA replication and increased expression of cyclin-dependent kinase inhibitors, including p27Kip1 , in an AhR-dependent manner. CRISPR/Cas9 knockout of individual genes revealed the requirement for both p53 and p27Kip1 for the AhR-mediated anti-proliferative effects. Our results identify 11-Cl-BBQ as a potential lung cancer therapeutic, highlight the feasibility of targeting AhR and provide important mechanistic insights into AhR-mediated-anticancer actions.

The Superfund Research Program Analytics Portal: linking environmental chemical exposure to biological phenotypes

The OSU/PNNL Superfund Research Program (SRP) represents a longstanding collaboration to quantify Polycyclic Aromatic Hydrocarbons (PAHs) at various superfund sites in the Pacific Northwest and assess their potential impact on human health. To link the chemical measurements to biological activity, we describe the use of the zebrafish as a high-throughput developmental toxicity model that provides quantitative measurements of the exposure to chemicals. Toward this end, we have linked over 150 PAHs found at Superfund sites to the effect of these same chemicals in zebrafish, creating a rich dataset that links environmental exposure to biological response. To quantify this response, we have implemented a dose-response modelling pipeline to calculate benchmark dose parameters which enable potency comparison across over 500 chemicals and 12 of the phenotypes measured in zebrafish. We provide a rich dataset for download and analysis as well as a web portal that provides public access to this dataset via an interactive web site designed to support exploration and re-use of these data by the scientific community at http://srp.pnnl.gov .

Coupling Environmental Whole Mixture Toxicity Screening with Unbiased RNA-Seq Reveals Site-Specific Biological Responses in Zebrafish

Passive sampling device (PSD) extracts paired with developmental toxicity assays in Danio Rerio (zebrafish) are excellent sensors for whole mixture toxicity associated with the bioavailable non-polar organics at environmental sites. We expand this concept by incorporating RNA-Seq in 48-h post fertilization zebrafish statically exposed to PSD extracts from two Portland Harbor Superfund Site locations: river mile 6.5W (RM 6.5W) and river mile 7W (RM 7W). RM 6.5W contained higher concentrations of polycyclic aromatic hydrocarbons (PAHs), but the diagnostic ratios of both extracts indicated similar PAH sourcing and composition. Developmental screens determined RM 6.5W to be more toxic with the most sensitive endpoint being a "wavy" notochord malformation. Differential gene expression from exposure to both extracts was largely parallel, although more pronounced for RM 6.5W. When compared to the gene expression associated with individual chemical exposures, PSD extracts produced some gene signatures parallel to PAHs but were more closely matched by oxygenated-PAHs. Additionally, differential expression, reminiscent of the wavy notochord phenotype, was not accounted for by either class of chemical, indicating the potential of other contaminants driving mixture toxicity. These techniques offer a compelling method for non-targeted hazard characterization of whole mixtures in an in vivo vertebrate system without requiring complete chemical characterization.

Utilizing a Population-Genetic Framework to Test for Gene-Environment Interactions between Zebrafish Behavior and Chemical Exposure

Individuals within genetically diverse populations display broad susceptibility differences upon chemical exposures. Understanding the role of gene-environment interactions (GxE) in differential susceptibility to an expanding exposome is key to protecting public health. However, a chemical's potential to elicit GxE is often not considered during risk assessment. Previously, we've leveraged high-throughput zebrafish (Danio rerio) morphology screening data to reveal patterns of potential GxE effects. Here, using a population genetics framework, we apportioned variation in larval behavior and gene expression in three different PFHxA environments via mixed-effect modeling to assess significance of GxE term. We estimated the intraclass correlation (ICC) between full siblings from different families using one-way random-effects model. We found a significant GxE effect upon PFHxA exposure in larval behavior, and the ICC of behavioral responses in the PFHxA exposed population at the lower concentration was 43.7%, while that of the control population was 14.6%. Considering global gene expression data, a total of 3746 genes showed statistically significant GxE. By showing evidence that heritable genetics are directly affecting gene expression and behavioral susceptibility of individuals to PFHxA exposure, we demonstrate how standing genetic variation in a heterogeneous population such as ours can be leveraged to test for potential GxE.

Leveraging Multiple Data Streams for Prioritization of Mixtures for Hazard Characterization

There is a growing need to establish alternative approaches for mixture safety assessment of polycyclic aromatic hydrocarbons (PAHs). Due to limitations with current component-based approaches, and the lack of established methods for using whole mixtures, a promising alternative is to use sufficiently similar mixtures; although, an established framework is lacking. In this study, several approaches are explored to form sufficiently similar mixtures. Multiple data streams including environmental concentrations and empirically and predicted toxicity data for cancer and non-cancer endpoints were used to prioritize chemical components for mixture formations. Air samplers were analyzed for unsubstituted and alkylated PAHs. A synthetic mixture of identified PAHs was created (Creosote-Fire Mix). Existing toxicity values and chemical concentrations were incorporated to identify hazardous components in the Creosote-Fire Mix. Sufficiently similar mixtures of the Creosote-Fire Mix were formed based on (1) relative abundance; (2) toxicity values; and (3) a combination approach incorporating toxicity and abundance. Hazard characterization of these mixtures was performed using high-throughput screening in primary normal human bronchial epithelium (NHBE) and zebrafish. Differences in chemical composition and potency were observed between mixture formation approaches. The toxicity-based approach (Tox Mix) was the most potent mixture in both models. The combination approach (Weighted-Tox Mix) was determined to be the ideal approach due its ability to prioritize chemicals with high exposure and hazard potential.

Background per- and polyfluoroalkyl substances (PFAS) in laboratory fish diet: Implications for zebrafish toxicological studies

Current attention is focused on determining the potential for per- and polyfluoroalkyl substances (PFAS) to adversely impact human health. Zebrafish are a popular biological model because they share early development pathways with humans. A dietary exposure paradigm is growing in popularity in the zebrafish model because the outcomes often translate to humans. To create a diet of known composition, it is crucial to understand background PFAS levels present in zebrafish diet. Background PFAS, if present, potentially confounds interpretation of toxicological data. To date, no studies document the PFAS background levels in laboratory fish diet and there is only limited information on some pet foods. The objective of this study was to develop and validate an analytical method for up to 50 target PFAS in high lipid and protein content laboratory fish diets and pet foods. Long-chain perfluoroalkyl carboxylic acids (C9-C13) and perfluorooctane sulfonate (PFOS) were quantified in 11 out of 16 laboratory fish diets and in three out of five pet fish foods. Foods for pet birds, lizards, and dogs were below the limit of detection for all PFAS. In two of the laboratory fish diets, PFOS concentrations were >1.3 ng/g and the total PFAS for the three laboratory fish diets exceeded 1.0 ng/g. Hundreds of biomedical laboratories across the world utilize these commercial laboratory fish diets, and these results indicate that numerous zebrafish colonies may be inadvertently receiving significant dietary PFAS exposures. In light of this new information, it is critical to design PFAS studies with appropriate controls with measured background PFAS concentrations in the diet and to urge caution when interpreting the results.

Oxidative stress, glutathione, and CYP2E1 in 1,4-dioxane liver cytotoxicity and genotoxicity: insights from animal models

1,4-Dioxane (DX) is an emerging drinking water contaminant worldwide, which poses a threat to public health due to its demonstrated liver carcinogenicity and potential for human exposure. The lack of drinking water standards for DX is attributed to undetermined mechanisms of DX carcinogenicity. This mini-review provides a brief discussion of a series of mechanistic studies, wherein unique mouse models were exposed to DX in drinking water to elucidate redox changes associated with DX cytotoxicity and genotoxicity. The overall conclusions from these studies support a direct genotoxic effect by high dose DX and imply that oxidative stress involving CYP2E1 activation may play a causal role in DX liver genotoxicity and potentially carcinogenicity. The mechanistic data derived from these studies can serve as important references to refine the assessment of carcinogenic pathways that may be triggered at environmentally relevant low doses of DX in future animal and human studies.

Experimental methods modestly impact interpretation of the effect of environmental exposures on the larval zebrafish gut microbiome

Rapidly growing fields, such as microbiome science, often lack standardization of procedures across research groups. This is especially the case for microbiome investigations in the zebrafish (Danio rerio) model system, which is quickly becoming a workhorse system for understanding the exposure-microbiome-physiology axis. To guide future investigations using this model system, we defined how various experimental decisions affect the outcomes of studies on the effects of exogenous exposure on the zebrafish gut microbiome. Using a model toxicant, benzo[a]pyrene (BaP), we assessed how each of two dissection methods (gut dissection vs. whole fish), three DNA extraction kits (Qiagen Blood & Tissue, Macherey-Nagel NucleoSpin, and Qiagen PowerSoil), and inclusion of PCR replicates (single vs. pooled triplicate reactions) affected our interpretation of how exposure influences the diversity and composition of the gut microbiome, as well as our ability to identify microbiome biomarkers of exposure. We found that inclusion of PCR replicates had the smallest effect on our final interpretations, and the effects of dissection method and DNA extraction kit had significant effects in specific contexts, primarily in the cases of identifying microbial biomarkers.

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.

Nitrate exposure reprograms hepatic amino acid and nutrient sensing pathways prior to exercise: A metabolomic and transcriptomic investigation in zebrafish (Danio rerio)

Scope: Nitrate supplementation is a popular ergogenic aid that improves exercise performance by reducing oxygen consumption during exercise. We investigated the effect of nitrate exposure and exercise on metabolic pathways in zebrafish liver.

Materials and methods: Fish were exposed to sodium nitrate (606.9 mg/L), or control water, for 21 days and analyzed at intervals during an exercise test. We utilized untargeted liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and measured gene expression of 24 genes central to energy metabolism and redox signaling.

Results: We observed a greater abundance of metabolites involved in endogenous nitric oxide (NO) metabolism and amino acid metabolism in nitrate-treated liver at rest, compared to rested controls. In the absence of exercise, nitrate treatment upregulated expression of genes central to nutrient sensing (pgc1a), protein synthesis (mtor) and purine metabolism (pnp5a and ampd1) and downregulated expression of genes involved in mitochondrial fat oxidation (acaca and cpt2).

Conclusion: Our data support a role for sub-chronic nitrate treatment in the improvement of exercise performance, in part, by improving NO bioavailability, sparing arginine, and modulating hepatic gluconeogenesis and glycolytic capacity in the liver.

Dietary Perfluorohexanoic Acid (PFHxA) Exposures in Juvenile Zebrafish Produce Subtle Behavioral Effects across Generations

Ubiquitous anthropogenic contaminants of concern, per- and polyfluoroalkyl substances (PFAS) are frequently detected in the environment and human populations around the world. Diet is a predominate route of human exposure, and PFAS are frequently measured in food. Manufacturing trends have shifted from legacy PFAS to shorter-chain alternatives that are suggested to be safer, such as perfluorohexanoic acid (PFHxA). However, the current amount of data to support safety assessments of these alternatives is not yet sufficient. The present study investigated the effects of a 42-day dietary exposure to 1, 10, or 100 ng/g PFHxA in juvenile zebrafish. The zebrafish model was leveraged to interrogate morphometrics, fecundity, and numerous behavior endpoints across multiple generations. Dietary PFHxA exposure did not result in measurable body burden and did not affect growth, fecundity, adult social perception behavior, or associative learning. PFHxA exposure did induce abnormal adult anxiety behaviors in the F0 generation that persisted transgenerationally in the F1 and F2. Abnormal larval and juvenile behavior was observed in the F1 generation, but not in the F2. PFHxA juvenile dietary exposure induced subtle and multigenerational behavior effects that warrant further investigation of this and other alternative short-chain PFAS.

Systematic developmental toxicity assessment of a structurally diverse library of PFAS in zebrafish

Per- and polyfluoroalkyl substances (PFAS) are a class of widely used chemicals with limited human health effects data relative to the diversity of structures manufactured. To help fill this data gap, an extensive in vivo developmental toxicity screen was performed on 139 PFAS provided by the US EPA. Dechorionated embryonic zebrafish were exposed to 10 nominal water concentrations of PFAS (0.015-100 µM) from 6 to 120 h post-fertilization (hpf). The embryos were assayed for embryonic photomotor response (EPR), larval photomotor response (LPR), and 13 morphological endpoints. A total of 49 PFAS (35%) were bioactive in one or more assays (11 altered EPR, 25 altered LPR, and 31 altered morphology). Perfluorooctanesulfonamide (FOSA) was the only structure that was bioactive in all 3 assays, while Perfluorodecanoic acid (PFDA) was the most potent teratogen. Low PFAS volatility was associated with developmental toxicity (p < 0.01), but no association was detected between bioactivity and five other physicochemical parameters. The bioactive PFAS were enriched for 6 supergroup chemotypes. The results illustrate the power of a multi-dimensional in vivo platform to assess the developmental (neuro)toxicity of diverse PFAS and in the acceleration of PFAS safety research.

Leveraging a High-Throughput Screening Method to Identify Mechanisms of Individual Susceptibility Differences in a Genetically Diverse Zebrafish Model

Understanding the mechanisms behind chemical susceptibility differences is key to protecting sensitive populations. However, elucidating gene-environment interactions (GxE) presents a daunting challenge. While mammalian models have proven useful, problems with scalability to an enormous chemical exposome and clinical translation faced by all models remain; therefore, alternatives are needed. Zebrafish (Danio rerio) have emerged as an excellent model for investigating GxE. This study used a combined bioinformatic and experimental approach to probe the mechanisms underlying chemical susceptibility differences in a genetically diverse zebrafish population. Starting from high-throughput screening (HTS) data, a genome-wide association study (GWAS) using embryonic fish exposed to 0.6 μM Abamectin revealed significantly different effects between individuals. A hypervariable region with two distinct alleles–one with G at the SNP locus (GG) and one with a T and the 16 bp deletion (TT)–associated with differential susceptibility was found. Sensitive fish had significantly lower sox7 expression. Due to their location and the observed expression differences, we hypothesized that these sequences differentially regulate sox7. A luciferase reporter gene assay was used to test if these sequences, alone, could lead to expression differences. The TT allele showed significantly lower expression than the GG allele in MCF-7 cells. To better understand the mechanism behind these expression differences, predicted transcription factor binding differences between individuals were compared in silico, and several putative binding differences were identified. EMSA was used to test for binding differences in whole embryo protein lysate to investigate these TF binding predictions. We confirmed that the GG sequence is bound to protein in zebrafish. Through a competition EMSA using an untagged oligo titration, we confirmed that the GG oligo had a higher binding affinity than the TT oligo, explaining the observed expression differences. This study identified differential susceptibility to chemical exposure in a genetically diverse population, then identified a plausible mechanism behind those differences from a genetic to molecular level. Thus, an HTS-compatible zebrafish model is valuable and adaptable in identifying GxE mechanisms behind susceptibility differences to chemical exposure.

Developmental, Behavioral and Transcriptomic Changes in Zebrafish Embryos after Smoke Dye Exposure

(1) Background: Disperse Blue 14, Disperse Red 9, Solvent Red 169 and Solvent Yellow 33 have been used to color smoke; however, they have not been comprehensively assessed for their potential health hazards. (2) Methods: To assess the effects of these dyes, zebrafish embryos were exposed from 6 to 120 h post fertilization (hpf) to 10-55 µM Disperse Red 9, 1-50 µM Solvent Red 169, 7.5-13.5 µM Solvent Yellow 33 or 133-314 µM Disperse Blue 14. Embryos were monitored for adverse effects on gene expression at 48 hpf as well as for mortality, development and behavior at 120 hpf. The dyes were examined for their potential to cross the blood-brain barrier. (3) Results: Solvent Yellow 33 and Disperse Blue 14 impaired development and behavior at all concentrations. Disperse Red 9 impaired behavior at all concentrations and development at all concentrations except for 10 µM. Solvent Red 169 caused no effects. Mortality was only seen in Disperse Blue 14 at 261.5 and 314 µM. Gene expression indicated impacts on neurodevelopment and folate and retinol metabolism as potential mechanisms of toxicity. (4) Conclusions: Smoke dyes have a high potential for causing developmental changes and neurotoxicity and should be examined more closely using comprehensive approaches as used here.

Transcriptomic and Long-Term Behavioral Deficits Associated with Developmental 3.5 GHz Radiofrequency Radiation Exposures in Zebrafish

The rapid deployment of the fifth-generation (5G) spectrum by the telecommunication industry is intended to promote better connectivity and data integration among various industries. However, concerns among the public about the safety and health effects of radiofrequency radiations (RFRs) emitted from the newer-generation cell phone frequencies remain, partly due to the lack of robust scientific data. Previously, we used developmental zebrafish to model the bioactivity of 3.5 GHz RFR, a frequency used by 5G-enabled cell phones, in a novel RFR exposure chamber. With RFR exposures from 6 h post-fertilization (hpf) to 48 hpf, we observed that, despite no teratogenic effects, embryos showed subtle hypoactivity in a startle response behavior assay, suggesting abnormal sensorimotor behavior. This study builds upon the previous one by investigating the transcriptomic basis of RFR-associated behavior effects and their persistence into adulthood. Using mRNA sequencing, we found a modest transcriptomic disruption at 48 hpf, with 28 differentially expressed genes. KEGG pathway analysis showed that biochemical pathways related to metabolism were significantly perturbed. Embryos were grown to adulthood, and then a battery of behavioral assays suggested subtle but significant abnormal responses in RFR-exposed fish across the different assays evaluated that suggest potential long-term behavioral effects. Overall, our study suggests the impacts of RFRs on the developing brain, behavior, and the metabolome should be further explored.

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.

A Comparative Multi-System Approach to Characterizing Bioactivity of Commonly Occurring Chemicals

A 2019 retrospective study analyzed wristband personal samplers from fourteen different communities across three different continents for over 1530 organic chemicals. Investigators identified fourteen chemicals (G14) detected in over 50% of personal samplers. The G14 represent a group of chemicals that individuals are commonly exposed to, and are mainly associated with consumer products including plasticizers, fragrances, flame retardants, and pesticides. The high frequency of exposure to these chemicals raises questions of their potential adverse human health effects. Additionally, the possibility of exposure to mixtures of these chemicals is likely due to their co-occurrence; thus, the potential for mixtures to induce differential bioactivity warrants further investigation. This study describes a novel approach to broadly evaluate the hazards of personal chemical exposures by coupling data from personal sampling devices with high-throughput bioactivity screenings using in vitro and non-mammalian in vivo models. To account for species and sensitivity differences, screening was conducted using primary normal human bronchial epithelial (NHBE) cells and early life-stage zebrafish. Mixtures of the G14 and most potent G14 chemicals were created to assess potential mixture effects. Chemical bioactivity was dependent on the model system, with five and eleven chemicals deemed bioactive in NHBE and zebrafish, respectively, supporting the use of a multi-system approach for bioactivity testing and highlighting sensitivity differences between the models. In both NHBE and zebrafish, mixture effects were observed when screening mixtures of the most potent chemicals. Observations of BMC-based mixtures in NHBE (NHBE BMC Mix) and zebrafish (ZF BMC Mix) suggested antagonistic effects. In this study, consumer product-related chemicals were prioritized for bioactivity screening using personal exposure data. High-throughput high-content screening was utilized to assess the chemical bioactivity and mixture effects of the most potent chemicals.

Size- and Oxidation-Dependent Toxicity of Graphene Oxide Nanomaterials in Embryonic Zebrafish

Graphene oxides (GOs) are a popular graphene alternative. The goal of this study was to compare the biocompatibility of a diversity of well-characterized GOs. Our previous work advanced developmental zebrafish as a model to interrogate the interactions and biological responses following exposures to engineered nanomaterials (ENMs). Here, we investigated GO 250 nm × 250 nm (sGO), 400 nm × 400 nm (mGO), and 1 μm × 1 μm (lGO), partially reduced GO (prGO) 400 nm × 400 nm, and reduced GO (rGO) 400 nm × 400 nm and 2 μm × 2 μm, which first underwent extensive characterization under the support of the Nanomaterials Health Implications Research (NHIR) Consortium. GOs were stabilized in water (GOs), while prGO and rGOs were dispersed in sodium cholate. Zebrafish were statically exposed to up to 50 μg/mL of each material from 6 h post-fertilization (hpf) until 120 hpf. Toxicity was dependent on GO properties. mGO was the most toxic material; its effects manifested in the yolk syncytial layer (YSL). Additionally, sodium cholate stabilization significantly increased GO toxicity. The observed effects were size- and oxidation-state-dependent, revealing the importance of identifying the structure-specific toxicity of GOs.

Sulfonamide functional head on short-chain perfluorinated substance drives developmental toxicity

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously detected in environmental and biological samples and cause adverse health effects. Studies have predominately focused on long-chain PFAS, with far fewer addressing short-chain alternatives. This study leveraged embryonic zebrafish to investigate developmental toxicity of a short-chain series: perfluorobutane sulfonate (PFBS), perfluoropentanoic acid (PFPeA), perfluorobutane sulfonamide (FBSA), and 4:2 fluorotelomer sulfonic acid (4:2 FTS). Following static exposures at 8 h postfertilization (hpf) to each chemical (1-100 μM), morphological and behavioral endpoints were assessed at 24 and 120 hpf. Only FBSA induced abnormal morphology, while exposure to all chemicals caused aberrant larval behavior. RNA sequencing at 48 hpf following 47 μM exposures revealed only FBSA significantly disrupted normal gene expression. Measured tissue concentrations were FBSA > PFBS > 4:2 FTS > PFPeA. This study demonstrates functional head groups impact bioactivity and bioconcentration.

Zebrafish Behavioral Assays in Toxicology

Zebrafish behavioral assays are commonly used to identify and study environmental stressors that elicit adverse effects on neurobehavior. Behavioral assay platforms are available for multiple life stages (embryonic, juvenile, and adults) and are robust in detecting stressor-induced acute effects on neurodevelopment as well as long term deficits in sensory mechanisms, social behavior, learning, memory, and neurodegenerative diseases. Within this chapter, we present an overview of zebrafish behavioral assays that are commonly used to study environmental neurotoxicants.

Developmental titanium dioxide nanoparticle exposure induces oxidative stress and neurobehavioral changes in zebrafish

The widespread commercial application of titanium dioxide nanoparticles (TiO₂ NPs) leads to ubiquitous presence of TiO₂ NPs in the aquatic environment, which highlights the necessity to determine their potential adverse effects on aquatic organisms. The developing nerve system is particularly susceptible to environment perturbation. However, few studies have explored the developmental neurobehavioral toxicity of TiO₂ NPs, especially at smaller particle size ranges (≤20 nm) that have relatively longer retention time in the water column. In this study, zebrafish embryos were exposed to non-teratogenic concentrations of 0.1 and 1 mg/L TiO₂ NPs (average size of 14-20 nm) from 8 to 108 h post-fertilization (hpf) followed by various assessments at different time points up to 12 days post-fertilization (dpf). Our findings revealed that 1 mg/L TiO₂ NPs perturbed the motor and social behaviors in larval zebrafish. These behavioral changes were characterized by decreased swimming speed in a locomotor response test at 5 dpf, increased travel distance in a flash stimulus test at 5 dpf, increased preference to the light zone in a light/dark preference test at 10 dpf, and increased mirror attack and percent time spent in the mirror zone in a mirror stimulus response assay at 12 dpf. Mechanistic examinations at 5 dpf revealed elevated cell apoptosis and oxidative stress. Cell apoptosis was characterized by increased acridine orange (AO) positive cells in the olfactory region and neuromasts of the lateral line system. Oxidative stress was characterized by increased lipid peroxidation, increased ROS production, and upregulated catalase (cat) gene expression. In addition, TiO₂ NP exposure also upregulated genes associated with the developmental nervous system such as the growth associated protein 43 (gap43) and proliferating cell nuclear antigen (pcna). Our results suggest that the neurobehavioral changes in larvae exposed to 1 mg/L TiO₂ NPs during early development may result from cell apoptosis and oxidative stress induced neuronal damages.

Gene Expression of CRAL_TRIO Family Proteins modulated by Vitamin E Deficiency in Zebrafish (Danio Rerio)

An evaluation of the impact of vitamin E deficiency on expression of the alpha-tocopherol transfer protein (α-TTP) and related CRAL_TRIO genes was undertaken using livers from adult zebrafish based on the hypothesis that increased lipid peroxidation would modulate gene expression. Zebrafish were fed either a vitamin E sufficient (E+) or deficient (E-) diet for 9 months, then fish were euthanized, and livers were harvested. Livers from the E+ relative to E- fish contained 40-times more α-tocopherol (P <0.0001) and one fourth the malondialdehyde (P = 0.0153). RNA was extracted from E+ and E- livers, then subject to evaluation of gene expression of ttpa and other genes of the CRAL_TRIO family, genes of antioxidant markers, and genes related to lipid metabolism. Ttpa expression was not altered by vitamin E status. However, one member of the CRAL_TRIO family, tyrosine-protein phosphatase non-receptor type 9 gene (ptpn9a), showed a 2.4-fold increase (P = 0.029) in E- relative to E+ livers. Further, we identified that the gene for choline kinase alpha (chka) showed a 3.0-fold increase (P = 0.010) in E- livers. These outcomes are consistent with our previous findings that show vitamin E deficiency increased lipid peroxidation causing increases in phospholipid turnover.

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.

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 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.

Early life stage transient aristolochic acid exposure induces behavioral hyperactivity but not nephrotoxicity in larval zebrafish

Aristolochic acids (AA) are nitrophenanthrene carboxylic acids found in plants of the Aristolochiaceae family. Humans are exposed to AA by deliberately taking herbal medicines or unintentionally as a result of environmental contamination. AA is notorious for its nephrotoxicity, however, fewer studies explore potential neurotoxicity associated with AA exposure. The developing nervous system is vulnerable to xenobiotics, and pregnant women exposed to AA may put their fetuses at risk. In the present study, we used the embryonic zebrafish model to evaluate the developmental neurotoxicity associated with AA exposure. At non-teratogenic concentrations (≤ 4 µM), continuous AA exposure from 8 to 120 hours post fertilization (hpf) resulted in larval hyperactivity that was characterized by increased moving distance, elevated activity and faster swimming speeds in several behavioral assays. Further analysis revealed that 8-24 hpf is the most sensitive exposure window for AA-induced hyperactivity. AA exposures specifically increased motor neuron proliferation, increased apoptosis in the eye, and resulted in cellular oxidative stress. In addition, AA exposures increased larval eye size and perturbed the expression of vision genes. Our study, for the first time, demonstrates that AA is neurotoxic to the developmental zebrafish with a sensitive window distinct from its well-documented nephrotoxicity.

Leveraging high-throughput screening data, deep neural networks, and conditional generative adversarial networks to advance predictive toxicology

There are currently 85,000 chemicals registered with the Environmental Protection Agency (EPA) under the Toxic Substances Control Act, but only a small fraction have measured toxicological data. To address this gap, high-throughput screening (HTS) and computational methods are vital. As part of one such HTS effort, embryonic zebrafish were used to examine a suite of morphological and mortality endpoints at six concentrations from over 1,000 unique chemicals found in the ToxCast library (phase 1 and 2). We hypothesized that by using a conditional generative adversarial network (cGAN) or deep neural networks (DNN), and leveraging this large set of toxicity data we could efficiently predict toxic outcomes of untested chemicals. Utilizing a novel method in this space, we converted the 3D structural information into a weighted set of points while retaining all information about the structure. In vivo toxicity and chemical data were used to train two neural network generators. The first was a DNN (Go-ZT) while the second utilized cGAN architecture (GAN-ZT) to train generators to produce toxicity data. Our results showed that Go-ZT significantly outperformed the cGAN, support vector machine, random forest and multilayer perceptron models in cross-validation, and when tested against an external test dataset. By combining both Go-ZT and GAN-ZT, our consensus model improved the SE, SP, PPV, and Kappa, to 71.4%, 95.9%, 71.4% and 0.673, respectively, resulting in an area under the receiver operating characteristic (AUROC) of 0.837. Considering their potential use as prescreening tools, these models could provide in vivo toxicity predictions and insight into the hundreds of thousands of untested chemicals to prioritize compounds for HT testing.

Concurrent Evaluation of Mortality and Behavioral Responses: A Fast and Efficient Testing Approach for High-Throughput Chemical Hazard Identification

The continual introduction of new chemicals into the market necessitates fast, efficient testing strategies for evaluating their toxicity. Ideally, these high-throughput screening (HTS) methods should capture the entirety of biological complexity while minimizing reliance on expensive resources that are required to assess diverse phenotypic endpoints. In recent years, the zebrafish (Danio rerio) has become a preferred vertebrate model to conduct rapid in vivo toxicity tests. Previously, using HTS data on 1060 chemicals tested as part of the ToxCast program, we showed that early, 24 h post-fertilization (hpf), behavioral responses of zebrafish embryos are predictive of later, 120 h post-fertilization, adverse developmental endpoints—indicating that embryonic behavior is a useful endpoint related to observable morphological effects. Here, our goal was to assess the contributions (i.e., information gain) from multiple phenotypic data streams and propose a framework for efficient identification of chemical hazards. We systematically swept through analysis parameters for data on 24 hpf behavior, 120 hpf behavior, and 120 hpf morphology to optimize settings for each of these assays. We evaluated the concordance of data from behavioral assays with that from morphology. We found that combining information from behavioral and mortality assessments captures early signals of potential chemical hazards, obviating the need to evaluate a comprehensive suite of morphological endpoints in initial screens for toxicity. We have demonstrated that such a screening strategy is useful for detecting compounds that elicit adverse morphological responses, in addition to identifying hazardous compounds that do not disrupt the underlying morphology. The application of this design for rapid preliminary toxicity screening will accelerate chemical testing and aid in prioritizing chemicals for risk assessment.

Developmental co-exposure of TBBPA and titanium dioxide nanoparticle induced behavioral deficits in larval zebrafish

Both tetrabromobisphenol A (TBBPA) and titanium dioxide nanoparticle (TiO₂ NP) have widespread commercial applications, resulting in their ubiquitous co-presence in the environment and biota. Although environmental chemicals exist as mixtures, toxicity studies are nearly always conducted with single chemicals. Few studies explore potential interactions of different chemical mixtures. In this study, we employ the sensitive developing nerve system in zebrafish to assess the neurotoxicity of TBBPA/TiO₂ NP mixtures. Specifically, zebrafish embryos were exposed to solvent control (0.1% DMSO), 2 μM TBBPA, 0.1 mg/L TiO₂ NP, and their mixture from 8 to 120 h post fertilization (hpf), and motor/social behavioral assessments were conducted on embryos/larvae at different developmental stages. Our results showed that TBBPA/TiO₂ NP single or co-exposures increased spontaneous movement, decreased touch response and swim speed, and affected social behaviors of light/dark preference, shoaling, mirror attack and social contact. In particular, many of these phenotypes were manifested with higher magnitude of changes from the mixture exposure. These behavioral deficits were also accompanied with increased cell death in olfactory region and neuromasts in the lateral line system, increased ROS in gallbladder, pancreas, liver, and intestine, as well as increased lipid peroxidation and decreased ATP levels in whole larval tissue homogenates. Further, genes coding for key cell apoptosis marker and antioxidant enzyme were significantly upregulated by these two chemicals, in particular to their mixture. Interestingly, the co-presence of TBBPA also increased the mean particle size of TiO₂ NP in the exposure solutions and the TiO₂ NP content in larval tissue. Together, our analysis suggests that TBBPA/TiO₂ NP induced behavioral changes may be due to physical accumulation of these two chemicals in the target organs, and TiO₂ NP may serve as carriers for increased accumulation of TBBPA. To conclude, we demonstrated that TBBPA/TiO₂ NP together cause increased bioaccumulation of TiO₂, and heightened responses in behavior, cell apoptosis and oxidative stress. Our findings also highlight the importance of toxicity assessment using chemical mixtures.

Nitrate-induced improvements in exercise performance are coincident with exuberant changes in metabolic genes and the metabolome in zebrafish (Danio rerio) skeletal muscle

Dietary nitrate supplementation improves exercise performance by reducing the oxygen cost of exercise and enhancing skeletal muscle function. However, the mechanisms underlying these effects are not well understood. The purpose of this study was to assess changes in skeletal muscle energy metabolism associated with exercise performance in a zebrafish model. Fish were exposed to sodium nitrate (60.7 mg/L, 303.5 mg/L, 606.9 mg/L), or control water, for 21 days and analyzed at intervals (5, 10, 20, 30, 40 cm/s) during a 2-h strenuous exercise test. We measured oxygen consumption during an exercise test and assessed muscle nitrate concentrations, gene expression, and the muscle metabolome before, during, and after exercise. Nitrate exposure reduced the oxygen cost of exercise and increased muscle nitrate concentrations at rest, which were reduced with increasing exercise duration. In skeletal muscle, nitrate treatment upregulated expression of genes central to nutrient sensing (mtor), redox signaling (nrf2a), and muscle differentiation (sox6). In rested muscle, nitrate treatment increased phosphocreatine (P = 0.002), creatine (P = 0.0005), ATP (P = 0.0008), ADP (P = 0.002), and AMP (P = 0.004) compared with rested-control muscle. Following the highest swimming speed, concentration of phosphocreatine (P = 8.0 × 10^-5), creatine (P = 6.0 × 10^-7), ATP (P = 2.0 × 10^-6), ADP (P = 0.0002), and AMP (P = 0.004) decreased compared with rested nitrate muscle. Our data suggest nitrate exposure in zebrafish lowers the oxygen cost of exercise by changing the metabolic programming of muscle prior to exercise and increasing availability of energy-rich metabolites required for exercise.NEW & NOTEWORTHY We show that skeletal muscle nitrate concentration is higher with supplementation at rest and was lower in groups with increasing exercise duration in a zebrafish model. The higher availability of nitrate at rest is associated with upregulation of key nutrient-sensing genes and greater availability of energy-producing metabolites (i.e., ATP, phosphocreatine, glycolytic intermediates). Overall, nitrate supplementation may lower oxygen cost of exercise through improved fuel availability resulting from metabolic programming of muscle prior to exercise.

Behavior Effects of Structurally Diverse Per- and Polyfluoroalkyl Substances in Zebrafish

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously detected in the environment, and some pose significant human and environmental health concerns globally. While some PFAS induce adverse health effects, relatively few toxicological studies adequately address the broad structural diversity of this chemical class. In the current study, we evaluated 58 individual PFAS spanning 14 structural subclasses and 2 mixtures at single concentrations for developmental toxicity in zebrafish using highly sensitive behavior endpoints. Following developmental exposure to PFAS, zebrafish were assessed for mortality and challenged with an embryonic photomotor response (EPR) assay at 24 h postfertilization (hpf) and with larval photomotor response (LPR) and larval startle response assays at 120 hpf. We found that none of the tested PFAS exposures elicited significant mortality or aberrant EPR; however, exposure to 21 individual PFAS from multiple structural subclasses and 1 mixture induced aberrant larval behavior. We then evaluated developmental toxicity across a concentration range of 0-100 μM for 10 perfluoroalkyl carboxylic acids (PFCAs; 4-carbon perfluorobutanoic acid through the 13-carbon perfluorotridecanoic acid). Exposure to the PFCAs did not cause significant mortality or morphological effects, with the exception of perfluorooctanoic acid and perfluorononanoic acid, and did not induce aberrant EPR. All PFCAs, except for longer-chain perfluorododecanoic acid caused abnormal LPR following exposure to at least one concentration. In this study, we evaluated a broad set of PFAS not previously assessed for in vivo sublethal behavior endpoints and confirmed previous findings that exposure to some PFAS induces abnormal behavior in developing zebrafish. The data from this study will guide the selection of PFAS for which to investigate modes of toxic action.

FutureTox IV Workshop Summary: Predictive Toxicology for Healthy Children

FutureTox IV, a Society of Toxicology Contemporary Concepts in Toxicology workshop, was held in November 2018. Building upon FutureTox I, II, and III, this conference focused on the latest science and technology for in vitro profiling and in silico modeling as it relates to predictive developmental and reproductive toxicity (DART). Publicly available high-throughput screening data sets are now available for broad in vitro profiling of bioactivities across large inventories of chemicals. Coupling this vast amount of mechanistic data with a deeper understanding of molecular embryology and post-natal development lays the groundwork for using new approach methodologies (NAMs) to evaluate chemical toxicity, drug efficacy, and safety assessment for embryo-fetal development. NAM is a term recently adopted in reference to any technology, methodology, approach, or combination thereof that can be used to provide information on chemical hazard and risk assessment to avoid the use of intact animals (U.S. Environmental Protection Agency [EPA], Strategic plan to promote the development and implementation of alternative test methods within the tsca program, 2018, https://www.epa.gov/sites/production/files/2018-06/documents/epa_alt_strat_plan_6-20-18_clean_final.pdf). There are challenges to implementing NAMs to evaluate chemicals for developmental toxicity compared with adult toxicity. This forum article reviews the 2018 workshop activities, highlighting challenges and opportunities for applying NAMs for adverse pregnancy outcomes (eg, preterm labor, malformations, low birth weight) as well as disorders manifesting postnatally (eg, neurodevelopmental impairment, breast cancer, cardiovascular disease, fertility). DART is an important concern for different regulatory statutes and test guidelines. Leveraging advancements in such approaches and the accompanying efficiencies to detecting potential hazards to human development are the unifying concepts toward implementing NAMs in DART testing. Although use of NAMs for higher level regulatory decision making is still on the horizon, the conference highlighted novel testing platforms and computational models that cover multiple levels of biological organization, with the unique temporal dynamics of embryonic development, and novel approaches for estimating toxicokinetic parameters essential in supporting in vitro to in vivo extrapolation.

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.

The Zebrafish Xenograft Models for Investigating Cancer and Cancer Therapeutics

Simple Summary

The identification and development of new anti-cancer drugs requires extensive testing in animal models to establish safety and efficacy of drug candidates. The transplantation of human tumor tissue into mouse (tumor xenografts) is commonly used to study cancer progression and to test potential drugs for their anti-cancer activity. Mouse models do not afford the ability to test a large number of drug candidates quickly as it takes several weeks to conduct these experiments. In contrast, tumor xenograft studies in zebrafish provide an efficient platform for rapid testing of safety and efficacy in less than two weeks.

Abstract

In order to develop new cancer therapeutics, rapid, reliable, and relevant biological models are required to screen and validate drug candidates for both efficacy and safety. In recent years, the zebrafish (Danio rerio) has emerged as an excellent model organism suited for these goals. Larval fish or immunocompromised adult fish are used to engraft human cancer cells and serve as a platform for screening potential drug candidates. With zebrafish sharing ~80% of disease-related orthologous genes with humans, they provide a low cost, high-throughput alternative to mouse xenografts that is relevant to human biology. In this review, we provide background on the methods and utility of zebrafish xenograft models in cancer research.

Developmental Hazard of Environmentally Persistent Free Radicals and Protective Effect of TEMPOL in Zebrafish Model

Environmentally persistent free radicals (EPFRs) can be detected in ambient PM_2.5, cigarette smoke, and soils and are formed through combustion and thermal processing of organic materials. The hazards of EPFRs are largely unknown. In this study, we assess the developmental toxicity of EPFRs and the ability of TEMPOL (4-Hydroxy-2,2,6,6-tetramethylpiperidine 1-oxyl) to protect against such hazards using zebrafish embryos. Particles containing EPFRs were acquired by dosing dichlorobenzene (DCB) vapor on the Cab-o-sil/5% CuO particles at 230 °C in vacuo (referred to as DCB-230). The particles were suspended in ultrapure water to make 1 mg/mL of stock solution from which series dilution was undertaken to obtain 10, 20, 30, 40, 50, 60, 80, and 100 µg/mL final test solutions, which were then placed in individual wells with a 4 h postfertilization (hpf) zebrafish embryo. Plates were run in duplicate to obtain a sample size of 24 animals per concentration; 12 embryos were exposed per concentration per plate. Statistical analysis of the morphology endpoints was performed. We investigated overt toxicity responses to DCB-230 in a 22-endpoint battery that included developing zebrafish from 24–120 hpf. Exposure to concentrations greater than 60 µg/mL of DCB-230 induced high mortality in the developmental zebrafish model. Exposure to EPFRs induced developmental hazards that were closely related to the concentrations of free radicals and EPFRs. The potential protective effects of TEMPOL against EPFRs’ toxicity in zebrafish were investigated. Exposure to EPFRs plus TEMPOL shifted the concentration to an induced 50% adverse effect (EC50), from 23.6 to 30.8 µg/mL, which verifies TEMPOL’s protective effect against EPFRs in the early phase of zebrafish development.

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.

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.

Vitamin E deficiency dysregulates thiols, amino acids and related molecules during zebrafish embryogenesis

Vitamin E (α-tocopherol, VitE) was discovered as a nutrient essential to protect fetuses, but its molecular role in embryogenesis remains undefined. We hypothesize that the increased lipid peroxidation due to VitE deficiency drives a complex mechanism of overlapping biochemical pathways needed to maintain glutathione (GSH) homeostasis that is dependent on betaine and its methyl group donation. We assess amino acids and thiol changes that occur during embryogenesis [12, 24 and 48 h post fertilization (hpf)] in VitE-sufficient (E+) and deficient (E-) embryos using two separate, novel protocols to quantitate changes using UPLC-MS/MS. Using partial least squares discriminant analysis, we found that betaine is a critical feature separating embryos by VitE status and is higher in E- embryos at all time points. Other important features include: glutamic acid, increased in E- embryos at 12 hpf; choline, decreased in E- embryos at 24 hpf; GSH, decreased in E- embryos at 48 hpf. By 48 hpf, GSH was significantly lower in E- embryos (P < 0.01), as were both S-adenosylmethionine (SAM, P < 0.05) and S-adenosylhomocysteine (SAH, P < 0.05), while glutamic acid was increased (P < 0.01). Since GSH synthesis requires cysteine (which was unchanged), these data suggest that both the conversion of homocysteine and the uptake of cystine via the X_c^- exchanger are dysregulated. Our data clearly demonstrates the highly inter-related dependence of methyl donors (choline, betaine, SAM) and the methionine cycle for maintenance of thiol homeostasis. Additional quantitative flux studies are needed to clarify the quantitative importance of these routes.