Determination of benchmark concentrations and their statistical uncertainty for cytotoxicity test data and functional in vitro assays

Integrating gene expression and splicing dynamics across dose-response oxidative modulators

Toxicological risk assessment increasingly utilizes transcriptomics to derive point of departure (POD) and modes of action (MOA) for chemicals. One essential biological process that allows a single gene to generate several different RNA isoforms is called alternative splicing. To comprehensively assess the role of splicing dysregulation in toxicological evaluation and elucidate its potential as a complementary endpoint, we performed RNA-seq on A549 cells treated with five oxidative stress modulators across a wide dose range. Differential gene expression (DGE) showed limited pathway enrichment except at high concentrations. However, alternative splicing analysis revealed variable intron retention events affecting diverse pathways for all chemicals in the absence of significant expression changes. For instance, diazinon elicited negligible gene expression changes but progressive increase in the number of intron retention events, suggesting splicing alterations precede expression responses. Benchmark dose modeling of intron retention data highlighted relevant pathways overlooked by expression analysis. Systematic integration of splicing datasets should be a useful addition to the toxicogenomic toolkit. Combining both modalities paint a more complete picture of transcriptomic dose-responses. Overall, evaluating intron retention dynamics afforded by toxicogenomics may provide biomarkers that can enhance chemical risk assessment and regulatory decision making. This work highlights splicing-aware toxicogenomics as a possible additional tool for examining cellular responses.

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.

Guidance for statistical design and analysis of toxicological dose-response experiments, based on a comprehensive literature review

The analysis of dose-response, concentration-response, and time-response relationships is a central component of toxicological research. A major decision with respect to the statistical analysis is whether to consider only the actually measured concentrations or to assume an underlying (parametric) model that allows extrapolation. Recent research suggests the application of modelling approaches for various types of toxicological assays. However, there is a discrepancy between the state of the art in statistical methodological research and published analyses in the toxicological literature. The extent of this gap is quantified in this work using an extensive literature review that considered all dose-response analyses published in three major toxicological journals in 2021. The aspects of the review include biological considerations (type of assay and of exposure), statistical design considerations (number of measured conditions, design, and sample sizes), and statistical analysis considerations (display, analysis goal, statistical testing or modelling method, and alert concentration). Based on the results of this review and the critical assessment of three selected issues in the context of statistical research, concrete guidance for planning, execution, and analysis of dose-response studies from a statistical viewpoint is proposed.

The Human Induced Pluripotent Stem Cell Test as an Alternative Method for Embryotoxicity Testing

The evaluation of substances for their potency to induce embryotoxicity is controlled by safety regulations. Test guidelines for reproductive and developmental toxicity rely mainly on animal studies, which make up the majority of animal usage in regulatory toxicology. Therefore, there is an urgent need for alternative in vitro methods to follow the 3R principles. To improve human safety, cell models based on human cells are of great interest to overcome species differences. Here, human induced pluripotent stem cells (hiPSCs) are an ideal cell source as they largely recapitulate embryonic stem cells without bearing ethical concerns and they are able to differentiate into most cell types of the human body. Here, we set up and characterized a fetal bovine serum (FBS)-free hiPSC-based in vitro test method, called the human induced pluripotent stem cell test (hiPS Test), to evaluate the embryotoxic potential of substances. After 10 days in culture, hiPSCs develop into beating cardiomyocytes. As terminal endpoint evaluations, cell viability, qPCR analyses as well as beating frequency and area of beating cardiomyocytes by video analyses are measured. The embryotoxic positive and non-embryotoxic negative controls, 5-Fluorouracil (5-FU) and Penicillin G (PenG), respectively, were correctly assessed in the hiPS Test. More compounds need to be screened in the future for defining the assay’s applicability domain, which will inform us of the suitability of the hiPS Test for detecting adverse effects of substances on embryonic development.

Acute effects of the imidacloprid metabolite desnitro-imidacloprid on human nACh receptors relevant for neuronal signaling

Several neonicotinoids have recently been shown to activate the nicotinic acetylcholine receptor (nAChR) on human neurons. Moreover, imidacloprid (IMI) and other members of this pesticide family form a set of diverse metabolites within crops. Among these, desnitro-imidacloprid (DN-IMI) is of special toxicological interest, as there is evidence (i) for human dietary exposure to this metabolite, (ii) and that DN-IMI is a strong trigger of mammalian nicotinic responses. We set out here to quantify responses of human nAChRs to DN-IMI and an alternative metabolite, IMI-olefin. To evaluate toxicological hazards, these data were then compared to those of IMI and nicotine. Ca²⁺-imaging experiments on human neurons showed that DN-IMI exhibits an agonistic effect on nAChRs at sub-micromolar concentrations (equipotent with nicotine) while IMI-olefin activated the receptors less potently (in a similar range as IMI). Direct experimental data on the interaction with defined receptor subtypes were obtained by heterologous expression of various human nAChR subtypes in Xenopus laevis oocytes and measurement of the transmembrane currents evoked by exposure to putative ligands. DN-IMI acted on the physiologically important human nAChR subtypes α7, α3β4, and α4β2 (high-sensitivity variant) with similar potency as nicotine. IMI and IMI-olefin were confirmed as nAChR agonists, although with 2–3 orders of magnitude lower potency. Molecular docking studies, using receptor models for the α7 and α4β2 nAChR subtypes supported an activity of DN-IMI similar to that of nicotine. In summary, these data suggest that DN-IMI functionally affects human neurons similar to the well-established neurotoxicant nicotine by triggering α7 and several non-α7 nAChRs.

Neurodevelopmental toxicity assessment of flame retardants using a human DNT in vitro testing battery

Due to their neurodevelopmental toxicity, flame retardants (FRs) like polybrominated diphenyl ethers are banned from the market and replaced by alternative FRs, like organophosphorus FRs, that have mostly unknown toxicological profiles. To study their neurodevelopmental toxicity, we evaluated the hazard of several FRs including phased-out polybrominated FRs and organophosphorus FRs: 2,2′,4,4′-tetrabromodiphenylether (BDE-47), 2,2′,4,4′,5-pentabromodiphenylether (BDE-99), tetrabromobisphenol A, triphenyl phosphate, tris(2-butoxyethyl) phosphate and its metabolite bis-(2-butoxyethyl) phosphate, isodecyl diphenyl phosphate, triphenyl isopropylated phosphate, tricresyl phosphate, tris(1,3-dichloro-2-propyl) phosphate, tert-butylphenyl diphenyl phosphate, 2-ethylhexyl diphenyl phosphate, tris(1-chloroisopropyl) phosphate, and tris(2-chloroethyl) phosphate. Therefore, we used a human cell–based developmental neurotoxicity (DNT) in vitro battery covering a large variety of neurodevelopmental endpoints. Potency according to the respective most sensitive benchmark concentration (BMC) across the battery ranked from <1 μM (5 FRs), 1<10 μM (7 FRs) to the >10 μM range (3 FRs). Evaluation of the data with the ToxPi tool revealed a distinct ranking (a) than with the BMC and (b) compared to the ToxCast data, suggesting that DNT hazard of these FRs is not well predicted by ToxCast assays. Extrapolating the DNT in vitro battery BMCs to human FR exposure via breast milk suggests low risk for individual compounds. However, it raises a potential concern for real-life mixture exposure, especially when different compounds converge through diverse modes-of-action on common endpoints, like oligodendrocyte differentiation in this study. This case study using FRs suggests that human cell–based DNT in vitro battery is a promising approach for neurodevelopmental hazard assessment and compound prioritization in risk assessment.

Functional alterations by a subgroup of neonicotinoid pesticides in human dopaminergic neurons

Neonicotinoid pesticides, originally developed to target the insect nervous system, have been reported to interact with human receptors and to activate rodent neurons. Therefore, we evaluated in how far these compounds may trigger signaling in human neurons, and thus, affect the human adult or developing nervous system. We used SH-SY5Y neuroblastoma cells as established model of nicotinic acetylcholine receptor (nAChR) signaling. In parallel, we profiled dopaminergic neurons, generated from LUHMES neuronal precursor cells, as novel system to study nAChR activation in human post-mitotic neurons. Changes of the free intracellular Ca²⁺ concentration ([Ca²⁺]_i) were used as readout, and key findings were confirmed by patch clamp recordings. Nicotine triggered typical neuronal signaling responses that were blocked by antagonists, such as tubocurarine and mecamylamine. Pharmacological approaches suggested a functional expression of α7 and non-α7 nAChRs on LUHMES cells. In this novel test system, the neonicotinoids acetamiprid, imidacloprid, clothianidin and thiacloprid, but not thiamethoxam and dinotefuran, triggered [Ca²⁺]_i signaling at 10-100 µM. Strong synergy of the active neonicotinoids (at low micromolar concentrations) with the α7 nAChR-positive allosteric modulator PNU-120596 was observed in LUHMES and SH-SY5Y cells, and specific antagonists fully inhibited such signaling. To provide a third line of evidence for neonicotinoid signaling via nAChR, we studied cross-desensitization: pretreatment of LUHMES and SH-SY5Y cells with active neonicotinoids (at 1-10 µM) blunted the signaling response of nicotine. The pesticides (at 3-30 µM) also blunted the response to the non-α7 agonist ABT 594 in LUHMES cells. These data show that human neuronal cells are functionally affected by low micromolar concentrations of several neonicotinoids. An effect of such signals on nervous system development is a toxicological concern.

Neurotoxicity and underlying cellular changes of 21 mitochondrial respiratory chain inhibitors

Inhibition of complex I of the mitochondrial respiratory chain (cI) by rotenone and methyl-phenylpyridinium (MPP +) leads to the degeneration of dopaminergic neurons in man and rodents. To formally describe this mechanism of toxicity, an adverse outcome pathway (AOP:3) has been developed that implies that any inhibitor of cI, or possibly of other parts of the respiratory chain, would have the potential to trigger parkinsonian motor deficits. We used here 21 pesticides, all of which are described in the literature as mitochondrial inhibitors, to study the general applicability of AOP:3 or of in vitro assays that are assessing its activation. Five cI, three complex II (cII), and five complex III (cIII) inhibitors were characterized in detail in human dopaminergic neuronal cell cultures. The NeuriTox assay, examining neurite damage in LUHMES cells, was used as in vitro proxy of the adverse outcome (AO), i.e., of dopaminergic neurodegeneration. This test provided data on whether test compounds were unspecific cytotoxicants or specifically neurotoxic, and it yielded potency data with respect to neurite degeneration. The pesticide panel was also examined in assays for the sequential key events (KE) leading to the AO, i.e., mitochondrial respiratory chain inhibition, mitochondrial dysfunction, and disturbed proteostasis. Data from KE assays were compared to the NeuriTox data (AO). The cII-inhibitory pesticides tested here did not appear to trigger the AOP:3 at all. Some of the cI/cIII inhibitors showed a consistent AOP activation response in all assays, while others did not. In general, there was a clear hierarchy of assay sensitivity: changes of gene expression (biomarker of neuronal stress) correlated well with NeuriTox data; mitochondrial failure (measured both by a mitochondrial membrane potential-sensitive dye and a respirometric assay) was about 10–260 times more sensitive than neurite damage (AO); cI/cIII activity was sometimes affected at > 1000 times lower concentrations than the neurites. These data suggest that the use of AOP:3 for hazard assessment has a number of caveats: (i) specific parkinsonian neurodegeneration cannot be easily predicted from assays of mitochondrial dysfunction; (ii) deriving a point-of-departure for risk assessment from early KE assays may overestimate toxicant potency.

Comparison of 21 data-rich mitochondrial toxicants for neurotoxicity

Quantitative comparison of key event triggering thresholds for AOP:3

Comparison of two cell models and two exposure times for neurotoxicity

Comparison of transcriptome changes and classical key event measures for sensitivity

Setting the stage for next-generation risk assessment with non-animal approaches: the EU-ToxRisk project experience

In 2016, the European Commission launched the EU-ToxRisk research project to develop and promote animal-free approaches in toxicology. The 36 partners of this consortium used in vitro and in silico methods in the context of case studies (CSs). These CSs included both compounds with a highly defined target (e.g. mitochondrial respiratory chain inhibitors) as well as compounds with poorly defined molecular initiation events (e.g. short-chain branched carboxylic acids). The initial project focus was on developing a science-based strategy for read-across (RAx) as an animal-free approach in chemical risk assessment. Moreover, seamless incorporation of new approach method (NAM) data into this process (= NAM-enhanced RAx) was explored. Here, the EU-ToxRisk consortium has collated its scientific and regulatory learnings from this particular project objective. For all CSs, a mechanistic hypothesis (in the form of an adverse outcome pathway) guided the safety evaluation. ADME data were generated from NAMs and used for comprehensive physiological-based kinetic modelling. Quality assurance and data management were optimized in parallel. Scientific and Regulatory Advisory Boards played a vital role in assessing the practical applicability of the new approaches. In a next step, external stakeholders evaluated the usefulness of NAMs in the context of RAx CSs for regulatory acceptance. For instance, the CSs were included in the OECD CS portfolio for the Integrated Approach to Testing and Assessment project. Feedback from regulators and other stakeholders was collected at several stages. Future chemical safety science projects can draw from this experience to implement systems toxicology-guided, animal-free next-generation risk assessment.

Multiparametric assessment of mitochondrial respiratory inhibition in HepG2 and RPTEC/TERT1 cells using a panel of mitochondrial targeting agrochemicals

Evidence is mounting for the central role of mitochondrial dysfunction in several pathologies including metabolic diseases, accelerated ageing, neurodegenerative diseases and in certain xenobiotic-induced organ toxicity. Assessing mitochondrial perturbations is not trivial and the outcomes of such investigations are dependent on the cell types used and assays employed. Here we systematically investigated the effect of electron transport chain (ETC) inhibitors on multiple mitochondrial-related parameters in two human cell types, HepG2 and RPTEC/TERT1. Cells were exposed to a broad range of concentrations of 20 ETC-inhibiting agrochemicals and capsaicin, consisting of inhibitors of NADH dehydrogenase (Complex I, CI), succinate dehydrogenase (Complex II, CII) and cytochrome bc1 complex (Complex III, CIII). A battery of tests was utilised, including viability assays, lactate production, mitochondrial membrane potential (MMP) and the Seahorse bioanalyser, which simultaneously measures extracellular acidification rate [ECAR] and oxygen consumption rate [OCR]. CI inhibitors caused a potent decrease in OCR, decreased mitochondrial membrane potential, increased ECAR and increased lactate production in both cell types. Twenty-fourhour exposure to CI inhibitors decreased viability of RPTEC/TERT1 cells and 3D spheroid-cultured HepG2 cells in the presence of glucose. CI inhibitors decreased 2D HepG2 viability only in the absence of glucose. CII inhibitors had no notable effects in intact cells up to 10 µM. CIII inhibitors had similar effects to the CI inhibitors. Antimycin A was the most potent CIII inhibitor, with activity in the nanomolar range. The proposed CIII inhibitor cyazofamid demonstrated a mitochondrial uncoupling signal in both cell types. The study presents a comprehensive example of a mitochondrial assessment workflow and establishes measurable key events of ETC inhibition.

The EU-ToxRisk method documentation, data processing and chemical testing pipeline for the regulatory use of new approach methods

Hazard assessment, based on new approach methods (NAM), requires the use of batteries of assays, where individual tests may be contributed by different laboratories. A unified strategy for such collaborative testing is presented. It details all procedures required to allow test information to be usable for integrated hazard assessment, strategic project decisions and/or for regulatory purposes. The EU-ToxRisk project developed a strategy to provide regulatorily valid data, and exemplified this using a panel of > 20 assays (with > 50 individual endpoints), each exposed to 19 well-known test compounds (e.g. rotenone, colchicine, mercury, paracetamol, rifampicine, paraquat, taxol). Examples of strategy implementation are provided for all aspects required to ensure data validity: (i) documentation of test methods in a publicly accessible database; (ii) deposition of standard operating procedures (SOP) at the European Union DB-ALM repository; (iii) test readiness scoring accoding to defined criteria; (iv) disclosure of the pipeline for data processing; (v) link of uncertainty measures and metadata to the data; (vi) definition of test chemicals, their handling and their behavior in test media; (vii) specification of the test purpose and overall evaluation plans. Moreover, data generation was exemplified by providing results from 25 reporter assays. A complete evaluation of the entire test battery will be described elsewhere. A major learning from the retrospective analysis of this large testing project was the need for thorough definitions of the above strategy aspects, ideally in form of a study pre-registration, to allow adequate interpretation of the data and to ensure overall scientific/toxicological validity.