Development and application of consensus in silico models for advancing high-throughput toxicological predictions

Computational toxicology models have been successfully implemented to prioritize and screen chemicals. There are numerous in silico (quantitative) structure-activity relationship ([Q]SAR) models for the prediction of a range of human-relevant toxicological endpoints, but for a given endpoint and chemical, not all predictions are identical due to differences in their training sets, algorithms, and methodology. This poses an issue for high-throughput screening of a large chemical inventory as it necessitates several models to cover diverse chemistries but will then generate data conflicts. To address this challenge, we developed a consensus modeling strategy to combine predictions obtained from different existing in silico (Q)SAR models into a single predictive value while also expanding chemical space coverage. This study developed consensus models for nine toxicological endpoints relating to estrogen receptor (ER) and androgen receptor (AR) interactions (i.e., binding, agonism, and antagonism) and genotoxicity (i.e., bacterial mutation, in vitro chromosomal aberration, and in vivo micronucleus). Consensus models were created by combining different (Q)SAR models using various weighting schemes. As a multi-objective optimization problem, there is no single best consensus model, and therefore, Pareto fronts were determined for each endpoint to identify the consensus models that optimize the multiple-criterion decisions simultaneously. Accordingly, this work presents sets of solutions for each endpoint that contain the optimal combination, regardless of the trade-off, with the results demonstrating that the consensus models improved both the predictive power and chemical space coverage. These solutions were further analyzed to find trends between the best consensus models and their components. Here, we demonstrate the development of a flexible and adaptable approach for in silico consensus modeling and its application across nine toxicological endpoints related to ER activity, AR activity, and genotoxicity. These consensus models are developed to be integrated into a larger multi-tier NAM-based framework to prioritize chemicals for further investigation and support the transition to a non-animal approach to risk assessment in Canada.

Interpretation of in vitro concentration-response data for risk assessment and regulatory decision-making: Report from the 2022 IWGT quantitative analysis expert working group meeting

Quantitative risk assessments of chemicals are routinely performed using in vivo data from rodents; however, there is growing recognition that non-animal approaches can be human-relevant alternatives. There is an urgent need to build confidence in non-animal alternatives given the international support to reduce the use of animals in toxicity testing where possible. In order for scientists and risk assessors to prepare for this paradigm shift in toxicity assessment, standardization and consensus on in vitro testing strategies and data interpretation will need to be established. To address this issue, an Expert Working Group (EWG) of the 8th International Workshop on Genotoxicity Testing (IWGT) evaluated the utility of quantitative in vitro genotoxicity concentration-response data for risk assessment. The EWG first evaluated available in vitro methodologies and then examined the variability and maximal response of in vitro tests to estimate biologically relevant values for the critical effect sizes considered adverse or unacceptable. Next, the EWG reviewed the approaches and computational models employed to provide human-relevant dose context to in vitro data. Lastly, the EWG evaluated risk assessment applications for which in vitro data are ready for use and applications where further work is required. The EWG concluded that in vitro genotoxicity concentration-response data can be interpreted in a risk assessment context. However, prior to routine use in regulatory settings, further research will be required to address the remaining uncertainties and limitations.

Quantitative in vitro to in vivo extrapolation of genotoxicity data provides protective estimates of in vivo dose

Genotoxicity assessment is a critical component in the development and evaluation of chemicals. Traditional genotoxicity assays (i.e., mutagenicity, clastogenicity, and aneugenicity) have been limited to dichotomous hazard classification, while other toxicity endpoints are assessed through quantitative determination of points-of-departures (PODs) for setting exposure limits. The more recent higher-throughput in vitro genotoxicity assays, many of which also provide mechanistic information, offer a powerful approach for determining defined PODs for potency ranking and risk assessment. In order to obtain relevant human dose context from the in vitro assays, in vitro to in vivo extrapolation (IVIVE) models are required to determine what dose would elicit a concentration in the body demonstrated to be genotoxic using in vitro assays. Previous work has demonstrated that application of IVIVE models to in vitro bioactivity data can provide PODs that are protective of human health, but there has been no evaluation of how these models perform with in vitro genotoxicity data. Thus, the Genetic Toxicology Technical Committee, under the Health and Environmental Sciences Institute, conducted a case study on 31 reference chemicals to evaluate the performance of IVIVE application to genotoxicity data. The results demonstrate that for most chemicals considered here (20/31), the PODs derived from in vitro data and IVIVE are health protective relative to in vivo PODs from animal studies. PODs were also protective by assay target: mutations (8/13 chemicals), micronuclei (9/12), and aneugenicity markers (4/4). It is envisioned that this novel testing strategy could enhance prioritization, rapid screening, and risk assessment of genotoxic chemicals.

Application of a new approach methodology (NAM)-based strategy for genotoxicity assessment of data-poor compounds

The conventional battery for genotoxicity testing is not well suited to assessing the large number of chemicals needing evaluation. Traditional in vitro tests lack throughput, provide little mechanistic information, and have poor specificity in predicting in vivo genotoxicity. New Approach Methodologies (NAMs) aim to accelerate the pace of hazard assessment and reduce reliance on in vivo tests that are time-consuming and resource-intensive. As such, high-throughput transcriptomic and flow cytometry-based assays have been developed for modernized in vitro genotoxicity assessment. This includes: the TGx-DDI transcriptomic biomarker (i.e., 64-gene expression signature to identify DNA damage-inducing (DDI) substances), the MicroFlow^® assay (i.e., a flow cytometry-based micronucleus (MN) test), and the MultiFlow^® assay (i.e., a multiplexed flow cytometry-based reporter assay that yields mode of action (MoA) information). The objective of this study was to investigate the utility of the TGx-DDI transcriptomic biomarker, multiplexed with the MicroFlow^® and MultiFlow^® assays, as an integrated NAM-based testing strategy for screening data-poor compounds prioritized by Health Canada's New Substances Assessment and Control Bureau. Human lymphoblastoid TK6 cells were exposed to 3 control and 10 data-poor substances, using a 6-point concentration range. Gene expression profiling was conducted using the targeted TempO-Seq™ assay, and the TGx-DDI classifier was applied to the dataset. Classifications were compared with those based on the MicroFlow^® and MultiFlow^® assays. Benchmark Concentration (BMC) modeling was used for potency ranking. The results of the integrated hazard calls indicate that five of the data-poor compounds were genotoxic in vitro, causing DNA damage via a clastogenic MoA, and one via a pan-genotoxic MoA. Two compounds were likely irrelevant positives in the MN test; two are considered possibly genotoxic causing DNA damage via an ambiguous MoA. BMC modeling revealed nearly identical potency rankings for each assay. This ranking was maintained when all endpoint BMCs were converted into a single score using the Toxicological Prioritization (ToxPi) approach. Overall, this study contributes to the establishment of a modernized approach for effective genotoxicity assessment and chemical prioritization for further regulatory scrutiny. We conclude that the integration of TGx-DDI, MicroFlow^®, and MultiFlow^® endpoints is an effective NAM-based strategy for genotoxicity assessment of data-poor compounds.

Establishing a quantitative framework for regulatory interpretation of genetic toxicity dose-response data: Margin of exposure case study of 48 compounds with both in vivo mutagenicity and carcinogenicity dose-response data

Quantitative relationships between carcinogenic potency and mutagenic potency have been previously examined using a benchmark dose (BMD)-based approach. We extended those analyses by using human exposure data for 48 compounds to calculate carcinogenicity-derived and genotoxicity-derived margin of exposure values (MOEs) that can be used to prioritize substances for risk management. MOEs for 16 of the 48 compounds were below 10,000, and consequently highlighted for regulatory concern. Of these, 15 were highlighted using genotoxicity-derived (micronucleus [MN] dose-response data) MOEs. A total of 13 compounds were highlighted using carcinogenicity-derived MOEs; 12 compounds were overlapping. MOEs were also calculated using transgenic rodent (TGR) mutagenicity data. For 10 of the 12 compounds examined using TGR data, the results similarly revealed that mutagenicity-derived MOEs yield regulatory decisions that correspond with those based on carcinogenicity-derived MOEs. The effect of benchmark response (BMR) on MOE determination was also examined. Reinterpretation of the analyses using a BMR of 50% indicated that four out of 15 compounds prioritized using MN-derived MOEs based on a default BMR of 5% would have been missed. The results indicate that regulatory decisions based on in vivo genotoxicity dose-response data would be consistent with those based on carcinogenicity dose-response data; in some cases, genotoxicity-based decisions would be more conservative. Going forward, and in the absence of carcinogenicity data, in vivo genotoxicity assays (MN and TGR) can be used to effectively prioritize substances for regulatory action. Routine use of the MOE approach necessitates the availability of reliable human exposure estimates, and consensus regarding appropriate BMRs for genotoxicity endpoints.

2,3,7,8-Tetrachlorodibenzo-p-Dioxin (TCDD)-Inducible Poly-ADP-Ribose Polymerase (TIPARP/PARP7) Catalytic Mutant Mice (TiparpH532A) Exhibit Increased Sensitivity to TCDD-Induced Hepatotoxicity and Lethality

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-inducible poly-adenosine diphosphate (ADP)-ribose polymerase (TIPARP/PARP7), an aryl hydrocarbon receptor (AHR) target gene and mono-ADP-ribosyltransferase, acts as part of a negative feedback loop to repress AHR signaling. This process is prevented by a single H532A mutation in TIPARP that destroys its catalytic activity. We hypothesized that the loss of TIPARP catalytic activity would increase sensitivity to TCDD-induced toxicity in vivo. To test this, we created a catalytically deficient mouse line (Tiparp^H532A) by introducing a single H532A mutation in TIPARP. Treatment of mouse embryonic fibroblasts or hepatocytes isolated from Tiparp^H532A mice confirmed the increased TCDD-induced expression of the AHR target genes Cyp1a1, Cyp1b1, and Tiparp. Tiparp^H532A mice given a single injection of 10 µg/kg TCDD, a nonlethal dose in Tiparp^+/+ mice, did not survive beyond day 10. All Tiparp^+/+ mice survived the 30-day treatment. TCDD-treated Tiparp^H532A mice displayed increased expression of AHR target genes, increased steatohepatitis and hepatotoxicity. Hepatic RNA-sequencing revealed 7-fold more differentially expressed genes in Tiparp^H532A mice than in Tiparp^+/+ mice (4542 vs 647 genes) 6 days after TCDD treatment. Differentially expressed genes included genes involved in xenobiotic metabolism, lipid homeostasis and inflammation. Taken together, these data further support TIPARP as a critical negative regulator of AHR activity and show that loss of its catalytic activity is sufficient to increase sensitivity to TCDD-induced steatohepatitis and lethality. Since TIPARP inhibition has recently emerged as a potential anticancer therapy, the impact on AHR signaling, TCDD and polycyclic aromatic hydrocarbon toxicity will need to be carefully considered under conditions of therapeutic TIPARP inhibition.

In vitro mutagenicity of selected environmental carcinogens and their metabolites in MutaMouse FE1 lung epithelial cells

Chemicals in commerce or under development must be assessed for genotoxicity; assessment is generally conducted using validated assays (e.g. Tk mouse lymphoma assay) as part of a regulatory process. Currently, the MutaMouse FE1 cell mutagenicity assay is undergoing validation for eventual use as a standard in vitro mammalian mutagenicity assay. FE1 cells have been shown to be metabolically competent with respect to some cytochrome P450 (CYP) isozymes; for instance, they can convert the human carcinogen benzo[a]pyrene into its proximate mutagenic metabolite. However, some contradictory results have been noted for other genotoxic carcinogens that require two-step metabolic activation (e.g. 2-acetylaminofluorene and 2-amino-3-methylimidazo[4,5-f]quinoxaline). Here, we examined three known or suspected human carcinogens, namely acrylamide, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) and 4-aminobiphenyl (4-ABP), together with their proximate metabolites (i.e. glycidamide, N-OH-PhIP and N-OH-4-ABP), to aid in the validation of the FE1 cell mutagenicity assay. Assessments of the parent compounds were conducted both in the presence and absence of an exogenous metabolic activation mixture S9; assessments of the metabolites were in the absence of S9. The most potent compound was N-OH-PhIP -S9, which elicited a mutant frequency (MF) level 5.3-fold over background at 5 µM. There was a 4.3-fold increase for PhIP +S9 at 5 µM, a 1.7-fold increase for glycidamide -S9 at 3.5 mM and a 1.5-fold increase for acrylamide +S9 at 4 mM. Acrylamide -S9 elicited a marginal 1.4-fold MF increase at 8 mM. Treatment with PhIP -S9, 4-ABP ±S9 and N-OH-4-ABP -S9 failed to elicit significant increases in lacZ MF with any of the treatment conditions tested. Gene expression of key CYP isozymes was quantified by RT-qPCR. Cyp1a1, 1a2 and 1b1 are required to metabolise PhIP and 4-ABP. Results showed that treatment with both compounds induced expression of Cyp1a1 and Cyp1b1 but not Cyp1a2. Cyp2e1, which catalyses the bioactivation of acrylamide to glycidamide, was not induced after acrylamide treatment. Overall, our results confirm that the FE1 cell mutagenicity assay has the potential for use alongside other, more traditional in vitro mutagenicity assays.

Quantitative Interpretation of Genetic Toxicity Dose-Response Data for Risk Assessment and Regulatory Decision-Making: Current Status and Emerging Priorities

The screen-and-bin approach for interpretation of genotoxicity data is predicated on three false assumptions: that genotoxicants are rare, that genotoxicity dose-response functions do not contain a low-dose region mechanistically characterized by zero-order kinetics, and that genotoxicity is not a bona fide toxicological endpoint. Consequently, there is a need to develop and implement quantitative methods to interpret genotoxicity dose-response data for risk assessment and regulatory decision-making. Standardized methods to analyze dose-response data, and determine point-of-departure (PoD) metrics, have been established; the most robust PoD is the benchmark dose (BMD). However, there are no standards for regulatory interpretation of mutagenicity BMDs. Although 5-10% is often used as a critical effect size (CES) for BMD determination, values for genotoxicity endpoints have not been established. The use of BMDs to determine health-based guidance values (HBGVs) requires assessment factors (AFs) to account for interspecies differences and variability in human sensitivity. Default AFs used for other endpoints may not be appropriate for interpretation of in vivo mutagenicity BMDs. Analyses of published dose-response data showing the effects of compensatory pathway deficiency indicate that AFs for sensitivity differences should be in the range of 2-20. Additional analyses indicate that the AF to compensate for short treatment durations should be in the range of 5-15. Future work should use available data to empirically determine endpoint-specific CES values; similarly, to determine AF values for BMD adjustment. Future work should also evaluate the ability to use in vitro dose-response data for risk assessment, and the utility of probabilistic methods for determination of mutagenicity HBGVs. Environ. Mol. Mutagen. 61:66-83, 2020. © 2019 Her Majesty the Queen in Right of Canada.

A comparison of transgenic rodent mutation and in vivo comet assay responses for 91 chemicals

A database of 91 chemicals with published data from both transgenic rodent mutation (TGR) and rodent comet assays has been compiled. The objective was to compare the sensitivity of the two assays for detecting genotoxicity. Critical aspects of study design and results were tabulated for each dataset. There were fewer datasets from rats than mice, particularly for the TGR assay, and therefore, results from both species were combined for further analysis. TGR and comet responses were compared in liver and bone marrow (the most commonly studied tissues), and in stomach and colon evaluated either separately or in combination with other GI tract segments. Overall positive, negative, or equivocal test results were assessed for each chemical across the tissues examined in the TGR and comet assays using two approaches: 1) overall calls based on weight of evidence (WoE) and expert judgement, and 2) curation of the data based on a priori acceptability criteria prior to deriving final tissue specific calls. Since the database contains a high prevalence of positive results, overall agreement between the assays was determined using statistics adjusted for prevalence (using AC1 and PABAK). These coefficients showed fair or moderate to good agreement for liver and the GI tract (predominantly stomach and colon data) using WoE, reduced agreement for stomach and colon evaluated separately using data curation, and poor or no agreement for bone marrow using both the WoE and data curation approaches. Confidence in these results is higher for liver than for the other tissues, for which there were less data. Our analysis finds that comet and TGR generally identify the same compounds (mainly potent mutagens) as genotoxic in liver, stomach and colon, but not in bone marrow. However, the current database content precluded drawing assay concordance conclusions for weak mutagens and non-DNA reactive chemicals.

Benchmark dose analyses of multiple genetic toxicity endpoints permit robust, cross-tissue comparisons of MutaMouse responses to orally delivered benzo[a]pyrene

Genetic damage is a key event in tumorigenesis, and chemically induced genotoxic effects are a human health concern. Although genetic toxicity data have historically been interpreted using a qualitative screen-and-bin approach, there is increasing interest in quantitative analysis of genetic toxicity dose-response data. We demonstrate an emerging use of the benchmark dose (BMD)-approach for empirically ranking cross-tissue sensitivity. Using a model environmental carcinogen, we quantitatively examined responses for four genetic damage endpoints over an extended dose range, and conducted cross-tissue sensitivity rankings using BMD100 values and their 90% confidence intervals (CIs). MutaMouse specimens were orally exposed to 11 doses of benzo[a]pyrene. DNA adduct frequency and lacZ mutant frequency (MF) were measured in up to 8 tissues, and Pig-a MF and micronuclei (MN) were assessed in immature (RETs) and mature red blood cells (RBCs). The cross-tissue BMD pattern for lacZ MF is similar to that observed for DNA adducts, and is consistent with an oral route-of-exposure and differences in tissue-specific metabolism and proliferation. The lacZ MF BMDs were significantly correlated with the tissue-matched adduct BMDs, demonstrating a consistent adduct conversion rate across tissues. The BMD CIs, for both the Pig-a and the MN endpoints, overlapped for RETs and RBCs, suggesting comparable utility of both cell populations for protracted exposures. Examination of endpoint-specific response maxima illustrates the difficulty of comparing BMD values for a fixed benchmark response across endpoints. Overall, the BMD-approach permitted robust comparisons of responses across tissues/endpoints, which is valuable to our mechanistic understanding of how benzo[a]pyrene induces genetic damage.

Performance of the in vitro transgene mutation assay in MutaMouse FE1 cells: Evaluation of nine misleading ("False") positive chemicals

The screening of chemicals for the protection of human health and the environment requires the assessment of genetic toxicity. However, existing, internationally-accepted in vitro mammalian genotoxicity tests have been criticized for their low specificity (i.e. high frequency of "false" or "misleading" positive results for compounds that are negative in vivo). An in vitro transgene mutation assay has been established that uses a metabolically competent cell line derived from MutaMouse lung (i.e. FE1 cells). Mutation scoring employs the well-characterized lacZ positive selection system, and the assay is proposed as an alternative in vitro assessment tool. In this study, the performance of the FE1 cell assay was evaluated by examining responses to nine non-DNA-reactive chemicals that previously elicited misleading positive results in other mammalian cell genotoxicity assays. FE1 cells were exposed to concentrations up to approximately 10 mM and/or concentrations that yielded approximately 80-90% cytotoxicity (as measured by relative increase in cell count). The assay demonstrated excellent specificity; exposures to the chemicals examined did not yield any positive responses even when tested in the presence of an exogenous metabolic activation system (i.e. S9) or with an extended sampling time. These results indicate that the FE1 cell mutagenicity assay is an effective and practical alternative to traditional mammalian cell gene mutation assays. The development and validation of effective in vitro tools such as the MutaMouse FE1 cell assay will contribute to international efforts to reduce, refine, and replace experimental animals for toxicity assessment. Environ. Mol. Mutagen. 58:582-591, 2017. © 2017 Wiley Periodicals, Inc.

Genetic Toxicity of Complex Mixtures of Polycyclic Aromatic Hydrocarbons: Evaluating Dose-Additivity in a Transgenic Mouse Model

This study evaluates the risk assessment approach currently employed for polycyclic aromatic hydrocarbon (PAH)-contaminated media, wherein carcinogenic hazards are evaluated using a dose-addition model that employs potency equivalency factors (PEFs) for targeted carcinogenic PAHs. Here, MutaMouse mice were subchronically exposed to PAH mixtures (p.o.), and mutagenic potency (MP) values were determined for five tissues. Predicted dose-additive mixture MPs were generated by summing the products of the concentrations and MPs of the individual targeted PAHs; values were compared to the experimental MPs of the mixtures to evaluate dose-additivity. Additionally, the PEF-determined BaP-equivalent concentrations were compared to those determined using a bioassay-derived method (BDM) (i.e., an additivity-independent approach). In bone marrow, mixture mutagenicity was less than dose-additive and the PEF-method provided higher estimates of BaP-equivalents than the BDM. Conversely, mixture mutagenicity in site-of-contact tissues (e.g., small intestine) was generally more than dose-additive and the PEF-method provided lower estimates of BaP-equivalents than the BDM. Overall, this study demonstrates that dose-additive predictions of mixture mutagenic potency based on the concentrations and potencies of a small number of targeted PAHs results in values that are surprisingly close to those determined experimentally, providing support for the dose-additive assumption employed for human health risk assessment of PAH mixtures.

Mitochondrial DNA exhibits resistance to induced point and deletion mutations

The accumulation of somatic mitochondrial DNA (mtDNA) mutations contributes to the pathogenesis of human disease. Currently, mitochondrial mutations are largely considered results of inaccurate processing of its heavily damaged genome. However, mainly from a lack of methods to monitor mtDNA mutations with sufficient sensitivity and accuracy, a link between mtDNA damage and mutation has not been established. To test the hypothesis that mtDNA-damaging agents induce mtDNA mutations, we exposed Muta^TMMouse mice to benzo[a]pyrene (B[a]P) or N-ethyl-N-nitrosourea (ENU), daily for 28 consecutive days, and quantified mtDNA point and deletion mutations in bone marrow and liver using our newly developed Digital Random Mutation Capture (dRMC) and Digital Deletion Detection (3D) assays. Surprisingly, our results demonstrate mutagen treatment did not increase mitochondrial point or deletion mutation frequencies, despite evidence both compounds increase nuclear DNA mutations and demonstrated B[a]P adduct formation in mtDNA. These findings contradict models of mtDNA mutagenesis that assert the elevated rate of mtDNA mutation stems from damage sensitivity and abridged repair capacity. Rather, our results demonstrate induced mtDNA damage does not readily convert into mutation. These findings suggest robust mitochondrial damage responses repress induced mutations after mutagen exposure.

Oral exposure to commercially available coal tar-based pavement sealcoat induces murine genetic damage and mutations

Coal tar (CT) is a thick black liquid produced as a by-product of coal carbonization to produce coke or manufactured gas. It is comprised a complex mixture of polycyclic aromatic compounds, including a wide range of polycyclic aromatic hydrocarbons (PAHs), many of which are genotoxic and carcinogenic. CT is used in some pavement sealants (also known as sealcoat), which are applied to pavement in order to seal and beautify the surface. Human exposure is known to occur not only during application, but also as a result of the weathering process, as elevated levels of PAHs have been found in settled house dust in residences adjacent to CT-sealed surfaces. In this study we examined the genotoxicity of an extract of a commercially available CT-based sealcoat in the transgenic Muta™Mouse model. Mice were orally exposed to 3 doses of sealcoat extract daily for 28 days. We evaluated genotoxicity by examining: (1) stable DNA adducts and (2) lacZ mutations in bone marrow, liver, lung, small intestine, and glandular stomach, as well as (3) micronucleated red blood cells. Significant increases were seen for each endpoint and in all tissues. The potency of the response differed across tissues, with the highest frequency of adducts occurring in liver and lung, and the highest frequency of mutations occurring in small intestine. The results of this study are the first demonstration of mammalian genotoxicity following exposure to CT-containing pavement sealcoat. This work provides in vivo evidence to support the contention that there may be adverse health effects in mammals, and potentially in humans, from exposure to coal tar. Environ. Mol. Mutagen. 57:535-545, 2016. © 2016 Her Majesty the Queen in Right of Canada.

Transcriptional profiling of the mouse hippocampus supports an NMDAR-mediated neurotoxic mode of action for benzo[a]pyrene

Benzo[a]pyrene (BaP) is a genotoxic carcinogen and a neurotoxicant. The neurotoxicity of BaP is proposed to arise from either genotoxicity leading to neuronal cell death, or perturbed expression of N-methyl-d-aspartate receptor (NMDAR) subunits. To explore these hypotheses, we profiled hippocampal gene expression of adult male Muta(™) Mouse administered 0, 1, 35, or 70 mg BaP/kg bw per day by oral gavage for 3 days. Transcriptional profiles were examined by RNA-sequencing (RNA-seq), DNA microarrays, and real-time quantitative reverse transcription polymerase chain reaction (RT-PCR). BaP-DNA adducts in the cerebellum were quantified by (32) P-post-labeling to measure genotoxicity. RNA-seq revealed altered expression of 0, 260, and 219 genes (P-value < 0.05, fold-change ≥ ± 1.5) following exposure to the low, medium, and high doses, respectively; 54 genes were confirmed by microarrays. Microarray and RT-PCR analysis showed increased expression of NMDAR subunits Grina and Grin2a. In contrast, no effects on DNA-damage response genes were observed despite comparable BaP-DNA adduct levels in the cerebellum and in the lungs and livers of mice at similar BaP doses in previous studies. The results suggest that DNA-damage response does not play a major role in BaP-induced adult neurotoxicity. Meta-analysis revealed that BaP-induced transcriptional profiles are highly correlated with those from the hippocampus of transgenic mice exhibiting similar neurotoxicity outcomes to BaP-exposed mice and rats (i.e., defects in learning and memory). Overall, we suggest that BaP-induced neurotoxicity is more likely to be a consequence of NMDAR perturbation than genotoxicity, and identify other important genes potentially mediating this adverse outcome. Environ. Mol. Mutagen. 57:350-363, 2016. © 2016 Her Majesty the Queen in Right of Canada. Environmental and Molecular Mutagenesis © 2016 Environmental Mutagen Society.

Empirical analysis of BMD metrics in genetic toxicology part II: in vivo potency comparisons to promote reductions in the use of experimental animals for genetic toxicity assessment

Genotoxicity tests have traditionally been used only for hazard identification, with qualitative dichotomous groupings being used to identify compounds that have the capacity to induce mutations and/or cytogenetic alterations. However, there is an increasing interest in employing quantitative analysis of in vivo dose-response data to derive point of departure (PoD) metrics that can be used to establish human exposure limits or margins of exposure (MOEs), thereby supporting human health risk assessments and regulatory decisions. This work is an extension of our companion article on in vitro dose-response analyses and outlines how the combined benchmark dose (BMD) approach across included covariates can be used to improve the analyses and interpretation of in vivo genetic toxicity dose-response data. Using the BMD-covariate approach, we show that empirical comparisons of micronucleus frequency dose-response data across multiple studies justifies dataset merging, with subsequent analyses improving the precision of BMD estimates and permitting attendant potency ranking of seven clastogens. Similarly, empirical comparisons of Pig-a mutant phenotype frequency data collected in males and females justified dataset merging across sex. This permitted more effective scrutiny regarding the effect of post-exposure sampling time on the mutagenicity of N-ethyl-N-nitrosourea observed in reticulocytes and erythrocytes in the Pig-a assay. The BMD-covariate approach revealed tissue-specific differences in the induction of lacZ transgene mutations in Muta™Mouse specimens exposed to benzo[a]pyrene (BaP), with the results permitting the formulation of mechanistic hypotheses regarding the observed potency ranking. Lastly, we illustrate how historical dose-response data for assessments that examined numerous doses (i.e. induced lacZ mutant frequency (MF) across 10 doses of BaP) can be used to improve the precision of BMDs derived from datasets with far fewer doses (i.e. lacZ MF for 3 doses of dibenz[a,h]anthracene). Collectively, the presented examples illustrate how innovative use of the BMD approach can permit refinement of the use of in vivo data; improving the efficacy of experimental animal use in genetic toxicology without sacrificing PoD precision.

Tissue-specific in vivo genetic toxicity of nine polycyclic aromatic hydrocarbons assessed using the Muta™Mouse transgenic rodent assay

Test batteries to screen chemicals for mutagenic hazard include several endpoints regarded as effective for detecting genotoxic carcinogens. Traditional in vivo methods primarily examine clastogenic endpoints in haematopoietic tissues. Although this approach is effective for identifying systemically distributed clastogens, some mutagens may not induce clastogenic effects; moreover, genotoxic effects may be restricted to the site of contact and/or related tissues. An OECD test guideline for transgenic rodent (TGR) gene mutation assays was released in 2011, and the TGR assays permit assessment of mutagenicity in any tissue. This study assessed the responses of two genotoxicity endpoints following sub-chronic oral exposures of male Muta™Mouse to 9 carcinogenic polycyclic aromatic hydrocarbons (PAHs). Clastogenicity was assessed via induction of micronuclei in peripheral blood, and mutagenicity via induction of lacZ transgene mutations in bone marrow, glandular stomach, small intestine, liver, and lung. Additionally, the presence of bulky PAH-DNA adducts was examined. Five of the 9 PAHs elicited positive results across all endpoints in at least one tissue, and no PAHs were negative or equivocal across all endpoints. All PAHs were positive for lacZ mutations in at least one tissue (sensitivity=100%), and for 8 PAHs, one or more initial sites of chemical contact (i.e., glandular stomach, liver, small intestine) yielded a greater response than bone marrow. Five PAHs were positive in the micronucleus assay (sensitivity=56%). Furthermore, all PAHs produced DNA adducts in at least one tissue. The results demonstrate the utility of the TGR assay for mutagenicity assessment, especially for compounds that may not be systemically distributed.

Transcriptional Profiling of Dibenzo[def,p]chrysene-induced Spleen Atrophy Provides Mechanistic Insights into its Immunotoxicity in MutaMouse

Dibenzo[def,p]chrysene (DBC) is the most carcinogenic polycyclic aromatic hydrocarbon (PAH) examined to date. We investigated the immunotoxicity of DBC, manifested as spleen atrophy, following acute exposure of adult MutaMouse males by oral gavage. Mice were exposed to 0, 2.0, 6.2, or 20.0 mg DBC /kg-bw per day, for 3 days. Genotoxic endpoints (DBC-DNA adducts and lacZ mutant frequency in spleen and bone marrow, and red blood cell micronucleus frequency) and global gene expression changes were measured. All of the genotoxicity measures increased in a dose-dependent manner in spleen and bone marrow. Gene expression analysis showed that DBC activates p53 signaling pathways related to cellular growth and proliferation, which was evident even at the low dose. Strikingly, the expression profiles of DBC exposed mouse spleens were highly inversely correlated with the expression profiles of the only published toxicogenomics dataset of enlarged mouse spleen. This analysis suggested a central role for Bnip3l, a pro-apoptotic protein involved in negative regulation of erythroid maturation. RT-PCR confirmed expression changes in several genes related to apoptosis, iron metabolism, and aryl hydrocarbon receptor signaling that are regulated in the opposite direction during spleen atrophy versus benzo[a]pyrene-mediated splenomegaly. In addition, benchmark dose modeling of toxicogenomics data yielded toxicity estimates that are very close to traditional toxicity endpoints. This work illustrates the power of toxicogenomics to reveal rich mechanistic information for immunotoxic compounds and its ability to provide information that is quantitatively similar to that derived from standard toxicity methods in health risk assessment.

HPLC Measurement of the DNA Oxidation Biomarker, 8-oxo-7,8-dihydro-2'-deoxyguanosine, in Cultured Cells and Animal Tissues

Oxidative stress is associated with many physiological and pathological processes, as well as xenobiotic metabolism, leading to the oxidation of biomacromolecules, including DNA. Therefore, efficient detection of DNA oxidation is important for a variety of research disciplines, including medicine and toxicology. A common biomarker of oxidatively damaged DNA is 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dGuo; often erroneously referred to as 8-hydroxy-2'-deoxyguanosine (8-OH-dGuo or 8-oxo-dG)). Several protocols for 8-oxo-dGuo measurement by high pressure liquid chromatography with electrochemical detection (HPLC-ED) have been described. However, these were mainly applied to purified DNA treated with pro-oxidants. In addition, due to methodological differences between laboratories, mainly due to differences in analytical equipment, the adoption of published methods for detection of 8-oxo-dGuo by HPLC-ED requires careful optimization by each laboratory. A comprehensive protocol, describing such an optimization process, is lacking. Here, a detailed protocol is described for the detection of 8-oxo-dGuo by HPLC-ED, in DNA from cultured cells or animal tissues. It illustrates how DNA sample preparation can be easily and rapidly optimized to minimize undesirable DNA oxidation that can occur during sample preparation. This protocol shows how to detect 8-oxo-dGuo in cultured human alveolar adenocarcinoma cells (i.e., A549 cells) treated with the oxidizing agent KBrO3, and from the spleen of mice exposed to the polycyclic aromatic hydrocarbon dibenzo(def,p)chrysene (DBC, formerly known as dibenzo(a,l)pyrene, DalP). Overall, this work illustrates how an HPLC-ED methodology can be readily optimized for the detection of 8-oxo-dGuo in biological samples.

Cancer risk assessment of polycyclic aromatic hydrocarbon contaminated soils determined using bioassay-derived levels of benzo[a]pyrene equivalents

Here we evaluate the excess lifetime cancer risk (ELCR) posed by 10 PAH-contaminated soils using (i) the currently advocated, targeted chemical-specific approach that assumes dose additivity for carcinogenic PAHs and (ii) a bioassay-based approach that employs the in vitro mutagenic activity of the soil fractions to determine levels of benzo[a]pyrene equivalents and, by extension, ELCR. Mutagenic activity results are presented in our companion paper.1 The results show that ELCR values for the PAH-containing fractions, determined using the chemical-specific approach, are generally (i.e., 8 out of 10) greater than those calculated using the bioassay-based approach; most are less than 5-fold greater. Only two chemical-specific ELCR estimates are less than their corresponding bioassay-derived values; differences are less than 10%. The bioassay-based approach, which permits estimation of ELCR without a priori knowledge of mixture composition, proved to be a useful tool to evaluate the chemical-specific approach. The results suggest that ELCR estimates for complex PAH mixtures determined using a targeted, chemical-specific approach are reasonable, albeit conservative. Calculated risk estimates still depend on contentious PEFs and cancer slope factors. Follow-up in vivo mutagenicity assessments will be required to validate the results and their relevance for human health risk assessment of PAH-contaminated soils.

In vitro mammalian mutagenicity of complex polycyclic aromatic hydrocarbon mixtures in contaminated soils

This study employed an in vitro version of the lacZ transgenic rodent mutation assay to assess the mutagenicity of nonpolar neutral and semipolar aromatic soil fractions from 10 PAH-contaminated sites, and evaluated the assumption of dose additivity that is routinely employed to calculate the risk posed by PAH mixtures. Significant mutagenic activity was detected in all nonpolar neutral fractions, and 8 of 10 semipolar aromatic fractions (nonpolar > semipolar). Mutagenic activity of synthetic PAH mixtures that mimic the PAH content of the soils (i.e., 5-PAH or 16-PAH mix) were greater than that of the PAH-containing soil fractions, with 5-PAH mix >16-PAH-mix. Predictions of mutagenic activity, calculated as the sum of the contributions from the mutagenic mixture components, were all within 2-fold of the observed activity of the nonpolar neutral fractions, with one exception. Observed differences in mutagenic activity are likely the result of dynamic metabolic processes, involving a complex interplay of AhR agonsim and saturation of metabolic machinery by competitive inhibition of mixture components. The presence of hitherto unidentified polar compounds present in PAH-contaminated soils may also contribute to overall hazard; however, these compounds are generally not included in current contaminated site risk assessment protocols.

Human urinary mutagenicity after wood smoke exposure during traditional temazcal use

In Central America, the traditional temazcales or wood-fired steam baths, commonly used by many Native American populations, are often heated by wood fires with little ventilation, and this use results in high wood smoke exposure. Urinary mutagenicity has been previously employed as a non-invasive biomarker of human exposure to combustion emissions. This study examined the urinary mutagenicity in 19 indigenous Mayan families from the highlands of Guatemala who regularly use temazcales (N = 32), as well as control (unexposed) individuals from the same population (N = 9). Urine samples collected before and after temazcal exposure were enzymatically deconjugated and extracted using solid-phase extraction. The creatinine-adjusted mutagenic potency of urine extracts was assessed using the plate-incorporation version of the Salmonella mutagenicity assay with strain YG1041 in the presence of exogenous metabolic activation. The post-exposure mutagenic potency of urine extracts were, on average, 1.7-fold higher than pre-exposure samples (P < 0.005) and also significantly more mutagenic than the control samples (P < 0.05). Exhaled carbon monoxide (CO) was ~10 times higher following temazcal use (P < 0.0001), and both CO level and time spent in temazcal were positively associated with urinary mutagenic potency (i.e. P < 0.0001 and P = 0.01, respectively). Thus, the wood smoke exposure associated with temazcal use contributes to increased excretion of conjugated mutagenic metabolites. Moreover, urinary mutagenic potency is correlated with other metrics of exposure (i.e. exhaled CO, duration of exposure). Since urinary mutagenicity is a biomarker associated with genetic damage, temazcal use may therefore be expected to contribute to an increased risk of DNA damage and mutation, effects associated with the initiation of cancer.

Hepatic genotoxicity and toxicogenomic responses in Muta™Mouse males treated with dibenz[a,h]anthracene

Dibenz[a,h]anthracene (DB[a,h]A) is a polycyclic aromatic hydrocarbon that is a by-product of combustion and a potent carcinogen. Few studies have investigated the effects of DB[a,h]A on mRNA and microRNA expression to dissect the mechanisms involved in carcinogenesis. In this study, mature male mice (Muta(™)Mouse) were exposed to 6.25, 12.5 and 25mg/kg/day DB[a,h]A by oral gavage for 28 consecutive days. Results were compared with mice similarly exposed to benzo[a]pyrene (B[a]P) in our previous work. Liver DNA adduct levels and lacZ mutant frequency increased dose dependently for both chemicals. Aryl hydrocarbon receptor (AhR) potency was greater for DB[a,h]A than B[a]P using the chemical-activated luciferase expression assay. Microarray analysis revealed 19 up-regulated and 22 down-regulated genes (false discovery rate-adjusted P ≤ 0.05; fold change ≥ 1.5) following treatment with 6.25 mg/kg/day DB[a,h]A. Thirteen transcripts were up-regulated and 32 down-regulated in the 12.5mg/kg/day group. The 25mg/kg/day dose had major effects on mRNA expression with 135 up-regulated and 104 down-regulated genes. Overall, perturbations were greater for DB[a,h]A than for B[a]P; in vitro chemical-activated luciferase expression supports that this may be driven by the AhR. Many of the DB[a,h]A-affected genes are implicated in cancer and are essential in vital biological functions including circadian rhythm, glucose metabolism, lipid metabolism, immune response, cell cycle and apoptosis. Although a number of functional groups were similarly affected by B[a]P and DB[a,h]A, in general the responses generated by each chemical were quite distinct. Commonalities included a DNA damage response leading to induction of cell cycle arrest and apoptosis in both Tp53-dependent and Tp53-independent manners. MicroRNA expression was identical for both chemicals, with only miR-34a showing a dose-dependent increase in treated mice.

Physical-chemical and microbiological characterization, and mutagenic activity of airborne PM sampled in a biomass-fueled electrical production facility

Biomass combustion is used in heating and electric power generation in many areas of the world. Airborne particulate matter (PM) is released when biomass is brought to a facility, stored, and combusted. Occupational exposure to airborne PM within biomass-fueled facilities may lead to health problems. In March and August of 2006, airborne PM was collected from a biomass-fueled facility located in Denmark. In addition, source-specific PM was generated from straw and wood pellets using a rotating drum. The PM was analyzed for polycyclic aromatic hydrocarbons (PAHs), metals, microbial components, mutagenic activity, and ability to generate highly reactive oxygen species (hROS) in cell-free aqueous suspensions. PM collected from the boiler room and the biomass storage hall had higher levels of mutagenic activity, PAHs and metals, and a higher hROS generating potential than the source specific PM. The mutagenic activity was generally more potent without S9 activation, and on the metabolically enhanced strain YG1041, relative to TA98. Significant correlations were found between mutagenicity on YG1041 (without S9) and PAH concentration and mutagenicity on YG1041 (with S9) and hROS generating ability. PM collected in March was more toxic than PM collected in August. Overall, airborne PM collected from the facility, especially that from the boiler room, were more toxic than PM generated from straw and wood chips. The results suggest that exposure to combustion PM in a biomass-fueled facility, which likely includes PM from biomass combustion as well as internal combustion vehicles, may contribute to an elevated risk of adverse health effects.