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.

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.

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.

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

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.

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.

Integration of sperm DNA damage assessment into OECD test guidelines for genotoxicity testing using the MutaMouse model.

The Organisation for Economic Co-operation and Development (OECD) endorses test guidelines (TG) for identifying chemicals that are genotoxic, such as the transgenic rodent gene mutation assay (TG 488). Current OECD TG do not include assays for sperm DNA damage resulting in a critical testing gap. We evaluated the performance of the Sperm Chromatin Structure Assay (SCSA) and the Terminal Deoxynucleotidyl Transferase-Mediated Deoxyuridine Triphosphate Nick end Labeling (TUNEL) assay to detect sperm DNA damage within the recommended TG 488 protocol. MutaMouse males received 0, 0.5, 1, or 2 mg/kg/day triethylenemelamine (TEM), a multifunctional alkylating agent, for 28 days orally and tissues were collected two (blood) and three (sperm and bone marrow) days later. TEM significantly increased the frequency of lacZ mutants in bone marrow, and of micronuclei (MN) in both reticulocytes (%MN-RET) and normochromatic erythrocytes (%MN-NCE) in a dose-dependent manner (P < 0.05). The percentage of DNA fragmentation index (%DFI) and %TUNEL positive cells demonstrated dose-related increases in sperm (P < 0.05), and the two assay results were strongly correlated (R = 0.9298). Within the same animal, a good correlation was observed between %MN-NCE and %DFI (R = 0.7189). Finally, benchmark dose modelling (BMD) showed comparable BMD10 values among the somatic and germ cell assays. Our results suggest that sperm DNA damage assays can be easily integrated into standard OECD designs investigating genotoxicity in somatic tissues to provide key information on whether a chemical is genotoxic in germ cells and impact its risk assessment.

Tris(1,3-dichloro-2-propyl) phosphate perturbs the expression of genes involved in immune response and lipid and steroid metabolism in chicken embryos.

We previously demonstrated that in ovo exposure to the flame retardant tris(1,3-dichloro-2-propyl) phosphate (TDCPP) decreased plasma thyroxine levels, reduced growth parameters, and decreased gallbladder size in chicken embryos. In the current study DNA microarrays were used to evaluate global mRNA expression in liver tissue of male chicken embryos that exhibited the above mentioned effects. Injected doses were dimethyl sulfoxide vehicle control, 7.6 or 45 µg TDCPP/g egg. TDCPP caused significant changes in the expression of five genes at the low dose and 47 genes at the high dose (False Discovery Rate p ≤ 0.1, fold change ≥ 1.5). The gene expression analysis suggested a compromised immune function, a state of cholestatic liver/biliary fibrosis, and disrupted lipid and steroid metabolism. Circulating bile acid levels were elevated, which is an indication of liver dysfunction, and plasma cholesterol levels were reduced; however, hepatic bile acid and cholesterol levels were unaltered. Interactome analyses identified apolipoprotein E, hepatocyte nuclear factor 4 alpha, and peroxisome proliferator-activated receptor alpha as key regulatory molecules involved in the effects of TDCPP. Our results demonstrate a targeted effect of TDCPP toxicity on lipid metabolism, including cholesterol, that helps explain the aforementioned phenotypic effects, as chicken embryos are highly dependent on yolk lipids for growth and maintenance throughout development. Finally, our results are in concordance with the literature that describes TDCPP as a cancer-causing agent, since the majority of dysregulated genes were involved in cancer pathways.

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.

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.

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.

No evidence for transgenerational genomic instability in the F1 or F2 descendants of Muta™Mouse males exposed to N-ethyl-N-nitrosourea.

Exposure of male mice to genotoxic agents can increase mutation frequencies in their unexposed descendants. This phenomenon, known as transgenerational genomic instability (TGI), can persist for several generations. However, little is known about the underlying mechanisms. Chemically-induced TGI has been demonstrated in non-coding unstable tandem repeat DNA regions, but it is unclear whether it extends to other genetic endpoints. We investigated whether exposure of Muta™Mouse males to a single dose of 75mg/kg N-ethyl-N-nitrosourea (ENU) increased the spontaneous frequency of gene mutations or chromosome damage in their offspring. Treated males were mated with untreated females 3 days, 6 weeks or 10 weeks post-exposure to produce the F1 generation. Offspring were thus conceived from germ cells exposed to ENU as mature spermatozoa, dividing spermatogonia, or spermatogonial stem cells, respectively. F2 mice were generated by mating F1 descendants with untreated partners. Mutations in the lacZ transgene were quantified in bone marrow and micronucleus frequencies were evaluated in red blood cells by flow-cytometry for all F0 and their descendants. LacZ mutant frequencies were also determined in sperm for all exposed males and their male descendants. In F0 males, lacZ mutant frequencies were significantly increased in bone marrow at least 10-fold at all three time points investigated. In sperm, lacZ mutant frequency was significantly increased 7-11-fold after exposure of dividing and stem cell spermatogonia, but not in replication-deficient haploid sperm. Micronucleus frequencies assessed two days after ENU treatment were increased 5-fold in F0 males, but returned to control levels after 10 weeks. Despite the strong mutagenic response in F0 males, pre- and post-meiotic ENU exposure did not significantly increase lacZ mutant or micronucleus frequencies in F1 or F2 offspring. These findings suggest that TGI may not extend to all genetic endpoints and that further investigation of this phenomenon and its health relevance will require multiple measures of genomic damage.

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.

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.

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.

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.

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.

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.

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.

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.