Contribution of emissions from the oil sands activities in Alberta, Canada to atmospheric concentration and deposition of nitrogen and sulfur species at a downwind site.

Oil sands activities in the Athabasca Oil Sands Region in Alberta, Canada, are large sources of atmospheric NOx and SO2. This study investigated the impact of oil sands emissions on the atmospheric deposition of nitrogen and sulfur species at a downwind site, about 350 km from the oil sands facilities. Measurement data are from the Canadian Air and Precipitation Monitoring Network (CAPMoN) from 2015 to 2019, including ambient concentrations of HNO3, pNO3-, NO2, pNH4+, NH3, SO2, pSO42- and base cations, as well as concentrations of NO3-, SO42-, NH4+, and base cations in precipitation. Sector analysis of air mass back trajectories was conducted to distinguish measurements with different air mass origins. Median atmospheric concentrations and dry deposition fluxes of HNO3, pNO3-, NO2, pNH4+, pSO42-, and SO2 on days when the air masses came from the oil sands sector were significantly greater than those with the "Clean" sector by 34-67%, whereas the difference in NH3 concentration was not significant. Contributions of the oil sands emissions to dry deposition fluxes of these species ranged from 3.8 to 13.1%. The precipitation-weighted mean concentrations of NO3-, SO42-, and NH4+ in samples with the oil sands sector were 76 %, 65 % and 81 % greater than those with the "Clean" sector, respectively. Contributions of the oil sands emissions to wet deposition of NO3-, SO42-, and NH4+ were 12.5 ± 8.9 %, 8.7 ± 4.4 %, and 6.0 ± 3.3 %, respectively. The annual total deposition of nitrogen and sulfur were 1.9 kg-N ha-1 and 0.74 kg-S ha-1, respectively, of which 8.0 ± 3.5 % and 8.7 ± 3.6 % were from oil sands emissions. The total deposition of sulfur and nitrogen did not exceed the critical loads (CL) of acidity, but nitrogen deposition exceeded the CLs of nutrient nitrogen in the region.

Transcriptomic point of departure determination: a comparison of distribution-based and gene set-based approaches.

A key step in assessing the potential human and environmental health risks of industrial and agricultural chemicals is to determine the toxicity point of departure (POD), which is the highest dose level that causes no adverse effect. Transcriptomic POD (tPOD) values have been suggested to accurately estimate toxicity POD values. One step in the most common approach for tPOD determination involves mapping genes to annotated gene sets, a process that might lead to substantial information loss particularly in species with poor gene annotation. Alternatively, methods that calculate tPOD values directly from the distribution of individual gene POD values omit this mapping step. Using rat transcriptome data for 79 molecules obtained from Open TG-GATEs (Toxicogenomics Project Genomics Assisted Toxicity Evaluation System), the hypothesis was tested that methods based on the distribution of all individual gene POD values will give a similar tPOD value to that obtained via the gene set-based method. Gene set-based tPOD values using four different gene set structures were compared to tPOD values from five different individual gene distribution methods. Results revealed a high tPOD concordance for all methods tested, especially for molecules with at least 300 dose-responsive probesets: for 90% of those molecules, the tPOD values from all methods were within 4-fold of each other. In addition, random gene sets based upon the structure of biological knowledge-derived gene sets produced tPOD values with a median absolute fold change of 1.3-1.4 when compared to the original biological knowledge-derived gene set counterparts, suggesting that little biological information is used in the gene set-based tPOD generation approach. These findings indicate using individual gene distributions to calculate a tPOD is a viable and parsimonious alternative to using gene sets. Importantly, individual gene distribution-based tPOD methods do not require knowledge of biological organization and can be applied to any species including those with poorly annotated gene sets.

Empirical Characterization of False Discovery Rates of Differentially Expressed Genes and Transcriptomic Benchmark Concentrations in Zebrafish Embryos.

High-throughput transcriptomics (HTTr) is increasingly applied to zebrafish embryos to survey the toxicological effects of environmental chemicals. Before the adoption of this approach in regulatory testing, it is essential to characterize background noise in order to guide experimental designs. We thus empirically quantified the HTTr false discovery rate (FDR) across different embryo pool sizes, sample sizes, and concentration groups for toxicology studies. We exposed zebrafish embryos to 0.1% dimethyl sulfoxide (DMSO) for 5 days. Pools of 1, 5, 10, and 20 embryos were created (n = 24 samples for each pool size). Samples were sequenced on the TempO-Seq platform and then randomly assigned to mock treatment groups before differentially expressed gene (DEG), pathway, and benchmark concentration (BMC) analyses. Given that all samples were treated with DMSO, any significant DEGs, pathways, or BMCs are false positives. As expected, we found decreasing FDRs for DEG and pathway analyses with increasing pool and sample sizes. Similarly, FDRs for BMC analyses decreased with increasing pool size and concentration groups, with more stringent BMC premodel filtering reducing BMC FDRs. Our study provides foundational data for determining appropriate experiment designs for regulatory toxicity testing with HTTr in zebrafish embryos.

Wildlife ecological risk assessment in the 21st century: Promising technologies to assess toxicological effects.

Despite advances in toxicity testing and the development of new approach methodologies (NAMs) for hazard assessment, the ecological risk assessment (ERA) framework for terrestrial wildlife (i.e., air-breathing amphibians, reptiles, birds, and mammals) has remained unchanged for decades. While survival, growth, and reproductive endpoints derived from whole-animal toxicity tests are central to hazard assessment, nonstandard measures of biological effects at multiple levels of biological organization (e.g., molecular, cellular, tissue, organ, organism, population, community, ecosystem) have the potential to enhance the relevance of prospective and retrospective wildlife ERAs. Other factors (e.g., indirect effects of contaminants on food supplies and infectious disease processes) are influenced by toxicants at individual, population, and community levels, and need to be factored into chemically based risk assessments to enhance the "eco" component of ERAs. Regulatory and logistical challenges often relegate such nonstandard endpoints and indirect effects to postregistration evaluations of pesticides and industrial chemicals and contaminated site evaluations. While NAMs are being developed, to date, their applications in ERAs focused on wildlife have been limited. No single magic tool or model will address all uncertainties in hazard assessment. Modernizing wildlife ERAs will likely entail combinations of laboratory- and field-derived data at multiple levels of biological organization, knowledge collection solutions (e.g., systematic review, adverse outcome pathway frameworks), and inferential methods that facilitate integrations and risk estimations focused on species, populations, interspecific extrapolations, and ecosystem services modeling, with less dependence on whole-animal data and simple hazard ratios. Integr Environ Assess Manag 2024;20:725-748. © 2023 His Majesty the King in Right of Canada and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). Reproduced with the permission of the Minister of Environment and Climate Change Canada. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Cross-species applicability of an adverse outcome pathway network for thyroid hormone system disruption.

Thyroid hormone system disrupting compounds are considered potential threats for human and environmental health. Multiple adverse outcome pathways (AOPs) for thyroid hormone system disruption (THSD) are being developed in different taxa. Combining these AOPs results in a cross-species AOP network for THSD which may provide an evidence-based foundation for extrapolating THSD data across vertebrate species and bridging the gap between human and environmental health. This review aimed to advance the description of the taxonomic domain of applicability (tDOA) in the network to improve its utility for cross-species extrapolation. We focused on the molecular initiating events (MIEs) and adverse outcomes (AOs) and evaluated both their plausible domain of applicability (taxa they are likely applicable to) and empirical domain of applicability (where evidence for applicability to various taxa exists) in a THSD context. The evaluation showed that all MIEs in the AOP network are applicable to mammals. With some exceptions, there was evidence of structural conservation across vertebrate taxa and especially for fish and amphibians, and to a lesser extent for birds, empirical evidence was found. Current evidence supports the applicability of impaired neurodevelopment, neurosensory development (eg, vision) and reproduction across vertebrate taxa. The results of this tDOA evaluation are summarized in a conceptual AOP network that helps prioritize (parts of) AOPs for a more detailed evaluation. In conclusion, this review advances the tDOA description of an existing THSD AOP network and serves as a catalog summarizing plausible and empirical evidence on which future cross-species AOP development and tDOA assessment could build.

Impact of sampling time on the detection of mutations in rapidly proliferating tissues using transgenic rodent gene mutation models: A review.

The OECD Test Guideline 488 (TG 488) for the Transgenic Rodent Gene Mutation Assay has undergone several revisions to update the recommended design for studying mutations in somatic tissues and male germ cells. The recently revised TG recommends a single sampling time of 28 days following 28 days of exposure (i.e., 28 + 28 days) for all tissues, irrespective of proliferation rates. An alternative design (i.e., 28 + 3 days) is appropriate when germ cell data is not required, nor considered. While the 28 + 28 days design is clearly preferable for slowly proliferating somatic tissues and germ cells, there is still uncertainty about the impact of extending the sampling time to 28 days for rapidly somatic tissues. Here, we searched the available literature for evidence supporting the applicability and utility of the 28 + 28 days design for rapidly proliferating tissues. A total of 79 tests were identified. When directly comparing results from both designs in the same study, there was no evidence that the 28 + 28 days regimen resulted in a qualitatively different outcome from the 28 + 3 days design. Studies with a diverse range of agents that employed only a 28 + 28 days protocol provide further evidence that this design is appropriate for rapidly proliferating tissues. Benchmark dose analyses demonstrate high quantitative concordance between the 28 + 3 and 28 + 28 days designs for rapidly proliferating tissues. Accordingly, our review confirms that the 28 + 28 days design is appropriate to assess mutagenicity in both slowly and rapidly proliferating somatic tissues, and germ cells, and provides further support for the recommended design in the recently adopted TG 488.

A Case Study on Integrating a New Key Event Into an Existing Adverse Outcome Pathway on Oxidative DNA Damage: Challenges and Approaches in a Data-Rich Area.

Adverse outcome pathways (AOPs) synthesize toxicological information to convey and weigh evidence in an accessible format. AOPs are constructed in modules that include key events (KEs) and key event relationships (KERs). This modular structure facilitates AOP expansion and network development. AOP development requires finding relevant information to evaluate the weight of evidence supporting each KER. To do this, the use of transparent/reproducible search methods, such as systematic review (SR), have been proposed. Applying SR to AOP development in a data-rich area is difficult as SR requires screening each article returned from a search. Here we describe a case study to integrate a single new KE into an existing AOP. We explored the use of SR concepts and software to conduct a transparent and documented literature search to identify empirical data supporting the incorporation of a new KE, increase in cellular reactive oxygen species (ROS), upstream of an existing AOP: "Oxidative DNA Damage Leading to Chromosomal Aberrations and Mutations". Connecting this KE to the AOP is supported by the development of five new KERs, the most important being the first adjacent KER (increase in ROS leading to oxidative DNA damage). We initially searched for evidence of all five KERs and screened 100 papers to develop a preliminary evidence map. After removing papers not containing relevant data based on our Population, Exposure, Comparator and Outcome statement, 39 articles supported one or more KERs; these primarily addressed temporal or dose concordance of the non-adjacent KERs with limited evidence supporting the first adjacent KER. We thus conducted a second focused set of searches using search terms for specific methodologies to measure these first two KEs. After screening, 12 articles were identified that contained quantitative evidence supporting the first adjacent KER. Given that integrating a new KE into an existing AOP requires the development of multiple KERs, this approach of building a preliminary evidence map, focusing evidence gathering on the first adjacent KER, and applying reproducible search strategies using specific methodologies for the first adjacent KER, enabled us to prioritize studies to support expansion of this data-rich AOP.

Microsatellite mutation frequencies in river otters (Lontra Canadensis) from the Athabasca Oil Sands region are correlated to polycyclic aromatic compound tissue burden.

Mining activities in the Athabasca oil sands region (AOSR) have contributed to an increase of polycyclic aromatic compounds (PACs) locally. However, many PACs found in the AOSR, and the combined effects of PAC mixtures have not been evaluated for genotoxicity in wildlife. Here, we examine whether mutation frequencies in AOSR river otters are correlated to PAC tissue burdens. We used single-molecule polymerase chain reaction (SM-PCR) to measure the mutant frequency of unstable DNA microsatellite loci in the bone marrow of wild river otters (n = 11) from the AOSR. Microsatellite mutation frequencies were regressed against liver PAC burden (total, low/high molecular weight [LMW/HMW], and parent/alkylated PACs), and to the distances from where the samples were collected to nearby bitumen upgraders. We found that microsatellite mutation frequency was positively correlated with total liver PAC burden. LMW and alkylated PACs were detected at higher levels and had a stronger positive relationship with mutation frequency than HMW (alkylated and parent) PACs. There were no significant relationships detected between mutation frequency and LMW parent PACs or the distance from bitumen upgraders. Furthermore, pyrogenic and petrogenic signatures suggest PACs in animals with high mutation frequencies were associated with combustion processes; although further investigation is warranted, due to limitations of diagnostic ratio determination with biotic models. Our findings support the hypothesis that PACs found in the AOSR increase mutation frequency in wildlife. Further investigation is required to determine if the elevated PAC levels associated with higher mutation frequency are due to natural exposure or elevated human activity.

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.

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.

A Pragmatic Approach to Adverse Outcome Pathway Development and Evaluation.

The adverse outcome pathway (AOP) framework provides a practical means for organizing scientific knowledge that can be used to infer cause-effect relationships between stressor events and toxicity outcomes in intact organisms. It has reached wide acceptance as a tool to aid chemical safety assessment and regulatory toxicology by supporting a systematic way of predicting adverse health outcomes based on accumulated mechanistic knowledge. A major challenge for broader application of the AOP concept in regulatory toxicology, however, has been developing robust AOPs to a level where they are peer reviewed and accepted. This is because the amount of work required to substantiate the modular units of a complete AOP is considerable, to the point where it can take years from start to finish. To help alleviate this bottleneck, we propose a more pragmatic approach to AOP development whereby the focus becomes on smaller blocks. First, we argue that the key event relationship (KER) should be formally recognized as the core building block of knowledge assembly within the AOP knowledge base (AOP-KB), albeit framing them within full AOPs to ensure regulatory utility. Second, we argue that KERs should be developed using systematic review approaches, but only in cases where the underlying concept does not build on what is considered canonical knowledge. In cases where knowledge is considered canonical, rigorous systematic review approaches should not be required. It is our hope that these approaches will contribute to increasing the pace at which the AOP-KB is populated with AOPs with utility for chemical safety assessors and regulators.

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.

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.

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

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

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.

Chemically induced mutations in a MutaMouse reporter gene inform mechanisms underlying human cancer mutational signatures.

Transgenic rodent (TGR) models use bacterial reporter genes to quantify in vivo mutagenesis. Pairing TGR assays with next-generation sequencing (NGS) enables comprehensive mutation pattern analysis to inform mutational mechanisms. We used this approach to identify 2751 independent lacZ mutations in the bone marrow of MutaMouse animals exposed to four chemical mutagens: benzo[a]pyrene, N-ethyl-N-nitrosourea, procarbazine, and triethylenemelamine. We also collected published data for 706 lacZ mutations from eight additional environmental mutagens. We report that lacZ gene sequencing generates chemical-specific mutation signatures observed in human cancers with established environmental causes. For example, the mutation signature of benzo[a]pyrene, a carcinogen present in tobacco smoke, matched the signature associated with tobacco-induced lung cancers. Our results suggest that the analysis of chemically induced mutations in the lacZ gene shortly after exposure provides an effective approach to characterize human-relevant mechanisms of carcinogenesis and propose novel environmental causes of mutation signatures observed in human cancers.

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.