Development of a toxicogenomics signature for genotoxicity using a dose-optimization and informatics strategy in human cells

Assessment of TGx-DDI genes for genotoxicity in a comprehensive panel of chemicals

BACKGROUND

The TGx-DDI biomarker identifies transcripts specifically induced by primary DNA damage. Profiling similarity of TGx-DDI signatures can allow clustering compounds by genotoxic mechanism. This transcriptomics-based approach complements conventional toxicology testing by enhancing mechanistic resolution.

METHODS

Unsupervised hierarchical clustering and t-distributed stochastic neighbor embedding (tSNE) were utilized to assess similarity of publicly-available per- and polyfluoroalkyl substances (PFAS) and ToxCast chemicals based on TGx-DDI modulation. TempO-seq transcriptomic data after highest chemical concentrations were analyzed.

RESULTS

Clustering discriminated between genotoxic and non-genotoxic compounds while drawing similarity among chemicals with shared mechanisms. PFAS largely clustered distinctly from classical mutagens. However, dynamic range across PFAS types and durations indicated variable potential for DNA damage. tSNE visualization reinforced phenotypic groupings, with genotoxins clustering separately from non-DNA damaging agents.

DISCUSSION

Unsupervised learning approaches applied to TGx-DDI profiles effectively categorizes chemical genotoxicity potential, aiding elucidation of biological response pathways. This transcriptomics-based strategy gives further insight into the role and effect of individual TGx-DDI biomarker genes and complements existing assays by enhancing mechanistic resolution. Overall, TGx-DDI biomarker profiling holds promise for predictive safety screening.

Non-clinical safety assessment of Annona atemoya leaf extract: evaluation of genotoxicity

The leaves, stems, and fruits of Annona atemoya (A. atemoya; AA), a fruit-bearing plant of the family Annonaceae, exhibit anti-angiogenic, anti-oxidative, anti-inflammatory, and neuroprotective activities. However, the safety of AA has not been comprehensively elucidated. In this study, we evaluated the potential genotoxicity of an AA leaf (AAL) ethanol extract using a standard three-test battery constituting in vitro mammalian chromosomal aberration, in vivo micronucleus, and bacterial reverse mutation (also known as the Ames test) tests, as recommended by the Ministry of Food and Drug Safety of Korea. In vitro chromosomal aberration assay revealed that AAL extract did not induce structural or numerical aberrations, with or without metabolic activation (S9). In vivo micronucleus assay revealed that the number of micronucleated polychromatic erythrocytes (PCEs) and the PCE/normochromatic erythrocyte ratio after AAL extract treatment were not substantially different from those in the negative control. Changes in body weight and mortality were not observed. However, AAL extract partially induced mutagenic activity in all three bacterial strains in the bacterial reverse mutation assay, indicating that it could potentially aid in determining the genotoxic safety of AAL. QuantSeq 3' mRNA sequencing analysis to elucidate the genotoxicity mechanisms of AAL extract using TK6 cells revealed that the genotoxic effects of AAL may be associated with cellular morphology-associated (cell development and keratinization), nucleotide metabolism, and electron transport chain functions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43188-024-00241-4.

FAILLA MEMORIAL LECTURE How We Got Here: One Laboratory's Odyssey in the Field of Radiation-Inducible Genes

This review focuses on early discoveries that contributed to our understanding and the scope of transcriptional responses after radiation damage. Before the development of modern approaches to assess overall global transcriptomic responses, the idea that mammalian cells could respond to DNA-damaging agents in a manner analogous to bacteria was not generally accepted. To investigate this possibility, the development of technology to identify differentially expressed low-abundance transcripts substantially facilitated our appreciation that DNA damaging agents like UV radiation and subsequently ionizing radiation did in fact produce robust transcriptional responses. Here we focus on our identification and characterization of radiation-inducible genes, and how even early studies on stress gene signaling highlighted the broad scope of transcriptional responses to radiation damage. Since then, the central role of transcriptional responses to radiation injury in maintaining genome integrity has been highlighted in many processes, including cell cycle checkpoint control, resistance to cancer by p53 and other key factors, cell senescence, and metabolism.

Application of HepaRG cells for genotoxicity assessment: a review

There has been growing interest in the use of human-derived metabolically competent cells for genotoxicity testing. The HepaRG cell line is considered one of the most promising cell models because it is TP53-proficient and retains many characteristics of primary human hepatocytes. In recent years, HepaRG cells, cultured in both a traditional two-dimensional (2D) format and as more advanced in-vivo-like 3D spheroids, have been employed in assays that measure different types of genetic toxicity endpoints, including DNA damage, mutations, and chromosomal damage. This review summarizes published studies that have used HepaRG cells for genotoxicity assessment, including cell model evaluation studies and risk assessment for various compounds. Both 2D and 3D HepaRG models can be adapted to several high-throughput genotoxicity assays, generating a large number of data points that facilitate quantitative benchmark concentration modeling. With further validation, HepaRG cells could serve as a unique, human-based new alternative methodology for in vitro genotoxicity testing.

Unlocking the Power of Transcriptomic Biomarkers in Qualitative and Quantitative Genotoxicity Assessment of Chemicals

To modernize genotoxicity assessment and reduce reliance on experimental animals, new approach methodologies (NAMs) that provide human-relevant dose-response data are needed. Two transcriptomic biomarkers, GENOMARK and TGx-DDI, have shown a high classification accuracy for genotoxicity. As these biomarkers were extracted from different training sets, we investigated whether combining the two biomarkers in a human-derived metabolically competent cell line (i.e., HepaRG) provides complementary information for the classification of genotoxic hazard identification and potency ranking. First, the applicability of GENOMARK to TempO-Seq, a high-throughput transcriptomic technology, was evaluated. HepaRG cells were exposed for 72 h to increasing concentrations of 10 chemicals (i.e., eight known in vivo genotoxicants and two in vivo nongenotoxicants). Gene expression data were generated using the TempO-Seq technology. We found a prediction performance of 100%, confirming the applicability of GENOMARK to TempO-Seq. Classification using TGx-DDI was then compared to GENOMARK. For the chemicals identified as genotoxic, benchmark concentration modeling was conducted to perform potency ranking. The high concordance observed for both hazard classification and potency ranking by GENOMARK and TGx-DDI highlights the value of integrating these NAMs in a weight of evidence evaluation of genotoxicity.

Evaluation of weak genotoxicity of hydroxychloroquine in human TK6 cells

Hydroxychloroquine (HCQ), a derivative of chloroquine (CQ), is an antimalarial and antirheumatic drug. Since there is limited data available on the genotoxicity of HCQ, in the current study, we used a battery of in vitro assays to systematically examine the genotoxicity of HCQ in human lymphoblastoid TK6 cells. We first showed that HCQ is not mutagenic in TK6 cells up to 80 μM with or without exogenous metabolic activation. Subsequently, we found that short-term (3-4 h) HCQ treatment did not cause DNA strand breakage as measured by the comet assay and the phosphorylation of histone H2A.X (γH2A.X), and did not induce chromosomal damage as determined by the micronucleus (MN) assay. However, after 24-h treatment, both CQ and HCQ induced comparable and weak DNA damage and MN formation in TK6 cells; upregulated p53 and p53-mediated DNA damage responsive genes; and triggered apoptosis and mitochondrial damage that may partially contribute to the observed MN formation. Using a benchmark dose (BMD) modeling analysis, the lower 95% confidence limit of BMD50 values (BMDL50) for MN induction in TK6 cells were about 19.7 μM for CQ and 16.3 μM for HCQ. These results provide additional information for quantitative genotoxic risk assessment of these drugs.

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.

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.

Integrated Genotoxicity Testing of three anti-infective drugs using the TGx-DDI transcriptomic biomarker and high-throughput CometChip® assay in TK6 cells

Genotoxicity testing relies on the detection of gene mutations and chromosome damage and has been used in the genetic safety assessment of drugs and chemicals for decades. However, the results of standard genotoxicity tests are often difficult to interpret due to lack of mode of action information. The TGx-DDI transcriptomic biomarker provides mechanistic information on the DNA damage-inducing (DDI) capability of chemicals to aid in the interpretation of positive in vitro genotoxicity data. The CometChip^® assay was developed to assess DNA strand breaks in a higher-throughput format. We paired the TGx-DDI biomarker with the CometChip^® assay in TK6 cells to evaluate three model agents: nitrofurantoin (NIT), metronidazole (MTZ), and novobiocin (NOV). TGx-DDI was analyzed by two independent labs and technologies (nCounter^® and TempO-Seq^®). Although these anti-infective drugs are, or have been, used in human and/or veterinary medicine, the standard genotoxicity testing battery showed significant genetic safety findings. Specifically, NIT is a mutagen and causes chromosome damage, and MTZ and NOV cause chromosome damage in conventional in vitro tests. Herein, the TGx-DDI biomarker classified NIT and MTZ as non-DDI at all concentrations tested, suggesting that NIT’s mutagenic activity is bacterial specific and that the observed chromosome damage by MTZ might be a consequence of in vitro test conditions. In contrast, NOV was classified as DDI at the second highest concentration tested, which is in line with the fact that NOV is a bacterial DNA-gyrase inhibitor that also affects topoisomerase II at high concentrations. The lack of DNA damage for NIT and MTZ was confirmed by the CometChip^® results, which were negative for all three drugs except at overtly cytotoxic concentrations. This case study demonstrates the utility of combining the TGx-DDI biomarker and CometChip^® to resolve conflicting genotoxicity data and provides further validation to support the reproducibility of the biomarker.

Search for the optimal genotoxicity assay for routine testing of chemicals: Sensitivity and specificity of conventional and new test systems

Many conventional in vitro tests that are currently widely used for routine screening of chemicals have a sensitivity/specificity in the range between 60 % and 80 % for the detection of carcinogens. Most procedures were developed 30-40 years ago. In the last decades several assays became available which are based on the use of metabolically competent cell lines, improvement of the cultivation conditions and development of new endpoints. Validation studies indicate that some of these models may be more reliable for the detection of genotoxicants (i.e. many of them have sensitivity and specificity values between 80 % and 95 %). Therefore, they could replace conventional tests in the future. The bone marrow micronucleus (MN) assay with rodents is at present the most widely used in vivo test. The majority of studies indicate that it detects only 5-6 out of 10 carcinogens while experiments with transgenic rodents and comet assays seem to have a higher predictive value and detect genotoxic carcinogens that are negative in MN experiments. Alternatives to rodent experiments could be MN experiments with hen eggs or their replacement by combinations of new in vitro tests. Examples for promising candidates are ToxTracker, TGx-DDI, multiplex flow cytometry, γH2AX experiments, measurement of p53 activation and MN experiments with metabolically competent human derived liver cells. However, the realization of multicentric collaborative validation studies is mandatory to identify the most reliable tests.

Genomic analysis of Fisher F344 rat kidneys from a reproductive study following dietary ochratoxin A exposure

Ochratoxin A (OTA) is a mycotoxin produced by species of Penicillium and Aspergillus, and is found in many commodities including cereal grains, nuts, and coffee. OTA is a renal carcinogen and nephrotoxin at high concentrations, targeting the proximal tubules. This study uses transcriptomics and the previously reported apical data (Bondy et al., 2021) to infer mode-of-action of OTA toxicity in male and female rats exposed to low doses of OTA in utero and throughout development. Our findings support a male-specific activation of the innate and adaptive immune responses in F1 pups to OTA exposure. This was not found in the female F1 pups, and may be due to female-specific increased p38 activity and VDR signaling. Differentially expressed genes related to karyomegaly, MAPK activity, and immune activation appears to develop from in utero exposure to OTA whereas those related to decreased kidney and liver function, and changes to reproductive pathways occur in both rat generations. Together, these transcriptional results confirm that dietary exposure to OTA causes renal toxicity as well as alterations to hepatic and reproductive pathways in rats. In utero exposure of rats to OTA results in sex-specific alterations in immune response pathways, VDR signaling, and p38 activity.

Comprehensive interpretation of in vitro micronucleus test results for 292 chemicals: from hazard identification to risk assessment application

Risk assessments are increasingly reliant on information from in vitro assays. The in vitro micronucleus test (MNvit) is a genotoxicity test that detects chromosomal abnormalities, including chromosome breakage (clastogenicity) and/or whole chromosome loss (aneugenicity). In this study, MNvit datasets for 292 chemicals, generated by the US EPA's ToxCast program, were evaluated using a decision tree-based pipeline for hazard identification. Chemicals were tested with 19 concentrations (n = 1) up to 200 µM, in the presence and absence of Aroclor 1254-induced rat liver S9. To identify clastogenic chemicals, %MN values at each concentration were compared to a distribution of batch-specific solvent controls; this was followed by cytotoxicity assessment and benchmark concentration (BMC) analyses. The approach classified 157 substances as positives, 25 as negatives, and 110 as inconclusive. Using the approach described in Bryce et al. (Environ Mol Mutagen 52:280-286, 2011), we identified 15 (5%) aneugens. IVIVE (in vitro to in vivo extrapolation) was employed to convert BMCs into administered equivalent doses (AEDs). Where possible, AEDs were compared to points of departure (PODs) for traditional genotoxicity endpoints; AEDs were generally lower than PODs based on in vivo endpoints. To facilitate interpretation of in vitro MN assay concentration-response data for risk assessment, exposure estimates were utilized to calculate bioactivity exposure ratio (BER) values. BERs for 50 clastogens and two aneugens had AEDs that approached exposure estimates (i.e., BER < 100); these chemicals might be considered priorities for additional testing. This work provides a framework for the use of high-throughput in vitro genotoxicity testing for priority setting and chemical risk assessment.

A Collaborative Initiative to Establish Genomic Biomarkers for Assessing Tumorigenic Potential to Reduce Reliance on Conventional Rodent Carcinogenicity Studies

There is growing recognition across broad sectors of the scientific community that use of genomic biomarkers has the potential to reduce the need for conventional rodent carcinogenicity studies of industrial chemicals, agrochemicals, and pharmaceuticals through a weight-of-evidence approach. These biomarkers fall into 2 major categories: (1) sets of gene transcripts that can identify distinct tumorigenic mechanisms of action; and (2) cancer driver gene mutations indicative of rapidly expanding growth-advantaged clonal cell populations. This call-to-action article describes a collaborative approach launched to develop and qualify biomarker gene expression panels that measure widely accepted molecular pathways linked to tumorigenesis and their activation levels to predict tumorigenic doses of chemicals from short-term exposures. Growing evidence suggests that application of such biomarker panels in short-term exposure rodent studies can identify both tumorigenic hazard and tumorigenic activation levels for chemical-induced carcinogenicity. In the future, this approach will be expanded to include methodologies examining mutations in key cancer driver gene mutation hotspots as biomarkers of both genotoxic and nongenotoxic chemical tumor risk. Analytical, technical, and biological validation studies of these complementary genomic tools are being undertaken by multisector and multidisciplinary collaborative teams within the Health and Environmental Sciences Institute. Success from these efforts will facilitate the transition from current heavy reliance on conventional 2-year rodent carcinogenicity studies to more rapid animal- and resource-sparing approaches for mechanism-based carcinogenicity evaluation supporting internal and regulatory decision-making.

A gene expression biomarker identifies inhibitors of two classes of epigenome effectors in a human microarray compendium

Biomarkers predictive of molecular and toxicological effects are needed to interpret emerging high-throughput transcriptomics (HTTr) data streams. To address the limited approaches available for identifying epigenotoxicants, we previously developed and validated an 81-gene biomarker that accurately predicts histone deacetylase inhibition (HDACi) in transcript profiles derived from chemically-treated TK6 cells. In the present study, we sought to determine if this biomarker (TGx-HDACi) could be used to identify HDACi chemicals in other cell lines using the Running Fisher correlation test. Using microarray comparisons derived from human cells exposed to HDACi, we found considerable heterogeneity in correlation with the TGx-HDACi biomarker dependent on chemical exposure conditions and tissue from which the cell line was derived. Using a defined set of conditions that overlapped with our earlier study, the biomarker was able to accurately identify HDACi chemicals (90-100% balanced accuracy). In an in silico screen of 2427 chemicals in 9660 chemical versus control comparisons, the biomarker coupled with the Running Fisher test was able to identify 14 additional HDACi chemicals as well as other chemicals not previously associated with HDACi. Most notable were 12 inhibitors of bromodomain (BRD) and extraterminal (BET) family proteins including BRD4 that bind to acetylated histones. The BET protein inhibitors could be distinguished from the HDACi based on differences in the expression of a small set of biomarker genes. Our results indicate that the TGx-HDACi biomarker will be useful for identifying inhibitors of two classes of epigenome effectors in HTTr screening studies.

Use of transcriptomics in hazard identification and next generation risk assessment: A case study with clothianidin

Toxicological risk assessment is essential in the evaluation and authorization of different classes of chemical substances. Genotoxicity and mutagenicity testing are of highest priority and rely on established in vitro systems with bacterial and mammalian cells, sometimes followed by in vivo testing using rodent animal models. Transcriptomic approaches have recently also shown their value to determine transcript signatures specific for genotoxicity. Here, we studied how transcriptomic data, in combination with in vitro tests with human cells, can be used for the identification of genotoxic properties of test compounds. To this end, we used liver samples from a 28-day oral toxicity study in rats with the pesticidal active substances imazalil, thiacloprid, and clothianidin, a neonicotinoid-type insecticide with, amongst others, known hepatotoxic properties. Transcriptomic results were bioinformatically evaluated and pointed towards a genotoxic potential of clothianidin. In vitro Comet and γH2AX assays in human HepaRG hepatoma cells, complemented by in silico analyses of mutagenicity, were conducted as follow-up experiments to check if the genotoxicity alert from the transcriptomic study is in line with results from a battery of guideline genotoxicity studies. Our results illustrate the combined use of toxicogenomics, classic toxicological data and new approach methods in risk assessment. By means of a weight-of-evidence decision, we conclude that clothianidin does most likely not pose genotoxic risks to humans.

•

Analysis of clothianidin genotoxicity in silico, in vitro and in vivo.

•

Application of a toxicogenomics approach to analyze genotoxicity.

•

Weight-of-evidence decision supports classification as “non-genotoxic”.

Current Therapeutic Landscape and Safety Roadmap for Targeting the Aryl Hydrocarbon Receptor in Inflammatory Gastrointestinal Indications

Target modulation of the AhR for inflammatory gastrointestinal (GI) conditions holds great promise but also the potential for safety liabilities both within and beyond the GI tract. The ubiquitous expression of the AhR across mammalian tissues coupled with its role in diverse signaling pathways makes development of a “clean” AhR therapeutically challenging. Ligand promiscuity and diversity in context-specific AhR activation further complicates targeting the AhR for drug development due to limitations surrounding clinical translatability. Despite these concerns, several approaches to target the AhR have been explored such as small molecules, microbials, PROTACs, and oligonucleotide-based approaches. These various chemical modalities are not without safety liabilities and require unique de-risking strategies to parse out toxicities. Collectively, these programs can benefit from in silico and in vitro methodologies that investigate specific AhR pathway activation and have the potential to implement thresholding parameters to categorize AhR ligands as “high” or “low” risk for sustained AhR activation. Exploration into transcriptomic signatures for AhR safety assessment, incorporation of physiologically-relevant in vitro model systems, and investigation into chronic activation of the AhR by structurally diverse ligands will help address gaps in our understanding regarding AhR-dependent toxicities. Here, we review the role of the AhR within the GI tract, novel therapeutic modality approaches to target the AhR, key AhR-dependent safety liabilities, and relevant strategies that can be implemented to address drug safety concerns. Together, this review discusses the emerging therapeutic landscape of modalities targeting the AhR for inflammatory GI indications and offers a safety roadmap for AhR drug development.

GADD45 in Stress Signaling, Cell Cycle Control, and Apoptosis

GADD45 is a gene family consisting of GADD45A, GADD45B, and GADD45G that is often induced by DNA damage and other stress signals associated with growth arrest and apoptosis. Many of these roles are carried out via signaling mediated by p38 mitogen-activated protein kinases (MAPKs). The GADD45 proteins can contribute to p38 activation either by activation of upstream kinase(s) or by direct interaction, as well as suppression of p38 activity in certain cases. In vivo, there are important tissue and cell type specific differences in the roles for GADD45 in MAPK signaling. In addition to being p53-regulated, GADD45A has also been found to contribute to p53 activation via p38. Like other stress and signaling proteins, GADD45 proteins show complex regulation and numerous effectors. More recently, aberrant GADD45 expression has been found in several human cancers, but the mechanisms behind these findings largely remain to be understood.

Comparative Analysis of Transcriptional Responses to Genotoxic and Non-Genotoxic Agents in the Blood Cell Model TK6 and the Liver Model HepaRG

Transcript signatures are a promising approach to identify and classify genotoxic and non-genotoxic compounds and are of interest as biomarkers or for future regulatory application. Not much data, however, is yet available about the concordance of transcriptional responses in different cell types or tissues. Here, we analyzed transcriptomic responses to selected genotoxic food contaminants in the human p53-competent lymphoblastoid cell line TK6 using RNA sequencing. Responses to treatment with five genotoxins, as well as with four non-genotoxic liver toxicants, were compared with previously published gene expression data from the human liver cell model HepaRG. A significant overlap of the transcriptomic changes upon genotoxic stress was detectable in TK6 cells, whereas the comparison with the HepaRG model revealed considerable differences, which was confirmed by bioinformatic data mining for cellular upstream regulators or pathways. Taken together, the study presents a transcriptomic signature for genotoxin exposure in the human TK6 blood cell model. The data demonstrate that responses in different cell models have considerable variations. Detection of a transcriptomic genotoxin signature in blood cells indicates that gene expression analyses of blood samples might be a valuable approach to also estimate responses to toxic exposure in target organs such as the liver.

Implementing in vitro bioactivity data to modernize priority setting of chemical inventories

Internationally, there are thousands of existing and newly introduced chemicals in commerce, highlighting the ongoing importance of innovative approaches to identify emerging chemicals of concern. For many chemicals, there is a paucity of hazard and exposure data. Thus, there is a crucial need for efficient and robust approaches to address data gaps and support risk-based prioritization. Several studies have demonstrated the utility of in vitro bioactivity data from the ToxCast program in deriving points of departure (PODs). ToxCast contains data for nearly 1,400 endpoints per chemical, and the bioactivity concentrations, indicative of potential adverse outcomes, can be converted to human-equivalent PODs using high-throughput toxicokinetics (HTTK) modeling. However, data gaps need to be addressed for broader application: the limited chemical space of HTTK and quantitative high-throughput screening data. Here we explore the applicability of in silico models to address these data needs. Specifically, we used ADMET predictor for HTTK predictions and a generalized read-across approach to predict ToxCast bioactivity potency. We applied these models to profile 5,801 chemicals on Canada’s Domestic Substances List (DSL). To evaluate the approach’s performance, bioactivity PODs were compared with in vivo results from the EPA Toxicity Values database for 1,042 DSL chemicals. Comparisons demonstrated that the bioac­tivity PODs, based on ToxCast data or read-across, were conservative for 95% of the chemicals. Comparing bioactivity PODs to human exposure estimates supports the identification of chemicals of potential interest for further work. The bioac­tivity workflow shows promise as a powerful screening tool to support effective triaging of chemical inventories.

CEBS update: curated toxicology database with enhanced tools for data integration

The Chemical Effects in Biological Systems database (CEBS) contains extensive toxicology study results and metadata from the Division of the National Toxicology Program (NTP) and other studies of environmental health interest. This resource grants public access to search and collate data from over 10 250 studies for 12 750 test articles (chemicals, environmental agents). CEBS has made considerable strides over the last 5 years to integrate growing internal data repositories into data warehouses and data marts to better serve the public with high quality curated datasets. This effort includes harmonizing legacy terms and metadata to current standards, mapping test articles to external identifiers, and aligning terms to OBO (Open Biological and Biomedical Ontology) Foundry ontologies. The data are made available through the CEBS Homepage (https://cebs.niehs.nih.gov/cebs/), guided search applications, flat files on FTP (file transfer protocol), and APIs (application programming interface) for user access and to provide a bridge for computational tools. The user interface is intuitive with a single search bar to query keywords related to study metadata, publications, and data availability. Results are consolidated to single pages for each test article with NTP conclusions, publications, individual studies, data collections, and links to related test articles and projects available together.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

A Modern Genotoxicity Testing Paradigm: Integration of the High-Throughput CometChip® and the TGx-DDI Transcriptomic Biomarker in Human HepaRG™ Cell Cultures

Higher-throughput, mode-of-action-based assays provide a valuable approach to expedite chemical evaluation for human health risk assessment. In this study, we combined the high-throughput alkaline DNA damage-sensing CometChip^® assay with the TGx-DDI transcriptomic biomarker (DDI = DNA damage-inducing) using high-throughput TempO-Seq^®, as an integrated genotoxicity testing approach. We used metabolically competent differentiated human HepaRG™ cell cultures to enable the identification of chemicals that require bioactivation to cause genotoxicity. We studied 12 chemicals (nine DDI, three non-DDI) in increasing concentrations to measure and classify chemicals based on their ability to damage DNA. The CometChip^® classified 10/12 test chemicals correctly, missing a positive DDI call for aflatoxin B1 and propyl gallate. The poor detection of aflatoxin B1 adducts is consistent with the insensitivity of the standard alkaline comet assay to bulky lesions (a shortcoming that can be overcome by trapping repair intermediates). The TGx-DDI biomarker accurately classified 10/12 agents. TGx-DDI correctly identified aflatoxin B1 as DDI, demonstrating efficacy for combined used of these complementary methodologies. Zidovudine, a known DDI chemical, was misclassified as it inhibits transcription, which prevents measurable changes in gene expression. Eugenol, a non-DDI chemical known to render misleading positive results at high concentrations, was classified as DDI at the highest concentration tested. When combined, the CometChip^® assay and the TGx-DDI biomarker were 100% accurate in identifying chemicals that induce DNA damage. Quantitative benchmark concentration (BMC) modeling was applied to evaluate chemical potencies for both assays. The BMCs for the CometChip^® assay and the TGx-DDI biomarker were highly concordant (within 4-fold) and resulted in identical potency rankings. These results demonstrate that these two assays can be integrated for efficient identification and potency ranking of DNA damaging agents in HepaRG™ cell cultures.

Comparing the Predictivity of Human Placental Gene, microRNA, and CpG Methylation Signatures in Relation to Perinatal Outcomes

Molecular signatures are being increasingly integrated into predictive biology applications. However, there are limited studies comparing the overall predictivity of transcriptomic vs. epigenomic signatures in relation to perinatal outcomes. This study set out to evaluate mRNA and microRNA (miRNA) expression and cytosine-guanine dinucleotide (CpG) methylation signatures in human placental tissues and relate these to perinatal outcomes known to influence maternal/fetal health; namely, birth weight, placenta weight, placental damage, and placental inflammation. The following hypotheses were tested: (1) different molecular signatures will demonstrate varying levels of predictivity towards perinatal outcomes, and (2) these signatures will show disruptions from an example exposure (i.e., cadmium) known to elicit perinatal toxicity. Multi-omic placental profiles from 390 infants in the Extremely Low Gestational Age Newborns cohort were used to develop molecular signatures that predict each perinatal outcome. Epigenomic signatures (i.e., miRNA and CpG methylation) consistently demonstrated the highest levels of predictivity, with model performance metrics including R^2 (predicted vs. observed) values of 0.36-0.57 for continuous outcomes and balanced accuracy values of 0.49-0.77 for categorical outcomes. Top-ranking predictors included miRNAs involved in injury and inflammation. To demonstrate the utility of these predictive signatures in screening of potentially harmful exogenous insults, top-ranking miRNA predictors were analyzed in a separate pregnancy cohort and related to cadmium. Key predictive miRNAs demonstrated altered expression in association with cadmium exposure, including miR-210, known to impact placental cell growth, blood vessel development, and fetal weight. These findings inform future predictive biology applications, where additional benefit will be gained by including epigenetic markers.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

Development and validation of the TGx-HDACi transcriptomic biomarker to detect histone deacetylase inhibitors in human TK6 cells

Transcriptomic biomarkers can be used to inform molecular initiating and key events involved in a toxicant's mode of action. To address the limited approaches available for identifying epigenotoxicants, we developed and assessed a transcriptomic biomarker of histone deacetylase inhibition (HDACi). First, we assembled a set of ten prototypical HDACi and ten non-HDACi reference compounds. Concentration-response experiments were performed for each chemical to collect TK6 human lymphoblastoid cell samples after 4 h of exposure and to assess cell viability following a 20-h recovery period in fresh media. One concentration was selected for each chemical for whole transcriptome profiling and transcriptomic signature derivation, based on cell viability at the 24-h time point and on maximal induction of HDACi-response genes (RGL1, NEU1, GPR183) or cellular stress-response genes (ATF3, CDKN1A, GADD45A) analyzed by TaqMan qPCR assays after 4 h of exposure. Whole transcriptomes were profiled after 4 h exposures by Templated Oligo-Sequencing (TempO-Seq). By applying the nearest shrunken centroid (NSC) method to the whole transcriptome profiles of the reference compounds, we derived an 81-gene toxicogenomic (TGx) signature, referred to as TGx-HDACi, that classified all 20 reference compounds correctly using NSC classification and the Running Fisher test. An additional 4 HDACi and 7 non-HDACi were profiled and analyzed using TGx-HDACi to further assess classification performance; the biomarker accurately classified all 11 compounds, including 3 non-HDACi epigenotoxicants, suggesting a promising specificity toward HDACi. The availability of TGx-HDACi increases the diversity of tools that can facilitate mode of action analysis of toxicants using gene expression profiling.

A comparison of classical and 21st century genotoxicity tools: A proof of concept study of 18 chemicals comparing in vitro micronucleus, ToxTracker and genomics-based methods (TGx-DDI, whole genome clustering and connectivity mapping)

A key step in the risk assessment process of a substance is the assessment of its genotoxic potential. Irrespective of the industry involved, current approaches rely on combinations of two or three in vitro tests and while highly sensitive, their specificity is thought to be limited. A refined in vitro genotoxicity testing strategy with improved predictive capacity would be beneficial and "3R" friendly as it helps to avoid unnecessary in vivo follow-up testing. Here, we describe a proof of concept study evaluating a balanced set of compounds that have in vivo negative or positive outcomes, but variable in vitro data, to determine if we could differentiate between direct and indirect acting genotoxicants. Compounds were examined in TK6 cells using an approach in which the same sample was used to evaluate both early genomic markers (Affymetrix analysis 4 hr post treatment), and the genotoxic outcome (micronuclei [MN] after 24 hr). The resulting genomic data was then analyzed using the TGx-DDI biomarker, Connectivity mapping and whole genome clustering. Chemicals were also tested in the ToxTracker assay, which uses six different biomarker genes. None of the methods correctly differentiated all direct from indirect acting genotoxicants when used alone, however, the ToxTracker assay, TGx-DDI biomarker and whole genome approaches provided high predictive capacity when used in combination with the MN assay (1/18, 2/18, 1/18 missed calls). Ultimately, a "fit for purpose" combination will depend on the specific tools available to the end user, as well as considerations of the unique benefits of the individual assays.

Image analysis of mechanistic protein biomarkers for the characterization of genotoxicants: Aneugens, clastogens, and reactive oxygen species inducers

The early detection of genotoxicity contributes to cutting-edge drug discovery and development, requiring effective identification of genotoxic hazards posed by drugs while providing mode of action (MoA) information in a high throughput manner. In other words, there is a need to complement standard genotoxicity testing according to the test battery given in ICH S2(R1) with new in vitro tools, thereby contributing to a more in-depth analysis of genotoxic effects. Here, we report on a proof-of-concept MoA approach based on post-translational modifications of proteins (PTMs) indicative of clastogenic and aneugenic effects in TK6 cells using imaging technology (with automated analysis). Cells were exposed in a 96-well plate format with a panel of reference (geno)toxic compounds and subsequently analyzed at 4 and 24 hr to detect dose-dependent changes in PTMs, relevant for mechanistic analysis. All tested compounds that interfere with the spindle apparatus yielded a BubR1 (S640) (3/3) and phospho-histone H3 (S28) (7/9) positive dose-response reflecting aneugenicity, whereas compounds inducing DNA double-strand-breaks were associated with positive FANCD2 (S1404) and 53BP1 (S1778) responses pointing to clastogenicity (2/3). The biomarker p53 (K373) was able to distinguish genotoxicants from non-genotoxicants (2/4), while the induction of reactive oxygen species (ROS), potentially causing DNA damage, was associated with a positive Nrf2 (S40) response (2/2). This work demonstrates that genotoxicants and non-genotoxicants induce different biomarker responses in TK6 cells which can be used for reliable classification into MoA groups (aneugens/clastogens/non-genotoxicants/ROS inducers), supporting a more in-depth safety assessment of drug candidates.

3Rs-friendly approach to exogenous metabolic activation that supports high-throughput genetic toxicology testing

MultiFlow® DNA Damage-p53, γH2AX, Phospho-Histone H3 is a miniaturized, flow cytometry-based assay that provides genotoxic mode of action information by distinguishing clastogens, aneugens, and nongenotoxicants. Work to date has focused on the p53-competent human cell line TK6. While mammalian cell genotoxicity assays typically supply exogenous metabolic activation in the form of concentrated rat liver S9, this is a less-than-ideal approach for several reasons, including 3Rs considerations. Here, we describe our experiences with low concentration S9 and saturating co-factors which were allowed to remain in contact with cells and test chemicals for 24 continuous hours. We exposed TK6 cells in 96-well plates to each of 15 reference chemicals over a range of concentrations, both in the presence and absence of 0.25% v/v phenobarbital/β-naphthoflavone-induced rat liver S9. After 4 and 24 hr of treatment cell aliquots were added to wells of a microtiter plate containing the working detergent/stain/antibody cocktail. After a brief incubation robotic sampling was employed for walk-away flow cytometric data acquisition. PROAST benchmark dose (BMD) modeling was used to characterize the resulting dose-response curves. For each of the 8 reference pro-genotoxicants studied, relative nuclei count, γH2AX, and/or p53 biomarker BMD values were order(s) of magnitude lower for 0.25% S9 conditions compared to 0% S9. Conversely, several of the direct-acting reference chemicals exhibited appreciably lower cytotoxicity and/or genotoxicity BMD values in the presence of S9 (eg, resorcinol). These results prove the efficacy of the low concentration S9 system, and indicate that an efficient and highly scalable multiplexed assay can effectively identify chemicals that require bioactivation to exert their genotoxic effects.

The promise of toxicogenomics for genetic toxicology: past, present and future

Toxicogenomics, the application of genomics to toxicology, was described as 'a new era' for toxicology. Standard toxicity tests typically involve a number of short-term bioassays that are costly, time consuming, require large numbers of animals and generally focus on a single end point. Toxicogenomics was heralded as a way to improve the efficiency of toxicity testing by assessing gene regulation across the genome, allowing rapid classification of compounds based on characteristic expression profiles. Gene expression microarrays could measure and characterise genome-wide gene expression changes in a single study and while transcriptomic profiles that can discriminate between genotoxic and non-genotoxic carcinogens have been identified, challenges with the approach limited its application. As such, toxicogenomics did not transform the field of genetic toxicology in the way it was predicted. More recently, next generation sequencing (NGS) technologies have revolutionised genomics owing to the fact that hundreds of billions of base pairs can be sequenced simultaneously cheaper and quicker than traditional Sanger methods. In relation to genetic toxicology, and thousands of cancer genomes have been sequenced with single-base substitution mutational signatures identified, and mutation signatures have been identified following treatment of cells with known or suspected environmental carcinogens. RNAseq has been applied to detect transcriptional changes following treatment with genotoxins; modified RNAseq protocols have been developed to identify adducts in the genome and Duplex sequencing is an example of a technique that has recently been developed to accurately detect mutation. Machine learning, including MutationSeq and SomaticSeq, has also been applied to somatic mutation detection and improvements in automation and/or the application of machine learning algorithms may allow high-throughput mutation sequencing in the future. This review will discuss the initial promise of transcriptomics for genetic toxicology, and how the development of NGS technologies and new machine learning algorithms may finally realise that promise.

Flow cytometric micronucleus assay and TGx-DDI transcriptomic biomarker analysis of ten genotoxic and non-genotoxic chemicals in human HepaRG™ cells

Background

Modern testing paradigms seek to apply human-relevant cell culture models and integrate data from multiple test systems to accurately inform potential hazards and modes of action for chemical toxicology. In genetic toxicology, the use of metabolically competent human hepatocyte cell culture models provides clear advantages over other more commonly used cell lines that require the use of external metabolic activation systems, such as rat liver S9. HepaRG™ cells are metabolically competent cells that express Phase I and II metabolic enzymes and differentiate into mature hepatocyte-like cells, making them ideal for toxicity testing. We assessed the performance of the flow cytometry in vitro micronucleus (MN) test and the TGx-DDI transcriptomic biomarker to detect DNA damage-inducing (DDI) chemicals in human HepaRG™ cells after a 3-day repeat exposure. The biomarker, developed for use in human TK6 cells, is a panel of 64 genes that accurately classifies chemicals as DDI or non-DDI. Herein, the TGx-DDI biomarker was analyzed by Ion AmpliSeq whole transcriptome sequencing to assess its classification accuracy using this more modern gene expression technology as a secondary objective.

Methods

HepaRG™ cells were exposed to increasing concentrations of 10 test chemicals (six genotoxic chemicals, including one aneugen, and four non-genotoxic chemicals). Cytotoxicity and genotoxicity were measured using the In Vitro MicroFlow® kit, which was run in parallel with the TGx-DDI biomarker.

Results

A concentration-related decrease in relative survival and a concomitant increase in MN frequency were observed for genotoxic chemicals in HepaRG™ cells. All five DDI and five non-DDI agents were correctly classified (as genotoxic/non-genotoxic and DDI/non-DDI) by pairing the test methods. The aneugenic agent (colchicine) yielded the expected positive result in the MN test and negative (non-DDI) result by TGx-DDI.

Conclusions

This next generation genotoxicity testing strategy is aligned with the paradigm shift occurring in the field of genetic toxicology. It provides mechanistic insight in a human-relevant cell-model, paired with measurement of a conventional endpoint, to inform the potential for adverse health effects. This work provides support for combining these assays in an integrated test strategy for accurate, higher throughput genetic toxicology testing in this metabolically competent human progenitor cell line.

A cross-sector call to improve carcinogenicity risk assessment through use of genomic methodologies

Robust genomic approaches are now available to realize improvements in efficiencies and translational relevance of cancer risk assessments for drugs and chemicals. Mechanistic and pathway data generated via genomics provide opportunities to advance beyond historical reliance on apical endpoints of uncertain human relevance. Published research and regulatory evaluations include many examples for which genomic data have been applied to address cancer risk assessment as a health protection endpoint. The alignment of mature, robust, reproducible, and affordable technologies with increasing demands for reduced animal testing sets the stage for this important transition. We present our shared vision for change from leading scientists from academic, government, nonprofit, and industrial sectors and chemical and pharmaceutical safety applications. This call to action builds upon a 2017 workshop on "Advances and Roadblocks for Use of Genomics in Cancer Risk Assessment." The authors propose a path for implementation of innovative cancer risk assessment including incorporating genomic signatures to assess mechanistic relevance of carcinogenicity and enhanced use of genomics in benchmark dose and point of departure evaluations. Novel opportunities for the chemical and pharmaceutical sectors to combine expertise, resources, and objectives to achieve a common goal of improved human health protection are identified.

TGx-DDI, a Transcriptomic Biomarker for Genotoxicity Hazard Assessment of Pharmaceuticals and Environmental Chemicals

Genotoxicity testing is an essential component of the safety assessment paradigm required by regulatory agencies world-wide for analysis of drug candidates, and environmental and industrial chemicals. Current genotoxicity testing batteries feature a high incidence of irrelevant positive findings—particularly for in vitro chromosomal damage (CD) assays. The risk management of compounds with positive in vitro findings is a major challenge and requires complex, time consuming, and costly follow-up strategies including animal testing. Thus, regulators are urgently in need of new testing approaches to meet legislated mandates. Using machine learning, we identified a set of transcripts that responds predictably to DNA-damage in human cells that we refer to as the TGx-DDI biomarker, which was originally referred to as TGx-28.65. We proposed to use this biomarker in conjunction with current genotoxicity testing batteries to differentiate compounds with irrelevant “false” positive findings in the in vitro CD assays from true DNA damaging agents (i.e., for de-risking agents that are clastogenic in vitro but not in vivo). We validated the performance of the TGx-DDI biomarker to identify true DNA damaging agents, assessed intra- and inter- laboratory reproducibility, and cross-platform performance. Recently, to augment the application of this biomarker, we developed a high-throughput cell-based genotoxicity testing system using the NanoString nCounter® technology. Here, we review the status of TGx-DDI development, its integration in the genotoxicity testing paradigm, and progress to date in its qualification at the US Food and Drug Administration (FDA) as a drug development tool. If successfully validated and implemented, the TGx-DDI biomarker assay is expected to significantly augment the current strategy for the assessment of genotoxic hazards for drugs and chemicals.

Genetic toxicology in silico protocol

In silico toxicology (IST) approaches to rapidly assess chemical hazard, and usage of such methods is increasing in all applications but especially for regulatory submissions, such as for assessing chemicals under REACH as well as the ICH M7 guideline for drug impurities. There are a number of obstacles to performing an IST assessment, including uncertainty in how such an assessment and associated expert review should be performed or what is fit for purpose, as well as a lack of confidence that the results will be accepted by colleagues, collaborators and regulatory authorities. To address this, a project to develop a series of IST protocols for different hazard endpoints has been initiated and this paper describes the genetic toxicity in silico (GIST) protocol. The protocol outlines a hazard assessment framework including key effects/mechanisms and their relationships to endpoints such as gene mutation and clastogenicity. IST models and data are reviewed that support the assessment of these effects/mechanisms along with defined approaches for combining the information and evaluating the confidence in the assessment. This protocol has been developed through a consortium of toxicologists, computational scientists, and regulatory scientists across several industries to support the implementation and acceptance of in silico approaches.

Identification of p53 Activators in a Human Microarray Compendium

Biomarkers predictive of molecular and toxicological effects are needed to interpret emerging high-throughput transcriptomic data streams. The previously characterized 63 gene TGx-DDI biomarker that includes 20 genes known to be regulated by p53 was previously shown to accurately predict DNA damage in chemically treated cells. We comprehensively evaluated whether the molecular basis of the DDI predictions was based on a p53-dependent response. The biomarker was compared to microarray data in a compendium derived from human cells using the Running Fisher test, a nonparametric correlation test. Using the biomarker, we identified conditions that led to p53 activation, including exposure to the chemical nutlin-3 which disrupts interactions between p53 and the negative regulator MDM2 or by knockdown of MDM2. The expression of most of the genes in the biomarker (75%) were found to depend on p53 activation status based on gene behavior after TP53 overexpression or knockdown. The biomarker identified DDI chemicals that were strong inducers of p53 in wild-type cells; these p53 responses were decreased or abolished in cells after p53 knockdown by siRNAs. Using the biomarker, we screened ∼1950 chemicals in ∼9800 human cell line chemical vs control comparisons and identified ∼100 chemicals that caused p53 activation. Among the positive chemicals were many that are known to activate p53 through direct and indirect DNA damaging mechanisms. These results contribute to the evidence that the TGx-DDI biomarker is useful for identifying chemicals that cause DDI and activate p53.

Considerations for Strategic Use of High-Throughput Transcriptomics Chemical Screening Data in Regulatory Decisions

Recently, numerous organizations, including governmental regulatory agencies in the U.S. and abroad, have proposed using data from New Approach Methodologies (NAMs) for augmenting and increasing the pace of chemical assessments. NAMs are broadly defined as any technology, methodology, approach or combination thereof that can be used to provide information on chemical hazard and risk assessment that avoids the use of intact animals. High-throughput transcriptomics (HTTr) is a type of NAM that uses gene expression profiling as an endpoint for rapidly evaluating the effects of large numbers of chemicals on in vitro cell culture systems. As compared to targeted high-throughput screening (HTS) approaches that measure the effect of chemical X on target Y, HTTr is a non-targeted approach that allows researchers to more broadly characterize the integrated response of an intact biological system to chemicals that may affect a specific biological target or many biological targets under a defined set of treatment conditions (time, concentration, etc.). HTTr screening performed in concentration-response mode can provide potency estimates for the concentrations of chemicals that produce perturbations in cellular response pathways. Here, we discuss study design considerations for HTTr concentration-response screening and present a framework for the use of HTTr-based biological pathway-altering concentrations (BPACs) in a screening-level, risk-based chemical prioritization approach. The framework involves concentration-response modeling of HTTr data, mapping gene level responses to biological pathways, determination of BPACs, in vitro-to-in vivo extrapolation (IVIVE) and comparison to human exposure predictions.

Toxicogenomics - What added Value Do These Approaches Provide for Carcinogen Risk Assessment?

It is still a major challenge to protect humans at workplaces and in the environment. To cope with this task, it is a prerequisite to obtain detailed information on the extent of chemical perturbations of biological pathways, in particular, adaptive vs. adverse effects and the dose-response relationships. This knowledge serves as the basis for the classification of non-carcinogens and carcinogens and for further distinguishing carcinogens in genotoxic (DNA damaging) or non-genotoxic compounds. Basing on quantitative dose-response relationships, points of departures can be derived for chemical risk assessment. In recent years, new methods have shown their capability to support the established rodent models of carcinogenicity testing. In vitro high throughput screening assays assess more comprehensively cell response. In addition, omics technologies were applied to study the mode of action of chemicals whereby the term "toxicogenomics" comprises various technologies such as transcriptomics, epigenomics, or metabolomics. This review aims to summarize the current state of toxicogenomic approaches in risk science and to compare them with established ones. For example, measurement of global transcriptional changes generates meaningful information for toxicological risk assessment such as accurate classification of genotoxic/non-genotoxic carcinogens. Alteration in mRNA expression offers previously unknown insights in the mode of action and enables the definition of key events. Based on these, benchmark doses can be calculated for the transition from an adaptive to an adverse state. In short, this review assesses the potential and challenges of transcriptomics and addresses the impact of other omics technologies on risk assessment in terms of hazard identification and dose-response assessment.

Assessment of the performance of the TGx-DDI biomarker to detect DNA damage-inducing agents using quantitative RT-PCR in TK6 cells

Gene expression biomarkers are now available for application in the identification of genotoxic hazards. The TGx-DDI transcriptomic biomarker can accurately distinguish DNA damage-inducing (DDI) from non-DDI exposures based on changes in the expression of 64 biomarker genes. The 64 genes were previously derived from whole transcriptome DNA microarray profiles of 28 reference agents (14 DDI and 14 non-DDI) after 4 h treatments of TK6 human lymphoblastoid cells. To broaden the applicability of TGx-DDI, we tested the biomarker using quantitative RT-PCR (qPCR), which is accessible to most molecular biology laboratories. First, we selectively profiled the expression of the 64 biomarker genes using TaqMan qPCR assays in 96-well arrays after exposing TK6 cells to the 28 reference agents for 4 h. To evaluate the classification capability of the qPCR profiles, we used the reference qPCR signature to classify 24 external validation chemicals using two different methods-a combination of three statistical analyses and an alternative, the Running Fisher test. The qPCR results for the reference set were comparable to the original microarray biomarker; 27 of the 28 reference agents (96%) were accurately classified. Moreover, the two classification approaches supported the conservation of TGx-DDI classification capability using qPCR; the combination of the two approaches accurately classified 21 of the 24 external validation chemicals, demonstrating 100% sensitivity, 81% specificity, and 91% balanced accuracy. This study demonstrates that qPCR can be used when applying the TGx-DDI biomarker and will improve the accessibility of TGx-DDI for genotoxicity screening. Environ. Mol. Mutagen. 60: 122-133, 2019. © 2018 Her Majesty the Queen in Right of Canada Environmental and Molecular Mutagenesis.

Assessment of in vivo genotoxicity of citrated-coated silver nanoparticles via transcriptomic analysis of rabbit liver tissue

Background

Silver nanoparticles (AgNPs) are widely used in industrial and household applications, arousing concern regarding their safety in humans. The risks posed by stabilizer-coated AgNPs continue to be unclear, and assessing their toxicity is for an understanding of the safety issues involved in their use in various applications.

Purpose

We aimed to investigated the long-term toxicity of citrate-coated silver nanoparticles (cAgNPs) in liver tissue using several toxicity tests and transcriptomic analysis at 7 and 28 days after a single intravenous injection into rabbit ear veins (n=4).

Materials and methods

The cAgNPs used in this study were in the form of a 20% (w/v) aqueous solution, and their size was 7.9±0.95 nm, measured using transmission electron microscopy. The animal experiments were performed based on the principles of good laboratory practice.

Results

Our results showed that the structure and function of liver tissue were disrupted due to a single exposure to cAgNPs. In addition, in vivo comet assay showed unrepaired genotoxicity in liver tissue until 4 weeks after a single injection, suggesting a potential carcinogenic effect of cAgNPs. In our transcriptomic analysis, a total of 244 genes were found to have differential expression at 28 days after a single cAgNP injection. Carefully curated pathway analysis of these genes using Pathway Studio and Ingenuity Pathway Analysis tools revealed major molecular networks responding to cAgNP exposure and indicated a high correlation of the genes with inflammation, hepatotoxicity, and cancer. Molecular validation suggested potential biomarkers for assessing the toxicity of accumulated cAgNPs.

Conclusion

Our investigation highlights the risk associated with a single cAgNP exposure with unrepaired damage persisting for at least a month.

Using a gene expression biomarker to identify DNA damage-inducing agents in microarray profiles

High-throughput transcriptomic technologies are increasingly being used to screen environmental chemicals in vitro to provide mechanistic context for regulatory testing. The TGx-DDI biomarker is a 64-gene expression profile generated from testing 28 model chemicals or treatments (13 that cause DNA damage and 15 that do not) in human TK6 cells. While the biomarker is very accurate at predicting DNA damage inducing (DDI) potential using the nearest shrunken centroid method, the broad utility of the biomarker using other computational methods is not fully known. Here, we determined the accuracy of the biomarker used with the Running Fisher test, a nonparametric correlation test. In TK6 cells, the methods could readily differentiate DDI and non-DDI compounds with balanced accuracies of 87-97%, depending on the threshold for determining DDI positives. The methods identified DDI agents in the metabolically competent hepatocyte cell line HepaRG (accuracy = 90%) but not in HepG2 cells or hepatocytes derived from embryonic stem cells (60 and 80%, respectively). DDI was also accurately classified when the gene expression changes were derived using the nCounter technology (accuracy = 89%). In addition, we found: (1) not all genes contributed equally to the correlations; (2) the minimal overlap in genes between the biomarker and the individual comparisons required for significant positive correlation was 10 genes, but usually was much higher; and (3) different sets of genes in the biomarker can by themselves contribute to the significant correlations. Overall, these results demonstrate the utility of the biomarker to accurately classify DDI agents. Environ. Mol. Mutagen. 59:772-784, 2018. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

A Qualitative Modeling Approach for Whole Genome Prediction Using High-Throughput Toxicogenomics Data and Pathway-Based Validation

Efficient high-throughput transcriptomics (HTT) tools promise inexpensive, rapid assessment of possible biological consequences of human and environmental exposures to tens of thousands of chemicals in commerce. HTT systems have used relatively small sets of gene expression measurements coupled with mathematical prediction methods to estimate genome-wide gene expression and are often trained and validated using pharmaceutical compounds. It is unclear whether these training sets are suitable for general toxicity testing applications and the more diverse chemical space represented by commercial chemicals and environmental contaminants. In this work, we built predictive computational models that inferred whole genome transcriptional profiles from a smaller sample of surrogate genes. The model was trained and validated using a large scale toxicogenomics database with gene expression data from exposure to heterogeneous chemicals from a wide range of classes (the Open TG-GATEs data base). The method of predictor selection was designed to allow high fidelity gene prediction from any pre-existing gene expression data set, regardless of animal species or data measurement platform. Predictive qualitative models were developed with this TG-GATES data that contained gene expression data of human primary hepatocytes with over 941 samples covering 158 compounds. A sequential forward search-based greedy algorithm, combining different fitting approaches and machine learning techniques, was used to find an optimal set of surrogate genes that predicted differential expression changes of the remaining genome. We then used pathway enrichment of up-regulated and down-regulated genes to assess the ability of a limited gene set to determine relevant patterns of tissue response. In addition, we compared prediction performance using the surrogate genes found from our greedy algorithm (referred to as the SV2000) with the landmark genes provided by existing technologies such as L1000 (Genometry) and S1500 (Tox21), finding better predictive performance for the SV2000. The ability of these predictive algorithms to predict pathway level responses is a positive step toward incorporating mode of action (MOA) analysis into the high throughput prioritization and testing of the large number of chemicals in need of safety evaluation.

Enniatin B induces expression changes in the electron transport chain pathway related genes in lymphoblastic T-cell line

Enniatin B is a ionophoric and lipophilic mycotoxin which reaches the bloodstream and has the ability to penetrate into cellular membranes. The purpose of this study was to reveal changes in the gene expression profile caused by enniatin B in human Jurkat lymphoblastic T-cells after 24 h of exposure at 1.5, 3 and 5 μM by next generation sequencing. It was found that up to 27% of human genome expression levels were significantly altered (5750 genes for both down-regulation and up-regulation). In the three enniatin B concentrations studied 245 differentially expressed genes were found to be overlapped, 83 were down and 162 up-regulated. ConsensusPathDB analysis of over-representation of differentially expressed genes provided a list of gene ontology terms in which several biological processes related to nucleoside monophosphate metabolic process, respiratory chain complex, electron transport chain, oxidative phosphorylation and cellular respiration were the most altered. Also, an interesting correlation was found between enniatin B toxicity and the up-regulation of the UCP protein complex. In summary, the transcriptomic analysis revealed that mitochondria are the organelles showing more related differentially expressed genes. Consequently, differentially expressed genes involved in biological processes, molecular functions and pathways related to mitochondrial metabolism and respiration were significantly changed.

Investigating the Generalizability of the MultiFlow ® DNA Damage Assay and Several Companion Machine Learning Models With a Set of 103 Diverse Test Chemicals

The in vitro MultiFlow DNA Damage assay multiplexes p53, γH2AX, phospho-histone H3, and polyploidization biomarkers into 1 flow cytometric analysis (Bryce, S. M., Bernacki, D. T., Bemis, J. C., and Dertinger, S. D. (2016). Genotoxic mode of action predictions from a multiplexed flow cytometric assay and a machine learning approach. Environ. Mol. Mutagen. 57, 171-189). The work reported herein evaluated the generalizability of the method, as well as several data analytics strategies, to a range of chemical classes not studied previously. TK6 cells were exposed to each of 103 diverse chemicals, 86 of which were supplied by the National Toxicology Program (NTP) and selected based upon responses in genetic damage assays conducted under the Tox21 program. Exposures occurred for 24 h over a range of concentrations, and cell aliquots were removed at 4 and 24 h for analysis. Multiplexed response data were evaluated using 3 machine learning models designed to predict genotoxic mode of action based on data from a training set of 85 previously studied chemicals. Of 54 chemicals with sufficient information to make an a priori call on genotoxic potential, the prediction models' accuracies were 79.6% (random forest), 88.9% (logistic regression), and 90.7% (artificial neural network). A majority vote ensemble of the 3 models provided 92.6% accuracy. Forty-nine NTP chemicals were not adequately tested (maximum concentration did not approach assay's cytotoxicity limit) and/or had insufficient conventional genotoxicity data to allow their genotoxic potential to be defined. For these chemicals MultiFlow data will be useful in future research and hypothesis testing. Collectively, the results suggest the MultiFlow assay and associated data analysis strategies are broadly generalizable, demonstrating high predictability when applied to new chemicals and classes of compounds.