Adverse outcome pathway-driven identification of rat liver tumorigens in short-term assays

Chemical-induced liver cancer: an adverse outcome pathway perspective

INTRODUCTION

The evaluation of the potential carcinogenicity is a key consideration in the risk assessment of chemicals. Predictive toxicology is currently switching toward non-animal approaches that rely on the mechanistic understanding of toxicity.

AREAS COVERED

Adverse outcome pathways (AOPs) present toxicological processes, including chemical-induced carcinogenicity, in a visual and comprehensive manner, which serve as the conceptual backbone for the development of non-animal approaches eligible for hazard identification. The current review provides an overview of the available AOPs leading to liver cancer and discusses their use in advanced testing of liver carcinogenic chemicals. Moreover, the challenges related to their use in risk assessment are outlined, including the exploitation of available data, the need for semantic ontologies, and the development of quantitative AOPs.

EXPERT OPINION

To exploit the potential of liver cancer AOPs in the field of risk assessment, 3 immediate prerequisites need to be fulfilled. These include developing human relevant AOPs for chemical-induced liver cancer, increasing the number of AOPs integrating quantitative toxicodynamic and toxicokinetic data, and developing a liver cancer AOP network. As AOPs and other areas in the field continue to evolve, liver cancer AOPs will progress into a reliable and robust tool serving future risk assessment and management.

Evaluating the mode of action of perfluorooctanoic acid-induced liver tumors in male Sprague-Dawley rats using a toxicogenomic approach

The mode of action (MOA) underlying perfluorooctanoic acid (PFOA)-induced liver tumors in rats is proposed to involve peroxisome proliferator-activated receptor α (PPARα) agonism. Despite clear PPARα activation evidence in rodent livers, the mechanisms driving cell growth remain elusive. Herein, we used dose-responsive apical endpoints and transcriptomic data to examine the proposed MOA. Male Sprague-Dawley rats were treated with 0, 1, 5, and 15 mg/kg PFOA for 7, 14, and 28 days via oral gavage. We showed PFOA induced hepatomegaly along with hepatocellular hypertrophy in rats. PPARα was activated in a dose-dependent manner. Toxicogenomic analysis revealed six early biomarkers (Cyp4a1, Nr1d1, Acot1, Acot2, Ehhadh, and Vnn1) in response to PPARα activation. A transient rise in hepatocellular DNA synthesis was demonstrated while Ki-67 labeling index showed no change. Transcriptomic analysis indicated no significant enrichment in pathways related to DNA synthesis, apoptosis, or the cell cycle. Key cyclins including Ccnd1, Ccnb1, Ccna2, and Ccne2 were dose-dependently suppressed by PFOA. Oxidative stress and the nuclear factor-κB signaling pathway were unaffected. Overall, evidence for PFOA-induced hepatocellular proliferation was transient within the studied timeframe. Our findings underscore the importance of considering inter-species differences and chemical-specific effects when evaluating the carcinogenic risk of PFOA in humans.

Applying genomics in regulatory toxicology: a report of the ECETOC workshop on omics threshold on non-adversity

In a joint effort involving scientists from academia, industry and regulatory agencies, ECETOC's activities in Omics have led to conceptual proposals for: (1) A framework that assures data quality for reporting and inclusion of Omics data in regulatory assessments; and (2) an approach to robustly quantify these data, prior to interpretation for regulatory use. In continuation of these activities this workshop explored and identified areas of need to facilitate robust interpretation of such data in the context of deriving points of departure (POD) for risk assessment and determining an adverse change from normal variation. ECETOC was amongst the first to systematically explore the application of Omics methods, now incorporated into the group of methods known as New Approach Methodologies (NAMs), to regulatory toxicology. This support has been in the form of both projects (primarily with CEFIC/LRI) and workshops. Outputs have led to projects included in the workplan of the Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) group of the Organisation for Economic Co-operation and Development (OECD) and to the drafting of OECD Guidance Documents for Omics data reporting, with potentially more to follow on data transformation and interpretation. The current workshop was the last in a series of technical methods development workshops, with a sub-focus on the derivation of a POD from Omics data. Workshop presentations demonstrated that Omics data developed within robust frameworks for both scientific data generation and analysis can be used to derive a POD. The issue of noise in the data was discussed as an important consideration for identifying robust Omics changes and deriving a POD. Such variability or "noise" can comprise technical or biological variation within a dataset and should clearly be distinguished from homeostatic responses. Adverse outcome pathways (AOPs) were considered a useful framework on which to assemble Omics methods, and a number of case examples were presented in illustration of this point. What is apparent is that high dimension data will always be subject to varying processing pipelines and hence interpretation, depending on the context they are used in. Yet, they can provide valuable input for regulatory toxicology, with the pre-condition being robust methods for the collection and processing of data together with a comprehensive description how the data were interpreted, and conclusions reached.

Short-term in vivo testing to discriminate genotoxic carcinogens from non-genotoxic carcinogens and non-carcinogens using next-generation RNA sequencing, DNA microarray, and qPCR

Next-generation RNA sequencing (RNA-Seq) has identified more differentially expressed protein-coding genes (DEGs) and provided a wider quantitative range of expression level changes than conventional DNA microarrays. JEMS·MMS·Toxicogenomics group studied DEGs with targeted RNA-Seq on freshly frozen rat liver tissues and on formalin-fixed paraffin-embedded (FFPE) rat liver tissues after 28 days of treatment with chemicals and quantitative real-time PCR (qPCR) on rat and mouse liver tissues after 4 to 48 h treatment with chemicals and analyzed by principal component analysis (PCA) as statics. Analysis of rat public DNA microarray data (Open TG-GATEs) was also performed. In total, 35 chemicals were analyzed [15 genotoxic hepatocarcinogens (GTHCs), 9 non-genotoxic hepatocarcinogens (NGTHCs), and 11 non-genotoxic non-hepatocarcinogens (NGTNHCs)]. As a result, 12 marker genes (Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2, and Tubb4b) were proposed to discriminate GTHCs from NGTHCs and NGTNHCs. U.S. Environmental Protection Agency studied DEGs induced by 4 known GTHCs in rat liver using DNA microarray and proposed 7 biomarker genes, Bax, Bcmp1, Btg2, Ccng1, Cdkn1a, Cgr19, and Mgmt for GTHCs. Studies involving the use of whole-transcriptome RNA-Seq upon exposure to chemical carcinogens in vivo have also been performed in rodent liver, kidney, lung, colon, and other organs, although discrimination of GTHCs from NGTHCs was not examined. Candidate genes published using RNA-Seq, qPCR, and DNA microarray will be useful for the future development of short-term in vivo studies of environmental carcinogens using RNA-Seq.

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.

Deriving time-concordant event cascades from gene expression data: A case study for Drug-Induced Liver Injury (DILI)

Adverse event pathogenesis is often a complex process which compromises multiple events ranging from the molecular to the phenotypic level. In toxicology, Adverse Outcome Pathways (AOPs) aim to formalize this as temporal sequences of events, in which event relationships should be supported by causal evidence according to the tailored Bradford-Hill criteria. One of the criteria is whether events are consistently observed in a certain temporal order and, in this work, we study this time concordance using the concept of “first activation” as data-driven means to generate hypotheses on potentially causal mechanisms. As a case study, we analysed liver data from repeat-dose studies in rats from the TG-GATEs database which comprises measurements across eight timepoints, ranging from 3 hours to 4 weeks post-treatment. We identified time-concordant gene expression-derived events preceding adverse histopathology, which serves as surrogate readout for Drug-Induced Liver Injury (DILI). We find known mechanisms in DILI to be time-concordant, and show further that significance, frequency and log fold change (logFC) of differential expression are metrics which can additionally prioritize events although not necessary to be mechanistically relevant. Moreover, we used the temporal order of transcription factor (TF) expression and regulon activity to identify transcriptionally regulated TFs and subsequently combined this with prior knowledge on functional interactions to derive detailed gene-regulatory mechanisms, such as reduced Hnf4a activity leading to decreased expression and activity of Cebpa. At the same time, also potentially novel events are identified such as Sox13 which is highly significantly time-concordant and shows sustained activation over time. Overall, we demonstrate how time-resolved transcriptomics can derive and support mechanistic hypotheses by quantifying time concordance and how this can be combined with prior causal knowledge, with the aim of both understanding mechanisms of toxicity, as well as potential applications to the AOP framework. We make our results available in the form of a Shiny app (https://anikaliu.shinyapps.io/dili_cascades), which allows users to query events of interest in more detail.

Understanding mechanisms from systems-scale biological data is of great relevance in toxicology as well as drug discovery; however how to generate causal hypotheses instead of correlations is by no means clear. In this work, we study the conserved temporal order of events and present an automatable framework to quantify and characterize time concordance across a large set of time-series. We apply this concept to events derived from time-resolved gene expression and histopathology from the TG-GATEs in vivo liver data as a case study. We were able to recover known events involved in the pathogenesis of Drug-Induced Liver Injury (DILI), and identify potentially novel pathway and transcription factors (TFs) which precede adverse histopathology. As complementary sources of evidence for causality, we additionally show how time concordance and prior knowledge on plausible interactions between TFs can be combined to derive causal hypotheses on the TFs’ mode of regulation and interaction partners. Overall, the results derived in our case study can serve as valuable hypothesis-free starting points for the development of Adverse Outcome Pathways for DILI, and demonstrate that our approach provides a novel angle to prioritize mechanistically relevant events.

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.

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.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Adverse Outcome Pathways as Versatile Tools in Liver Toxicity Testing

Adverse outcome pathways (AOPs) are tools to capture and visualize mechanisms driving toxicological effects. They share a common structure consisting of a molecular initiating event, a series of key events connected by key event relationships and an adverse outcome. Development and evaluation of AOPs ideally comply with guidelines issued by the Organization for Economic Cooperation and Development. AOPs have been introduced for major types of hepatotoxicity, which is not a surprise, as the liver is a frequent target for systemic adversity. Various applications for AOPs have been proposed in the areas of toxicology and chemical risk assessment, in particular in relation to the establishment of quantitative structure-activity relationships, the elaboration of prioritization strategies, and the development of novel in vitro toxicity screening tests and testing strategies.

Application of humanized mice to toxicology studies: Evaluation of the human relevance of the mode of action for rodent liver tumor formation by activators of the constitutive androstane receptor (CAR)

The constitutive androstane receptor (CAR)-mediated mode of action (MOA) for phenobarbital (PB)-induced rodent liver tumor formation has been established, with increased hepatocyte proliferation, which is a key event in tumor formation. Previous studies have demonstrated that PB and other CAR-activators stimulate proliferation in cultured rodent hepatocytes, but not in cultured human hepatocytes. However, in the genetically humanized CAR and pregnane X receptor (PXR) mouse (hCAR/hPXR mouse, downstream genes are still mouse), PB increased hepatocyte proliferation and tumor production in vivo. In contrast to the hCAR/hPXR mouse, studies with chimeric mice with human hepatocytes (PXB-mouse, both receptor and downstream genes are human) demonstrated that PB did not increase human hepatocyte proliferation in vivo. PB increased hepatocyte proliferation in a chimeric mouse model with rat hepatocytes, indicating that the lack of human hepatocyte proliferation is not due to any functional defect in the chimeric mouse liver environment. Gene expression analysis demonstrated that the downstream genes of CAR/PXR activation were similar in hCAR/hPXR and CD-1 mice, but differed from those observed in chimeric mice with human hepatocytes. These findings strongly support the conclusion that the MOA for CAR-mediated rodent liver tumor formation is qualitatively implausible for humans. Indeed, epidemiological studies have found no causal link between PB and human liver tumors. There are many similarities with respect to hepatic effects and species differences between rodent CAR and peroxisome proliferator-activated receptor α activators. Based on our research, the chimeric mouse with human hepatocytes (PXB-mouse) is reliable for human cancer risk assessment of test chemicals.

An evaluation of carcinogenicity predictors from short-term and sub chronic repeat-dose studies of agrochemicals in rats: Opportunities to refine and reduce animal use

The aim of this study was to examine whether short term, repeat dose, rat studies provide sufficient information about potential carcinogenicity to enable predictions about the carcinogenic potential of agrochemicals to be made earlier in compound development. This study aimed to identify any correlations between toxicity findings obtained for short term rat studies (28 day and 90 day) and neoplastic findings obtained from 24 month rat carcinogenicity studies for agrochemical compounds (18 compounds) tested in Han Wistar and Sprague Dawley rats. The macroscopic pathology, microscopic pathology, hematology, biochemistry, organ weights, estrogen receptor activation and genotoxicity results were examined. Seven out of 18 non genotoxic compounds developed tumors in treated rats in the carcinogenicity study and of these, two compounds showed no preneoplastic findings in the affected tissues (false positives). Of the remaining five true positives, correlations were noted between corneal opacity and keratitis (90 day study) as early indicators of squamous cell carcinoma and papilloma of the cornea of the eye (compound 1, a hydroxyphenylpyruvate dioxygenase inhibitor) and inflammation of the stomach and kidney (90 day study) and gastric squamous cell papilloma and squamous cell carcinoma and renal tubular adenoma and carcinoma, respectively (compound 12, a fungicide with multisite activity). Minor decreases in uterine weight and increases in estradiol hydroxylation activity at 28 days were associated with endometrial adenocarcinoma (compound 18, a mitochondrial complex II electron transport inhibitor). Early liver weight increases and hepatocellular centrilobular hypertrophy (28 day study) were associated with thyroid follicular adenomas (compound 11, a succinate dehydrogenase inhibitor) in female animals only. Hepatic centrilobular hypertrophy (28 day studies) correlated with thyroid adenomas in males in carcinogenicity studies (compound 2, a hydroxyphenylpyruvate dioxygenase inhibitor). In contrast, treatment related, nasopharynx tumors (compound 3, an elongase inhibitor) and uterine adenocarcinoma (compound 9, a succinate dehydrogenase inhibitor) could not be correlated with findings from the short term studies examined. Eleven compounds displayed preneoplastic findings with no tumors (false negatives) and there were no compounds with no preneoplastic findings and no tumors (true negatives). This work indicates the value of examining historical, short term studies for specific, nonneoplastic findings which correlate with tumors in carcinogenicity studies, which may obviate the need for further animal carcinogenicity studies.

Quantitative Transcriptional Biomarkers of Xenobiotic Receptor Activation in Rat Liver for the Early Assessment of Drug Safety Liabilities

The robust transcriptional plasticity of liver mediated through xenobiotic receptors underlies its ability to respond rapidly and effectively to diverse chemical stressors. Thus, drug-induced gene expression changes in liver serve not only as biomarkers of liver injury, but also as mechanistic sentinels of adaptation in metabolism, detoxification, and tissue protection from chemicals. Modern RNA sequencing methods offer an unmatched opportunity to quantitatively monitor these processes in parallel and to contextualize the spectrum of dose-dependent stress, adaptation, protection, and injury responses induced in liver by drug treatments. Using this approach, we profiled the transcriptional changes in rat liver following daily oral administration of 120 different compounds, many of which are known to be associated with clinical risk for drug-induced liver injury by diverse mechanisms. Clustering, correlation, and linear modeling analyses were used to identify and optimize coexpressed gene signatures modulated by drug treatment. Here, we specifically focused on prioritizing 9 key signatures for their pragmatic utility for routine monitoring in initial rat tolerability studies just prior to entering drug development. These signatures are associated with 5 canonical xenobiotic nuclear receptors (AHR, CAR, PXR, PPARα, ER), 3 mediators of reactive metabolite-mediated stress responses (NRF2, NRF1, P53), and 1 liver response following activation of the innate immune response. Comparing paradigm chemical inducers of each receptor to the other compounds surveyed enabled us to identify sets of optimized gene expression panels and associated scoring algorithms proposed as quantitative mechanistic biomarkers with high sensitivity, specificity, and quantitative accuracy. These findings were further qualified using public datasets, Open TG-GATEs and DrugMatrix, and internal development compounds. With broader collaboration and additional qualification, the quantitative toxicogenomic framework described here could inform candidate selection prior to committing to drug development, as well as complement and provide a deeper understanding of the conventional toxicology study endpoints used later in drug development.

Critical evaluation of the human relevance of the mode of action for rodent liver tumor formation by activators of the constitutive androstane receptor (CAR)

Many nongenotoxic chemicals have been shown to produce liver tumors in mice and/or rats by a mode of action (MOA) involving activation of the constitutive androstane receptor (CAR). Studies with phenobarbital (PB) and other compounds have identified the key events for this MOA: CAR activation; increased hepatocellular proliferation; altered foci formation; and ultimately the development of adenomas/carcinomas. In terms of human relevance, the pivotal species difference is that CAR activators are mitogenic agents in mouse and rat hepatocytes, but they do not stimulate increased hepatocellular proliferation in humans. This conclusion is supported by substantial in vitro studies with cultured rodent and human hepatocytes and also by in vivo studies with chimeric mice with human hepatocytes. Examination of the literature reveals many similarities in the hepatic effects and species differences between activators of rodent CAR and the peroxisome proliferator-activated receptor alpha (PPARα), with PPARα activators also not being mitogenic agents in human hepatocytes. Overall, a critical analysis of the available data demonstrates that the established MOA for rodent liver tumor formation by PB and other CAR activators is qualitatively not plausible for humans. This conclusion is supported by data from several human epidemiology studies.

A Gene Expression Biomarker Identifies Chemical Modulators of Estrogen Receptor α in an MCF-7 Microarray Compendium

Identification of chemicals that affect hormone-regulated systems will help to predict endocrine disruption. In our previous study, a 46 gene biomarker was found to be an accurate predictor of estrogen receptor (ER) α modulation in chemically treated MCF-7 cells. Here, potential ERα modulators were identified using the biomarker by screening a microarray compendium consisting of ∼1600 gene expression comparisons representing exposure to ∼1200 chemicals. A total of ∼170 chemicals were identified as potential ERα modulators. In the Connectivity Map 2.0 collection, 75 and 39 chemicals were predicted to activate or suppress ERα, and they included 12 and six known ERα agonists and antagonists/selective ERα modulators, respectively. Nineteen and eight of the total number were also identified as active in an ERα transactivation assay carried out in an MCF-7-derived cell line used to screen the Tox21 10K chemical library in agonist or antagonist modes, respectively. Chemicals predicted to modulate ERα in MCF-7 cells were examined further using global and targeted gene expression in wild-type and ERα-null cells, transactivation assays, and cell-free ERα coregulator interaction assays. Environmental chemicals classified as weak and very weak agonists were confirmed to activate ERα including apigenin, kaempferol, and oxybenzone. Novel activators included digoxin, nabumetone, ivermectin, and six progestins. Novel suppressors included emetine, mifepristone, niclosamide, and proscillaridin. Our strategy will be useful to identify environmentally relevant ERα modulators in future high-throughput transcriptomic screens.

Thresholds Derived From Common Measures in Rat Studies Are Predictive of Liver Tumorigenic Chemicals

We hypothesized that typical tissue and clinical chemistry (ClinChem) end points measured in rat toxicity studies exhibit chemical-independent biological thresholds beyond which cancer occurs. Using the rat in vivo TG-GATES study, 75 chemicals were examined across chemical-dose-time comparisons that could be linked to liver tumor outcomes. Thresholds for liver weight to body weight (LW/BW) and 21 serum ClinChem end points were defined as the maximum and minimum values for those exposures that did not lead to liver tumors in rats. Upper thresholds were identified for LW/BW (117%), aspartate aminotransferase (195%), alanine aminotransferase (141%), alkaline phosphatase (152%), and total bilirubin (115%), and lower thresholds were identified for phospholipids (82%), relative albumin (93%), total cholesterol (82%), and total protein (94%). Thresholds derived from the TG-GATES data set were consistent across other acute and subchronic rat studies. A training set of ClinChem and LW/BW thresholds derived from a 38 chemical training set from TG-GATES was predictive of liver tumor outcomes for a test set of 37 independent TG-GATES chemicals (91%). The thresholds were most predictive when applied to 7d treatments (98%). These findings provide support that biological thresholds for common end points in rodent studies can be used to predict chemical tumorigenic potential.

A set of six Gene expression biomarkers and their thresholds identify rat liver tumorigens in short-term assays

Traditional methods for cancer risk assessment are retrospective, resource-intensive, and not feasible for the vast majority of environmental chemicals. In earlier studies, we used a set of six biomarkers to accurately identify liver tumorigens in transcript profiles derived from chemically-treated rats using either a Toxicological Priority Index (ToxPi) approach or using derived biomarker thresholds for cancer. The biomarkers consisting of 7-113 genes are used to predict the most common liver cancer molecular initiating events: genotoxicity, cytotoxicity and activation of the xenobiotic receptors AhR, CAR, ER, and PPARα. In the present study, we apply and evaluate the performance of these methods for cancer prediction in an independent rat liver study of 44 chemicals (6 h-7d exposures) examined by Affymetrix arrays. In the first approach, ToxPi ranking of biomarker scores consistently gave the highest scores to tumorigenic chemical-dose pairs; balanced accuracies for identification of liver tumorigenic chemicals were up to 89 %. The second approach used tumorigenic thresholds derived in the present study or from our earlier study that were set at the maximum value for chemical-dose exposures without detectable liver tumor outcomes. Using these thresholds, balanced accuracies were up to 90 %. Both approaches identified all tumorigenic chemicals. Almost all of the tumorigenic chemicals activated more than one MIE. We also compared biomarker responses between two types of profiling platforms (Affymetrix full-genome array, TempO-Seq 1500+ array containing ∼2600 genes) and found that the lack of the full set of biomarker genes on the 1500+ array resulted in decreased ability to identify chemicals that activate the MIEs. Overall, these results demonstrate that predictive approaches based on the 6 biomarkers could be used in short-term assays to identify chemicals and their doses that induce liver tumors, the most common endpoint in rodent bioassays.

Mining a human transcriptome database for chemical modulators of NRF2

Nuclear factor erythroid-2 related factor 2 (NRF2) encoded by the NFE2L2 gene is a transcription factor critical for protecting cells from chemically-induced oxidative stress. We developed computational procedures to identify chemical modulators of NRF2 in a large database of human microarray data. A gene expression biomarker was built from statistically-filtered gene lists derived from microarray experiments in primary human hepatocytes and cancer cell lines exposed to NRF2-activating chemicals (oltipraz, sulforaphane, CDDO-Im) or in which the NRF2 suppressor Keap1 was knocked down by siRNA. Directionally consistent biomarker genes were further filtered for those dependent on NRF2 using a microarray dataset from cells after NFE2L2 siRNA knockdown. The resulting 143-gene biomarker was evaluated as a predictive tool using the correlation-based Running Fisher algorithm. Using 59 gene expression comparisons from chemically-treated cells with known NRF2 activating potential, the biomarker gave a balanced accuracy of 93%. The biomarker was comprised of many well-known NRF2 target genes (AKR1B10, AKR1C1, NQO1, TXNRD1, SRXN1, GCLC, GCLM), 69% of which were found to be bound directly by NRF2 using ChIP-Seq. NRF2 activity was assessed across ~9840 microarray comparisons from ~1460 studies examining the effects of ~2260 chemicals in human cell lines. A total of 260 and 43 chemicals were found to activate or suppress NRF2, respectively, most of which have not been previously reported to modulate NRF2 activity. Using a NRF2-responsive reporter gene in HepG2 cells, we confirmed the activity of a set of chemicals predicted using the biomarker. The biomarker will be useful for future gene expression screening studies of environmentally-relevant chemicals.

A Set of Six Gene Expression Biomarkers Identify Rat Liver Tumorigens in Short-term Assays

Chemical-induced liver cancer occurs in rodents through well-characterized adverse outcome pathways. We hypothesized that measurement of the 6 most common molecular initiating events (MIEs) in liver cancer adverse outcome pathways in short-term assays using only gene expression will allow early identification of chemicals and their associated doses that are likely to be tumorigenic in the liver in 2-year bioassays. We tested this hypothesis using transcript data from a rat liver microarray compendium consisting of 2013 comparisons of 146 chemicals administered at doses with previously established effects on rat liver tumor induction. Five MIEs were measured using previously characterized gene expression biomarkers composed of gene sets predictive for genotoxicity and activation of 1 or more xenobiotic receptors (aryl hydrocarbon receptor, constitutive activated receptor, estrogen receptor, and peroxisome proliferator-activated receptor α). Because chronic injury can be important in tumorigenesis, we also developed a biomarker for cytotoxicity that had a 96% balanced accuracy. Characterization of the genes in each biomarker set using the unsupervised TXG-MAP network model demonstrated that the genes were associated with distinct functional coexpression modules. Using the Toxicological Priority Index to rank chemicals based on their ability to activate the MIEs showed that chemicals administered at tumorigenic doses clearly gave the highest ranked scores. Balanced accuracies using thresholds derived from either TG-GATES or DrugMatrix data sets to predict tumorigenicity in independent sets of chemicals were up to 93%. These results show that a MIE-directed approach using only gene expression biomarkers could be used in short-term assays to identify chemicals and their doses that cause tumors.

Gene Expression Thresholds Derived From Short-term Exposures Identify Rat Liver Tumorigens

Traditional methods for cancer risk assessment are resource-intensive, retrospective, and not feasible for the vast majority of environmental chemicals. In this study, we investigated whether quantitative genomic data from short-term studies may be used to set protective thresholds for potential tumorigenic effects. We hypothesized that gene expression biomarkers measuring activation of the key early events in established pathways for rodent liver cancer exhibit cross-chemical thresholds for tumorigenesis predictive for liver cancer risk. We defined biomarker thresholds for 6 major liver cancer pathways using training sets of chemicals with short-term genomic data (3-29 days of exposure) from the TG-GATES (n = 77 chemicals) and DrugMatrix (n = 86 chemicals) databases and then tested these thresholds within and between datasets. The 6 pathway biomarkers represented genotoxicity, cytotoxicity, and activation of xenobiotic, steroid, and lipid receptors (aryl hydrocarbon receptor, constitutive activated receptor, estrogen receptor, and peroxisome proliferator-activated receptor α). Thresholds were calculated as the maximum values derived from exposures without detectable liver tumor outcomes. We identified clear response values that were consistent across training and test sets. Thresholds derived from the TG-GATES training set were highly predictive (97%) in a test set of independent chemicals, whereas thresholds derived from the DrugMatrix study were 96%-97% predictive for the TG-GATES study. Threshold values derived from an abridged gene list (2/biomarker) also exhibited high predictive accuracy (91%-94%). These findings support the idea that early genomic changes can be used to establish threshold estimates or "molecular tipping points" that are predictive of later-life health outcomes.

Early microRNA indicators of PPARα pathway activation in the liver

MicroRNAs (miRNAs) are short non-coding RNA species that play key roles in post-transcriptional regulation of gene expression. MiRNAs also serve as a promising source of early biomarkers for different environmental exposures and health effects, although there is limited information linking miRNA changes to specific target pathways. In this study, we measured liver miRNAs in male B6C3F1 mice exposed to a known chemical activator of the peroxisome proliferator-activated receptor alpha (PPARα) pathway, di(2-ethylhexyl) phthalate (DEHP), for 7 and 28 days at concentrations of 0, 750, 1500, 3000, or 6000 ppm in feed. At the highest dose tested, DEHP altered 61 miRNAs after 7 days and 171 miRNAs after 28 days of exposure, with 48 overlapping miRNAs between timepoints. Analysis of these 48 common miRNAs indicated enrichment in PPARα–related targets and other pathways related to liver injury and cancer. Four of the 10 miRNAs exhibiting a clear dose trend were linked to the PPARα pathway: mmu-miRs-125a-5p, -182−5p, -20a−5p, and -378a−3p. mmu-miRs-182−5p and -378a−3p were subsequently measured using digital drop PCR across a dose range for DEHP and two related phthalates with weaker PPARα activity, di-n-octyl phthalate and n-butyl benzyl phthalate, following 7-day exposures. Analysis of mmu-miRs-182−5p and -378a−3p by transcriptional benchmark dose analysis correctly identified DEHP as having the greatest potency. However, benchmark dose estimates for DEHP based on these miRNAs (average 163; range 126−202 mg/kg-day) were higher on average than values for PPARα target genes (average 74; range 29−183 mg/kg-day). These findings identify putative miRNA biomarkers of PPARα pathway activity and suggest that early miRNA changes may be used to stratify chemical potency.

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.

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.

AhR Activation in Pharmaceutical Development: Applying Liver Gene Expression Biomarker Thresholds to Identify Doses Associated with Tumorigenic Risks in Rats

Aryl hydrocarbon receptor (AhR) activation is associated with carcinogenicity of non-genotoxic AhR-activating carcinogens such as 2, 3, 7, 8-tetrachlorodibenzodioxin (TCDD), and is often observed with drug candidate molecules in development and raises safety concerns. As downstream effectors of AhR signaling, the expression and activity of Cyp1a1 and Cyp1a2 genes are commonly monitored as evidence of AhR activation to inform carcinogenic risk of compounds in question. However, many marketed drugs and phytochemicals are reported to induce these Cyps modestly and are not associated with dioxin-like toxicity or carcinogenicity. We hypothesized that a threshold of AhR activation needs to be surpassed in a sustained manner in order for the dioxin-like toxicity to manifest, and a simple liver gene expression signature based on Cyp1a1 and Cyp1a2 from a short-term rat study could be used to assess AhR activation strength and differentiate tumorigenic dose levels from non-tumorigenic ones. To test this hypothesis, short term studies were conducted in Wistar Han rats with two AhR-activating carcinogens (TCDD and PCB126) at minimally carcinogenic and non-carcinogenic dose levels, and three AhR-activating non-carcinogens (omeprazole, mexiletine, and canagliflozin) at the top doses used in their reported 2-year rat carcinogenicity studies. A threshold of AhR activation was identified in rat liver that separated a meaningful "tumorigenic-strength AhR signal" from a statistically-significant AhR activation signal that was not associated with dioxin-like carcinogenicity. These studies also confirmed the importance of the sustainability of AhR activation for carcinogenic potential. A sustained activation of AhR above the threshold could thus be used in early pharmaceutical development to identify dose levels of drug candidates expected to exhibit dioxin-like carcinogenic potential.

Frequent Modulation of the Sterol Regulatory Element Binding Protein (SREBP) by Chemical Exposure in the Livers of Rats

Inappropriate activation of sterol regulatory element-binding proteins (SREBPs) can lead to non-alcoholic fatty liver disease (NAFLD). To link chemical exposure to SREBP activity, a previously characterized gene expression biomarker (Rooney et al., 2019) was used to identify microarray comparisons from rat liver that exhibited significant positive or negative correlation to the biomarker. The effects of 620 chemicals on SREBP activity were examined across 9305 chemical-dose-time microarray comparisons. SREBP was found to be frequently modulated by chemical exposure with 54% of the chemicals affecting SREBP activity. Activators included inhibitors of cholesterogenesis that act to inhibit HMG-CoA reductase (statins) or inhibit Cyp51 (conazoles). Most chemical effects were transient, lasting usually no more than 2-4 days. Modulation of SREBP in most cases led to coordinated increases or decreases in lipogenic and cholesterogenic genes. However, 570 chemical exposure conditions were identified in which regulation was uncoupled. Most of these conditions affected cholesterogenic genes in the absence of parallel effects on lipogenic genes. Together, these findings show that SREBP is a frequent target of chemical exposure and expression of lipogenic and cholesterogenic genes can be uncoupled.