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.

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.

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

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.

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.

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.

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 (EMM) between each congener and each individual amino acid (AA) residue in AHR1-LBD was calculated. These EMM 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.

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, r2 = 0.57), and the average ratios of transcriptional potency for each SCCP in RZT to that in RHT were ∼100,000. C10H14Cl8 was the most potent SCCP, while C10H17Cl5 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.

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.

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.

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

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.

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.

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.

Transcriptomic points-of-departure from short-term exposure studies are protective of chronic effects for fish exposed to estrogenic chemicals.

Resource limitations often require risk assessors to extrapolate chronic toxicity from acute tests using assessment factors. Transcriptomic dose-response analysis following short-term exposures may provide a more reliable and biologically-based alternative for estimating chronic toxicity. Here, we demonstrate that transcriptomic dose-response analysis in fish following short-term exposure to endocrine disrupting chemicals (EDCs) provides estimates of chronic toxicity that may be used as protective points-of-departure (POD) for risk assessment. The benchmark dose (BMD) method was used on publicly available datasets (n = 5) to determine transcriptomic PODs in fish exposed to three EDCs (bisphenol A, ethinylestradiol, and diethylstilbestrol). To test for potential bias related to data processing, our analysis compared the effect of different normalization, filtering, and BMD-grouping methods on the transcriptomic PODs. The resulting PODs were then compared to the empirically-derived chronic LOEC of each substance. Normalization and filtering methods had limited impact on the final PODs. However, we found that PODs derived from ontology- or pathway-based gene grouping methods were highly variable, whereas PODs from grouping methods that focused on the most responsive genes were more stable and provided POD estimates that were most similar to the chronic LOEC. Overall, 72% of transcriptomic PODs were within 1 order of magnitude of the chronic LOEC, regardless of data analysis method. When our recommended analysis approach was applied, the concordance improved to 100%. These results suggest that toxicogenomic dose-response analysis has the potential to be a protective decision-support tool for compounds with chronic toxicity, such as EDCs.

Extracts of Passive Samplers Deployed in Variably Contaminated Wetlands in the Athabasca Oil Sands Region Elicit Biochemical and Transcriptomic Effects in Avian Hepatocytes.

Recent contaminant monitoring in boreal wetlands situated in Alberta's Athabasca oil sands region revealed increased concentrations of polycyclic aromatic compounds (PACs) in passive sampling devices deployed in wetlands close to bitumen surface mining operations. In this study, graded concentrations of semipermeable membrane device (SPMD) extracts, collected from 4 wetlands with variable burdens of PACs, were administered to chicken and double-crested cormorant (DCCO) embryonic hepatocytes to determine effects on 7-ethoxyresorufin-O-deethylase (EROD) activity and mRNA expression. Concentrations and composition of PACs detected in SPMDs varied among sites, and the proportion of alkyl PACs was greater than parent compounds at all sites. ΣPACs was the highest in SPMDs deployed within 10 km of mining activity (near-site wetland [5930 ng SPMD-1]) compared to those ∼50 km south (far-site wetland [689 ng SPMD-1]). Measures of EROD activity and Cyp1a4 mRNA expression allowed the ranking of wetland sites based on aryl hydrocarbon receptor-mediated end points; EROD activity and Cyp1a4 mRNA induction were the highest at the near-site wetland. ToxChip PCR arrays (one chicken and one DCCO) provided a more exhaustive transcriptomic evaluation across multiple toxicological pathways following exposure to the SPMD extracts. Study sites with the greatest PAC concentrations had the most genes altered on the chicken ToxChip (12-15/43 genes). Exposure of avian hepatocytes to SPMD extracts from variably contaminated wetlands highlighted traditional PAC-related toxicity pathways as well as other novel mechanisms of action. A novel combination of passive sampling techniques and high-throughput toxicity evaluation techniques shows promise in terms of identifying hotspots of chemical concern in the natural environment.

Organophosphate Ester, 2-Ethylhexyl Diphenyl Phosphate (EHDPP), Elicits Cytotoxic and Transcriptomic Effects in Chicken Embryonic Hepatocytes and Its Biotransformation Profile Compared to Humans.

The effects of 2-ethylhexyl diphenyl phosphate (EHDPP) on cytotoxicity and mRNA expression, as well as its metabolism, were investigated using a chicken embryonic hepatocyte (CEH) assay. After incubation for 36 h, the lethal concentration 50 (LC50) was 50 ± 11 µM, suggesting that EHDPP is one of a small cohort of highly toxic organophosphate esters (OPEs). By use of a ToxChip polymerase chain reaction (PCR) array, we report modulation of 6, 11, or 16/43 genes in CEH following exposure to 0.1, 1, or 10 µM EHDPP, respectively. The altered genes were from all nine biological pathways represented on the ToxChip including bile acids/cholesterol regulation, glucose metabolism, lipid homeostasis, and the thyroid hormone pathway. After incubation for 36 h, 92.5% of EHDPP was transformed, and one of its presumed metabolites, diphenyl phosphate (DPHP), only accounted for 12% of the original EHDPP concentration. Further screening by use of high-resolution mass spectrometry revealed a novel EHDPP metabolite, hydroxylated 2-ethylhexyl monophenyl phosphate (OH-EHMPP), which was also detected in a human blood pool. Additional EHDPP metabolites detected in the human blood pool included EHMPP and DPHP. Overall, this study provided novel information regarding the toxicity of EHDPP and identified a potential EHDPP metabolite, OH-EHMPP, in both avian species and humans.

Down-Regulation of hspb9 and hspb11 Contributes to Wavy Notochord in Zebrafish Embryos Following Exposure to Polychlorinated Diphenylsulfides.

It is hypothesized that key genes, other than ahr2, are present and associated with the development of a unique type of notochord malformation known as wavy notochord in early life stages of zebrafish following exposure to polychlorinated diphenylsulfides (PCDPSs). To investigate the potential mechanism(s), time-dependent developmental morphologies of zebrafish embryos following exposure to 2500 nM 2,4,4',5-tetra-CDPS, 2,2',4-tri-CDPS or 4,4'-di-CDPS were observed to determine the developmental time point when notochord twists began to occur (i.e., 21 h-postfertilization (hpf)). Simultaneously, morphometric measurements suggested that PCDPS exposure did not affect notochord growth at 21 or 120 hpf; however, elongation of the body axis was significantly inhibited at 120 hpf. Transcriptome analysis revealed that the retardation of body growth was potentially related with dysregulation of transcripts predominantly associated with the insulin-associated Irs-Akt-FoxO cascade. Moreover, knockdown and gain-of-function experiments in vivo on codifferentially expressed genes demonstrated that reduced expression of hspb9 and hspb11 contributed to the occurrence of wavy notochord. The results of this study strongly support the hypothesis that the notochord kinks and twists are triggered by the down-regulation of hspb9 and hspb11, and intensified by body growth retardation along with normal notochord length in PCDPS-exposed zebrafish embryos.

Athabasca Oil Sands Petcoke Extract Elicits Biochemical and Transcriptomic Effects in Avian Hepatocytes.

Petroleum coke or "petcoke" is a granular carbonaceous material produced during the upgrading of heavy crude oils, including bitumen. Petcoke dust was recently reported as an environmental contaminant in the Athabasca oil sands region, but the ecotoxicological hazards posed by this complex bitumen-derived material-including those to avian species-have not been characterized. In this study, solvent extracts (x) of delayed and fluid petcoke (xDP and xFP) were prepared and dissolved in dimethyl sulfoxide. A water-accommodated fraction of delayed petcoke (waDP) was also prepared. Graded concentrations of xDP, xFP, and waDP were administered to chicken and double-crested cormorant hepatocytes to determine effects on 7-ethoxyresorufin-O-deethylase (EROD) activity, porphyrin accumulation, and mRNA expression. Polycyclic aromatic compounds (PACs) were characterized, and xDP, xFP, and waDP had total PAC concentrations of 93 000, 270, and 5.3 ng/mL. The rank order of biochemical and transcriptomic responses was xDP > xFP > waDP (e.g., EROD EC50s were lower for xDP compared to xFP and waDP). A total of 22, 18, and 4 genes were altered following exposure to the highest concentrations of xDP, xFP, and waDP, respectively, using a chicken PCR array comprising 27 AhR-related genes. To provide more exhaustive coverage of potential toxicity pathways being impacted, two avian ToxChip PCR arrays-chicken and double-crested cormorant-were utilized, and xDP altered the expression of more genes than xFP. Traditional PAC-related toxicity pathways and novel mechanisms of action were identified in two avian species following petcoke extract exposure. Extrapolation to real-world exposure scenarios must consider the bioavailability of the extracted PACs compared to those in exposed organisms.