The Power of Resolution: Contextualized Understanding of Biological Responses to Liver Injury Chemicals Using High-throughput Transcriptomics and Benchmark Concentration Modeling

Empirical Characterization of False Discovery Rates of Differentially Expressed Genes and Transcriptomic Benchmark Concentrations in Zebrafish Embryos

High-throughput transcriptomics (HTTr) is increasingly applied to zebrafish embryos to survey the toxicological effects of environmental chemicals. Before the adoption of this approach in regulatory testing, it is essential to characterize background noise in order to guide experimental designs. We thus empirically quantified the HTTr false discovery rate (FDR) across different embryo pool sizes, sample sizes, and concentration groups for toxicology studies. We exposed zebrafish embryos to 0.1% dimethyl sulfoxide (DMSO) for 5 days. Pools of 1, 5, 10, and 20 embryos were created (n = 24 samples for each pool size). Samples were sequenced on the TempO-Seq platform and then randomly assigned to mock treatment groups before differentially expressed gene (DEG), pathway, and benchmark concentration (BMC) analyses. Given that all samples were treated with DMSO, any significant DEGs, pathways, or BMCs are false positives. As expected, we found decreasing FDRs for DEG and pathway analyses with increasing pool and sample sizes. Similarly, FDRs for BMC analyses decreased with increasing pool size and concentration groups, with more stringent BMC premodel filtering reducing BMC FDRs. Our study provides foundational data for determining appropriate experiment designs for regulatory toxicity testing with HTTr in zebrafish embryos.

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.

Investigating the dependency of in vitro benchmark concentrations on exposure time in transcriptomics experiments

There is increasing interest to employ in vitro transcriptomics experiments in toxicological testing, for example to determine a point-of-departure (PoD) for chemical safety assessment. However current practices to derive PoD tend to utilise a single exposure time despite the importance of exposure time on the manifestation of toxicity caused by a chemical. Therefore it is important to investigate both concentration and exposure time to determine how these factors affect biological responses, and as a consequence, the derivation of PoDs. In this study, metabolically competent HepaRG cells were exposed to five known toxicants over a range of concentrations and time points for subsequent gene expression analysis, using a targeted RNA expression assay (TempO-Seq). A non-parametric factor-modelling approach was used to model the collective response of all significant genes, which exploited the interdependence of differentially expressed gene responses. This in turn allowed the determination of an isobenchmark response (isoBMR) curve for each chemical in a reproducible manner. For 2 of the 5 chemicals tested, the PoD was observed to vary by 0.5-1 log-order within the 48-h timeframe of the experiment. The approach and findings presented here clearly demonstrate the need to take both concentration and exposure time into account when designing in vitro toxicogenomics experiments to determine PoD. Doing so also provides a means to use concentration-time-response modelling as a basis to extrapolate a PoD from shorter to longer exposure durations, and to identify chemicals of concern that can cause cumulative effects over time.

Definition of the Neurotoxicity-Associated Metabolic Signature Triggered by Berberine and Other Respiratory Chain Inhibitors

To characterize the hits from a phenotypic neurotoxicity screen, we obtained transcriptomics data for valinomycin, diethylstilbestrol, colchicine, rotenone, 1-methyl-4-phenylpyridinium (MPP), carbaryl and berberine (Ber). For all compounds, the concentration triggering neurite degeneration correlated with the onset of gene expression changes. The mechanistically diverse toxicants caused similar patterns of gene regulation: the responses were dominated by cell de-differentiation and a triggering of canonical stress response pathways driven by ATF4 and NRF2. To obtain more detailed and specific information on the modes-of-action, the effects on energy metabolism (respiration and glycolysis) were measured. Ber, rotenone and MPP inhibited the mitochondrial respiratory chain and they shared complex I as the target. This group of toxicants was further evaluated by metabolomics under experimental conditions that did not deplete ATP. Ber (204 changed metabolites) showed similar effects as MPP and rotenone. The overall metabolic situation was characterized by oxidative stress, an over-abundance of NADH (>1000% increase) and a re-routing of metabolism in order to dispose of the nitrogen resulting from increased amino acid turnover. This unique overall pattern led to the accumulation of metabolites known as biomarkers of neurodegeneration (saccharopine, aminoadipate and branched-chain ketoacids). These findings suggest that neurotoxicity of mitochondrial inhibitors may result from an ensemble of metabolic changes rather than from a simple ATP depletion. The combi-omics approach used here provided richer and more specific MoA data than the more common transcriptomics analysis alone. As Ber, a human drug and food supplement, mimicked closely the mode-of-action of known neurotoxicants, its potential hazard requires further investigation.

High-Throughput Transcriptomics Differentiates Toxic versus Non-Toxic Chemical Exposures Using a Rat Liver Model

To address the challenge of limited throughput with traditional toxicity testing, a newly developed high-throughput transcriptomics (HTT) platform, together with a 5-day in vivo rat model, offers an alternative approach to estimate chemical exposures and provide reasonable estimates of toxicological endpoints. This study contains an HTT analysis of 18 environmental chemicals with known liver toxicity. They were evaluated using male Sprague Dawley rats exposed to various concentrations daily for five consecutive days via oral gavage, with data collected on the sixth day. Here, we further explored the 5-day rat model to identify potential gene signatures that can differentiate between toxic and non-toxic liver responses and provide us with a potential histopathological endpoint of chemical exposure. We identified a distinct gene expression pattern that differentiated non-hepatotoxic compounds from hepatotoxic compounds in a dose-dependent manner, and an analysis of the significantly altered common genes indicated that toxic chemicals predominantly upregulated most of the genes and several pathways in amino acid and lipid metabolism. Finally, our liver injury module analysis revealed that several liver-toxic compounds showed similarities in the key injury phenotypes of cellular inflammation and proliferation, indicating potential molecular initiating processes that may lead to a specific end-stage liver disease.

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.

From vision toward best practices: Evaluating in vitro transcriptomic points of departure for application in risk assessment using a uniform workflow

The growing number of chemicals in the current consumer and industrial markets presents a major challenge for regulatory programs faced with the need to assess the potential risks they pose to human and ecological health. The increasing demand for hazard and risk assessment of chemicals currently exceeds the capacity to produce the toxicity data necessary for regulatory decision making, and the applied data is commonly generated using traditional approaches with animal models that have limited context in terms of human relevance. This scenario provides the opportunity to implement novel, more efficient strategies for risk assessment purposes. This study aims to increase confidence in the implementation of new approach methods in a risk assessment context by using a parallel analysis to identify data gaps in current experimental designs, reveal the limitations of common approaches deriving transcriptomic points of departure, and demonstrate the strengths in using high-throughput transcriptomics (HTTr) to derive practical endpoints. A uniform workflow was applied across six curated gene expression datasets from concentration-response studies containing 117 diverse chemicals, three cell types, and a range of exposure durations, to determine tPODs based on gene expression profiles. After benchmark concentration modeling, a range of approaches was used to determine consistent and reliable tPODs. High-throughput toxicokinetics were employed to translate in vitro tPODs (µM) to human-relevant administered equivalent doses (AEDs, mg/kg-bw/day). The tPODs from most chemicals had AEDs that were lower (i.e., more conservative) than apical PODs in the US EPA CompTox chemical dashboard, suggesting in vitro tPODs would be protective of potential effects on human health. An assessment of multiple data points for single chemicals revealed that longer exposure duration and varied cell culture systems (e.g., 3D vs. 2D) lead to a decreased tPOD value that indicated increased chemical potency. Seven chemicals were flagged as outliers when comparing the ratio of tPOD to traditional POD, thus indicating they require further assessment to better understand their hazard potential. Our findings build confidence in the use of tPODs but also reveal data gaps that must be addressed prior to their adoption to support risk assessment applications.

Determination of in vitro hepatotoxic potencies of a series of perfluoroalkyl substances (PFASs) based on gene expression changes in HepaRG liver cells

Per- and polyfluoroalkyl substances (PFASs) are omnipresent and have been shown to induce a wide range of adverse health effects, including hepatotoxicity, developmental toxicity, and immunotoxicity. The aim of the present work was to assess whether human HepaRG liver cells can be used to obtain insight into differences in hepatotoxic potencies of a series of PFASs. Therefore, the effects of 18 PFASs on cellular triglyceride accumulation (AdipoRed assay) and gene expression (DNA microarray for PFOS and RT-qPCR for all 18 PFASs) were studied in HepaRG cells. BMDExpress analysis of the PFOS microarray data indicated that various cellular processes were affected at the gene expression level. From these data, ten genes were selected to assess the concentration-effect relationship of all 18 PFASs using RT-qPCR analysis. The AdipoRed data and the RT-qPCR data were used for the derivation of in vitro relative potencies using PROAST analysis. In vitro relative potency factors (RPFs) could be obtained for 8 PFASs (including index chemical PFOA) based on the AdipoRed data, whereas for the selected genes, in vitro RPFs could be obtained for 11-18 PFASs (including index chemical PFOA). For the readout OAT5 expression, in vitro RPFs were obtained for all PFASs. In vitro RPFs were found to correlate in general well with each other (Spearman correlation) except for the PPAR target genes ANGPTL4 and PDK4. Comparison of in vitro RPFs with RPFs obtained from in vivo studies in rats indicate that best correlations (Spearman correlation) were obtained for in vitro RPFs based on OAT5 and CXCL10 expression changes and external in vivo RPFs. HFPO-TA was found to be the most potent PFAS tested, being around tenfold more potent than PFOA. Altogether, it may be concluded that the HepaRG model may provide relevant data to provide insight into which PFASs are relevant regarding their hepatotoxic effects and that it can be applied as a screening tool to prioritize other PFASs for further hazard and risk assessment.

Establishment of a point of departure for CBD hepatotoxicity employing human HepaRG spheroids

The United States Food and Drug Administration recently approved the use of Cannabis sativa derived cannabidiol (CBD) in the treatment of Dravet Syndrome and Lennox-Gastaut Syndrome, under the trade name, Epidiolex. In double-blinded, placebo-controlled clinical trials, elevated ALT levels were observed in some patients, but these findings could not be uncoupled from the confounds of potential drug-drug interactions with co-administration of valproate and clobazam. Given the uncertainty of the potential hepatatoxic effects of CBD, the objective of the present study was to determine a point of departure for CBD, using human HepaRG spheroid cultures, followed by transcriptomic benchmark dose analysis. Treatment of HepaRG spheroids with CBD for 24 and 72 h, resulted in EC50 concentrations for cytotoxicity of 86.27 µM and 58.04 µM, respectively. Subsequent transcriptomic analysis at these timepoints demonstrated little alteration of gene and pathway data sets at a CBD concentration at or below 10 µM. Although this current analysis was conducted using liver cells, interestingly the findings at 72 h post CBD treatment showed suppression of many genes more commonly associated with immune regulation. Indeed, the immune system is a well-established target for CBD based on immune function assays. Collectively, in the present studies a point of departure was derived using transcriptomic changes produced by CBD in a human cell-based model system, which has been shown to accurately translate to human hepatotoxicity modeling.

Considerations for application of benchmark dose modeling in radiation research: workshop highlights

BACKGROUND

Exposure to different forms of ionizing radiation occurs in diverse occupational, medical, and environmental settings. Improving the accuracy of the estimated health risks associated with exposure is therefore, essential for protecting the public, particularly as it relates to chronic low dose exposures. A key aspect to understanding health risks is precise and accurate modeling of the dose-response relationship. Toward this vision, benchmark dose (BMD) modeling may be a suitable approach for consideration in the radiation field. BMD modeling is already extensively used for chemical hazard assessments and is considered statistically preferable to identifying low and no observed adverse effects levels. BMD modeling involves fitting mathematical models to dose-response data for a relevant biological endpoint and identifying a point of departure (the BMD, or its lower bound). Recent examples in chemical toxicology show that when applied to molecular endpoints (e.g. genotoxic and transcriptional endpoints), BMDs correlate to points of departure for more apical endpoints such as phenotypic changes (e.g. adverse effects) of interest to regulatory decisions. This use of BMD modeling may be valuable to explore in the radiation field, specifically in combination with adverse outcome pathways, and may facilitate better interpretation of relevant in vivo and in vitro dose-response data. To advance this application, a workshop was organized on June 3rd, 2022, in Ottawa, Ontario that brought together BMD experts in chemical toxicology and the radiation scientific community of researchers, regulators, and policy-makers. The workshop's objective was to introduce radiation scientists to BMD modeling and its practical application using case examples from the chemical toxicity field and demonstrate the BMDExpress software using a radiation dataset. Discussions focused on the BMD approach, the importance of experimental design, regulatory applications, its use in supporting the development of adverse outcome pathways, and specific radiation-relevant examples.

CONCLUSIONS

Although further deliberations are needed to advance the use of BMD modeling in the radiation field, these initial discussions and partnerships highlight some key steps to guide future undertakings related to new experimental work.

ToxicR: A computational platform in R for computational toxicology and dose-response analyses

The need to analyze the complex relationships observed in high-throughput toxicogenomic and other omic platforms has resulted in an explosion of methodological advances in computational toxicology. However, advancements in the literature often outpace the development of software researchers can implement in their pipelines, and existing software is frequently based on pre-specified workflows built from well-vetted assumptions that may not be optimal for novel research questions. Accordingly, there is a need for a stable platform and open-source codebase attached to a programming language that allows users to program new algorithms. To fill this gap, the Biostatistics and Computational Biology Branch of the National Institute of Environmental Health Sciences, in cooperation with the National Toxicology Program (NTP) and US Environmental Protection Agency (EPA), developed ToxicR, an open-source R programming package. The ToxicR platform implements many of the standard analyses used by the NTP and EPA, including dose-response analyses for continuous and dichotomous data that employ Bayesian, maximum likelihood, and model averaging methods, as well as many standard tests the NTP uses in rodent toxicology and carcinogenicity studies, such as the poly-K and Jonckheere trend tests. ToxicR is built on the same codebase as current versions of the EPA's Benchmark Dose software and NTP's BMDExpress software but has increased flexibility because it directly accesses this software. To demonstrate ToxicR, we developed a custom workflow to illustrate its capabilities for analyzing toxicogenomic data. The unique features of ToxicR will allow researchers in other fields to add modules, increasing its functionality in the future.

Substantiate a read-across hypothesis by using transcriptome data-A case study on volatile diketones

This case study explores the applicability of transcriptome data to characterize a common mechanism of action within groups of short-chain aliphatic α-, β-, and γ-diketones. Human reference in vivo data indicate that the α-diketone diacetyl induces bronchiolitis obliterans in workers involved in the preparation of microwave popcorn. The other three α-diketones induced inflammatory responses in preclinical in vivo animal studies, whereas beta and gamma diketones in addition caused neuronal effects. We investigated early transcriptional responses in primary human bronchiolar (PBEC) cell cultures after 24 h and 72 h of air-liquid exposure. Differentially expressed genes (DEGs) were assessed based on transcriptome data generated with the EUToxRisk gene panel of Temp-O-Seq^®. For each individual substance, genes were identified displaying a consistent differential expression across dose and exposure duration. The log fold change values of the DEG profiles indicate that α- and β-diketones are more active compared to γ-diketones. α-diketones in particular showed a highly concordant expression pattern, which may serve as a first indication of the shared mode of action. In order to gain a better mechanistic understanding, the resultant DEGs were submitted to a pathway analysis using ConsensusPathDB. The four α-diketones showed very similar results with regard to the number of activated and shared pathways. Overall, the number of signaling pathways decreased from α-to β-to γ-diketones. Additionally, we reconstructed networks of genes that interact with one another and are associated with different adverse outcomes such as fibrosis, inflammation or apoptosis using the TRANSPATH-database. Transcription factor enrichment and upstream analyses with the geneXplain platform revealed highly interacting gene products (called master regulators, MRs) per case study compound. The mapping of the resultant MRs on the reconstructed networks, visualized similar gene regulation with regard to fibrosis, inflammation and apoptosis. This analysis showed that transcriptome data can strengthen the similarity assessment of compounds, which is of particular importance, e.g., in read-across approaches. It is one important step towards grouping of compounds based on biological profiles.

Derivation of metabolic point of departure using high-throughput in vitro metabolomics: investigating the importance of sampling time points on benchmark concentration values in the HepaRG cell line

Amongst omics technologies, metabolomics should have particular value in regulatory toxicology as the measurement of the molecular phenotype is the closest to traditional apical endpoints, whilst offering mechanistic insights into the biological perturbations. Despite this, the application of untargeted metabolomics for point-of-departure (POD) derivation via benchmark concentration (BMC) modelling is still a relatively unexplored area. In this study, a high-throughput workflow was applied to derive PODs associated with a chemical exposure by measuring the intracellular metabolome of the HepaRG cell line following treatment with one of four chemicals (aflatoxin B₁, benzo[a]pyrene, cyclosporin A, or rotenone), each at seven concentrations (aflatoxin B₁, benzo[a]pyrene, cyclosporin A: from 0.2048 μM to 50 μM; rotenone: from 0.04096 to 10 μM) and five sampling time points (2, 6, 12, 24 and 48 h). The study explored three approaches to derive PODs using benchmark concentration modelling applied to single features in the metabolomics datasets or annotated metabolites or lipids: (1) the 1st rank-ordered unannotated feature, (2) the 1st rank-ordered putatively annotated feature (using a recently developed HepaRG-specific library of polar metabolites and lipids), and (3) 25th rank-ordered feature, demonstrating that for three out of four chemical datasets all of these approaches led to relatively consistent BMC values, varying less than tenfold across the methods. In addition, using the 1st rank-ordered unannotated feature it was possible to investigate temporal trends in the datasets, which were shown to be chemical specific. Furthermore, a possible integration of metabolomics-driven POD derivation with the liver steatosis adverse outcome pathway (AOP) was demonstrated. The study highlights that advances in technologies enable application of in vitro metabolomics at scale; however, greater confidence in metabolite identification is required to ensure PODs are mechanistically anchored.

Analyses of Transcriptomics Cell Signalling for Pre-Screening Applications in the Integrated Approach for Testing and Assessment of Non-Genotoxic Carcinogens

With recent rapid advancement of methodological tools, mechanistic understanding of biological processes leading to carcinogenesis is expanding. New approach methodologies such as transcriptomics can inform on non-genotoxic mechanisms of chemical carcinogens and can be developed for regulatory applications. The Organisation for the Economic Cooperation and Development (OECD) expert group developing an Integrated Approach to the Testing and Assessment (IATA) of Non-Genotoxic Carcinogens (NGTxC) is reviewing the possible assays to be integrated therein. In this context, we review the application of transcriptomics approaches suitable for pre-screening gene expression changes associated with phenotypic alterations that underlie the carcinogenic processes for subsequent prioritisation of downstream test methods appropriate to specific key events of non-genotoxic carcinogenesis. Using case studies, we evaluate the potential of gene expression analyses especially in relation to breast cancer, to identify the most relevant approaches that could be utilised as (pre-) screening tools, for example Gene Set Enrichment Analysis (GSEA). We also consider how to address the challenges to integrate gene panels and transcriptomic assays into the IATA, highlighting the pivotal omics markers identified for assay measurement in the IATA key events of inflammation, immune response, mitogenic signalling and cell injury.

Evaluating the Influences of Confounding Variables on Benchmark Dose using a Case Study in the Field of Ionizing Radiation

Purpose A vast amount of data regarding the effects of radiation stressors on transcriptional changes has been produced over the past few decades. These data have shown remarkable consistency across platforms and experimental design, enabling increased understanding of early molecular effects of radiation exposure. However, the value of transcriptomic data in the context of risk assessment is not clear and represents a gap that is worthy of further consideration. Recently, benchmark dose (BMD) modeling has shown promise in correlating a transcriptional point of departure (POD) to that derived using phenotypic outcomes relevant to human health risk assessment. Although frequently applied in chemical toxicity evaluation, our group has recently demonstrated application within the field of radiation research. This approach allows the possibility to quantitatively compare radiation-induced gene and pathway alterations across various datasets using BMD values and derive meaningful biological effects. However, before BMD modeling can confidently be used, an understanding of the impact of confounding variables on BMD outputs is needed. Methods: To this end, BMD modeling was applied to a publicly available microarray dataset (Gene Expression Omnibus #GSE23515) that used peripheral blood ex-vivo gamma-irradiated at 0.82 Gy/min, at doses of 0, 0.1, 0.5 or 2 Gy, and assessed 6 hours post-exposure. The dataset comprised six female smokers (F-S), six female non-smokers (F-NS), six male smokers (M-S), and six male non-smokers (M-NS). Results: A combined total of 412 genes were fit to models and the BMD distribution was noted to be bi-modal across the four groups. A total of 74, 41, 62 and 62 genes were unique to the F-NS, M-NS, F-S and M-S groups. Sixty-two BMD modeled genes and nine pathways were common across all four groups. There were no differential sensitivity of responses in the robust common genes and pathways. Conclusion: For radiation-responsive genes and pathways common across the study groups, the BMD distribution of transcriptional activity was unaltered by sex and smoking status. Although further validation of the data is needed, these initial findings suggest BMD values for radiation relevant genes and pathways are robust and could be explored further in future studies.

Benchmark dose modeling of transcriptional data: a systematic approach to identify best practices for study designs used in radiation research

PURPOSE

Benchmark dose (BMD) modeling is a method commonly used in chemical toxicology to identify the point of departure (POD) from a dose-response curve linked to a health-related outcome. Recently, it is being explored on transcriptional data and in adverse outcome pathways (AOPs). As AOPs are informed by diverse data types, it is important to understand the impact of study parameters such as dose selection, number of replicates and dose range on BMD outputs for radiation induced genes and pathways.

MATERIALS AND METHODS

Data were selected from the Gene Expression Omnibus (GSE52403) that featured gene expression profiles of peripheral blood samples from C57BL/6 mice 6 hours post-exposure to ¹³⁷Cs gamma-radiation at 0, 1, 2, 3, 4.5, 6, 8 and 10.5 Gy. The dataset comprised a broad dose-range over multiple dose-points with consistent dose spacing and multiple biological replicates. This dataset was ideal for systematically transforming across three categories: (1) dose-range, (2) dose-spacing and (3) number of controls/replicates. Across these categories, 29 transformed datasets were compared to the original dataset to determine the impact of each transformation on the BMD outputs.

RESULTS

Most of the experimental changes did not impact the BMD outputs. The transformed datasets were largely consistent with the original dataset in terms of number of reproduced genes modeled and absolute BMD values for genes and pathways. Variations in dose selection identified the importance of the absolute value of the lowest and second dose. It was determined that dose selection should include at least two doses <1 Gy and two >5 Gy to achieve meaningful BMD outputs. Changes to the number of biological replicates in the control and non-zero dose groups impacted the overall accuracy and precision of the BMD outputs as well as the ability to fit dose-response models consistent with the original dataset.

CONCLUSION

Successful application of transcriptomic BMD modeling for radiation datasets requires considerations of the exposure dose and the number of biological replicates. Most important is the selection of the lowest doses and dose spacing. Reflections on these parameters in experimental design will provide meaningful BMD outputs that could correlate well to apical endpoints of relevance to radiation exposure assessment.

ALOHA: Aggregated local extrema splines for high-throughput dose-response analysis

Computational methods for genomic dose-response integrate dose-response modeling with bioinformatics tools to evaluate changes in molecular and cellular functions related to pathogenic processes. These methods use parametric models to describe each gene's dose-response, but such models may not adequately capture expression changes. Additionally, current approaches do not consider gene co-expression networks. When assessing co-expression networks, one typically does not consider the dose-response relationship, resulting in 'co-regulated' gene sets containing genes having different dose-response patterns. To avoid these limitations, we develop an analysis pipeline called Aggregated Local Extrema Splines for High-throughput Analysis (ALOHA), which computes individual genomic dose-response functions using a flexible class Bayesian shape constrained splines and clusters gene co-regulation based upon these fits. Using splines, we reduce information loss due to parametric lack-of-fit issues, and because we cluster on dose-response relationships, we better identify co-regulation clusters for genes that have co-expressed dose-response patterns from chemical exposure. The clustered pathways can then be used to estimate a dose associated with a pre-specified biological response, i.e., the benchmark dose (BMD), and approximate a point of departure dose corresponding to minimal adverse response in the whole tissue/organism. We compare our approach to current parametric methods and our biologically enriched gene sets to cluster on normalized expression data. Using this methodology, we can more effectively extract the underlying structure leading to more cohesive estimates of gene set potency.

Deep Learning Image Analysis of High-Throughput Toxicology Assay Images

High-throughput chemical screening approaches often employ microscopy to capture photomicrographs from multi-well cell culture plates, generating thousands of images that require time-consuming human analysis. To automate this subjective and time-consuming manual process, we have developed a method that uses deep learning to automatically classify digital assay images. We have trained a convolutional neural network (CNN) to perform binary and multi-class classification. The binary classifier binned assay images into healthy (comparable to untreated controls) and altered (not comparable to untreated-control) classes with >98% accuracy; the multi-class classifier assigned "Healthy," "Intermediate" and "Altered" labels to assay images with >95% accuracy. Our dataset comprised high-resolution assay images from primary human hepatocytes and undifferentiated (proliferating) and differentiated 2D cultures of HepaRG cells. In this study we have focused on testing and fine-tuning various CNN architectures, including ResNet 34, 50 and 101. To visualize regions in the images that the CNN model used for classification, we employed Class Activation Maps (CAM). This allowed us to better understand the inner workings of the neural network and led to additional optimizations of the algorithm. The results indicate a strong correspondence between dosage and classifier-predicted scores, suggesting that these scores might be useful in further characterizing benchmark dose. Together, these results clearly demonstrate that deep-learning based automated image classification of cell morphology changes upon chemical-induced stress can yield highly accurate and reproducible assessments of cytotoxicity across a variety of cell types.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

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.

Reproducibility and robustness of high-throughput S1500+ transcriptomics on primary rat hepatocytes for chemical-induced hepatotoxicity assessment

•

TempO-Seq assays of rat hepatocytes in collagen sandwich are highly reproducible.

•

Gene expression analysis shows S1500+ is representative of the whole transcriptome.

•

Connectivity mapping shows consistency between TempO-Seq and Affymetrix data.

•

Gene set enrichment shows consistency between S1500+ and the whole transcriptome.

•

Gene set enrichment using hallmark gene sets informs hepatotoxicity.

Cell-based in vitro models coupled with high-throughput transcriptomics (HTTr) are increasingly utilized as alternative methods to animal-based toxicity testing. Here, using a panel of 14 chemicals with different risks of human drug-induced liver injury (DILI) and two dosing concentrations, we evaluated an HTTr platform comprised of collagen sandwich primary rat hepatocyte culture and the TempO-Seq surrogate S1500+ (ST) assay. First, the HTTr platform was found to exhibit high reproducibility between technical and biological replicates (r greater than 0.85). Connectivity mapping analysis further demonstrated a high level of inter-platform reproducibility between TempO-Seq data and Affymetrix GeneChip data from the Open TG-GATES project. Second, the TempO-Seq ST assay was shown to be a robust surrogate to the whole transcriptome (WT) assay in capturing chemical-induced changes in gene expression, as evident from correlation analysis, PCA and unsupervised hierarchical clustering. Gene set enrichment analysis (GSEA) using the Hallmark gene set collection also demonstrated consistency in enrichment scores between ST and WT assays. Lastly, unsupervised hierarchical clustering of hallmark enrichment scores broadly divided the samples into hepatotoxic, intermediate, and non-hepatotoxic groups. Xenobiotic metabolism, bile acid metabolism, apoptosis, p53 pathway, and coagulation were found to be the key hallmarks driving the clustering. Taken together, our results established the reproducibility and performance of collagen sandwich culture in combination with TempO-Seq S1500+ assay, and demonstrated the utility of GSEA using the hallmark gene set collection to identify potential hepatotoxicants for further validation.

3D cell culture models: Drug pharmacokinetics, safety assessment, and regulatory consideration

Nonclinical testing has served as a foundation for evaluating potential risks and effectiveness of investigational new drugs in humans. However, the current two-dimensional (2D) in vitro cell culture systems cannot accurately depict and simulate the rich environment and complex processes observed in vivo, whereas animal studies present significant drawbacks with inherited species-specific differences and low throughput for increased demands. To improve the nonclinical prediction of drug safety and efficacy, researchers continue to develop novel models to evaluate and promote the use of improved cell- and organ-based assays for more accurate representation of human susceptibility to drug response. Among others, the three-dimensional (3D) cell culture models present physiologically relevant cellular microenvironment and offer great promise for assessing drug disposition and pharmacokinetics (PKs) that influence drug safety and efficacy from an early stage of drug development. Currently, there are numerous different types of 3D culture systems, from simple spheroids to more complicated organoids and organs-on-chips, and from single-cell type static 3D models to cell co-culture 3D models equipped with microfluidic flow control as well as hybrid 3D systems that combine 2D culture with biomedical microelectromechanical systems. This article reviews the current application and challenges of 3D culture systems in drug PKs, safety, and efficacy assessment, and provides a focused discussion and regulatory perspectives on the liver-, intestine-, kidney-, and neuron-based 3D cellular models.

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

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

Estimating Hepatotoxic Doses Using High-Content Imaging in Primary Hepatocytes

Using in vitro data to estimate point of departure (POD) values is an essential component of new approach methodologies (NAM)-based chemical risk assessments. In this case study, we evaluated a NAM for hepatotoxicity based on rat primary hepatocytes, high-content imaging (HCI), and toxicokinetic modeling. First, we treated rat primary hepatocytes with 10 concentrations (0.2 to 100 µM) of 51 chemicals that produced hepatotoxicity in repeat-dose subchronic and chronic exposures. Second, we used HCI to measure endoplasmic reticulum stress, mitochondrial function, lysosomal mass, steatosis, apoptosis, DNA texture, nuclear size, and cell number at 24, 48, and 72 h and calculated concentrations at 50% maximal activity (AC50). Third, we estimated administered equivalent doses (AEDs) from AC50 values using toxicokinetic modeling. AEDs using physiologically-based toxicokinetic models were 4.1-fold (SD 6.3) and 8.1-fold (SD 15.5) lower than subchronic and chronic lowest observed adverse effect levels (LOAELs), respectively. In contrast, AEDs from ToxCast and Tox21 assays were 89.8-fold (SD 149.5) and 168-fold (SD 323.7) lower than subchronic and chronic LOAELs. Individual HCI end-points also estimated AEDs for specific hepatic lesions that were lower than in vivo PODs. Lastly, AEDs were similar for different in vitro exposure durations, but steady-state toxicokinetic models produced 7.6-fold lower estimates than dynamic physiologically-based ones. Our findings suggest that NAMs from diverse cell types provide conservative estimates of PODs. In contrast, NAMs based on the same species and cell type as the adverse outcome may produce estimates closer to the traditional in vivo PODs.

Benchmark Concentrations for Untargeted Metabolomics Versus Transcriptomics for Liver Injury Compounds in In Vitro Liver Models

Interpretation of untargeted metabolomics data from both in vivo and physiologically relevant in vitro model systems continues to be a significant challenge for toxicology research. Potency-based modeling of toxicological responses has served as a pillar of interpretive context and translation of testing data. In this study, we leverage the resolving power of concentration-response modeling through benchmark concentration (BMC) analysis to interpret untargeted metabolomics data from differentiated cultures of HepaRG cells exposed to a panel of reference compounds and integrate data in a potency-aligned framework with matched transcriptomic data. For this work, we characterized biological responses to classical human liver injury compounds and comparator compounds, known to not cause liver injury in humans, at 10 exposure concentrations in spent culture media by untargeted liquid chromatography-mass spectrometry analysis. The analyte features observed (with limited metabolites identified) were analyzed using BMC modeling to derive compound-induced points of departure. The results revealed liver injury compounds produced concentration-related increases in metabolomic response compared to those rarely associated with liver injury (ie, sucrose, potassium chloride). Moreover, the distributions of altered metabolomic features were largely comparable with those observed using high throughput transcriptomics, which were further extended to investigate the potential for in vitro observed biological responses to be observed in humans with exposures at therapeutic doses. These results demonstrate the utility of BMC modeling of untargeted metabolomics data as a sensitive and quantitative indicator of human liver injury potential.

High-Throughput Transcriptomics Platform for Screening Environmental Chemicals

New approach methodologies (NAMs) that efficiently provide information about chemical hazard without using whole animals are needed to accelerate the pace of chemical risk assessments. Technological advancements in gene expression assays have made in vitro high-throughput transcriptomics (HTTr) a feasible option for NAMs-based hazard characterization of environmental chemicals. In this study, we evaluated the Templated Oligo with Sequencing Readout (TempO-Seq) assay for HTTr concentration-response screening of a small set of chemicals in the human-derived MCF7 cell model. Our experimental design included a variety of reference samples and reference chemical treatments in order to objectively evaluate TempO-Seq assay performance. To facilitate analysis of these data, we developed a robust and scalable bioinformatics pipeline using open-source tools. We also developed a novel gene expression signature-based concentration-response modeling approach and compared the results to a previously implemented workflow for concentration-response analysis of transcriptomics data using BMDExpress. Analysis of reference samples and reference chemical treatments demonstrated highly reproducible differential gene expression signatures. In addition, we found that aggregating signals from individual genes into gene signatures prior to concentration-response modeling yielded in vitro transcriptional biological pathway altering concentrations (BPACs) that were closely aligned with previous ToxCast high-throughput screening assays. Often these identified signatures were associated with the known molecular target of the chemicals in our test set as the most sensitive components of the overall transcriptional response. This work has resulted in a novel and scalable in vitro HTTr workflow that is suitable for high-throughput hazard evaluation of environmental chemicals.

Systematic transcriptome-based comparison of cellular adaptive stress response activation networks in hepatic stem cell-derived progeny and primary human hepatocytes

Various adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte-like cells (HLCs) provide a promising tool to study cellular stress response pathways, but in this context there is limited insight on how HLCs compare to other in vitro liver models. Here, we systematically compared the transcriptomic profiles upon chemical activation in HLCs, hiPSC, primary human hepatocytes (PHH) and HepG2 liver cancer cells. We used targeted RNA-sequencing to map concentration transcriptional response using benchmark concentration modeling for the various stress responses in the different test systems. We found that HLCs are very sensitive towards oxidative stress and inflammation conditions as corresponding genes were activated at over 3 fold lower concentrations of the corresponding pathway inducing compounds as compared to PHH. PHH were the most sensitive model when studying UPR related effects. Due to the non-proliferative nature of PHH and HLCs, these do not pose a good/sensitive model to pick up DNA damage responses, while hiPSC and HepG2 were more sensitive in these conditions. We envision that this study contributes to a better understanding on how HLCs can contribute to the assessment of cell physiological stress response activation to predict hepatotoxic events.

High-throughput toxicogenomic screening of chemicals in the environment using metabolically competent hepatic cell cultures

The ToxCast in vitro screening program has provided concentration-response bioactivity data across more than a thousand assay endpoints for thousands of chemicals found in our environment and commerce. However, most ToxCast screening assays have evaluated individual biological targets in cancer cell lines lacking integrated physiological functionality (such as receptor signaling, metabolism). We evaluated differentiated HepaRGTM cells, a human liver-derived cell model understood to effectively model physiologically relevant hepatic signaling. Expression of 93 gene transcripts was measured by quantitative polymerase chain reaction using Fluidigm 96.96 dynamic arrays in response to 1060 chemicals tested in eight-point concentration-response. A Bayesian framework quantitatively modeled chemical-induced changes in gene expression via six transcription factors including: aryl hydrocarbon receptor, constitutive androstane receptor, pregnane X receptor, farnesoid X receptor, androgen receptor, and peroxisome proliferator-activated receptor alpha. For these chemicals the network model translates transcriptomic data into Bayesian inferences about molecular targets known to activate toxicological adverse outcome pathways. These data also provide new insights into the molecular signaling network of HepaRGTM cell cultures.

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the “measured only” and “whole transcriptome” (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the “measured only” genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original “measured only” dataset. Furthermore, the inclusion of the extrapolated genes raised “tunicamycin” from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

Identification of novel activators of the metal responsive transcription factor (MTF-1) using a gene expression biomarker in a microarray compendium

Environmental exposure to metals is known to cause a number of human toxicities including cancer. Metal-responsive transcription factor 1 (MTF-1) is an important component of metal regulation systems in mammalian cells. Here, we describe a novel method to identify chemicals that activate MTF-1 based on microarray profiling data. MTF-1 biomarker genes were identified that exhibited consistent, robust expression across 10 microarray comparisons examining the effects of metals (zinc, nickel, lead, arsenic, mercury, and silver) on gene expression in human cells. A subset of the resulting 81 biomarker genes was shown to be altered by knockdown of the MTF1 gene including metallothionein family members and a zinc transporter. The ability to correctly identify treatment conditions that activate MTF-1 was determined by comparing the biomarker to microarray comparisons from cells exposed to reference metal activators of MTF-1 using the rank-based Running Fisher algorithm. The balanced accuracy for prediction was 93%. The biomarker was then used to identify organic chemicals that activate MTF-1 from a compendium of 11 725 human gene expression comparisons representing 2582 chemicals. There were 700 chemicals identified that included those known to interact with cellular metals, such as clioquinol and disulfiram, as well as a set of novel chemicals. All nine of the novel chemicals selected for validation were confirmed to activate MTF-1 biomarker genes in MCF-7 cells and to lesser extents in MTF1-null cells by qPCR and targeted RNA-Seq. Overall, our work demonstrates that the biomarker for MTF-1 coupled with the Running Fisher test is a reliable strategy to identify novel chemical modulators of metal homeostasis using gene expression profiling.

Comparison of Approaches for Determining Bioactivity Hits from High-Dimensional Profiling Data

Phenotypic profiling assays are untargeted screening assays that measure a large number (hundreds to thousands) of cellular features in response to a stimulus and often yield diverse and unanticipated profiles of phenotypic effects, leading to challenges in distinguishing active from inactive treatments. Here, we compare a variety of different strategies for hit identification in imaging-based phenotypic profiling assays using a previously published Cell Painting data set. Hit identification strategies based on multiconcentration analysis involve curve fitting at several levels of data aggregation (e.g., individual feature level, aggregation of similarly derived features into categories, and global modeling of all features) and on computed metrics (e.g., Euclidean and Mahalanobis distance metrics and eigenfeatures). Hit identification strategies based on single-concentration analysis included measurement of signal strength (e.g., total effect magnitude) and correlation of profiles among biological replicates. Modeling parameters for each approach were optimized to retain the ability to detect a reference chemical with subtle phenotypic effects while limiting the false-positive rate to 10%. The percentage of test chemicals identified as hits was highest for feature-level and category-based approaches, followed by global fitting, whereas signal strength and profile correlation approaches detected the fewest number of active hits at the fixed false-positive rate. Approaches involving fitting of distance metrics had the lowest likelihood for identifying high-potency false-positive hits that may be associated with assay noise. Most of the methods achieved a 100% hit rate for the reference chemical and high concordance for 82% of test chemicals, indicating that hit calls are robust across different analysis approaches.

Evaluation of 5-day In Vivo Rat Liver and Kidney With High-throughput Transcriptomics for Estimating Benchmark Doses of Apical Outcomes

A 5-day in vivo rat model was evaluated as an approach to estimate chemical exposures that may pose minimal risk by comparing benchmark dose (BMD) values for transcriptional changes in the liver and kidney to BMD values for toxicological endpoints from traditional toxicity studies. Eighteen chemicals, most having been tested by the National Toxicology Program in 2-year bioassays, were evaluated. Some of these chemicals are potent hepatotoxicants (eg, DE71, PFOA, and furan) in rodents, some exhibit toxicity but have minimal hepatic effects (eg, acrylamide and α,β-thujone), and some exhibit little overt toxicity (eg, ginseng and milk thistle extract) based on traditional toxicological evaluations. Male Sprague Dawley rats were exposed once daily for 5 consecutive days by oral gavage to 8-10 dose levels for each chemical. Liver and kidney were collected 24 h after the final exposure and total RNA was assayed using high-throughput transcriptomics (HTT) with the rat S1500+ platform. HTT data were analyzed using BMD Express 2 to determine transcriptional gene set BMD values. BMDS was used to determine BMD values for histopathological effects from chronic or subchronic toxicity studies. For many of the chemicals, the lowest transcriptional BMDs from the 5-day assays were within a factor of 5 of the lowest histopathological BMDs from the toxicity studies. These data suggest that using HTT in a 5-day in vivo rat model provides reasonable estimates of BMD values for traditional apical endpoints. This approach may be useful to prioritize chemicals for further testing while providing actionable data in a timely and cost-effective manner.

Identifying and Characterizing Stress Pathways of Concern for Consumer Safety in Next-Generation Risk Assessment

Many substances for which consumer safety risk assessments need to be conducted are not associated with specific toxicity modes of action, but rather exhibit nonspecific toxicity leading to cell stress. In this work, a cellular stress panel is described, consisting of 36 biomarkers representing mitochondrial toxicity, cell stress, and cell health, measured predominantly using high content imaging. To evaluate the panel, data were generated for 13 substances at exposures consistent with typical use-case scenarios. These included some that have been shown to cause adverse effects in a proportion of exposed humans and have a toxicological mode-of-action associated with cellular stress (eg, doxorubicin, troglitazone, and diclofenac), and some that are not associated with adverse effects due to cellular stress at human-relevant exposures (eg, caffeine, niacinamide, and phenoxyethanol). For each substance, concentration response data were generated for each biomarker at 3 timepoints. A Bayesian model was then developed to quantify the evidence for a biological response, and if present, a credibility range for the estimated point of departure (PoD) was determined. PoDs were compared with the plasma Cmax associated with the typical substance exposures, and indicated a clear differentiation between "low" risk and "high" risk chemical exposure scenarios. Developing robust methods to characterize the in vitro bioactivity of xenobiotics is an important part of non-animal safety assessment. The results presented in this work show that the cellular stress panel can be used, together with other new approach methodologies, to identify chemical exposures that are protective of consumer health.

A Next-Generation Risk Assessment Case Study for Coumarin in Cosmetic Products

Next-Generation Risk Assessment is defined as an exposure-led, hypothesis-driven risk assessment approach that integrates new approach methodologies (NAMs) to assure safety without the use of animal testing. These principles were applied to a hypothetical safety assessment of 0.1% coumarin in face cream and body lotion. For the purpose of evaluating the use of NAMs, existing animal and human data on coumarin were excluded. Internal concentrations (plasma Cmax) were estimated using a physiologically based kinetic model for dermally applied coumarin. Systemic toxicity was assessed using a battery of in vitro NAMs to identify points of departure (PoDs) for a variety of biological effects such as receptor-mediated and immunomodulatory effects (Eurofins SafetyScreen44 and BioMap Diversity 8 Panel, respectively), and general bioactivity (ToxCast data, an in vitro cell stress panel and high-throughput transcriptomics). In addition, in silico alerts for genotoxicity were followed up with the ToxTracker tool. The PoDs from the in vitro assays were plotted against the calculated in vivo exposure to calculate a margin of safety with associated uncertainty. The predicted Cmax values for face cream and body lotion were lower than all PoDs with margin of safety higher than 100. Furthermore, coumarin was not genotoxic, did not bind to any of the 44 receptors tested and did not show any immunomodulatory effects at consumer-relevant exposures. In conclusion, this case study demonstrated the value of integrating exposure science, computational modeling and in vitro bioactivity data, to reach a safety decision without animal data.

Comparison of Normalization Methods for Analysis of TempO-Seq Targeted RNA Sequencing Data

Analysis of bulk RNA sequencing (RNA-Seq) data is a valuable tool to understand transcription at the genome scale. Targeted sequencing of RNA has emerged as a practical means of assessing the majority of the transcriptomic space with less reliance on large resources for consumables and bioinformatics. TempO-Seq is a templated, multiplexed RNA-Seq platform that interrogates a panel of sentinel genes representative of genome-wide transcription. Nuances of the technology require proper preprocessing of the data. Various methods have been proposed and compared for normalizing bulk RNA-Seq data, but there has been little to no investigation of how the methods perform on TempO-Seq data. We simulated count data into two groups (treated vs. untreated) at seven-fold change (FC) levels (including no change) using control samples from human HepaRG cells run on TempO-Seq and normalized the data using seven normalization methods. Upper Quartile (UQ) performed the best with regard to maintaining FC levels as detected by a limma contrast between treated vs. untreated groups. For all FC levels, specificity of the UQ normalization was greater than 0.84 and sensitivity greater than 0.90 except for the no change and +1.5 levels. Furthermore, K-means clustering of the simulated genes normalized by UQ agreed the most with the FC assignments [adjusted Rand index (ARI) = 0.67]. Despite having an assumption of the majority of genes being unchanged, the DESeq2 scaling factors normalization method performed reasonably well as did simple normalization procedures counts per million (CPM) and total counts (TCs). These results suggest that for two class comparisons of TempO-Seq data, UQ, CPM, TC, or DESeq2 normalization should provide reasonably reliable results at absolute FC levels ≥2.0. These findings will help guide researchers to normalize TempO-Seq gene expression data for more reliable results.

Phenotypic Profiling of Reference Chemicals across Biologically Diverse Cell Types Using the Cell Painting Assay

Cell Painting is a high-throughput phenotypic profiling assay that uses fluorescent cytochemistry to visualize a variety of organelles and high-content imaging to derive a large number of morphological features at the single-cell level. Most Cell Painting studies have used the U-2 OS cell line for chemical or functional genomics screening. The Cell Painting assay can be used with many other human-derived cell types, given that the assay is based on the use of fluoroprobes that label organelles that are present in most (if not all) human cells. Questions remain, however, regarding the optimization steps required and overall ease of deployment of the Cell Painting assay to novel cell types. Here, we used the Cell Painting assay to characterize the phenotypic effects of 14 phenotypic reference chemicals in concentration-response screening mode across six biologically diverse human-derived cell lines (U-2 OS, MCF7, HepG2, A549, HTB-9 and ARPE-19). All cell lines were labeled using the same cytochemistry protocol, and the same set of phenotypic features was calculated. We found it necessary to optimize image acquisition settings and cell segmentation parameters for each cell type, but did not adjust the cytochemistry protocol. For some reference chemicals, similar subsets of phenotypic features corresponding to a particular organelle were associated with the highest-effect magnitudes in each affected cell type. Overall, for certain chemicals, the Cell Painting assay yielded qualitatively similar biological activity profiles among a group of diverse, morphologically distinct human-derived cell lines without the requirement for cell type-specific optimization of cytochemistry protocols.

Transcriptomics in Toxicogenomics, Part I: Experimental Design, Technologies, Publicly Available Data, and Regulatory Aspects

The starting point of successful hazard assessment is the generation of unbiased and trustworthy data. Conventional toxicity testing deals with extensive observations of phenotypic endpoints in vivo and complementing in vitro models. The increasing development of novel materials and chemical compounds dictates the need for a better understanding of the molecular changes occurring in exposed biological systems. Transcriptomics enables the exploration of organisms' responses to environmental, chemical, and physical agents by observing the molecular alterations in more detail. Toxicogenomics integrates classical toxicology with omics assays, thus allowing the characterization of the mechanism of action (MOA) of chemical compounds, novel small molecules, and engineered nanomaterials (ENMs). Lack of standardization in data generation and analysis currently hampers the full exploitation of toxicogenomics-based evidence in risk assessment. To fill this gap, TGx methods need to take into account appropriate experimental design and possible pitfalls in the transcriptomic analyses as well as data generation and sharing that adhere to the FAIR (Findable, Accessible, Interoperable, and Reusable) principles. In this review, we summarize the recent advancements in the design and analysis of DNA microarray, RNA sequencing (RNA-Seq), and single-cell RNA-Seq (scRNA-Seq) data. We provide guidelines on exposure time, dose and complex endpoint selection, sample quality considerations and sample randomization. Furthermore, we summarize publicly available data resources and highlight applications of TGx data to understand and predict chemical toxicity potential. Additionally, we discuss the efforts to implement TGx into regulatory decision making to promote alternative methods for risk assessment and to support the 3R (reduction, refinement, and replacement) concept. This review is the first part of a three-article series on Transcriptomics in Toxicogenomics. These initial considerations on Experimental Design, Technologies, Publicly Available Data, Regulatory Aspects, are the starting point for further rigorous and reliable data preprocessing and modeling, described in the second and third part of the review series.

NanoSolveIT Project: Driving nanoinformatics research to develop innovative and integrated tools for in silico nanosafety assessment

Nanotechnology has enabled the discovery of a multitude of novel materials exhibiting unique physicochemical (PChem) properties compared to their bulk analogues. These properties have led to a rapidly increasing range of commercial applications; this, however, may come at a cost, if an association to long-term health and environmental risks is discovered or even just perceived. Many nanomaterials (NMs) have not yet had their potential adverse biological effects fully assessed, due to costs and time constraints associated with the experimental assessment, frequently involving animals. Here, the available NM libraries are analyzed for their suitability for integration with novel nanoinformatics approaches and for the development of NM specific Integrated Approaches to Testing and Assessment (IATA) for human and environmental risk assessment, all within the NanoSolveIT cloud-platform. These established and well-characterized NM libraries (e.g. NanoMILE, NanoSolutions, NANoREG, NanoFASE, caLIBRAte, NanoTEST and the Nanomaterial Registry (>2000 NMs)) contain physicochemical characterization data as well as data for several relevant biological endpoints, assessed in part using harmonized Organisation for Economic Co-operation and Development (OECD) methods and test guidelines. Integration of such extensive NM information sources with the latest nanoinformatics methods will allow NanoSolveIT to model the relationships between NM structure (morphology), properties and their adverse effects and to predict the effects of other NMs for which less data is available. The project specifically addresses the needs of regulatory agencies and industry to effectively and rapidly evaluate the exposure, NM hazard and risk from nanomaterials and nano-enabled products, enabling implementation of computational 'safe-by-design' approaches to facilitate NM commercialization.

Pathway-based assessment of single chemicals and mixtures by a high-throughput transcriptomics approach

The ever-increasing number of chemicals and complex mixtures demands a high-throughput and cost-effective approach for chemical safety assessment. High-throughput transcriptomics (HTT) is promising in investigating genome-scale perturbation of chemical exposure in concentration-dependent manner. However, the application of HTT has been limited due to lack of methodology for single chemicals and mixture assessment. This study aimed to evaluate the ability of a newly-developed human reduced transcriptomics (RHT) approach to assess pathway-based profiles of single chemicals, and to develop a biological pathway-based approach for benchmarking mixture potency using single chemical-based prediction model. First, concentration-dependent RHT were used to qualitatively and quantitatively differentiate pathway-based patterns of different chemicals, using three model toxicants, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), triclosan (TCS) and 5-Chloro-6-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (5-Cl-6-OH-BDE-47). AHR-regulated genes and pathways were most sensitively induced by TCDD, while TCS and 5-Cl-6-OH-BDE-47 were much less potent in AHR-associated activation, which was concordant with known MoA of each single chemical. Second, two artificial mixtures and their components of twelve individual chemicals were performed with concentration-dependent RHT. Concentration addition (CA) and independent action (IA) models were used to predict transcriptional potency of mixtures from transcriptomics of individual chemicals. For overall bioactivity, CA and IA models can both predict potency of observed responses within 95% confidence interval. For specific biological processes, multiple biological processes such as hormone signaling and DNA damage can be predicted using CA models for mixtures. The concentration-dependent RHT can provide a powerful approach for qualitative and quantitative assessment of biological pathway perturbated by environment chemical and mixtures.

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.

Improving natural product research translation: From source to clinical trial

While great interest in health effects of natural product (NP) including dietary supplements and foods persists, promising preclinical NP research is not consistently translating into actionable clinical trial (CT) outcomes. Generally considered the gold standard for assessing safety and efficacy, CTs, especially phase III CTs, are costly and require rigorous planning to optimize the value of the information obtained. More effective bridging from NP research to CT was the goal of a September, 2018 transdisciplinary workshop. Participants emphasized that replicability and likelihood of successful translation depend on rigor in experimental design, interpretation, and reporting across the continuum of NP research. Discussions spanned good practices for NP characterization and quality control; use and interpretation of models (computational through in vivo) with strong clinical predictive validity; controls for experimental artefacts, especially for in vitro interrogation of bioactivity and mechanisms of action; rigorous assessment and interpretation of prior research; transparency in all reporting; and prioritization of research questions. Natural product clinical trials prioritized based on rigorous, convergent supporting data and current public health needs are most likely to be informative and ultimately affect public health. Thoughtful, coordinated implementation of these practices should enhance the knowledge gained from future NP research.

Structure-activity relationship study: Mechanism of cyto-genotoxicity of Nitropyrazole-derived high energy density materials family

Stringent toxicological tests have to be performed prior to the industrial development of alternative chemicals particularly high energy dense materials (HEDMs) such as explosives. The properties (e.g., power, stability) of these compounds are constantly being improved, the current axis of research being the nitration of nitrogen heterocycles leading to HEDMs such as nitropyrazole-derived molecules. However, except for 3,4,5-trinitropyrazole (3,4,5-TNP), which was shown to be highly toxic in mice, the toxicological impact of these HEDMs has so far not been investigated. Furthermore, as industrials are strongly advised to develop alternative safety testing assays to in vivo experiments, we herein focused on determining the cytotoxic and genotoxic effects of seven Nitropyrazole-derived HEDMs on three rodent cell lines (mouse embryonic BALB/3T3 clone A31 cells, Chinese hamster ovary cells CHO-K1 and mouse lymphoma L5178Y TK +/- clone (3.7.2C) cells), two human fibroblast lines (CRC05, PFS04062) and on the human hepatic HepaRG model (both in proliferative and differentiated cells). A stronger cytotoxic effect was observed for 1,3-dinitropyrazole (1, 3-DNP) and 3,4,5-TNP in all cell lines, though differentiated HepaRG cells clearly displayed fewer likely due to the metabolism and elimination of these molecules by their functional biotransformation pathways. At the mechanistic level, the sub-chronic cytotoxic and genotoxic effects were linked to ROS/RNS production (experimental assays), HA2.X and to transcriptomic data highlighting the increase in DNA repair mechanisms.

Gene expression and cytosine DNA methylation alterations in induced pluripotent stem-cell-derived human hepatocytes treated with low doses of chemical carcinogens

The increasing number of man-made chemicals in the environment that may pose a carcinogenic risk emphasizes the need to develop reliable time- and cost-effective approaches for carcinogen detection. To address this issue, we have investigated the utility of human hepatocytes for the in vitro identification of genotoxic and non-genotoxic carcinogens. Induced pluripotent stem-cell (iPSC)-derived human hepatocytes were treated with the genotoxic carcinogens aflatoxin B₁ (AFB1) and benzo[a]pyrene (B[a]P), the non-genotoxic liver carcinogen methapyrilene, and the non-carcinogens aflatoxin B₂ (AFB2) and benzo[e]pyrene (B[e]P) at non-cytotoxic concentrations for 7 days, and transcriptomic and DNA methylation profiles were examined. 1569, 1693, and 2061 differentially expressed genes (DEGs) were detected in cells treated with AFB1, B[a]P, and methapyrilene, respectively, whereas no DEGs were found in cells treated with AFB2 or B[e]P. In contrast to the profound cellular transcriptomic responses, exposure of iPSC-derived hepatocytes to the test chemicals resulted in minor random alterations in global DNA methylome, most of which were not associated with changes in gene expression. Overall, our results demonstrate that the major non-genotoxic effect of exposure to carcinogens, regardless of their mode of action, is a profound global transcriptomic response rather than global DNA methylome alterations, indicating the significance of transcriptomic alterations as an informative endpoint in short-term in vitro carcinogen testing.

Development of a Zebrafish S1500+ Sentinel Gene Set for High-Throughput Transcriptomics

Sentinel gene sets have been developed with the purpose of maximizing the information from targeted transcriptomic platforms. We recently described the development of an S1500+ sentinel gene set, which was built for the human transcriptome, utilizing a data- and knowledge-driven hybrid approach to select a small subset of genes that optimally capture transcriptional diversity, correlation with other genes based on large-scale expression profiling, and known pathway annotation within the human genome. While this detailed bioinformatics approach for gene selection can in principle be applied to other species, the reliability of the resulting gene set depends on availability of a large body of transcriptomics data. For the model organism zebrafish, we aimed to create a similar sentinel gene set (Zf S1500+ gene set); however, there is insufficient standardized expression data in the public domain to train the gene correlation model. Therefore, our strategy was to use human-zebrafish ortholog mapping of the human S1500+ genes and nominations from experts in the zebrafish scientific community. In this study, we present the bioinformatics curation and refinement process to produce the final Zf S1500+ gene set, explore whole transcriptome extrapolation using this gene set, and assess pathway-level inference. This gene set will add value to targeted high-throughput transcriptomics in zebrafish for toxicogenomic screening and other research domains.