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 Transformative Vision for an Omics-Based Regulatory Chemical Testing Paradigm.

Use of molecular data in human and ecological health risk assessments of industrial chemicals and agrochemicals has been anticipated by the scientific community for many years; however, these data are rarely used for risk assessment. Here, a logic framework is proposed to explore the feasibility and future development of transcriptomic methods to refine and replace the current apical endpoint-based regulatory toxicity testing paradigm. Four foundational principles are outlined and discussed that would need to be accepted by stakeholders prior to this transformative vision being realized. Well-supported by current knowledge, the first principle is that transcriptomics is a reliable tool for detecting alterations in gene expression that result from endogenous or exogenous influences on the test organism. The second principle states that alterations in gene expression are indicators of adverse or adaptive biological responses to stressors in an organism. Principle 3 is that transcriptomics can be employed to establish a benchmark dose-based point of departure (POD) from short-term, in vivo studies at a dose level below which a concerted molecular change (CMC) is not expected. Finally, Principle 4 states that the use of a transcriptomic POD (set at the CMC dose level) will support a human health-protective risk assessment. If all four principles are substantiated, this vision is expected to transform aspects of the industrial chemical and agrochemical risk assessment process that are focused on establishing safe exposure levels for mammals across numerous toxicological contexts resulting in a significant reduction in animal use while providing equal or greater protection of human health. Importantly, these principles and approaches are also generally applicable for ecological safety assessment.

Reactive Oxygen Species in the Adverse Outcome Pathway Framework: Toward Creation of Harmonized Consensus Key Events.

Reactive oxygen species (ROS) and reactive nitrogen species (RNS) are formed as a result of natural cellular processes, intracellular signaling, or as adverse responses associated with diseases or exposure to oxidizing chemical and non-chemical stressors. The action of ROS and RNS, collectively referred to as reactive oxygen and nitrogen species (RONS), has recently become highly relevant in a number of adverse outcome pathways (AOPs) that capture, organize, evaluate and portray causal relationships pertinent to adversity or disease progression. RONS can potentially act as a key event (KE) in the cascade of responses leading to an adverse outcome (AO) within such AOPs, but are also known to modulate responses of events along the AOP continuum without being an AOP event itself. A substantial discussion has therefore been undertaken in a series of workshops named "Mystery or ROS" to elucidate the role of RONS in disease and adverse effects associated with exposure to stressors such as nanoparticles, chemical, and ionizing and non-ionizing radiation. This review introduces the background for RONS production, reflects on the direct and indirect effects of RONS, addresses the diversity of terminology used in different fields of research, and provides guidance for developing a harmonized approach for defining a common event terminology within the AOP developer community.

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.

Temporal change and the influence of climate and weather factors on mercury concentrations in Hudson Bay polar bears, caribou, and seabird eggs.

Temporal trends of mercury in Arctic wildlife are inconsistent within and between species and are often insignificant, which limits data interpretation. Recent multivariate analyses have shown that weather and climate factors (e.g. temperatures, sea ice conditions) are related to total Hg (THg) concentrations in wildlife tissues, though relatively few studies have explored these relationships. The present study compared time series of THg concentrations in liver of polar bear (Ursus maritimus, 2007/08-2015/16), eggs of thick-billed murres (Uria lomvia, 1993-2015) and kidney of caribou (Rangifer tarandus groenlandicus, 2006-2015) from the Hudson Bay region of Canada and statistically modelled THg over time with available climate and weather data. Significant temporal trends of THg concentrations were not detected in any species. However, in multivariate models that included time-lagged sea ice freeze up dates, THg concentrations increased 4.4% yr-1 in Qamanirjuaq caribou. Sea ice conditions were also related to THg levels in polar bear liver but not those in eggs of murres, though year was not a signifcant factor. Greater precipitation levels one to two years prior to sampling were associated with greater THg concentrations in polar bears and caribou, likely due to greater deposition, flooding and discharge from nearby wetlands and rivers. Time-lagged Arctic and/or North Atlantic Oscillation (AO/NAO) indices also generated significant, inverse models for all three species, agreeing with relationships in other time series of similar length. The magnitude and direction of many relationships were affected by season, duration of time-lags, and the length of the time series. Our findings support recent observations suggesting that temporal studies monitoring Hg in Arctic wildlife should consider including key climatic or weather factors to help identify consistent variables of influence and to improve temporal analyses of THg time series.

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.

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

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

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

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

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.

Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology.

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

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.