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

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.

Carcinogenicity testing in drug development: Getting it right

For a pharmaceutical drug, carcinogenicity testing occurs in rodents to identify its tumorigenic potential to allow assessment of the risk from its use in humans. Testing takes the form of 2-year studies in mice and rats and/or more recently, a 6-month study in transgenic mice. This paper examines the process of regulatory interaction regarding carcinogenicity testing, notably through the United States (US) Food and Drug Administration (FDA) Special Protocol Assessment (SPA) process to seek Executive Carcinogenicity Assessment Committee (ECAC) approval. The content of 37 submissions to CAC were examined. The paper also examines the outcome from such agency engagement, notably around study dose level selection as well as looking at the design of proposed carcinogenicity study protocols used in submissions (including numbers of animals, control group aspects and toxicokinetic [TK] evaluation). Overall, it was shown that the current process of regulatory interaction allows for studies acceptable to support eventual drug approval and marketing. However, it was established that areas exist to improve the content of submission documents and study design aspects.

Increased Cell Proliferation as a Key Event in Chemical Carcinogenesis: Application in an Integrated Approach for the Testing and Assessment of Non-Genotoxic Carcinogenesis

In contrast to genotoxic carcinogens, there are currently no internationally agreed upon regulatory tools for identifying non-genotoxic carcinogens of human relevance. The rodent cancer bioassay is only used in certain regulatory sectors and is criticized for its limited predictive power for human cancer risk. Cancer is due to genetic errors occurring in single cells. The risk of cancer is higher when there is an increase in the number of errors per replication (genotoxic agents) or in the number of replications (cell proliferation-inducing agents). The default regulatory approach for genotoxic agents whereby no threshold is set is reasonably conservative. However, non-genotoxic carcinogens cannot be regulated in the same way since increased cell proliferation has a clear threshold. An integrated approach for the testing and assessment (IATA) of non-genotoxic carcinogens is under development at the OECD, considering learnings from the regulatory assessment of data-rich substances such as agrochemicals. The aim is to achieve an endorsed IATA that predicts human cancer better than the rodent cancer bioassay, using methodologies that equally or better protect human health and are superior from the view of animal welfare/efficiency. This paper describes the technical opportunities available to assess cell proliferation as the central gateway of an IATA for non-genotoxic carcinogenicity.

A tiered testing strategy based on in vitro phenotypic and transcriptomic data for selecting representative petroleum UVCBs for toxicity evaluation in vivo

Hazard evaluation of substances of "unknown or variable composition, complex reaction products and biological materials" (UVCBs) remains a major challenge in regulatory science because their chemical composition is difficult to ascertain. Petroleum substances are representative UVCBs and human cell-based data have been previously used to substantiate their groupings for regulatory submissions. We hypothesized that a combination of phenotypic and transcriptomic data could be integrated to make decisions as to selection of group-representative worst-case petroleum UVCBs for subsequent toxicity evaluation in vivo. We used data obtained from 141 substances from 16 manufacturing categories previously tested in 6 human cell types (induced pluripotent stem cell [iPSC]-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, and MCF7 and A375 cell lines). Benchmark doses for gene-substance combinations were calculated, and both transcriptomic and phenotype-derived points of departure (PODs) were obtained. Correlation analysis and machine learning were used to assess associations between phenotypic and transcriptional PODs and to determine the most informative cell types and assays, thus representing a cost-effective integrated testing strategy. We found that 2 cell types-iPSC-derived-hepatocytes and -cardiomyocytes-contributed the most informative and protective PODs and may be used to inform selection of representative petroleum UVCBs for further toxicity evaluation in vivo. Overall, although the use of new approach methodologies to prioritize UVCBs has not been widely adopted, our study proposes a tiered testing strategy based on iPSC-derived hepatocytes and cardiomyocytes to inform selection of representative worst-case petroleum UVCBs from each manufacturing category for further toxicity evaluation in vivo.

The cell transformation assay to assess potential carcinogenic properties of nanoparticles

Nanoparticles (NPs) are present in many daily life products with particular physical-chemical properties (size, density, porosity, geometry …) giving very interesting technological properties. Their use is continuously growing and NPs represent a new challenge in terms of risk assessment, consumers being multi-exposed. Toxic effects have already been identified such as oxidative stress, genotoxicity, inflammatory effects, and immune reactions, some of which are leading to carcinogenesis. Cancer is a complex phenomenon implying multiple modes of action and key events, and prevention strategies in cancer include a proper assessment of the properties of NPs. Therefore, introduction of new agents like NPs into the market creates fresh regulatory challenges for an adequate safety evaluation and requires new tools. The Cell Transformation Assay (CTA) is an in vitro test able of highlighting key events of characteristic phases in the cancer process, initiation and promotion. This review presents the development of this test and its use with NPs. The article underlines also the critical issues to address for assessing NPs carcinogenic properties and approaches for improving its relevance.

Towards achieving a modern science-based paradigm for agrochemical carcinogenicity assessment

The rodent cancer bioassay has been the standard approach to fulfill regulatory requirements for assessing human carcinogenic potential of agrochemicals, food additives, industrial chemicals, and pharmaceuticals. Decades of research have described the limitations of the rodent cancer bioassay leading to international initiatives to seek alternatives and establish approaches that modernize carcinogenicity assessment. Biologically relevant approaches can provide mechanistic information and increased efficiency for evaluating hazard and risk of chemical carcinogenicity to humans. The application of human-relevant mechanistic understanding to support new approaches to carcinogenicity assessment will be invaluable for regulatory decision-making. The present work outlines the challenges and opportunities that authorities should consider as they come together to build a roadmap that leads to global acceptance and incorporation of fit-for-purpose, scientifically defensible new approaches for human-relevant carcinogenicity assessment of agrochemicals.

Predicting Chemical Carcinogens Using a Hybrid Neural Network Deep Learning Method

Determining environmental chemical carcinogenicity is urgently needed as humans are increasingly exposed to these chemicals. In this study, we developed a hybrid neural network (HNN) method called HNN-Cancer to predict potential carcinogens of real-life chemicals. The HNN-Cancer included a new SMILES feature representation method by modifying our previous 3D array representation of 1D SMILES simulated by the convolutional neural network (CNN). We developed binary classification, multiclass classification, and regression models based on diverse non-congeneric chemicals. Along with the HNN-Cancer model, we developed models based on the random forest (RF), bootstrap aggregating (Bagging), and adaptive boosting (AdaBoost) methods for binary and multiclass classification. We developed regression models using HNN-Cancer, RF, support vector regressor (SVR), gradient boosting (GB), kernel ridge (KR), decision tree with AdaBoost (DT), KNeighbors (KN), and a consensus method. The performance of the models for all classifications was assessed using various statistical metrics. The accuracy of the HNN-Cancer, RF, and Bagging models were 74%, and their AUC was ~0.81 for binary classification models developed with 7994 chemicals. The sensitivity was 79.5% and the specificity was 67.3% for the HNN-Cancer, which outperforms the other methods. In the case of multiclass classification models with 1618 chemicals, we obtained the optimal accuracy of 70% with an AUC 0.7 for HNN-Cancer, RF, Bagging, and AdaBoost, respectively. In the case of regression models, the correlation coefficient (R) was around 0.62 for HNN-Cancer and RF higher than the SVM, GB, KR, DTBoost, and NN machine learning methods. Overall, the HNN-Cancer performed better for the majority of the known carcinogen experimental datasets. Further, the predictive performance of HNN-Cancer on diverse chemicals is comparable to the literature-reported models that included similar and less diverse molecules. Our HNN-Cancer could be used in identifying potentially carcinogenic chemicals for a wide variety of chemical classes.

Grouping of UVCB substances with dose-response transcriptomics data from human cell-based assays

The application of in vitro biological assays as new approach methodologies (NAMs) to support grouping of UVCB (unknown or variable composition, complex reaction products, and biological materials) substances has recently been demonstrated. In addition to cell-based phenotyping as NAMs, in vitro transcriptomic profiling is used to gain deeper mechanistic understanding of biological responses to chemicals and to support grouping and read-across. However, the value of gene expression profiling for characterizing complex substances like UVCBs has not been explored. Using 141 petroleum substance extracts, we performed dose-response transcriptomic profiling in human induced pluripotent stem cell (iPSC)-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, as well as cell lines MCF7 and A375. The goal was to determine whether transcriptomic data can be used to group these UVCBs and to further characterize the molecular basis for in vitro biological responses. We found distinct transcriptional responses for petroleum substances by manufacturing class. Pathway enrichment informed interpretation of effects of substances and UVCB petroleum-class. Transcriptional activity was strongly correlated with concentration of polycyclic aromatic compounds (PAC), especially in iPSC-derived hepatocytes. Supervised analysis using transcriptomics, alone or in combination with bioactivity data collected on these same substances/cells, suggest that transcriptomics data provide useful mechanistic information, but only modest additional value for grouping. Overall, these results further demonstrate the value of NAMs for grouping of UVCBs, identify informative cell lines, and provide data that could be used for justifying selection of substances for further testing that may be required for registration.

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.

Rethinking chronic toxicity and carcinogenicity assessment for agrochemicals project (ReCAAP): A reporting framework to support a weight of evidence safety assessment without long-term rodent bioassays

Rodent cancer bioassays have been long-required studies for regulatory assessment of human cancer hazard and risk. These studies use hundreds of animals, are resource intensive, and certain aspects of these studies have limited human relevance. The past 10 years have seen an exponential growth of new technologies with the potential to effectively evaluate human cancer hazard and risk while reducing, refining, or replacing animal use. To streamline and facilitate uptake of new technologies, a workgroup comprised of scientists from government, academia, non-governmental organizations, and industry stakeholders developed a framework for waiver rationales of rodent cancer bioassays for consideration in agrochemical safety assessment. The workgroup used an iterative approach, incorporating regulatory agency feedback, and identifying critical information to be considered in a risk assessment-based weight of evidence determination of the need for rodent cancer bioassays. The reporting framework described herein was developed to support a chronic toxicity and carcinogenicity study waiver rationale, which includes information on use pattern(s), exposure scenario(s), pesticidal mode-of-action, physicochemical properties, metabolism, toxicokinetics, toxicological data including mechanistic data, and chemical read-across from similar registered pesticides. The framework could also be applied to endpoints other than chronic toxicity and carcinogenicity, and for chemicals other than agrochemicals.

DeepCarc: Deep Learning-Powered Carcinogenicity Prediction Using Model-Level Representation

Carcinogenicity testing plays an essential role in identifying carcinogens in environmental chemistry and drug development. However, it is a time-consuming and label-intensive process to evaluate the carcinogenic potency with conventional 2-years rodent animal studies. Thus, there is an urgent need for alternative approaches to providing reliable and robust assessments on carcinogenicity. In this study, we proposed a DeepCarc model to predict carcinogenicity for small molecules using deep learning-based model-level representations. The DeepCarc Model was developed using a data set of 692 compounds and evaluated on a test set containing 171 compounds in the National Center for Toxicological Research liver cancer database (NCTRlcdb). As a result, the proposed DeepCarc model yielded a Matthews correlation coefficient (MCC) of 0.432 for the test set, outperforming four advanced deep learning (DL) powered quantitative structure-activity relationship (QSAR) models with an average improvement rate of 37%. Furthermore, the DeepCarc model was also employed to screen the carcinogenicity potential of the compounds from both DrugBank and Tox21. Altogether, the proposed DeepCarc model could serve as an early detection tool (https://github.com/TingLi2016/DeepCarc) for carcinogenicity assessment.

CarcSeq Measurement of Rat Mammary Cancer Driver Mutations and Relation to Spontaneous Mammary Neoplasia

The ability to deduce carcinogenic potential from subchronic, repeat dose rodent studies would constitute a major advance in chemical safety assessment and drug development. This study investigated an error-corrected NGS method (CarcSeq) for quantifying cancer driver mutations (CDMs) and deriving a metric of clonal expansion predictive of future neoplastic potential. CarcSeq was designed to interrogate subsets of amplicons encompassing hotspot CDMs applicable to a variety of cancers. Previously, normal human breast DNA was analyzed by CarcSeq and metrics based on mammary-specific CDMs were correlated with tissue donor age, a surrogate of breast cancer risk. Here we report development of parallel methodologies for rat. The utility of the rat CarcSeq method for predicting neoplastic potential was investigated by analyzing mammary tissue of 16-week-old untreated rats with known differences in spontaneous mammary neoplasia (Fischer 344, Wistar Han, and Sprague Dawley). Hundreds of mutants with mutant fractions ≥ 10-4 were quantified in each strain, most were recurrent mutations, and 42.5% of the nonsynonymous mutations have human homologs. Mutants in the mammary-specific target of the most tumor-sensitive strain (Sprague Dawley) showed the greatest nonsynonymous/synonymous mutation ratio, indicative of positive selection consistent with clonal expansion. For the mammary-specific target (Hras, Pik3ca, and Tp53 amplicons), median absolute deviation correlated with percentages of rats that develop spontaneous mammary neoplasia at 104 weeks (Pearson r = 1.0000, 1-tailed p = .0010). Therefore, this study produced evidence CarcSeq analysis of spontaneously occurring CDMs can be used to derive an early metric of clonal expansion relatable to long-term neoplastic outcome.

Assessment of clonal expansion using CarcSeq measurement of lung cancer driver mutations and correlation with mouse strain- and sex-related incidence of spontaneous lung neoplasia

Quantification of variation in levels of spontaneously occurring cancer driver mutations (CDMs) was developed to assess clonal expansion and predict future risk of neoplasm development. Specifically, an error-corrected next generation sequencing method, CarcSeq, and a mouse CarcSeq panel (analogous to human and rat panels) were developed and used to quantify low-frequency mutations in a panel of amplicons enriched in hotspot CDMs. Mutations in a subset of panel amplicons, Braf, Egfr, Kras, Stk11 and Tp53, were related to incidence of lung neoplasms at two years. This was achieved by correlating median absolute deviation (MAD) from the overall median mutant fraction (MF) measured in the lung DNA of 16-week-old male and female, B6C3F1 and CD-1 mice (10 mice/sex/strain) with percentages of spontaneous alveolar/bronchioloalveolar adenomas and carcinomas reported in bioassay control groups. 1,586 mouse lung mutants with MFs >1 x 10-4 were recovered. The ratio of non-synonymous to synonymous mutations was used to assess the proportion of recovered mutations conferring a positive selective advantage. The greatest ratio was observed in what is considered the most lung tumor-sensitive model examined, male B6C3F1 mice. Of the recurrent, non-synonymous mouse mutations recovered, 55.5% have been reported in human tumors, with many located in or around the mouse equivalent of human cancer hotspot codons. MAD for the same subset of amplicons measured in normal human lung DNA samples showed a correlation of moderate strength and borderline significance) with age (a cancer risk factor), as well as age-related cumulative lung cancer risk, suggesting MAD may inform species extrapolation.

Repeat-dose toxicity prediction with Generalized Read-Across (GenRA) using targeted transcriptomic data: A proof-of-concept case study

Read-across is a data gap filling technique utilized to predict the toxicity of a target chemical using data from similar analogues. Recent efforts such as Generalized Read-Across (GenRA) facilitate automated read-across predictions for untested chemicals. GenRA makes predictions of toxicity outcomes based on "neighboring" chemicals characterized by chemical and bioactivity fingerprints. Here we investigated the impact of biological similarities on neighborhood formation and read-across performance in predicting hazard (based on repeat-dose testing outcomes from US EPA ToxRefDB v2.0). We used targeted transcriptomic data on 93 genes for 1060 chemicals in HepaRG^™ cells that measure nuclear receptor activation, xenobiotic metabolism, cellular stress, cell cycle progression, and apoptosis. Transcriptomic similarity between chemicals was calculated using binary hit-calls from concentration-response data for each gene. We evaluated GenRA performance in predicting ToxRefDB v2.0 hazard outcomes using the area under the Receiver Operating Characteristic (ROC) curve (AUC) for the baseline approach (chemical fingerprints) versus transcriptomic fingerprints and a combination of both (hybrid). For all endpoints, there were significant but only modest improvements in ROC AUC scores of 0.01 (2.1%) and 0.04 (7.3%) with transcriptomic and hybrid descriptors, respectively. However, for liver-specific toxicity endpoints, ROC AUC scores improved by 10% and 17% for transcriptomic and hybrid descriptors, respectively. Our findings suggest that using hybrid descriptors formed by combining chemical and targeted transcriptomic information can improve in vivo toxicity predictions in the right context.

Quantification of cancer driver mutations in human breast and lung DNA using targeted, error-corrected CarcSeq

There is a need for scientifically-sound, practical approaches to improve carcinogenicity testing. Advances in DNA sequencing technology and knowledge of events underlying cancer development have created an opportunity for progress in this area. The long-term goal of this work is to develop variation in cancer driver mutation (CDM) levels as a metric of clonal expansion of cells carrying CDMs because these important early events could inform carcinogenicity testing. The first step toward this goal was to develop and validate an error-corrected next-generation sequencing method to analyze panels of hotspot cancer driver mutations (hCDMs). The "CarcSeq" method that was developed uses unique molecular identifier sequences to construct single-strand consensus sequences for error correction. CarcSeq was used for mutational analysis of 13 amplicons encompassing >20 hotspot CDMs in normal breast, normal lung, ductal carcinomas, and lung adenocarcinomas. The approach was validated by detecting expected differences related to tissue type (normal vs. tumor and breast vs. lung) and mutation spectra. CarcSeq mutant fractions (MFs) correlated strongly with previously obtained ACB-PCR mutant fraction (MF) measurements from the same samples. A reconstruction experiment, in conjunction with other analyses, showed CarcSeq accurately quantifies MFs ≥10-4 . CarcSeq MF measurements were correlated with tissue donor age and breast cancer risk. CarcSeq MF measurements were correlated with deviation from median MFs analyzed to assess clonal expansion. Thus, CarcSeq is a promising approach to advance cancer risk assessment and carcinogenicity testing practices. Paradigms that should be investigated to advance this strategy for carcinogenicity testing are proposed.

A Rat Liver Transcriptomic Point of Departure Predicts a Prospective Liver or Non-liver Apical Point of Departure

Identifying a toxicity point of departure (POD) is a required step in human health risk characterization of crop protection molecules, and this POD has historically been derived from apical endpoints across a battery of animal-based toxicology studies. Using rat transcriptome and apical data for 79 molecules obtained from Open TG-GATES (Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System) (632 datasets), the hypothesis was tested that a short-term exposure, transcriptome-based liver biological effect POD (BEPOD) could estimate a longer-term exposure "systemic" apical endpoint POD. Apical endpoints considered were body weight, clinical observation, kidney weight and histopathology and liver weight and histopathology. A BMDExpress algorithm using Gene Ontology Biological Process gene sets was optimized to derive a liver BEPOD most predictive of a systemic apical POD. Liver BEPODs were stable from 3 h to 29 days of exposure; the median fold difference of the 29-day BEPOD to BEPODs from earlier time points was approximately 1 (range: 0.7-1.1). Strong positive correlation (Pearson R = 0.86) and predictive accuracy (root mean square difference = 0.41) were observed between a concurrent (29 days) liver BEPOD and the systemic apical POD. Similar Pearson R and root mean square difference values were observed for comparisons between a 29-day systemic apical POD and liver BEPODs derived from 3 h to 15 days of exposure. These data across 79 molecules suggest that a longer-term exposure study apical POD from liver and non-liver compartments can be estimated using a liver BEPOD derived from an acute or subacute exposure study.