Description of Solvent-Extractable Chemicals in Thermal Receipts and Toxicological Assessment of Bisphenol S and Diphenyl Sulfone

Research on thermal receipts has previously focused on the toxic effects of dermal exposure from the most publicized developers (e.g., bisphenol A (BPA) and bisphenol S (BPS)), while no studies have reported on the other solvent-extractable compounds therein. Diphenyl sulfone (DPS) is a sensitizer added to thermal receipts, but little is known about DPS concentrations in receipts or potential toxicity. Here, we quantified BPA, BPS, and DPS concentrations and tentatively identified the solvent-extractable compounds of thermal receipts collected from three South Dakota (USA) cities during 2016-2017. An immortalized chicken hepatic cell line, cultured as 3D spheroids, was used to screen effects of DPS, BPS, and 17ß estradiol (E2; 0.1-1000 µM) on cell viability and gene expression changes. These chemicals elicited limited cytotoxicity with LC50 values ranging from 113 to 143 µM, and induced dysregulation in genes associated with lipid and bile acid homeostasis. Taken together, this study generated novel information on solvent-extractable chemicals from thermal receipts and toxicity data for DPS.

Hepatic Transcriptomic Responses to Ethinylestradiol in Embryonic Japanese Quail and Double-Crested Cormorant

Understanding species differences in sensitivity to toxicants is a critical issue in ecotoxicology. We recently established that double-crested cormorant (DCCO) embryos are more sensitive than Japanese quail (JQ) to the developmental effects of ethinylestradiol (EE2). We explored how this difference in sensitivity between species is reflected at a transcriptomic level. The EE2 was dissolved in dimethyl sulfoxide and injected into the air cell of eggs prior to incubation at nominal concentrations of 0, 3.33, and 33.3 µg/g egg weight. At midincubation (JQ 9 days; DCCO 16 days), livers were collected from five embryos/treatment group for RNA sequencing. Data were processed and analyzed using EcoOmicsAnalyst and ExpressAnalyst. The EE2 exposure dysregulated 238 and 1,987 genes in JQ and DCCO, respectively, with 78 genes in common between the two species. These included classic biomarkers of estrogen exposure such as vitellogenin and apovitellenin. We also report DCCO-specific dysregulation of Phase I/II enzyme-coding genes and species-specific transcriptional ontogeny of vitellogenin-2. Twelve Kyoto Encyclopedia of Genes and Genomes pathways and two EcoToxModules were dysregulated in common in both species including the peroxisome proliferator-activated receptor (PPAR) signaling pathway and fatty acid metabolism. Similar to previously reported differences at the organismal level, DCCO were more responsive to EE2 exposure than JQ at the gene expression level. Our description of differences in transcriptional responses to EE2 in early life stage birds may contribute to a better understanding of the molecular basis for species differences. Environ Toxicol Chem 2024;00:1-12. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Contaminant biomonitoring augmented with a qPCR array indicates hepatic mRNA gene expression effects in wild-collected seabird embryos

Birds can bioaccumulate persistent contaminants, and maternal transfer to eggs may expose embryos to concentrations sufficient to cause adverse effects during sensitive early-life stages. However, using tissue residue concentrations alone to infer whether contaminant effects are occurring suffers from uncertainty, and efficient, sensitive biomarkers remain limited in wildlife. We studied relationships between whole embryo contaminant concentrations (total mercury, organochlorine pesticides, perfluoroalkyl substances, polychlorinated biphenyls, and halogenated flame retardants) together with mRNA expression in embryonic liver tissue from a Pacific Ocean seabird, the rhinoceros auklet (Cerorhinca monocerata). Fresh eggs were collected, incubated under controlled conditions, and from the pre-hatch embryo, hepatic RNA was extracted for qPCR array analysis to measure gene expression (2-∆Cq), while the remaining embryo was analyzed for contaminant residues. Contaminant and gene expression data were assessed with a combination of multivariate approaches and linear models. Results indicated correlations between embryonic total mercury and several genes such as sepp1, which encodes selenoprotein P. Correlation between the biotransformation gene cyp1a4 and the C7 perfluoroalkyl carboxylic acid PFHpA was also evident. This study demonstrates that egg collection from free-living populations for contaminant biomonitoring programs can relate chemical residues to in ovo mRNA gene expression effects in embryo hepatic tissue.

Comparing Transcriptomic Responses to Chemicals Across Six Species Using the EcoToxChip RNASeq Database

The EcoToxChip project includes RNA-sequencing data from experiments involving model (Japanese quail, fathead minnow, African clawed frog) and ecological (double-crested cormorant, rainbow trout, northern leopard frog) species at multiple life stages (whole embryo and adult) exposed to eight chemicals of environmental concern known to perturb a wide range of biological systems (ethinyl estradiol, hexabromocyclododecane, lead, selenomethionine, 17β trenbolone, chlorpyrifos, fluoxetine, and benzo[a]pyrene). The objectives of this short communication were to (1) present and make available this RNA-sequencing database (i.e., 724 samples from 49 experiments) under the FAIR principles (FAIR data are data which meet principles of findability, accessibility, interoperability, and reusability), while also summarizing key meta-data attributes and (2) use ExpressAnalyst (including the Seq2Fun algorithm and EcoOmicsDB) to perform a comparative transcriptomics analysis of this database focusing on baseline and differential transcriptomic changes across species-life stage-chemical combinations. The database is available in NCBI GEO under accession number GSE239776. Across all species, the number of raw reads per sample ranged between 13 and 58 million, with 30% to 79% of clean reads mapped to the "vertebrate" subgroup database in EcoOmicsDB. Principal component analyses of the reads illustrated separation across the three taxonomic groups as well as some between tissue types. The most common differentially expressed gene was CYP1A1 followed by CTSE, FAM20CL, MYC, ST1S3, RIPK4, VTG1, and VIT2. The most common enriched pathways were metabolic pathways, biosynthesis of cofactors and biosynthesis of secondary metabolites, and chemical carcinogenesis, drug metabolism, and metabolism of xenobiotics by cytochrome P450. The RNA-sequencing database in the present study may be used by the research community for multiple purposes, including, for example, cross-species investigations, in-depth analyses of a particular test compound, and transcriptomic meta-analyses. Environ Toxicol Chem 2024;00:1-6. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Development and initial testing of EcoToxChip, a novel toxicogenomics tool for environmental management and chemical risk assessment

New approach methods (NAMs) are increasingly important to help accelerate the pace of ecological risk assessment and offer more ethical, affordable, and efficient alternatives to traditional toxicity tests. In the present study, we describe the development, technical characterization and initial testing of a toxicogenomics tool, EcoToxChip (384-well qPCR array), to support chemical management and environmental monitoring for three laboratory model species - fathead minnow (Pimephales promelas), African clawed frog (Xenopus laevis), and Japanese quail (Coturnix japonica). Chip design, including gene selection, was informed by a diverse end-user group and quality control metrics (e.g. primer assay, reverse transcription and PCR efficiency) performed well based on a priori established criteria. Correlation with RNAseq data provided additional confidence in this novel toxicogenomics tool. While the present study represents an initial testing of only 24 EcoToxChips for each of the model species, the results provide increased confidence in the robustness/reproducibility of EcoToxChips for evaluating perturbations in gene expression associated with chemical exposure and thus, this NAM, combined with early-life stage toxicity testing, could augment current efforts for chemical prioritization and environmental management.

Bioaccumulation, Trophic Transfer, and Biotransformation of Polychlorinated Diphenyl Ethers in a Simulated Aquatic Food Chain

Polychlorinated diphenyl ethers (PCDEs) are detected in aquatic environments and demonstrate adverse effects in aquatic organisms. However, data regarding the environmental behavior of PCDEs in aquatic ecosystems are lacking. In the present study, a simulated aquatic food chain (Scenedesmus obliquus-Daphnia magna-Danio rerio) was constructed in a lab setting, and the bioaccumulation, trophic transfer, and biotransformation of 12 PCDE congeners were quantitatively investigated for the first time. The log-transformed bioaccumulation factors (BCFs) of PCDEs in S. obliquus, D. magna, and D. rerio were in the range of 2.94-3.77, 3.29-4.03, and 2.42-2.89 L/kg w.w., respectively, indicating the species-specific bioaccumulation of PCDE congeners. The BCF values increased significantly with the increasing number of substituted Cl atoms, with the exception of CDE 209. The number of Cl atoms at the para and meta positions were found to be the major positive contributing factors for BCFs in the case of the same number of substituted Cl. The lipid-normalized biomagnification factors (BMFs) of S. obliquus to D. magna, D. magna to D. rerio, and the whole food chain for the 12 PCDE congeners ranged at 1.08-2.27, 0.81-1.64, and 0.88-3.64, respectively, suggesting that some congeners had BMFs comparable to PBDEs and PCBs. Dechlorination was the only metabolic pathway observed for S. obliquus and D. magna. For D. rerio, dechlorination, methoxylation, and hydroxylation metabolic pathways were observed. ¹H nuclear magnetic resonance (NMR) experiments and theoretical calculations confirmed that methoxylation and hydroxylation occurred at the ortho position of the benzene rings. In addition, reliable quantitative structure-property relationship (QSPR) models were constructed to qualitatively describe the relationships between molecular structure descriptors and BCFs for PCDEs. These findings provide insights into the movement and transformation of PCDEs in aquatic ecosystems.

Maternal Transfer and Apical and Physiological Effects of Dietary Hexabromocyclododecane Exposure in Parental Fathead Minnows (Pimephales promelas)

Hexabromocyclododecane (HBCD) is a persistent organic pollutant that has been characterized as an endocrine disruptor, undergoes maternal transfer, and hinders development and growth in oviparous organisms. The present study examined the apical effects of dietary HBCD (11.5, 36.4, 106 mg/kg, wet wt) on adult fathead minnow exposed for 49 days and the subsequent accumulation and maternal transfer kinetics in adult tissue and eggs, respectively. Exposed adults displayed a significant increase in egg production in the medium treatment group, but no other significant effects were noted. Maternal transfer of dietary HBCD had a similar egg-to-muscle ratios (EMR) in the low and medium treatment groups (1.65 and 1.27 [wet wt], respectively). However, the high treatment group deviated from other treatments with an EMR of 4.2 (wet wt), potentially due to differences in total lipid content in food and/or reaching diffusion/lipid saturation limits in adult tissue, resulting in lower accumulation in the adult muscle tissue. A positive correlation was observed between egg HBCD concentration and time of exposure, which indicates that maternal transfer of HBCD is of concern in fish, and further studies should be conducted to fully elucidate the potential adverse effects that may be observed in the early life stage of oviparous organisms. Environ Toxicol Chem 2023;42:143-153. © 2022 SETAC.

Characterizing Variability and Uncertainty Associated with Transcriptomic Dose-Response Modeling

Transcriptomics dose-response analysis (TDRA) has emerged as a promising approach for integrating toxicogenomics data into a risk assessment context; however, variability and uncertainty associated with experimental design are not well understood. Here, we evaluated n = 55 RNA-seq profiles derived from Japanese quail liver tissue following exposure to chlorpyrifos (0, 0.04, 0.1, 0.2, 0.4, 1, 2, 4, 10, 20, and 40 μg/g; n = 5 replicates per group) via egg injection. The full dataset was subsampled 637 times to generate smaller datasets with different dose ranges and spacing (designs A-E) and number of replicates (n = 2-5). TDRA of the 637 datasets revealed substantial variability in the gene and pathway benchmark doses, but relative stability in overall transcriptomic point-of-departure (tPOD) values when tPODs were calculated with the "pathway" and "mode" methods. Further, we found that tPOD values were more dependent on the dose range and spacing than on the number of replicates, suggesting that optimal experimental designs should use fewer replicates (n = 2 or 3) and more dose groups to reduce uncertainty in the results. Finally, tPOD values ranged by over ten times for all surveyed experimental designs and tPOD types, suggesting that tPODs should be interpreted as order-of-magnitude estimates.

Hepatic Transcriptomic Responses to Ethinylestradiol in Two Life Stages of Japanese Quail

Chemical risk assessment for avian species typically depends on information from toxicity tests performed in adult birds. Early-life stage (ELS) toxicity tests have been proposed as an alternative, but incorporation of these data into existing frameworks will require knowledge about the similarities/differences between ELS and adult responses. The present study uses transcriptomics to assess hepatic gene expression in ELS and adult Japanese quail following exposure to ethinylestradiol (EE2). Prior to incubation, ELS quail were dosed with measured EE2 concentrations of 0.54, 6.3, and 54.2 µg/g egg weight via air cell injection. Adult quail were fed a single dose of EE2 at nominal concentrations of 0, 0.5, and 5 mg/kg body weight by gavage. Liver tissue was collected from five to six individuals per dose group at mid-incubation for ELS quail and 4 days after dosing for adults. A total of 283 and 111 differentially expressed genes (DEGs) were detected in ELS and adult quail, respectively, 16 of which were shared across life stages. Shared DEGs included estrogenic biomarkers such as vitellogenin genes and apovitellenin-1. For the dose groups that resulted in the highest number of DEGs (ELS, 6.3 µg/g; adult, 5 mg/kg), 21 and 35 Kyoto Encyclopedia of Genes and Genomes pathways were enriched, respectively. Ten of these pathways were shared between life stages, including pathways involved with signaling molecules and interaction and the endocrine system. Taken together, our results suggest conserved mechanisms of action following estrogenic exposure across two life stages, with evidence from differential expression of key biomarker genes and enriched pathways. The present study contributes to the development and evaluation of ELS tests and toxicogenomic approaches and highlights their combined potential for screening estrogenic chemicals. Environ Toxicol Chem 2022;41:2769-2781. © 2022 SETAC.

Consideration of metabolomics and transcriptomics data in the context of using avian embryos for toxicity testing

Early-life stage (ELS) avian toxicity tests have been proposed as a more ethical alternative to traditional standardized tests with adult birds. At the same time, 'omics approaches are gaining traction in the field of avian toxicology, but little has been done to characterize the metabolome and transcriptome at different life stages. The present study uses 'omics data from toxicity tests of 8 environmental chemicals in ELS and adult Japanese quail (Coturnix japonica) to address this data gap. Previous analyses of these data focused on responses to each of the individual chemicals. Here, we consider data from all studies to describe variation in the metabolome and transcriptome between life stages and across independent experiments, irrespective of chemical treatment. Of the 230 metabolites detected in liver, 163 were shared between the two life stages. However, many of the targeted bile acids that were present in the adult liver were absent from ELS samples. For the transcriptome, >90% of the 18,364 detected transcripts were common to both life stages. Based on the 213 genes solely detected in ELS liver, the neuroactive ligand-receptor interaction pathway was significantly enriched. Multivariate and hierarchical clustering analyses revealed that variability among independent experiments was higher for the adult than the ELS studies at both the metabolomic and transcriptomic levels. Our results indicate concordance of the two approaches, with less variation between independent experiments in the ELS metabolome and transcriptome than in adults, lending support for the use of ELS as an alternative toxicity testing strategy.

Chemical risk governance: Exploring stakeholder participation in Canada, the USA, and the EU

Hazardous chemicals are one of the greatest environmental challenges facing our planet, testing governments in the face of economic and social development. Chemical risks are often complex systemic risks, which require particular governance processes, stakeholder participation mechanisms, and communication procedures to manage. In this article we explore how such processes, mechanisms, and procedures could offer more effective pathways for international chemical risk assessment and management. Combining document analysis with 69 key informant interviews and four focus groups with experts from across government, non-governmental organizations (NGOs) and industry in Canada, the USA and Europe (n = 24), the results highlight the importance of increasing stakeholder participation and risk communication in existing chemical risk assessment and management frameworks. Similarities and differences between the three jurisdictions are discussed with a view to inform future risk governance strategies. The contemporary insights are of relevance to decision-makers seeking to improve the effectiveness of chemical risk governance in practice.

Occurrence, partitioning, and bioaccumulation of an emerging class of PBT substances (polychlorinated diphenyl sulfides) in Chaohu Lake, Southeast China

Polychlorinated diphenyl sulfides (PCDPSs) represent an emerging group of constituents that are persistent, bioaccumulative and toxic (PBT) substances of great concern in terms of human health and ecological integrity. However, little is known about the occurrence, environmental behaviour and ecological risks of PCDPSs in lake environments. In this study, the concentrations of 21 PCDPSs were determined in surface water, suspended particulate matter (SPM), sediments, and 8 fish species from Chaohu Lake, China. Eighteen PCDPS congeners were prevalently detected in the samples, with concentrations ranging from 0.272-1.69 ng/L (water), 0.477-2.03 ng/g d.w. (SPM), 0.719-4.07 ng/g d.w. (sediment) and 0-0.131 ng/g w.w. (fish), respectively. Medium- and high-chlorinated PCDPSs in SPM and sediment were significantly higher than those in water samples. Increased PCDPS concentrations were found in higher trophic level fishes and those with a demersal habitat preference, indicating their bioaccumulation and biomagnification potential. The logBCFs, BSSAFs, and BSAFs of PCDPS congeners in fishes were determined to be 3.91-5.18, 0.0500-2.33, and 0.0360-4.94 L/kg, respectively. The organic carbon normalized partition coefficients (logK_oc) of PCDPSs in surface water-SPM (4.61-5.54 L/g) and surface water-sediment (4.38-5.69 L/g) systems were determined, and it was found that highly chlorinated PCDPSs were more prone to migrate from water to sediment and SPM. The toxic equivalent (TEQ) values of PCDPSs in the samples (lower than 10^-1 pg/g or pg/L) and daily intake via fish consumption (0.180-0.340 μg/kg/day) were estimated for humans, and cumulative risk quotients (RQs) after correction at ten sampling sites (0.065-0.66) were calculated for green algae. The findings elucidated the environmental behaviour of PCDPSs in Chaohu Lake.

Toxicity screening of bisphenol A replacement compounds: cytotoxicity and mRNA expression in LMH 3D spheroids

Previously, we showed that the chicken LMH cell line cultured as 3D spheroids may be a suitable animal free alternative to primary chicken embryonic hepatocytes (CEH) for avian in vitro chemical screening. In this study, cytotoxicity and mRNA expression were determined in LMH 3D spheroids following exposure to bisphenol A (BPA), five BPA replacement compounds (BPF, TGSH, DD-70, BPAF, BPSIP), and 17β estradiol (E2). Results were compared to an earlier study that evaluated the same endpoints for these chemicals in CEH. BPA and the replacement compounds had LC50 values ranging from 16.6 to 81.8 μM; DD-70 and BPAF were the most cytotoxic replacements (LC50 = 17.23 ± 4.51 and 16.6 ± 4.78 μM). TGSH and DD-70 modulated the greatest number of genes, although fewer than observed in CEH. Based on the expression of apovitellenin and vitellogenin, BPAF was the most estrogenic compound followed by BPF, BPSIP, and BPA. More estrogen-responsive genes were modulated in LMH spheroids compared to CEH. Concentration-dependent gene expression revealed that DD-70 and BPAF altered genes related to lipid and bile acid regulation. Overall, cytotoxicity and clustering of replacements based on gene expression profiles were similar between LMH spheroids and CEH. In addition to generating novel gene expression data for five BPA replacement compounds in an in vitro avian model, this research demonstrates that LMH spheroids may represent a useful animal free alternative for avian toxicity testing.

Characterizing toxicity pathways of fluoxetine to predict adverse outcomes in adult fathead minnows (Pimephales promelas)

Current ecotoxicity testing programs are impeded as they predominantly rely on slow and expensive animal tests measuring adverse outcomes. Therefore, new approach methodologies (NAMs) increasingly involve short-term mechanistic assays that employ molecular endpoints to predict adverse outcomes of regulatory relevance. This study aimed to elucidate the application of NAMs in adult fathead minnows using fluoxetine (FLX) as a model compound. Fish were exposed to three FLX concentrations (measured: 2.42, 10.7, and 56.7 μgL^-1) and a control. After 96 h, molecular toxicity signatures were characterized using proteomics and transcriptomics analyses in livers and brains of a sub-set of fish. The remaining fish were sampled at 21 days and assessed for liver histopathology and morphometric measurements. Fecundity was monitored throughout the study. In the livers, 56.7 μgL^-1 FLX caused enrichment of PPAR signaling in the proteome and fatty acid-related pathways in the transcriptome, potential upstream responses that led to lipid-type vacuolation of hepatocytes, observed via histopathology. Upregulated genes in the brain suggested alterations in serotonin-related signaling processes and reproductive behaviour, which may explain the observed significant decrease in fecundity. While the relationships between molecular responses and adverse outcomes remain complex, this research provided important insights into the mechanistic toxicity of FLX.

Developmental and Hepatic Gene Expression Changes in Chicken Embryos Exposed to p-Tert-Butylphenyl Diphenyl Phosphate and Isopropylphenyl Phosphate via Egg Injection

Organophosphate flame retardants (OPFRs) are used in a variety of products such as clear coats, resins, and plastics; however, research into their toxicological effects is limited. p-Tert-butylphenyl diphenyl phosphate (BPDP) and isopropylphenyl phosphate (IPPP) are two OPFRs that were prioritized for whole-animal toxicological studies based on observed effects in cultured avian hepatocytes in a previous study. The present study investigates the toxicity of BPDP and IPPP in chicken embryos at different developmental stages by evaluating morphological and gene expression endpoints. Chicken eggs were exposed via air cell injection to 0-250 μg/g (nominal) of either compound and then artificially incubated. At day 11 (midincubation), liver samples were collected for mRNA expression analysis; and at day 20 (1 day prehatch), morphological measurements and liver samples for transcriptomic evaluation were collected. At 250 μg/g, gallbladder size was significantly reduced for both compounds, head/bill length and tarsus length were significantly decreased, and liver somatic index was significantly increased following IPPP exposure only. No effects on mortality were observed up to the highest administered concentration for either chemical. Using a ToxChip polymerase chain reaction array, we report significant differences in hepatic gene expression for both compounds and time points; the most pronounced transcriptomic effects occurred at midincubation. Genes related to xenobiotic metabolism, bile acid/cholesterol regulation, and oxidative stress were significantly dysregulated. Given these changes observed throughout avian embryonic development, further research into the long-term effects of BPDP and IPPP are warranted, especially as they pertain to liver cholestasis. Environ Toxicol Chem 2022;41:739-747. © 2021 Her Majesty the Queen in Right of Canada. Environmental Toxicology and Chemistry © 2021 SETAC. Reproduced with the permission of the Minister of Environment and Climate Change Canada.

Metabolomic profiles in relation to benchmark polycyclic aromatic compounds (PACs) and trace elements in two seabird species from Arctic Canada

While exposure of birds to oil-related contaminants has been documented, the related adverse effects this exposure has on Arctic marine birds remain unexplored. Metabolomics can play an important role to explore biologically relevant metabolite biomarkers in relation to different stressors, even at benchmark levels of contamination. The aim of this study was to characterize the metabolomics profiles in relation to polycyclic aromatic compounds (PACs) and trace elements in the liver of two seabird species in the Canadian Arctic. In July 2018, black guillemots (Cepphus grylle) and thick-billed murres (Uria lomvia) were collected by hunters from a region where natural oil seeps occur in the seabed near Qikiqtarjuaq, Nunavut, Canada. A total of 121 metabolites were identified in liver tissue samples using reversed phase and hydrophilic interaction liquid chromatography coupled to high resolution mass spectrometry platforms to detect non-polar and polar metabolites, respectively. Sixty-nine metabolites showed excellent repeatability and linearity and were used to examine possible effects of oil-related contaminants exposure (PACs and trace elements). Metabolites including 3-hydroxy anthranilic acid, adenine, adenosine, adenosine mono-phosphate, ascorbic acid, butyrylcarnitine, cholic acid, guanosine, guanosine mono-phosphate, inosine, norepinephrine and threonine showed significant differences (more than two fold) between the two species. Elevated adenine and adenosine, along with decreased reduced/oxidized glutathione ratio, highlighted the potential for oxidative stress in murres. Lipid peroxidation and superoxide dismutase activity assays also confirmed these metabolomic findings. These results will help to characterize the baseline metabolomic profiles of Arctic seabird species with different foraging behaviour and trace element burden.

Comparative analysis of transcriptomic points-of-departure (tPODs) and apical responses in embryo-larval fathead minnows exposed to fluoxetine

Current approaches in chemical hazard assessment face significant challenges because they rely on live animal testing, which is time-consuming, expensive, and ethically questionable. These concerns serve as an impetus to develop new approach methodologies (NAMs) that do not rely on live animal tests. This study explored a molecular benchmark dose (BMD) approach using a 7-day embryo-larval fathead minnow (FHM) assay to derive transcriptomic points-of-departure (tPODs) to predict apical BMDs of fluoxetine (FLX), a highly prescribed and potent selective serotonin reuptake inhibitor frequently detected in surface waters. Fertilized FHM embryos were exposed to graded concentrations of FLX (confirmed at < LOD, 0.19, 0.74, 3.38, 10.2, 47.5 μg/L) for 32 days. Subsets of fish were subjected to omics and locomotor analyses at 7 days post-fertilization (dpf) and to histological and biometric measurements at 32 dpf. Enrichment analyses of transcriptomics and proteomics data revealed significant perturbations in gene sets associated with serotonergic and axonal functions. BMD analysis resulted in tPOD values of 0.56 μg/L (median of the 20 most sensitive gene-level BMDs), 5.0 μg/L (tenth percentile of all gene-level BMDs), 7.51 μg/L (mode of the first peak of all gene-level BMDs), and 5.66 μg/L (pathway-level BMD). These tPODs were protective of locomotor and reduced body weight effects (LOEC of 10.2 μg/L) observed in this study and were reflective of chronic apical BMDs of FLX reported in the literature. Furthermore, the distribution of gene-level BMDs followed a bimodal pattern, revealing disruption of sensitive neurotoxic pathways at low concentrations and metabolic pathway perturbations at higher concentrations. This is one of the first studies to derive protective tPODs for FLX using a short-term embryo assay at a life stage not considered to be a live animal under current legislations.

Acute toxicity of polychlorinated diphenyl ethers (PCDEs) in three model aquatic organisms (Scenedesmus obliquus, Daphnia magna, and Danio rerio) of different trophic levels

The frequent detection of polychlorinated diphenyl ethers (PCDEs) in aquatic systems has aroused widespread concerns, however, their potential hazard to aquatic ecosystems has been poorly understood. Here the acute toxicity of 12 PCDE congeners was evaluated in three model aquatic organisms representing different trophic levels following OECD test guidelines, including green algae (Scenedesmus obliquus), water flea (Daphnia magna), and zebrafish (Danio rerio). Dose-dependent increases in growth inhibition and mortality were observed for all tested PCDE congeners. Most of the PCDE congeners, in particular 3,3',4,4'-tetra-CDE, were highly toxic to the three aquatic organisms with EC₅₀ or LC₅₀ values below 1 mg L^-1. Their toxicities were generally comparable with those of certain polychlorinated biphenyls (PCBs) and polybrominated diphenyl ethers (PBDEs). Moreover, D. magna was the most sensitive species among the three aquatic organisms. In addition, the EC₅₀ or LC₅₀ values had an extremely significant correlation with the n-octanol-water partition coefficient (logK_ow) of the PCDE congeners. The established quantitative structure-property relationship (QSPR) models indicated that the molecular polarizability (α) could significantly influence the acute toxicity of PCDEs on Daphnia magna and Danio rerio, and the energy of the lowest unoccupied molecular orbital (E_LUMO) is the key factor of the acute toxicity of PCDEs in Scenedesmus obliquus. In addition, even at environmental levels, 3,3',4,4'-tetra-CDE could induced seveve oxidative damages in the three aquactic species. These findings would contribute to the understanding of adverse effects of PCDEs in aquatic organisms.

Effects of Avian Eggshell Oiling With Diluted Bitumen Show Sublethal Embryonic Polycyclic Aromatic Compound Exposure

Breeding birds that become oiled may contaminate the shells of their eggs, and studies of conventional crude oil suggest that even small quantities can be absorbed through the eggshell and cause embryotoxicity. Unconventional crude oils remain untested, so we evaluated whether a major Canadian oil sands product, diluted bitumen (dilbit), would be absorbed and cause toxicity when applied to eggshells of two species, domestic chicken (Gallus gallus domesticus) and double-crested cormorant (Nannopterum auritum). We artificially incubated eggs and applied lightly weathered dilbit (Cold Lake blend) to the eggshells (0.015-0.15 mg g^-1 egg in chicken; 0.1-0.4 mg g^-1 egg in cormorant) at various points during incubation before sampling prehatch embryos. Polycyclic aromatic compound (PAC) residue in cormorant embryos was elevated only at the highest dilbit application (0.4 mg g^-1 egg) closest (day 16) to sampling on day 22. In contrast, cormorant liver cytochrome P450 1a4 (Cyp1a4) mRNA expression (quantitative polymerase chain reaction assay) was elevated only in embryos treated with the earliest and lowest dilbit application (0.1 mg g^-1 egg on day 4). These results confirm that dilbit can cross through the eggshell and be absorbed by embryos, and they imply rapid biotransformation of PACs and a nonmonotonic Cyp1a4 response. Despite evidence of exposure in cormorant, we found no detectable effects on the frequency of survival, deformity, and gross lesions, nor did we find effects on physiological endpoints indicative of growth and cardiovascular function in either chicken or cormorant. In ovo dilbit exposure may be less toxic than well-studied conventional crude oils. The effects of an oil spill scenario involving dilbit to bird embryos might be subtle, and PACs may be rapidly metabolized. Environ Toxicol Chem 2022;41:159-174. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Targeted Metabolomics to Assess Exposure to Environmental Chemicals of Concern in Japanese Quail at Two Life Stages

This proof-of-concept study characterizes the Japanese quail (Coturnix japonica) hepatic metabolome following exposure to benzo[a]pyrene, chlorpyrifos, ethinylestradiol, fluoxetine hydrochloride, hexabromocyclododecane, lead(II)nitrate, seleno-L-methionine, and trenbolone in embryos and adults. The analysis revealed effects on lipid metabolism following exposure to several chemicals at both life stages. The most pronounced effects were observed in embryos exposed to 41.1 μg/g chlorpyrifos. This work highlighted challenges and the need for further avian metabolomics studies.

Using In Vitro and Machine Learning Approaches to Determine Species-Specific Dioxin-like Potency and Congener-Specific Relative Sensitivity among Birds for Brominated Dioxin Analogues

There is a paucity of experimental data regarding dioxin-like toxicity of polybrominated dibenzo-p-dioxins/dibenzofurans (PBDD/Fs) and non-ortho polybrominated biphenyls (PBBs). In this study, avian aryl hydrocarbon receptor 1 (AHR1)-luciferase reporter gene assays were used to determine their species-specific dioxin-like potencies (DLPs) and congener-specific interspecies relative sensitivities in birds. The results suggested that DLPs of the brominated congeners for chicken-like (Ile324_Ser380) species did not always follow World Health Organization toxicity equivalency factors of their chlorinated analogues. For ring-necked pheasant-like (Ile324_Ala380) and Japanese quail-like (Val324_Ala380) species, the difference in DLP for several congeners was 1 or even 2 orders of magnitude. Moreover, molecular docking and molecular dynamics simulation were performed to explore the interactions between the brominated congeners and AHR1-ligand-binding domain (LBD). The molecular mechanics energy (E_MM) between each congener and each individual amino acid (AA) residue in AHR1-LBD was calculated. These E_MM values could finely characterize the final conformation of species-specific AHR1-LBD for each brominated congener. Based on this, mechanism-driven generalized linear models were successfully built using machine learning algorithms and the spline approximation method, and these models could qualitatively predict the complex relationships between AHR1 conformations and DLPs or avian interspecies relative sensitivity to brominated dioxin-like compounds (DLCs). In addition, several AAs conserved among birds were found to potentially interact with species-specific AAs, thereby inducing species-specific interactions between AHR1 and brominated DLCs. The present study provides a novel strategy to facilitate the development of mechanism-driven computational prediction models for supporting safety assessment of DLCs, as well as a basis for the ecotoxicological risk assessment of brominated congeners in birds.

Cytotoxic and Transcriptomic Effects in Avian Hepatocytes Exposed to a Complex Mixture from Air Samples, and Their Relation to the Organic Flame Retardant Signature

Assessing complex environmental mixtures and their effects is challenging. In this study, we evaluate the utility of an avian in vitro screening approach to determine the effects of passive air sampler extracts collected from different global megacities on cytotoxicity and gene expression. Concentrations of a suite of organic flame retardants (OFRs) were quantified in extracts from a total of 19 megacities/major cities in an earlier study, and levels were highly variable across sites. Chicken embryonic hepatocytes were exposed to serial dilutions of extracts from the 19 cities for 24 h. Cell viability results indicate a high level of variability in cytotoxicity, with extracts from Toronto, Canada, having the lowest LC50 value. Partial least squares (PLS) regression analysis was used to estimate LC50 values from OFR concentrations. PLS modeling of OFRs was moderately predictive of LC50 (p-value = 0.0003, r² = 0.66, slope = 0.76, when comparing predicted LC50 to actual values), although only after one outlier city was removed from the analysis. A chicken ToxChip PCR array, comprising 43 target genes, was used to determine effects on gene expression, and similar to results for cell viability, gene expression profiles were highly variable among the megacities. PLS modeling was used to determine if gene expression was related to the OFR profiles of the extracts. Weak relationships to the ToxChip expression profiles could be detected for only three of the 35 OFRs (indicated by regression slopes between 0.6 and 0.5 when comparing predicted to actual OFR concentrations). While this in vitro approach shows promise in terms of evaluating effects of complex mixtures, we also identified several limitations that, if addressed in future studies, might improve its performance.

Using Transcriptomics and Metabolomics to Understand Species Differences in Sensitivity to Chlorpyrifos in Japanese Quail and Double-Crested Cormorant Embryos

Modern 21st-century toxicity testing makes use of omics technologies to address critical questions in toxicology and chemical management. Of interest are questions relating to chemical mechanisms of toxicity, differences in species sensitivity, and translation of molecular effects to observable apical endpoints. Our study addressed these questions by comparing apical outcomes and multiple omics responses in early-life stage exposure studies with Japanese quail (Coturnix japonica) and double-crested cormorant (Phalacrocorax auritus), representing a model and ecological species, respectively. Specifically, we investigated the dose-dependent response of apical outcomes as well as transcriptomics and metabolomics in the liver of each species exposed to chlorpyrifos, a widely used organophosphate pesticide. Our results revealed a clear pattern of dose-dependent disruption of gene expression and metabolic profiles in Japanese quail but not double-crested cormorant at similar chlorpyrifos exposure concentrations. The difference in sensitivity between species was likely due to higher metabolic transformation of chlorpyrifos in Japanese quail compared to double-crested cormorant. The most impacted biological pathways after chlorpyrifos exposure in Japanese quail included hepatic metabolism, oxidative stress, endocrine disruption (steroid and nonsteroid hormones), and metabolic disease (lipid and fatty acid metabolism). Importantly, we show consistent responses across biological scales, suggesting that significant disruption at the level of gene expression and metabolite profiles leads to observable apical responses at the organism level. Our study demonstrates the utility of evaluating effects at multiple biological levels of organization to understand how modern toxicity testing relates to outcomes of regulatory relevance, while also highlighting important, yet poorly understood, species differences in sensitivity to chemical exposure. Environ Toxicol Chem 2021;40:3019-3033. © 2021 SETAC.

Polychlorinated Diphenyl Sulfides: An Emerging Class of Persistent, Bioaccumulative, and Toxic Substances in the Environment

Polychlorinated diphenyl sulfides (PCDPSs) have recently attracted increasing attention due to their potential adverse effects on human and ecosystem health. We present a review regarding their environmental occurrence, persistence, bioaccumulation, toxicity, and biotransformation. The existing literature demonstrates that PCDPSs are ubiquitous in various environmental matrices, are persistent in the environment, and have long-range transport potential. In addition, the high bioaccumulation potential of these emerging pollutants, especially the low chlorinated PCDPS congeners, has been confirmed based on both theoretical calculations and experimental investigations. Moreover, a spectrum of adverse effects, such as acute liver injury, retardation of development, reproductive disorders, and increased mortality have been widely reported in vertebrates. These adverse outcomes were associated with oxidative stress and activation of aryl hydrocarbon receptors. Given these findings, PCDPSs represent candidate persistent, bioaccumulative, and toxic substances and thus deserve further research to fully elucidate their environmental behavior and fate, and evaluate the risks to human and ecosystem health. Environ Toxicol Chem 2021;40:2657-2666. © 2021 SETAC.

Assessing the Toxicity of 17α-Ethinylestradiol in Rainbow Trout Using a 4-Day Transcriptomics Benchmark Dose (BMD) Embryo Assay

There is an urgent demand for more efficient and ethical approaches in ecological risk assessment. Using 17α-ethinylestradiol (EE2) as a model compound, this study established an embryo benchmark dose (BMD) assay for rainbow trout (RBT; Oncorhynchus mykiss) to derive transcriptomic points-of-departure (tPODs) as an alternative to live-animal tests. Embryos were exposed to graded concentrations of EE2 (measured: 0, 1.13, 1.57, 6.22, 16.3, 55.1, and 169 ng/L) from hatch to 4 and up to 60 days post-hatch (dph) to assess molecular and apical responses, respectively. Whole proteome analyses of alevins did not show clear estrogenic effects. In contrast, transcriptomics revealed responses that were in agreement with apical effects, including excessive accumulation of intravascular and hepatic proteinaceous fluid and significant increases in mortality at 55.1 and 169 ng/L EE2 at later time points. Transcriptomic BMD analysis estimated the median of the 20th lowest geneBMD to be 0.18 ng/L, the most sensitive tPOD. Other estimates (0.78, 3.64, and 1.63 ng/L for the 10th percentile geneBMD, first peak geneBMD distribution, and median geneBMD of the most sensitive over-represented pathway, respectively) were within the same order of magnitude as empirically derived apical PODs for EE2 in the literature. This 4-day alternative RBT embryonic assay was effective in deriving tPODs that are protective of chronic effects of EE2.

Extracts from Dated Lake Sediment Cores in the Athabasca Oil Sands Region Alter Ethoxyresorufin-O-deethylase Activity and Gene Expression in Avian Hepatocytes

Increases in oil sands mining operations in the Athabasca oil sands region have resulted in increased concentrations of polycyclic aromatic compounds (PACs) and heavy metals in aquatic systems located near surface mining operations. In the present study, sediment cores were collected from 3 lakes with varying proximity to surface mining operations to determine the differences in PAC concentrations. Sediment cores were separated into 2 sections-current mining (top; 2000-2017) and premining (bottom; pre-1945)-and extracts were prepared for in vitro screening using a well-established chicken embryonic hepatocyte (CEH) assay. Concentrations and composition of PACs varied between sites, with the highest ∑PACs in Saline Lake, 5 km from an active oil sands mine site. The proportion of alkylated PACs was greater than that of parent PACs in the top sediment sections compared with the bottom. Ethoxyresorufin-O-deethylase activity in CEH permitted the ranking of lake sites/core sections based on an aryl hydrocarbon receptor-mediated end point; mean median effect concentration values were lowest for the top cores from Saline Lake and another near-mining operations lake, referred to as WF1. A ToxChip polymerase chain reaction (PCR) array was used to evaluate gene expression changes across 43 target genes associated with numerous toxicological pathways following exposure to top and bottom sediment core extracts. The 2 study sites with the greatest ∑PAC concentrations (Saline Lake and WF1) had the highest gene expression alterations on the ToxChip PCR array (19 [top] and 17 [bottom]/43), compared with a reference site (13 [top] and 7 [bottom]/43). The avian in vitro bioassay was useful for identifying the toxicity of complex PAC extracts associated with variably contaminated sediment cores, supporting its potential use for hotspot identification and complex mixture screening. EnvironToxicol Chem 2021;40:1883-1893. © 2021 SETAC.

Polycyclic Aromatic Hydrocarbons Alter the Hepatic Expression of Genes Involved in Sanderling (Calidris alba) Pre-migratory Fueling

Polycyclic aromatic hydrocarbons (PAHs) impaired pre-migratory fueling in 49 orally dosed Sanderling (Calidris alba). In the present study, 8 genes related to fat deposition and PAH exposure were measured in liver subsamples from these same shorebirds. At the highest dose (1260 µg total PAH [tPAH]/kg body wt/day), PAH exposure decreased liver basic fatty acid binding protein 1 (Lbfabp) and hepatic lipase (Lipc) expression. The present study reveals candidate molecular-level pathways for observed avian pre-migratory refueling impairment. Environ Toxicol Chem 2021;40:1983-1991. © 2021 SETAC.

In Vitro Screening of 21 Bisphenol A Replacement Alternatives: Compared with Bisphenol A, the Majority of Alternatives Are More Cytotoxic and Dysregulate More Genes in Avian Hepatocytes

An avian in vitro screening approach was used to determine the effects of 21 bisphenol A (BPA) alternatives. Cytotoxicity and dysregulation of genes associated with estrogen response and other toxicologically relevant pathways evoked by these alternatives were compared with BPA. Most of the BPA alternatives (15/21) were equally or more cytotoxic than BPA in chicken embryonic hepatocytes; variability in cell viability was associated with chemical structure and the log octanol-water partition coefficient (logP) values. A negative linear relationship (r ²  = 0.745; p = 0.49^-07 ; n = 18) was observed between logP and the log median lethal concentration (logLC50) values. The least cytotoxic BPA alternatives elicited the greatest gene dysregulation and, overall, most of the alternatives altered more genes than BPA (measured with a custom polymerase chain reaction array). This overall approach shows promise for use as a screen for hazard-based prioritization of BPA replacement alternatives and to ideally identify those that may be less harmful and/or require additional toxicity testing. Environ Toxicol Chem 2021;40:2026-2033. © 2021 Her Majesty the Queen in Right of Canada Environmental Toxicology and Chemistry © 2021 SETAC. Reproduced with the permission of the Minister of Environment and Climate Change Canada.

Cross-Model Comparison of Transcriptomic Dose-Response of Short-Chain Chlorinated Paraffins

Short-chain chlorinated paraffins (SCCPs) have attracted attention because of their toxicological potential in humans and wildlife at environmentally relevant doses. However, limited information is available regarding mechanistic differences across species in terms of the biological pathways that are impacted by SCCP exposure. Here, a concentration-dependent reduced human transcriptome (RHT) approach was conducted to evaluate 15 SCCPs in HepG2 cells and compared with our previous results using a reduced zebrafish transcriptome (RZT) approach in zebrafish embryos (ZFEs). Generally, SCCPs induced a broader suite of biological pathways in ZFEs than HepG2 cells, and all of the 15 SCCPs were more potent in HepG2 cells compared to ZFEs. Despite these general differences, the transcriptional potency of SCCPs in both model systems showed a significant linear relationship (p = 0.0017, r² = 0.57), and the average ratios of transcriptional potency for each SCCP in RZT to that in RHT were ∼100,000. C₁₀H₁₄Cl₈ was the most potent SCCP, while C₁₀H₁₇Cl₅ was the least potent in both ZFEs and HepG2 cells. An adverse outcome pathway network-based analysis demonstrated model-specific responses, such as xenobiotic metabolism that may be mediated by different nuclear receptor-mediated pathways between HepG2 cells (e.g., CAR and AhR activation) and ZFEs (e.g., PXR activation). Moreover, induced transcriptional changes in ZFEs associated with pathways and molecular initiating events (e.g., activation of nicotinic acetylcholine receptor) suggest that SCCPs may disrupt neural development processes. The cross-model comparison of concentration-dependent transcriptomics represents a promising approach to assess and prioritize SCCPs.

ToxChip PCR Arrays for Two Arctic-Breeding Seabirds: Applications for Regional Environmental Assessments

Increasing pollution in the Arctic poses challenges in terms of geographical and ecological monitoring. The Baffin Bay-Davis Strait (BBDS) region in the Canadian Arctic Archipelago is of particular concern due to the potential for increased shipping traffic and oil exploration. However, data on background contaminants associated with oil exploration/spills/natural seeps (e.g., polycyclic aromatic compounds [PAC]) and measures of potential effects for Arctic birds are limited. We developed a toxicogenomics approach to investigate the background gene expression profiles for two Arctic-breeding seabirds, the thick-billed murre (Uria lomvia) and the black guillemot (Cepphus grylle), which will aid effects-based monitoring efforts. Chemical burdens (53 PACs and 5 trace elements) and transcriptomic profiles (31 genes using a ToxChip PCR array) were examined in liver tissues (n = 30) of each species collected from the Qaqulluit and Akpait National Wildlife Areas in the BBDS region. While chemical and transcriptomic profiles demonstrated low variability across individuals for each species, gene expression signatures were able to distinguish guillemots collected from two distinct colonies. This toxicogenomics approach provides benchmark data for two Arctic seabirds and is promising for future monitoring efforts and strategic environmental assessments in this sensitive ecosystem and areas elsewhere in the circumpolar Arctic that are undergoing change.

Effects of two Bisphenol A replacement compounds, 1,7-bis (4-hydroxyphenylthio)-3,5-dioxaheptane and Bisphenol AF, on development and mRNA expression in chicken embryos

Concerns about the estrogenic properties of Bisphenol A (BPA) have led to increased efforts to find BPA replacements. 1,7-bis(4-Hydroxyphenylthio)-3,5-dioxaheptane (DD-70) and 4,4'-(hexafluoroisopropylidene) diphenol (bisphenol AF, BPAF) are two potential chemical substitutes for BPA; however, toxicity data for these chemicals in avian species are limited. To determine effects on avian embryonic viability, development, and hepatic mRNA expression at two distinct developmental periods (mid-incubation [day 11] and term [day 20]), two egg injection studies were performed. Test chemicals were injected into the air cell of unincubated, fertilized chicken eggs at concentrations ranging from 0-88.2 µg/g for DD-70 and 0-114 µg/g egg for BPAF. Embryonic concentrations of DD-70 and BPAF decreased at mid-incubation and term compared to injected concentrations suggesting embryonic metabolism. Exposure to DD-70 (40.9 and 88.2 µg/g) and BPAF (114 µg/g) significantly decreased embryonic viability at mid-incubation. Exposure to DD-70 (88.2 µg/g) decreased embryo mass and increased gallbladder mass, while 114 µg/g BPAF resulted in increased gallbladder mass in term embryos. Expression of hepatic genes related to xenobiotic metabolism, lipid homeostasis, and response to estrogen were altered at both developmental stages. Given the importance of identifying suitable BPA replacements, the present study provides novel, whole animal avian toxicological data for two replacement compounds, DD-70 and BPAF. DATA AVAILABILITY: Data, associated metadata, and calculation tools are available from the corresponding author (doug.crump@canada.ca).

Concentration- and time-dependent induction of Cyp1a and DNA damage response by benzo(a)pyrene in LMH three-dimensional spheroids

In vitro liver toxicity tests performed using cell lines cultured as two-dimensional (2D) monolayer have limited CYP450 activity and may be inadequate for screening chemicals that require activation to exert toxicity. Metabolic competence is greatly improved using three-dimensional (3D) cell culture. In this study, Cyp1a induction, and subsequent DNA damage response induced by benzo(a)pyrene (BaP) were compared in 2D monolayer cells and 3D spheroids of the chicken hepatic cell line, LMH. Cells were exposed to BaP (0.1-100 μM) for different durations: 8, 24, 35, or 48 hr. Cyp1a activity, mRNA expression of Cyp1a and DNA damage response (DDR) genes, and phosphorylation of H2AX (γH2AX) were determined using the EROD assay, a customized PCR array, and flow cytometry, respectively. EROD activity was induced at 8 hr and achieved maximal induction at 24 hr in spheroids; earlier time points than for monolayer cells. In spheroids, BaP exposure resulted in a concentration-dependent increase in Cyp1a4 mRNA expression at 8 hr followed by upregulation of DDR genes at 24 hr, whereas Cyp1a4 mRNA induction was only observed at 48 hr in monolayer cells. Cyp1a5 mRNA was induced at 8 hr in monolayer cells but maximum induction was greater in spheroids. An increase in γH2AX was observed at 24 hr in spheroids; this endpoint was not evaluated in monolayer cells. These results suggest that BaP metabolism precedes the DNA damage response and occurs earlier in 3D spheroids. This study demonstrates that LMH 3D spheroids could be a suitable metabolically-competent in vitro model to measure genotoxicity of chemicals that require metabolic activation by Cyp1a.

Toxicity Screening of Bisphenol A Replacement Compounds: Cytotoxicity and mRNA Expression in Primary Hepatocytes of Chicken and Double-Crested Cormorant

A market for bisphenol A (BPA) replacement compounds has emerged as a result of restrictions on the use of BPA. Some of these compounds have been detected in the environment; however, little is known about their toxicological properties. In the present study, an avian in vitro toxicogenomic approach was used to compare the effects of 5 BPA alternatives. Cell viability and mRNA expression were compared in primary embryonic hepatocytes of chicken (CEH) and double-crested cormorant (DCEH) exposed to 4,4'-(propane-2,2-diyl) diphenol (BPA), bis (4-hydroxyphenyl) methane (BPF), bis (3-allyl-4-hydroxyphenyl) sulfone (TGSH), 7-bis (4-hydroxyphenylthio)-3,5-dioxaheptane (DD-70), 2,2-bis (4-hydroxyphenyl) hexafluoropropane (BPAF), and 4-hydroxyphenyl 4-isoprooxyphenylsulfone (BPSIP). Changes in gene expression were determined using 2 polymerase chain reaction (PCR) arrays: 1) species-specific ToxChips that contain genes representing toxicologically relevant pathways, and 2) chicken-specific AestroChip that measures estrogen responsive genes. In CEH and DCEH, BPA alternatives TGSH, DD-70, and BPAF were most cytotoxic. Some of the replacement compounds changed the expression of genes related to xenobiotic metabolism, bile acid, and cholesterol regulation. The rank order based on the number of genes altered on the chicken ToxChip array was TGSH > DD-70 > BPAF = BPF > 17β estradiol (E2) > BPSIP > BPA. On the cormorant ToxChip array, BPSIP altered the greatest number of genes. Based on the chicken AestroChip data, BPSIP and BPF were slightly estrogenic. These results suggest that the replacement compounds have comparable or even greater toxicity than BPA and act via different mechanisms. Environ Toxicol Chem 2021;40:1368-1378. © 2021 Her Majesty the Queen in Right of Canada. Reproduced with the permission of the Minister of Environment and Climate Change Canada.

Development of a Comprehensive Toxicity Pathway Model for 17α-Ethinylestradiol in Early Life Stage Fathead Minnows (Pimephales promelas)

There is increasing pressure to develop alternative ecotoxicological risk assessment approaches that do not rely on expensive, time-consuming, and ethically questionable live animal testing. This study aimed to develop a comprehensive early life stage toxicity pathway model for the exposure of fish to estrogenic chemicals that is rooted in mechanistic toxicology. Embryo-larval fathead minnows (FHM; Pimephales promelas) were exposed to graded concentrations of 17α-ethinylestradiol (water control, 0.01% DMSO, 4, 20, and 100 ng/L) for 32 days. Fish were assessed for transcriptomic and proteomic responses at 4 days post-hatch (dph), and for histological and apical end points at 28 dph. Molecular analyses revealed core responses that were indicative of observed apical outcomes, including biological processes resulting in overproduction of vitellogenin and impairment of visual development. Histological observations indicated accumulation of proteinaceous fluid in liver and kidney tissues, energy depletion, and delayed or suppressed gonad development. Additionally, fish in the 100 ng/L treatment group were smaller than controls. Integration of omics data improved the interpretation of perturbations in early life stage FHM, providing evidence of conservation of toxicity pathways across levels of biological organization. Overall, the mechanism-based embryo-larval FHM model showed promise as a replacement for standard adult live animal tests.

Ultrafast functional profiling of RNA-seq data for nonmodel organisms

Computational time and cost remain a major bottleneck for RNA-seq data analysis of nonmodel organisms without reference genomes. To address this challenge, we have developed Seq2Fun, a novel, all-in-one, ultrafast tool to directly perform functional quantification of RNA-seq reads without transcriptome de novo assembly. The pipeline starts with raw read quality control: sequencing error correction, removing poly(A) tails, and joining overlapped paired-end reads. It then conducts a DNA-to-protein search by translating each read into all possible amino acid fragments and subsequently identifies possible homologous sequences in a well-curated protein database. Finally, the pipeline generates several informative outputs including gene abundance tables, pathway and species hit tables, an HTML report to visualize the results, and an output of clean reads annotated with mapped genes ready for downstream analysis. Seq2Fun does not have any intermediate steps of file writing and loading, making I/O very efficient. Seq2Fun is written in C++ and can run on a personal computer with a limited number of CPUs and memory. It can process >2,000,000 reads/min and is >120 times faster than conventional workflows based on de novo assembly, while maintaining high accuracy in our various test data sets.

Effects on Apical Outcomes of Regulatory Relevance of Early-Life Stage Exposure of Double-Crested Cormorant Embryos to 4 Environmental Chemicals

Environmental risk assessment is often challenged by a lack of toxicity data for ecological species. The overall goal of the present study was to employ an avian early-life stage toxicity test to determine the effects of 4 chemicals (benzo[a]pyrene [BaP], chlorpyrifos, fluoxetine hydrochloride [FLX], and ethinyl estradiol [EE2]) on an ecologically relevant avian species, the double-crested cormorant (Phalacrocorax auritus), and to compare our results with those we previously reported for a laboratory model species, Japanese quail. Chemicals were dissolved in dimethyl sulfoxide and administered via air cell injection to fertilized, unincubated double-crested cormorant eggs at 3 nominal concentrations, the highest selected to approximate the 20% lethal dose. Of the 4 chemicals, only chlorpyrifos and FLX were detected in liver tissue of embryos at midincubation (day 14) and termination (day 26; 1-2 d prior to hatch); EE2 and BaP were not detectable, suggesting embryonic clearance/metabolism. No apical effects were observed in double-crested cormorant embryos up to the highest concentrations of chlorpyrifos (no-observed-effect level [NOEL] = 25 µg/g) or FLX (NOEL = 18 µg/g). Exposure to EE2 reduced embryonic viability and increased deformities at a concentration of 2.3 µg/g (NOEL = 0.18 µg/g), and BaP decreased embryonic viability (median lethal dose = 0.015 µg/g; NOEL = 0.0027 µg/g). Compared with Japanese quail, double-crested cormorant were more sensitive with regard to embryolethality and deformities for EE2 and embryolethality for BaP, whereas they were less sensitive to embryonic deformities associated with chlorpyrifos exposure. These data reinforce the idea that standardized toxicity tests using a laboratory model species may not always be protective of wild birds, and thus they stress the importance of developing such alternative testing strategies (e.g., the EcoToxChip Project) for ecologically relevant species to augment risk assessment efforts. Environ Toxicol Chem 2021;40:390-401. © 2020 SETAC.

Polycyclic aromatic compounds (PACs) and trace elements in four marine bird species from northern Canada in a region of natural marine oil and gas seeps

There is a growing understanding of how oil pollution can affect aquatic ecosystems, including physical and chemical effects. One of the biggest challenges with detecting the effects of oil-related contaminants on biota from resource development is understanding the background levels and potential effects of the exposure of biota to contaminants from various natural and anthropogenic sources prior to large scale oil and gas operations. Seabirds are effective indicators of pollution, and can be useful for tracking oil-related contaminants in the marine environment. We sampled four seabird species (black guillemot, Cepphus grylle; thick-billed murre, Uria lomvia; black-legged kittiwake, Rissa tridactyla; and northern fulmar, Fulmarus glacialis) in the Baffin Bay-Davis Strait region of the Northwest Atlantic and Arctic oceans, an area where natural oil and gas seeps are present but lacking any large-scale oil and gas projects. We found detectable levels of PACs and several trace elements in all species examined. Alkylated PAC levels were higher than parent compounds in all four seabird species examined, with fulmars and murres having the highest levels detected; mean hepatic concentrations of ∑₁₆PAC were 99.05, 46.42, 12.78 and 9.57 ng/g lw, respectively, for guillemots, murres, fulmars and kittiwakes. Overall, PAC concentrations in the seabird species examined were similar to PAC concentrations measured in other bird species in regions with more industrialization. These findings provide data which can be used to assess the current oil-related contaminant exposure of biota in the region. As well, they provide background levels for the region at a time when shipping activity is relatively low, which can used for future comparisons following expected increases in shipping and oil and gas activities in the region.

Drivers of and Obstacles to the Adoption of Toxicogenomics for Chemical Risk Assessment: Insights from Social Science Perspectives

BACKGROUND

Some 20 y ago, scientific and regulatory communities identified the potential of omics sciences (genomics, transcriptomics, proteomics, metabolomics) to improve chemical risk assessment through development of toxicogenomics. Recognizing that regulators adopt new scientific methods cautiously given accountability to diverse stakeholders, the scope and pace of adoption of toxicogenomics tools and data have nonetheless not met the ambitious, early expectations of omics proponents.

OBJECTIVE

Our objective was, therefore, to inventory, investigate, and derive insights into drivers of and obstacles to adoption of toxicogenomics in chemical risk assessment. By invoking established social science frameworks conceptualizing innovation adoption, we also aimed to develop recommendations for proponents of toxicogenomics and other new approach methodologies (NAMs).

METHODS

We report findings from an analysis of 56 scientific and regulatory publications from 1998 through 2017 that address the adoption of toxicogenomics for chemical risk assessment. From this purposeful sample of toxicogenomics discourse, we identified major categories of drivers of and obstacles to adoption of toxicogenomics tools and data sets. We then mapped these categories onto social science frameworks for conceptualizing innovation adoption to generate actionable insights for proponents of toxicogenomics.

DISCUSSION

We identify the most salient drivers and obstacles. From 1998 through 2017, adoption of toxicogenomics was understood to be helped by drivers such as those we labeled Superior scientific understanding, New applications, and Reduced cost & increased efficiency but hindered by obstacles such as those we labeled Insufficient validation, Complexity of interpretation, and Lack of standardization. Leveraging social science frameworks, we find that arguments for adoption that draw on the most salient drivers, which emphasize superior and novel functionality of omics as rationales, overlook potential adopters' key concerns: simplicity of use and compatibility with existing practices. We also identify two perspectives-innovation-centric and adopter-centric-on omics adoption and explain how overreliance on the former may be undermining efforts to promote toxicogenomics. https://doi.org/10.1289/EHP6500.

Evaluation of the Aryl Hydrocarbon Receptor Response in LMH 3D Spheroids

In the present study, we investigated whether the immortalized chicken hepatocellular carcinoma cell line, leghorn male hepatoma (LMH), had a comparable aryl hydrocarbon receptor (AhR) response to primary chicken embryonic hepatocytes (CEHs) when used in a well-established assay for chemical screening and prioritization. The LMH cells were grown as 2-dimensional (2D) confluent cells and 3D spheroids to determine the optimal cell culture states for chemical screening. Cytochrome P450 1A4 and 1A5 (CYP1A) activity and gene expression were compared between CEHs and LMH cells grown in 2 culture states following exposure to the dioxin-like compound 3,3',4,4',5-pentachlorobiphenyl (PCB-126). The CYP1A activity was measured using the ethoxyresorufin-O-deethylase (EROD) assay, and changes in mRNA expression associated with the AhR pathway were determined using a custom-designed polymerase chain reaction array. Among LMH cell culture states (i.e., 2D vs 3D), EROD induction was observed only in 3D LMH spheroids. Similarly, 3D spheroids had the greatest number of changes in AhR-related genes compared with confluent cells. Overall, these results suggest that LMH cells grown as 3D spheroids have a metabolic and gene expression profile that is comparable to that of CEH, and may represent a suitable animal-free alternative for in vitro screening of chemicals. Environ Toxicol Chem 2020;39:1693-1701. © 2020 SETAC.

Toxicogenomic Assessment of Complex Chemical Signatures in Double-Crested Cormorant Embryos from Variably Contaminated Great Lakes Sites

Using omics approaches to monitor complex environmental mixtures is challenging. Previously, we evaluated in vitro transcriptomic effects of complex organic extracts derived from avian eggs. However, there is a lack of studies using wild species that are naturally exposed to contaminant mixtures. Here, we examined polychlorinated biphenyl (PCB) and polybrominated diphenyl ether (PBDE) residues and gene expression in embryonic liver tissue of double-crested cormorants (Phalacrocorax auritus) collected from six variably contaminated colonies. Colonies near industrialized areas were distinguished from less contaminated sites based on their PCB and PBDE concentrations. The most variably expressed genes between sites were involved in pathways including, xenobiotic metabolism (e.g., Cyp1a4), lipid/bile acid homeostasis (e.g., Lbfabp), and oxidative stress (e.g., Mt4). Hierarchical clustering, based on relative gene expression, revealed a grouping pattern similar to chemical residue concentrations. Further, partial least squares regression analysis was used to estimate chemical concentrations from transcriptomics data. PCB 155 and BDE 47 showed the highest slopes (0.77 and 0.69, respectively) fitted by linear regression of measured and estimated chemical concentrations. The application of transcriptomics to a wild avian species, naturally exposed to complex chemical mixtures and other stressors, represents a promising means to distinguish and prioritize variably contaminated sites.

EcoToxModules: Custom Gene Sets to Organize and Analyze Toxicogenomics Data from Ecological Species

Traditional results from toxicogenomics studies are complex lists of significantly impacted genes or gene sets, which are challenging to synthesize down to actionable results with a clear interpretation. Here, we defined two sets of 21 custom gene sets, called the functional and statistical EcoToxModules, in fathead minnow (Pimephales promelas) to (1) re-cast predefined molecular pathways into a toxicological framework and (2) provide a data-driven, unsupervised grouping of genes impacted by exposure to environmental contaminants. The functional EcoToxModules were identified by re-organizing KEGG pathways into biological processes that are more relevant to ecotoxicology based on the input from expert scientists and regulators. The statistical EcoToxModules were identified using co-expression analysis of publicly available microarray data (n = 303 profiles) measured in livers of fathead minnows after exposure to 38 different conditions. Potential applications of the EcoToxModules were demonstrated with two case studies that represent exposure to a pure chemical and to environmental wastewater samples. In comparisons to differential expression and gene set analysis, we found that EcoToxModule responses were consistent with these traditional results. Additionally, they were easier to visualize and quantitatively compare across different conditions, which facilitated drawing conclusions about the relative toxicity of the exposures within each case study.

Computational evaluation of interactions between organophosphate esters and nuclear hormone receptors

Organophosphate esters (OPEs) have gained considerable interest from many environmental chemists and toxicologists due to their frequent detection in the environment and potential adverse effects on health. Nuclear hormone receptors (NHRs) were found to mediate many of their adverse effects. However, our knowledge regarding the direct binding and interaction between OPEs and NHRs is limited. In this study, Endocrine Disruptome, an online computational tool based on the technique of inverse docking, was used to calculate the binding affinity score of 25 individual OPEs with 12 different human NHRs. Results showed that 20% of potential binding interactions between the OPEs and NHRs had medium-to-high probabilities. The accuracy, sensitivity and specificity of the predictions were 78.8, 60.0 and 80.9%, respectively. OPEs with a benzene ring were more active than those without, among which, tri-o-tolyl phosphate and tri-m-tolyl phosphate displayed the highest activities, suggesting that they might pose the greatest potential risks for interference with endocrine functions. In addition, the antagonistic conformations of androgen receptor and estrogen receptor β were found to be the two most vulnerable NHR conformations. Our findings can further the understanding about the health risk(s) of OPEs.

Factors Affecting the Perception of New Approach Methodologies (NAMs) in the Ecotoxicology Community

Given current legislative mandates to assess the safety of thousands of chemicals and the slow pace at which conventional testing proceeds, there is a need to accelerate chemical risk assessment. Governments and businesses are increasingly interested in new approach methodologies (NAMs) that promise to reduce costs and delays. We explore 5 sociological factors within the ecotoxicology community that can influence the perception of NAMs: 1) professional profile (educational cohort, employer), 2) internal science communication within professional forums, 3) concern for "error cost," 4) collaboration across stakeholders, and 5) fundamental beliefs regarding toxicology. We conducted an online survey (n = 171; 2018) asking participants about their experiences and perspectives at events of the Society of Environmental Toxicology and Chemistry (SETAC) to assess 1) how NAMs are discussed compared to conventional testing and 2) how respondents perceive their viability. We developed ordered logistic regression (OLR) models to understand the influence of exploratory variables (cohort, core views on toxicology, frequency of collaboration) on respondents' evaluation of the viability of different NAMs. Our results showed that 1) NAMs were more likely than conventional methods to be challenged in forum discussions, which may be fueled by concerns for error costs in regulatory decision making; 2) perceptions of the viability of NAMs tended to follow a "pattern of familiarity," whereby respondents that were more knowledgeable about a test method tended to find it more viable; 3) respondents who agreed with the Paracelsus maxim had a greater likelihood of finding conventional testing viable; and 4) the more a respondent reported collaborating with industry on alternative testing strategies, the more likely she or he was to report that NAMs were less viable. These results suggest that there are professional and organizational barriers to greater acceptance of NAMs that can be addressed through a social learning process within the professional community. Integr Environ Assess Manag 2020;16:269-281. © 2020 SETAC.

An Early-Life Stage Alternative Testing Strategy for Assessing the Impacts of Environmental Chemicals in Birds

Early-life stage (ELS) toxicity tests are recognized as an advancement over current testing methodologies in terms of cost, animal use, and biological relevance. However, standardized ELS tests are not presently available for some vertebrate taxa, including birds. The present study describes a Japanese quail (Coturnix japonica) ELS test that is a promising candidate for standardization and applies it to test 8 environmental chemicals (ethinylestradiol, benzo[a]pyrene, chlorpyrifos, fluoxetine, lead(II)nitrate, trenbolone, seleno-L-methionine, hexabromocyclododecane). Individual chemicals were injected into the air cell of unincubated Japanese quail eggs at 3 concentrations, all predicted to cause ≤20% mortality. Survival to embryonic day 16 was consistently high (>90%) among the vehicle-injected controls. All chemicals, except ethinylestradiol, were detected in liver tissue, most at concentrations suggestive of embryonic clearance. Adverse effects were observed for 5 of the 8 chemicals; chlorpyrifos (41.1 µg/g) significantly increased developmental abnormalities and decreased embryo and gallbladder mass. Ethinylestradiol (54.2 µg/g) and hexabromocyclododecane (0.02 µg/g) decreased embryo mass and tarsus length, respectively. Benzo[a]pyrene (0.83 µg/g) and fluoxetine hydrochloride (32.7 µg/g) exceeded the 20% mortality cutoff. No effects were observed following lead(II)nitrate, seleno-L-methionine, or trenbolone exposure up to 10.7, 0.07, and 4.4 µg/g, respectively. Overall, our ELS approach was time- and cost-effective, caused minimal mortality in controls, effectively delivered diverse chemicals to the embryo, and permitted identification of apical outcomes, all of which provide support toward standardization. Environ Toxicol Chem 2019;39:141-154. © 2019 SETAC.

Persistent, bioaccumulative, and toxic properties of liquid crystal monomers and their detection in indoor residential dust

Liquid crystal monomers (LCMs) are used widely in liquid crystal displays (LCDs), which are dramatically changing the world due to the provision of convenient communication. However, there are essentially no published reports on the fate and/or effects of LCMs in the environment. Of 362 currently produced LCMs, 87 were identified as persistent and bioaccumulative (P&B) chemicals, which indicated that these chemicals would exhibit resistance to degradation and exhibit mobility after entering the environment. Following exposure to mixtures of LCM collected from 6 LCD devices, significant modulation of 5 genes, CYP1A4, PDK4, FGF19, LBFABP, and THRSP, was observed in vitro. Modulation of expressions of mRNAs coding for these genes has frequently been reported for toxic (T) persistent organic pollutants (POPs). In LCM mixtures, 33 individual LCMs were identified by use of mass spectrometry and screened for in 53 samples of dust from indoor environments. LCMs were detectable in 47% of analyzed samples, and 17 of the 33 LCMs were detectable in at least 1 sample of dust. Based on chemical properties, including P&B&T of LCMs and their ubiquitous detection in dust samples, the initial screening information suggests a need for studies to determine status and trends in concentrations of LCMs in various environmental matrices as well as tissues of humans and wildlife. There is also a need for more comprehensive in vivo studies to determine toxic effects and potencies of LCMs during chronic, sublethal exposures.

A comparative study of 3 alternative avian toxicity testing methods: Effects on hepatic gene expression in the chicken embryo

There is growing interest in developing alternative methods to screen and prioritize chemical hazards, although few studies have compared responses across different methods. The objective of the present study was to compare 3 alternative liver methods derived from white Leghorn chicken (Gallus gallus domesticus): primary hepatocyte culture, liver slices, and liver from in ovo injected embryos. We examined hepatic gene expression changes after exposure to 3 chemicals (17β-trenbolone [17βT], 17β-estradiol [E2], and 2,3,7,8-tetrachlorodibenzo-p-dioxin [TCDD]) using a custom quantitative polymerase chain reaction (qPCR) array with 7 genes (vitellogenin [VTG], apolipoprotein [Apo], cytochrome P450 1A4 [CYP1A4], liver basic fatty acid binding protein [LBFABP], 3β hydroxysteroid dehydrogenase [HSD3β1], stearoyl coenzyme A desaturase [SCD], and estrogen sulfotransferase [SULT1E1]). Gene expression across the 3 methods was examined using hierarchical clustering. Up-regulation of CYP1A4 in response to TCDD was consistent across all methods, and the magnitude was higher in hepatocytes (>150-fold) compared with slices (>31-fold) and in ovo liver (>27-fold). In hepatocytes, SCD and VTG up-regulation in response to 17βT and E2 was >4-fold and 16-fold, respectively. The rank order of cases with significant changes in gene expression among the 3 methods was: hepatocytes (22) > in ovo liver (11) > liver slices (6). Hierarchical clustering grouped liver slices and in ovo liver as more similar, whereas hepatocytes were grouped separately from in ovo liver. More introspective comparisons are needed to understand how and why alternative methods differ and to aid in their integration into toxicity testing. Environ Toxicol Chem 2019;38:2546-2555. © 2019 SETAC.

T1000: a reduced gene set prioritized for toxicogenomic studies

There is growing interest within regulatory agencies and toxicological research communities to develop, test, and apply new approaches, such as toxicogenomics, to more efficiently evaluate chemical hazards. Given the complexity of analyzing thousands of genes simultaneously, there is a need to identify reduced gene sets. Though several gene sets have been defined for toxicological applications, few of these were purposefully derived using toxicogenomics data. Here, we developed and applied a systematic approach to identify 1,000 genes (called Toxicogenomics-1000 or T1000) highly responsive to chemical exposures. First, a co-expression network of 11,210 genes was built by leveraging microarray data from the Open TG-GATEs program. This network was then re-weighted based on prior knowledge of their biological (KEGG, MSigDB) and toxicological (CTD) relevance. Finally, weighted correlation network analysis was applied to identify 258 gene clusters. T1000 was defined by selecting genes from each cluster that were most associated with outcome measures. For model evaluation, we compared the performance of T1000 to that of other gene sets (L1000, S1500, Genes selected by Limma, and random set) using two external datasets based on the rat model. Additionally, a smaller (T384) and a larger version (T1500) of T1000 were used for dose-response modeling to test the effect of gene set size. Our findings demonstrated that the T1000 gene set is predictive of apical outcomes across a range of conditions (e.g., in vitro and in vivo, dose-response, multiple species, tissues, and chemicals), and generally performs as well, or better than other gene sets available.