Vol.:(0123456789)1 3 Archives of Toxicology (2023) 97:1701–1721 https://doi.org/10.1007/s00204-023-03490-8 ORGAN TOXICITY AND MECHANISMS In vitro to in vivo extrapolation and high‑content imaging for simultaneous characterization of chemically induced liver steatosis and markers of hepatotoxicity Fabrice A. Müller1  · Marianna Stamou1  · Felix H. Englert1 · Ole Frenzel1 · Sabine Diedrich1 · Laura Suter‑Dick3,4  · John F. Wambaugh2  · Shana J. Sturla1 Received: 21 January 2023 / Accepted: 21 March 2023 / Published online: 12 April 2023 © The Author(s) 2023 Abstract Chemically induced steatosis is characterized by lipid accumulation associated with mitochondrial dysfunction, oxidative stress and nucleus distortion. New approach methods integrating in vitro and in silico models are needed to identify chemicals that may induce these cellular events as potential risk factors for steatosis and associated hepatotoxicity. In this study we used high-content imaging for the simultaneous quantification of four cellular markers as sentinels for hepatotoxicity and steatosis in chemically exposed human liver cells in vitro. Furthermore, we evaluated the results with a computational model for the extrapolation of human oral equivalent doses (OED). First, we tested 16 reference chemicals with known capacities to induce cellular alterations in nuclear morphology, lipid accumulation, mitochondrial membrane potential and oxidative stress. Then, using physiologically based pharmacokinetic modeling and reverse dosimetry, OEDs were extrapolated from data of any stimulated individual sentinel response. The extrapolated OEDs were confirmed to be within biologically relevant exposure ranges for the reference chemicals. Next, we tested 14 chemicals found in food, selected from thousands of puta- tive chemicals on the basis of structure-based prediction for nuclear receptor activation. Amongst these, orotic acid had an extrapolated OED overlapping with realistic exposure ranges. Thus, we were able to characterize known steatosis-inducing chemicals as well as data-scarce food-related chemicals, amongst which we confirmed orotic acid to induce hepatotoxicity. This strategy addresses needs of next generation risk assessment and can be used as a first chemical prioritization hazard screening step in a tiered approach to identify chemical risk factors for steatosis and hepatotoxicity-associated events. Keywords New approach method · High-content imaging · Hepatotoxicity · In vitro to in vivo extrapolation Introduction Steatosis is a hepatic accumulation of fatty acids in greater than 5% of hepatocytes (Chalasani et al. 2012; Musso et al. 2016). This cellular phenotype is observed in nonalco- holic fatty liver disease (NAFLD), a complex spectrum of diseases that affects around 25% of the adult population (Chalasani et al. 2012; Younossi et al. 2018) and ranges from benign liver steatosis to nonalcoholic steatohepatitis (NASH), leading to cirrhosis and eventually hepatocel- lular carcinoma. In addition to metabolic disorders (obe- sity, type 2 diabetes) and genetic factors, environmental exposure to chemicals can lead to hepatic steatosis (Kaiser et al. 2012; Wahlang et al. 2013). Also, drugs, such as ami- odarone, valproic acid, and cancer chemotherapeutics can cause steatosis in some patients, particularly after chronic therapy (Begriche et al. 2011; Kaiser et al. 2012; Wahlang * Shana J. Sturla sturlas@ethz.ch 1 Department of Health Sciences and Technology, ETH Zurich, Schmelzbergstrasse 9, 8092 Zurich, Switzerland 2 Center for Computational Toxicology and Exposure, Office of Research and Development, United States Environmental Protection Agency, Research Triangle Park, Durham, NC 27711, USA 3 School of Life Sciences, University of Applied Sciences and Arts Northwestern Switzerland, 4132 Muttenz, Switzerland 4 Swiss Centre for Applied Human Toxicology (SCAHT), 4001 Basel, Switzerland http://crossmark.crossref.org/dialog/?doi=10.1007/s00204-023-03490-8&domain=pdf http://orcid.org/0000-0002-5111-5583 http://orcid.org/0000-0003-0496-3229 http://orcid.org/0000-0002-1449-3913 http://orcid.org/0000-0002-4024-534X http://orcid.org/0000-0001-6808-5950 1702 Archives of Toxicology (2023) 97:1701–1721 1 3 et al. 2013; Jennings et al. 2014; Willebrords et al. 2015; Schumacher and Guo 2015). Chemically induced liver steatosis appears to induce fatty acid synthesis, decrease fatty acid β-oxidation, decrease lipoprotein export and increase uptake of free fatty acids (Pavlik et al. 2019). Yet, data derived from animal models are poorly trans- lated to humans due to interspecies differences concerning molecular mechanism of steatosis (Soret et al. 2020). The adverse outcome pathway (AOP) for chemically induced liver steatosis involves nuclear receptor binding, microsomal triglyceride transfer protein (MTP) or carni- tine palmitoyltransferase I (CPT-1) inhibition, coenzyme A binding and carnitine depletion as the molecular initi- ating events that lead to triglyceride accumulation (Allen et al. 2014; Vinken 2015). Activation of several nuclear receptors, such as ER, LXR, PXR or FXR, can lead to mitochondrial dysfunction via inhibition of β-oxidation and an increase of de novo fatty acid synthesis, creating an imbalance in fatty acid metabolism leading to the accu- mulation of liver triglycerides and eventually steatosis. Moreover, disruption of β-oxidation induces increased production of reactive oxygen species leading to oxida- tive stress (Masarone et al. 2018). Furthermore, ongoing accumulation of liver triglycerides can distort the nucleus, exacerbate mitochondrial dysfunction, oxidative stress and cause endoplasmic reticulum stress. (Moya et al. 2010; Mellor et al. 2015). Chemicals that give rise to fatty liver disease may invoke several of these cellular events as a basis for identification. There are several recent studies concerning chemically induced steatosis in vitro (Anthérieu et al. 2011, 2012; Rogue et al. 2014; Klein et al. 2016; Tolosa et al. 2016; Angrish et al. 2017; Cuykx et al. 2018; Luckert et al. 2018; Bucher et al. 2018; Allard et al. 2020; Lichtenstein et al. 2020; Lasch et al. 2021). These include at least four studies in which key events along the liver steatosis AOP were quan- tified in vitro (Supplementary table 1) (Donato et al. 2012; Tolosa et al. 2016; Luckert et al. 2018; Shah et al. 2021). Most recently, Shah et al. quantified mitochondrial function, lipid accumulation, endoplasmic reticulum stress, lysosomal mass, DNA texture, nuclear size, apoptosis, and cell number using a high-content imaging approach with rat hepatocytes and 51 hepatotoxicants. Moreover, these in vitro data were extrapolated to in vivo administered equivalent doses and compared to in vivo (rat) data (Shah et al. 2021). Relat- ing in vitro bioactivity data to exposure values by means of physiologically based pharmacokinetic modeling (PBPK) and reverse dosimetry is critical in order that in vitro bioas- says can be accepted for chemical hazard and eventually risk assessment by regulatory agencies (Wetmore 2015; Hartung 2018). To our knowledge, there are as yet no reports of com- bining an in vitro approach for simultaneous quantification of fat accumulation with other hepatotoxicity-associated relevant events in human cells with extrapolation to relevant human exposures. In this study, we used high-content imaging to simultane- ously quantify lipid accumulation, as a marker of steatosis, along with several measures of hepatotoxicity that are asso- ciated with chemically induced liver steatosis. These include mitochondrial membrane potential disruption, oxidative stress, and nuclear morphology changes. Furthermore, these were evaluated on single cell and population levels. Sixteen reference chemicals known to induce steatosis by diverse mechanisms were tested in HepaRG cells. To evaluate the physiological relevance of the resulting dose–response rela- tionships, we integrated the in vitro data for the reference compounds with physiologically based pharmacokinetic modeling (PBPK) and reverse dosimetry to derive human dose equivalents. These levels were then compared to in vivo exposure values. Finally, on the basis of predicted nuclear receptor agonist activity, we curated and tested a panel of food-related compounds and pesticides for their potential to induce steatosis and associated hepatotoxicity. Materials and methods Chemicals and reagents Chemicals (Supplementary table 2) were purchased from Sigma-Aldrich (St. Louis, Missouri). Stock solutions (200x) were prepared in dimethylsulfoxide (DMSO). Chemicals soluble in water were directly dissolved in cell culture medium. A mixture of oleic acid/palmitic acid (OA/PA) was combined with bovine serum albumin (BSA) in a ratio of 5.5:1 (OA/PA: BSA). OA/PA was first dissolved in DMSO and then added to treatment medium supplemented with the corresponding BSA concentration. The OA/PA medium mix was then warmed at 60 °C for 1 h (to aid dissolution), allowed to return to room temperature, and subsequently used for experiments. Cell culture and chemical exposure Undifferentiated HepaRG cells were purchased from Bio- predic international (Saint Grégoire, France). HepaRG cells have unlimited growth capacity and phenotypic stability, stable enzymatic activities over 4 weeks of cultivation, and similar metabolic function to primary human hepatocytes, and also avoid limitations of primary cells such as donor variability and phenotypic instability (Aninat et al. 2006; Guillouzo et al. 2007; Lambert et al. 2009; Lübberstedt et al. 2011; Andersson et al. 2012). The line contains two cell types, hepatocyte-like and cholangiocyte-like cells (Parent et al. 2004), which allows for elucidation of cell-specific toxicity (Jossé et al. 2008). Cells were seeded at a density 1703Archives of Toxicology (2023) 97:1701–1721 1 3 of approximately 27′000 cells/cm2 into a 150 × 25 mm petri dish and were cultured for 2 weeks in William’s E medium (cat. nr.: 32,551,020, ThermoFisher Scientific, Massachu- setts, USA) supplemented with 2 mM glutamine, 10% (v/v) fetal bovine serum (FBS) good forte (cat. nr.: P40-47,500, PAN-Biotech, Aidenbach, Germany), 100 U/ml penicillin and 100 µg/ml streptomycin (cat. nr.: 15,140,122, Ther- moFisher Scientific, Massachusetts, USA), 5 µg/ml Gibco recombinant AOF insulin (cat. nr.: A11382II, ThermoFisher Scientific, Massachusetts, USA) and 50 µM hydrocortisone hemisuccinate (cat. nr.: H4881, Sigma-Aldrich, St. Louis, USA) at 37 °C and 5% CO2. Differentiation was induced by culturing the cells for another 2 weeks in the above mentioned medium supplemented with 1.7% DMSO. After 4 weeks in culture, differentiated HepaRG cells were seeded at a high density of approximately 227′000 cells/cm2 (72,000 cells/well with 100 µl per well) into black 96-well plates with a clear bottom (cat.nr.: 3603, Corning, New York, USA) and maintained for 10 days. After 10 days, cells were first adapted to treatment medium (culture medium contain- ing 0.5% DMSO and 2% FBS) for 48 h and then used for cell viability and high-content imaging analysis. Cell viability assessment Cell viability was analyzed using the WST-1 assay which is based on the measurement of mitochondrial dehydroge- nase activity (cat. nr: ab155902, Abcam, Cambridge, UK). HepaRG cells were exposed to chemicals at defined con- centrations for 24 h. Triton X-100 (0.01%) was a positive control. One hour before the end of incubation period, 10 µl of WST-1 reagent was added to each well, quickly mixed by shaking the plate, and returned to the cell culture incuba- tor (37 °C, 5% CO2). After 1 h, absorbance was measured (450 nm, reference wavelength 620 nm) using a plate reader (Infinite M200 Pro, Tecan group, Männerdorf, Switzer- land). Absorbance was then subtracted from the reference wavelength, corrected for background by subtracting values measured without cells, and then normalized to the solvent control, which was set to 100%. Dose–response data were fitted using drc package (Version 3.0.1) (Ritz et al. 2015) in R 3.6.1. Data were from at least 3 independent experiments with 6 technical replicates per condition. High‑content imaging assay Cellular responses were quantified using a high-content imaging approach (Fig. 1). Based on cell viability assess- ment, concentrations below the EC50 level were tested in the high-content imaging assay. After 24 h chemical exposure, live cells were incubated with CellROX, which reflects hydroxyl radical and superoxide anion levels (ThermoFisher Scientific 2011) (5 µM, cat. nr.: C10422, ThermoFisher Scientific, Massachusetts, USA) by adding 10 µl of a 55 µM CellROX solution in treatment medium directly to the well. After incubation (30 min, 37 °C, 5% CO2), cells were rinsed twice with treatment medium with- out serum. After removal of the treatment medium, cells were incubated with Mitotracker Orange CM-H2TMROS which was used to measure changes in mitochondrial mem- brane potential (50 nM, 10 min, 37 °C, 5% CO2 cat. nr.: M7511, ThermoFisher Scientific, Massachusetts, USA). Subsequently, cells were washed three times with PBS and fixed with 4% PFA at room temperature for 20 min. Both dyes were retained after aldehyde fixation. After washing 3 × 5 min, cells were stained with Bodipy (50 nM, cat. nr.: D3921, ThermoFisher Scientific, Massachusetts, USA) and Hoechst (2 µM, room temperature, 15 min, cat. nr.: 62,249, ThermoFisher Scientific, Massachusetts, USA), rinsed once with PBS (cat. nr.: 14,190,250, ThermoFisher Scientific, Massachusetts, USA), and then kept in PBS. Images were acquired using an ImageXpress micro (IXM) High-Content Imaging System from Molecular Devices with a 20 × 0.75 NA S Fluor objective. A photometrics CoolSNAP HQ high resolution camera designed for quantitative fluorescence microscopy applications was set to 16 bits with a binning of one and a pixel size of 6.45 µm × 6.45 µm, and 25 images were captured per well. Automated quantitative image analysis Images were quantified using the open-source software CellProfiler 3.1.8 (McQuin et al. 2018) and ImageJ 1.52 g (Schindelin et  al. 2012). The image analysis pipelines are available as supplementary materials and are briefly described here. For cell population-level data, nuclei were segmented using Otsu thresholding (Sankur 2004) and eve- rything outside the diameter range of 10–100 pixel was dis- carded. For lipid droplet quantification, GFP images were first deconvoluted using ImageJ 1.52 g and then imported into the CellProfiler pipeline. A binary mask of the lipid droplets then was created, segmented, and quantified using the standard settings in the “IdentifyPrimaryObjects” mod- ule (Supplementary Fig. 1). Next, mitotracker and Cell- ROX image intensity were measured. All parameters were exported to Excel spreadsheets and analyzed in R 3.6.1, as further explained in the paragraph “Dose response data and model fitting” below. For the single cell analysis pipeline, nuclei images were used to identify individual cells. Cell borders were determined based on the mitotracker staining with the propagation method (Jones et al. 2005). All quanti- fications from the cell population-level pipelines were also performed at the single cell level. In addition, the follow- ing parameters were quantified per cell: size and shape of identified objects, object intensity, and object textures. The data were then exported to an SQLite database and analyzed 1704 Archives of Toxicology (2023) 97:1701–1721 1 3 using R 3.6.1 or GraphPad Prism 9. These pipelines can be used for a rapid population-based quantification which is amenable for high-throughput assay or for a detailed single cell analysis to identify heterogenous responses and distinct cellular phenotypes. Data analysis of single cell data The single cell SQLite database file was read using RSQLite package (Müller et al. 2020). Data across three independ- ent experiments with 6 technical replicates per concentra- tion were summarized in a data matrix. This data matrix consisted of over 106 cells with over 800 features per cell. It was imported into R 3.6.1. Constant and NA values were removed from the data matrix to perform a principle com- ponent analysis. Subsequently, the parameter cytoplasm area size was used to separate hepatocyte-like cells (< 5000  px2) and cholangiocyte-like cells (> 5000  px2) with the assump- tion that cells with a cytoplasm area below 5000 px2 are likely to be hepatocytes and those above likely to be chol- angiocytes. Separation by size was further validated by anti-ASGPR1 immunostaining, a basolaterally expressed surface protein highly specific to hepatocytes (Supplemen- tary Fig. 2). Its main function is the clearance of desialylated Fig. 1 Basis of high-content imaging assay. A Fluorescent probes for the selected endpoints including corresponding positive controls. B Description of the workflow including rapid population-based quanti- fication and detailed single cell analysis. During single cell level anal- ysis, over 800 different features (intensity, size/shape, and texture) per cell were quantified and analyzed using principal component analysis (PCA). Plotted data were then color labeled using the values from fol- lowing parameters: lipid droplet size, mitotracker intensity and Cell- ROX intensity 1705Archives of Toxicology (2023) 97:1701–1721 1 3 glycoproteins (D’Souza and Devarajan 2015). Principle component analysis was used for dimensionality reduction and performed on the hepatocyte data set and plotted using ggplot2. Machine learning approach to quantify apoptotic cells Apoptotic cells were quantified using a supervised machine learning approach. Nuclei images were analyzed with Cell- Profiler and intensity, size, shape and texture features were extracted. These data were loaded into CellProfiler Analyst (Jones et al. 2008) together with the raw nuclei images. A random forest classifier was trained to distinguish between healthy and apoptotic cells using a training set consisting of 477 different nuclei (Supplementary Fig. 3A). To take imbal- ance in the dataset into account, random undersampling was performed in the training set to achieve a 67/33 ratio (323 healthy, 154 apoptotic). The training set was assembled using nuclei from the menadione exposure across three inde- pendent experiments. Apoptotic nuclei were defined based on morphological nuclear fragmentation. The model was trained until a classification accuracy between 80 and 90% was reached. Classification accuracy was determined based on user-defined assessment. In addition, images from the positive control were visually inspected to check the overall classification performance of the model and were corrected if numerous misclassifications were identified. If the clas- sification was satisfying, the entire dataset was then scored. The resulting hit table was imported into R and the number of apoptotic cells were normalized to total cell number in the same well and further to untreated cells (Supplemen- tary Fig. 3). Image analysis and supervised machine learning were performed on a terminal server based on Windows 10 equipped with 40 cores Intel Xeon Gold 6150 @ 2.7 Ghz and 512 GB RAM. Dose–response data and model fitting Dose–response data were collected on a population level (i.e., per well averages) and combined with responses from all HepaRG cells per well. For every biological endpoint, the following raw parameters were selected per well and normalized to cell number per well: • Lipid accumulation: Total number of lipid droplets divided by total number of cells • Mitochondrial membrane potential: (Mean intensity– background intensity)/total number of cells • Oxidative stress: (Mean intensity–background intensity)/ total number of cells • Nuclear morphology/cell death: number of apoptotic cells/total number of cells These data were normalized to untreated cells to derive a fold change; unexposed cells had a fold change of 1. All data were normalized per plate to overcome possible plate effects. Non-linear regression of dose–response data was performed by fitting a log-logistic model with 3 or 4 parameters. The best fitting model was selected based on lowest Akaike Information Criterion (AIC) (Akaike 1973) and visual inspection. A positive hit, i.e., an induction of a biological endpoint by exposure, was identified when the non-linear regression exceeded the calculated threshold of ± 2 standard deviations of the response at the two lowest concentrations. Dose–response data were fitted using the drc package (Ritz et al. 2015) and plotted with ggplot2 (Hadley Wickham 2016). The R script of dose–response data and model fitting is available as supplementary material. Benchmark concentration modeling To determine a point of departure of the high-content imag- ing concentration–response data, benchmark concentration (BMC) was modeled and calculated using the benchmark dose software (BMDS) version 3.1.2 (United States Environ- mental Protection Agency 2023). The benchmark dose tech- nical guidance was used as a basis to perform the analysis (United States Environmental Protection Agency (USEPA) 2012). Normalized data including tested concentrations were entered as individual data points into the BMDS Excel appli- cation. The BMC was then defined as a benchmark response (BMR) of a change in the mean of one standard deviation from the control mean. Concentration–response data were then fitted to several models including exponential, hill, lin- ear, polynomial and power models. Normal and log-normal distribution with constant or non-constant variance were also tested to find the best fit. BMC and its corresponding 95% confidence interval (95% CI) were then calculated. Selec- tion criteria for the model with the best fit were the Akaike Information Criterion (AIC), goodness of fit p-value, scaled residuals for dose group near BMD, and for control dose group and BMDS recommendation. If there were several models which received a BMDS recommendation, BMCs from all models were averaged. BMCs were calculated for all chemicals tested (Supplementary table 3). In vitro to in vivo extrapolation (IVIVE) IVIVE was performed based on the previously published method by Rotroff and Wetmore et al. (Rotroff et al. 2010; Wetmore et al. 2015) and using the high-throughput toxi- cokinetic (httk) R package (version 1.10.1) developed by the US-EPA. The httk package contains four toxi- cokinetic models which can be parameterized using high throughput-derived in vitro data on plasma protein binding and hepatic clearance. Moreover, it has a Monte Carlo 1706 Archives of Toxicology (2023) 97:1701–1721 1 3 sampler to simulate population variability and includes tools for reverse dosimetry together with functions for the analysis of concentration versus time simulations (Pearce et al. 2017). In our study, a three-compartment steady state pharmacokinetic (PK) model was parameterized and used for the simulations. For all simulations, the follow- ing assumptions were made: oral route of exposure, daily dosing with constant dose rate, and 100% bioavailability. Human data were used except where otherwise stated. Steady-state concentration (Css) in the blood was calcu- lated using formula 1 (Pearce et al. 2017) where ko = constant dose rate (mg/kg BW/day), fup = in vitro measured chemical fraction unbound in plasma, Clmetabolism = metabolic clearance scaled from in  vitro intrinsic hepatic clearance, Rblood2plasma = ratio of blood concentration of a chemical to the plasma concentration, GFR = glomerular filtration rate (mean 5.17 ml/min/kg^3/4), Qliver = blood flow to the liver (mean 59.9 ml/min/kg^3/4), and Qgut = blood flow to the gut (mean 47.5 ml/min/kg^3/4). Fraction unbound, metabolic clearance, and Rblood2plasma values were obtained from the httk package or from litera- ture. Recognizing that formula 1 is linear in dose rate (ko), oral equivalent doses (OED) were then calculated using the Css predicted for a 1 mg/kg/day dose rate as in formula 2 (Rotroff et al. 2010). The calculated OED corresponds to an oral daily dose which would lead to a plasma concentration equal to the in vitro-derived BMC. The OED is linearly related to the BMC and inversely related to Css. This equation is only valid for first-order metabolism (Rotroff et al. 2010). A Monte Carlo analysis was used to simulate population variability in, GFR, Qliver, Qgut, and Clmetabolism across 1000 healthy individuals. Extrapolations were performed for each chemical and endpoint which was positive in the imaging assay. Every predicted parameter per chemical was then combined into one boxplot per endpoint, which contained all predicted OEDs across all positive endpoints. For the model evaluation, Css of selected reference chemicals were predicted using formula 1 with the same dose rate from the in vivo Css measurements. For chemi- cals with no in vitro intrinsic hepatic clearance and frac- tion unbound values available in the database included in the httk package, in silico predictions from (Sipes et al. (1) Css (mg l ) = ko∕(GFR × fup + ( Qliver + Qgut ) × fup × Clmetabolism ( Qliver + Qgut ) + fup×Clmetabolism Rblood2plasma ) (2)OED ( mg kg day ) = BMC(uM) × 1 mg kg day Css(uM) 2017) were loaded using load_sipes2017(). In addition to that, in vitro intrinsic hepatic clearance values were cal- culated for a selection of chemicals (Lomitapide, Fialu- ridine, Metformin, beta-naphthoflavone and Menadione) using in vivo clearance values (Supplementary table 4) and formula 3: where Clint = in vitro intrinsic hepatic clearance in ul/ min/10^6 cells, Clmetabolism = in vivo clearance in L/h/kg (if available hepatic otherwise total clearance), Fub.corr = assay correction factor which is assumed to be 1 (i.e., no correc- tions), ρliver = liver density of 1.05 g/ml, Ncells = 1.1 × 10^8 hepatocytes per gram of liver, Vliver = liver volume: 0.0245 L/kg (human), 0.0349  L/kg (rat), 0.04 L/kg (rab- bit), fmc2c = conversion factor for millions of cell to one cell, fml2l = conversion factor from ml to l, fmin2h = conversion fac- tor from minute to hour, and fl2ul = conversion factor from L to ul. The formula to calculate Clint was derived from the function calc_hepatic_clearance () within the httk pack- age. Predicted Css values were then compared to published in vivo Css values. If there was no published in vivo Css available, Css was then calculated based on published in vivo pharmacokinetic parameters using the following formula (Supplementary table 3) Vice versa, predicted OED derived from BMC were com- pared to published exposure values of the corresponding chemicals. Simulations were performed in humans where pharmacokinetic parameters for humans were available; otherwise, estimations were calculated in rats or rabbits, as stated in Fig. 7. The R script written to perform the IVIVE analysis is available as supplementary material. Results Multiparametric analysis of reference chemicals To validate the assay, positive and negative control chem- icals with previously established mechanisms of action were characterized by the high-content imaging approach described above (Fig. 1A). Thus, we tested a mixture of oleic acid and palmitic acid, which are often used to induce steatosis in vitro (Gómez-Lechón et al. 2007; Sharma et al. 2011; Graffmann et al. 2016; Michaut et al. 2016); rote- none, a selective complex I inhibitor (Siddiqui et al. 2013) and carbonylcyanide-p-trifluoromethoxyphenylhydrazone (FCCP), a weak acid that selectively increases proton (3) Clint = (Clmetabolism × Fub.corr) �liver × Ncells × Vliver × fmc2c × fml2l × fmin2h × fl2ul) (4)Css( mg L ) = (Dose × F)∕(Cl × dosinginterval) 1707Archives of Toxicology (2023) 97:1701–1721 1 3 permeability in lipid membranes (Benz and McLaughlin 1983; Sakamuru et al. 2016), were therefore used as posi- tive controls for inducing mitochondrial dysfunction via a decrease in mitochondrial membrane potential; amiodar- one and arsenite, which induce oxidative stress (Anthérieu et al. 2011; Tolosa et al. 2016; Zhao et al. 2019; Lv et al. 2020); menadione, which induces DNA damage and cell death (Loor et al. 2010); and caffeine as negative control (Persson et al. 2013; Saito et al. 2016) (Fig. 2). Concen- tration ranges were selected on the basis of preliminary cell viability screen, and each chemical was used at con- centrations ranging up to the corresponding EC50 value (Supplementary Fig. 4). Response profiles of the positive and negative control compounds were consistent with expected outcomes. Lipid accumulation increased in a dose-dependent manner upon exposure to oleic acid/palmitic acid, rotenone, amiodarone and menadione, with oleic acid/palmitic acid (1000 µM) being the most efficacious (2.3-fold increase). Caffeine induced a decrease of lipid accumulation that exceeded the threshold, but only at the highest concentration (5000 µM). In the case of mitochondrial membrane potential, oleic acid/ palmitic acid (62.5 µM), rotenone (2.5 µM), fccp (2.5 µM) and amiodarone (25 µM) all induced a dose-dependent decrease exceeding the threshold at the concentrations indi- cated. In the case of oxidative stress, amiodarone (25 µM) Fig. 2 Representative fluorescent images from chemically exposed HepaRG cells. A DMSO control B Oleic acid/palmitic acid (1000 µM) C Rotenone (20 µM) D Amiodarone (50 µM) E Menadi- one (50 µM) F) Caffeine (5 mM). Uncropped images were acquired at 20x. Scalebar = 50 µm. Blue nuclei, green lipid accumulation, red mitochondrial membrane potential, yellow oxidative stress 1708 Archives of Toxicology (2023) 97:1701–1721 1 3 and arsenite (25 µM) exposure led to a dose-dependent increase and exceeded the threshold at the concentrations indicated. The number of apoptotic cells (nuclear mor- phology) increased in oleic acid/palmitic acid-(250 µM), amiodarone- (25 µM), arsenite- (25 µM) and menadione- (25 µM) treated HepaRG cells in a dose-dependent fash- ion, exceeding the threshold at the concentrations indicated (Fig. 3). These observations are consistent with previous response patterns and sensitivities (Anthérieu et al. 2011; Tolosa et al. 2016; Angrish et al. 2017; Allard et al. 2020) (Gómez-Lechón et al. 2007; Anthérieu et al. 2011; Sharma et al. 2011; Siddiqui et al. 2013; Graffmann et al. 2016; Michaut et al. 2016; Tolosa et al. 2016), suggesting the accuracy of the in vitro assay and data analysis pipelines. Following the effective systematic characterization of positive and negative control chemicals for each endpoint, we tested four Food and Drug Administration (FDA)- approved drugs with known mechanisms of action. They included metformin, an anti-diabetic drug that inhibits mito- chondrial complex I (Pernicova and Korbonits 2014), tetra- cycline, an antibiotic drug that inhibits β-oxidation, micro- somal triglyceride transfer protein (MTP) and upregulates PPARγ and SREBP1-c (Anthérieu et al. 2011; Schumacher and Guo 2015), lomitapide, an anti-hypercholesterolemic drug that inhibits MTP (Lin et al. 2014), and valproic acid, an anti-convulsant drug that inhibits β-oxidation and upregu- lates CD36 and DGAT2 (Schumacher and Guo 2015; Bai et al. 2017). Lipids accumulated in a dose-dependent man- ner above the threshold for metformin, tetracycline and valproic acid, but not lomitapide. Mitochondrial membrane potential decreased when cells were exposed to metformin (50,000 µM) or lomitapide (12 µM), but only at the highest concentrations. An unexpected hormetic effect was observed in cells exposed to tetracycline (between 31.25 and 125 µM); a hyperpolarization, which to our knowledge has not been reported previously, was observed, such that only at the highest concentration (250 µM) did the mitochondrial mem- brane potential decrease. Oxidative stress dose-dependently increased after tetracycline (125 and 250 µM) exposure and exceeded the threshold at the concentrations indicated, but no changes were observed for lomitapide or valbroic acid. A unique decrease of oxidative stress below the threshold was observed after metformin exposure (50,000 µM). Finally, the number of apoptotic cells increased above the threshold in metformin, tetracycline and lomitapide treated cells (Fig. 3). These data thus confirmed the anticipated activation profiles of the cellular markers for these reference compounds. As a next step, we were interested to test etomoxir and fialuridine, both chemicals that had their clinical develop- ment terminated due to severe hepatotoxicity in clinical tri- als, but with limited knowledge of the underlying mecha- nisms (Manning and Swartz 1995; Holubarsch et al. 2007). Etomoxir was developed as an anti-diabetic drug and is an irreversible carnitine palmitoyltransferase 1 (CPT1) inhibitor (Merrill et al. 2002). Fialuridine is a nucleoside analog that was developed to treat chronic hepatitis and inhibits mito- chondrial DNA polymerase γ (McKenzie et al. 1995; Lewis et al. 2003). However, neither induced lipid accumulation, had any consistent impact on oxidative stress, nor induced apoptosis (Fig. 3). Nonetheless, fialuridine (25 µM) and eto- moxir (36 µM) induced an increase in MMP and exceeded the threshold the concentrations indicated (Fig. 3). The profiles for these known steatotic compound, i.e., lack of lipid accumulation with MMP increase, illustrate a concept behind using a hallmark-driven assay and potential diffi- culties in modeling lipid accumulation after short in vitro exposures. Thus, chemicals that induce related measures of hepatotoxicity, which may be easier to detect, may be flagged for more in-depth evaluation with a focus on lipid accumulation potential, such as longer duration exposures. We next were interested in testing chemicals known to bind to nuclear receptors, as nuclear receptor activa- tion is a molecular initiating event of steatosis. Thus, we characterized three experimental drugs: T0901317, a syn- thetic liver X receptor agonist that increases the expression of SREBP-1c and CHREBP (Mitro et al. 2007; Cha and Repa 2007); WY14643, a dual PPARα/PPARγ agonist that increases oxidative stress in mouse liver and affects fatty acid metabolism (Woods et al. 2007; Jennings et al. 2014); and β-naphthoflavone (BNF), an aryl hydrocarbon receptor agonist that upregulates CD36 (Lee et al. 2010; Jennings et al. 2014). Lipid accumulation dose-dependently increased and exceeded the threshold for cells exposed to T0901317 and WY14643, but not BNF. In addition to inducing lipid accumulation, T0901317 caused a decrease in MMP below the threshold, and oxidative stress and apoptosis exceed- ing the threshold at the highest T0901317 concentration (40 µM). For WY14643, which also induced lipid accumu- lation, we observed MMP hyperpolarization and it induced apoptosis Yet, none of the hallmarks were induced by BNF (Fig. 3). Thus, 2 out of 3 nuclear receptor agonists induced several steatosis and hepatotoxicity related cellular events. Analysis of food‑related chemicals and pesticides With nuclear receptor activation as a potential initial pre- dictor of steatotic potential, as further supported by the T0901317 and WY14643 results, we wanted to character- ize more broadly environmentally relevant compounds with structural alerts for nuclear receptor binding. Thus, about 6,000 Smiles (Simplified Molecular Input Line Entry Sys- tem, a chemical language system in which chemical struc- tures can be written using ASCII characters (Weininger 1988)) were extracted from the regulated food-related use/ occurrence chemical database from the Swiss Federal Food Safety and Veterinary Office and about 4,000 from the 1709Archives of Toxicology (2023) 97:1701–1721 1 3 Fig. 3 Cellular responses for four key steatosis-relevant endpoints. Gray area around non-linear regression represents 95% CI. Black dot- ted line represents a fold change of 1 (DMSO control). Red dotted lines represent threshold (± 2 standard deviations of the response at the two lowest concentrations). All data were first normalized to cell number and then to untreated cells. At least 3 independent experi- ments with 6 technical replicates per concentration per experiment were conducted for all chemicals 1710 Archives of Toxicology (2023) 97:1701–1721 1 3 Toxcast database (Richard et al. 2016). We assessed poten- tial for nuclear receptor binding on the basis of molecular fragments and other relevant chemical features, by a previ- ously reported approach (Mellor et al. 2016), and identified about 80 chemicals that might bind one or more receptors. Of these 80 compounds, 14 had prior evidence of hepato- toxicity and published toxicokinetic data (Supplementary Table 6), and were selected for further testing in the assay established in this study. Of these 14 food-related chemicals and pesticides, six of them (uric acid, tartrazine, bisphenol A, atrazine, metazachlor, and vinclozolin) induced signifi- cant lipid accumulation. Orotic acid, fructose and carbo- furan did not stimulate lipid accumulation, but did stimulate a dose-dependent increase in MMP exceeding the threshold. No response was observed in oxidative stress or the number of apoptotic cells after 24 h exposure (Fig. 4). Thus, 9 out of 14 food-related chemicals and pesticides induced steatosis and/or an increase in MMP, whereas none induced oxidative stress or an increase in the number of apoptotic cells. Chemical exposures induce distinct cellular endophenotypes within cell populations To evaluate the potential heterogeneity of single cell responses and interrogate potential endophenotypes within and/or between cell populations in the HepaRG cells, we further evaluated alterations of lipid accumulation, mito- chondrial membrane potential and oxidative stress following tetracycline exposure. Tetracycline was selected because it induced responses above the threshold in all four endpoints. First of all, we observed that in general, there were more lipid droplets in hepatocyte-like cells compared to cholangi- ocyte-like cells even in unexposed cells (Fig. 5), and, moreo- ver, at 125 µM and 250 µM, the median number of lipid droplets per hepatocyte-like cell did not drastically increase (~ 8 to ~ 12 lipid droplets per cell), but the proportion of hepatocyte-like cells with more lipid droplets increased. The opposite response was observed for cholangiocyte-like cells: the median number of lipid droplets per cholangiocyte- like cells increased, but the proportion of cholangiocyte-like cells with more lipid droplets remained unchanged (125 µM or 250 µM tetracycline, Fig. 5). Similar to lipid droplets, the mitochondrial membrane potential also increased more in hepatocyte- compared to cholangiocyte-like cells for all tet- racycline concentrations tested. Notably, at the highest con- centration of tetracycline (250 µM), a decrease of mitochon- drial membrane potential was observed in both cell types. Finally, for oxidative stress, similar patterns of changes were observed in hepatocyte-like cells which had increased lev- els of oxidative stress compared to cholangiocyte-like cells throughout all tested concentrations. At the highest concen- tration at which tetracycline was tested (250 µM), oxidative stress levels were elevated in both cell types. By evaluating increased oxidative stress on a single cell basis, it became apparent that the increased average values for oxidative stress on a cell population level arose from a small num- ber of cells with extremely high levels of oxidative stress, rather than constitutive moderate increase across the whole cell population (Fig. 5). In summary, hepatocyte-like cells appear to be more sensitive than cholangiocyte-like cells and Fig. 4 Cellular responses for four key steatosis-relevant endpoints of food-related chemicals. Gray area around non-linear regression repre- sents 95% CI. Black dotted line represents a fold change of 1 (DMSO control). Red dotted lines represent threshold (± 2 standard deviations of the response at the two lowest concentrations). All data were first normalized to cell number and then to untreated cells. 2-3 independ- ent experiments with 6 technical replicates per concentration per experiment were conducted for all chemicals 1711Archives of Toxicology (2023) 97:1701–1721 1 3 that population-level responses appear to be driven by the evolution of small populations with high responses rather than an equal distribution of responses amongst cells. Given the heterogeneity of responses in hepatocyte- and cholangiocyte-like cells, we investigated whether hepatocyte-like cells with a strong response in one end- point (e.g., large lipid droplet size) responded similarly in other endpoints (mitochondrial membrane potential and oxidative stress) at the same concentration. For this Fig. 5 Quantification of selected endpoints in individual hepatocytes versus cholangiocyte-like cells after tetracycline exposure. The same raw data used in the population-level analysis were also used for this single cell analysis. 30,000 to 150,000 cells per concentration and cell type were analyzed. Cells from 3 independent experiment were included in the analysis. Notched boxplots were overlaid in the violin plots to display the confidence interval around the median 1712 Archives of Toxicology (2023) 97:1701–1721 1 3 purpose, over 800 features per cell were extracted by high-content imaging and resulting PCA plots were color labeled with lipid droplet size, mitochondrial membrane potential, and oxidative stress level to visualize their distribution (Fig. 6). Heterogeneity in lipid droplet size and mitochondrial membrane potential was observed in both chemically and mock-exposed hepatocyte-like cells, whereas it was only observed for oxidative stress in chemi- cally exposed cells. Oxidative stress levels in non-exposed hepatocytes were similar across all measured hepatocytes (Fig. 6A). While mitochondrial membrane potential did not decrease in hepatocytes with larger lipid droplets it did decrease in hepatocytes with smaller lipid droplets (Fig. 6B). Similarly, oxidative stress levels in chemically exposed hepatocytes with larger lipid droplets were lower than in hepatocytes with medium and small sized lipid droplets. Interestingly, a subgroup of hepatocytes with medium and small sized lipid droplets also had lower oxi- dative stress levels compared to other hepatocytes with the same lipid droplet size. Thus, tetracycline exposure led to heterogenous response patterns in hepatocytes involving a subpopulation that was protected from mitochondrial dysfunction and oxidative stress. IVIVE We used pharmacokinetic modeling and reverse dosimetry to evaluate the relevance of the in vitro concentrations that stimulated one or more of the endpoints defined as exceed- ing a calculated threshold of ± two times the standard devia- tions of the response to the two lowest concentrations. Pre- dicted Css values were compared with in vivo Css values (Fig. 7A), which were derived from pharmacokinetic studies in human, rat or rabbit (Supplementary table 3). For the ten chemicals, the coefficient of variation (R2) on a logarithmic scale was 0.78 and the root mean squared error was 0.69 (a factor of 4.9x). The chemicals with the greatest discrepancy between actual and predicted Css were beta-naphthoflavone and WY-14643, where the predicted Css were about 13 to 26 times higher than the actual Css values. Benchmark concentrations (BMC) were then calculated for every positive response indication in the imaging assay to determine a point of departure (Supplementary table 3 and Table 1). A direct comparison between the in vitro-derived BMC and in vivo Css was then made to identify a potential in vivo hazard. For menadione, the BMC was about 1.5–4 times lower than the median predicted Css. For amiodarone and valproic acid, it was about 1.5–4 times higher and for Fig. 6 Heterogenous response pattern after tetracycline exposure in hepatocytes. A untreated B tetracycline (250 µM). A single point in the PCA plot represents a hepatocyte and its color code represents lipid droplet size (px^2), mitochondrial membrane potential (ΔΨm) (AU) or oxidative stress (AU). 80,000 to 125,000 cells from 3 inde- pendent experiments were analyzed per concentration 1713Archives of Toxicology (2023) 97:1701–1721 1 3 artrazine about double the predicted Css.. For orotic acid, uric acid, carbofuran, and carbosulfan, the calculated BMC was about 10–550 times higher than the predicted Css. For lomitapide, fialuridine, WY14643, metformin, and tetra- cycline, the BMC was about 10–10,000 times higher than the median PBPK-predicted Css. This direct comparison of in vitro-derived BMC with in vivo Css was consistent with the indication of menadione, but only this compound, which is a known hepatotoxicant, was identified as a hazard. As the Css metrics were not available for all the chemi- cals, we used a three-compartment steady state PK model to derive oral equivalent doses (OED) from respective BMC values. OEDs were predicted for chemicals with available exposure data such as human therapeutic doses, doses from clinical trials, food and drinking amounts, or experimental doses in animal models. OED was then predicted for every endpoint for which a significant response was observed in the imaging assay (Fig. 7 and Tables 1 and 2). Predicted OED values overlapped with realistic doses for the drugs fialuridine, tetracycline, amiodarone, valproic acid and menadione, whereas for Rotenone, WY14643, Lomitapide and Metformin, the OED values were 10–10,000 times higher than the doses. The OED values for orotic acid and carbosulfan overlapped with their reported no observable adverse effect level (NOAEL), lowest observed adverse effect level (LOAEL) concentrations, or the exposure estimate. Discussion In this study we combined a high-content imaging assay for chemically exposed metabolically competent human liver cells with PBPK modeling to derive relevant human exposure predictions for the potential of chemicals to induce hepatic steatosis or other hepatotoxic processes. Validation studies were carried out using 16 reference chemicals with an overall sensitivity of 66% and specificity of 69% (Fig. 8). Fig. 7 Pharmacokinetic modeling of selected refer- ence chemicals. A Forward dosimetry of chemicals with published Css values. Yellow boxplot represents predicted Css using IVIVE. B Reverse dosimetry of chemicals with published human exposure val- ues from different sources. Blue boxplot represents predicted oral equivalent doses using httk with a monte carlo population sampler. OED predictions were performed for every endpoint from the HCI assay which was above the calculated threshold 1714 Archives of Toxicology (2023) 97:1701–1721 1 3 Notably, 6/9 predicted OED values for selected reference chemicals overlapped with estimated exposure values. Orotic acid disrupted MMP, but not other markers, at a concentra- tion corresponding to an OED that overlaps with estimated human exposures. These outcomes are notable for the high accuracy of predicted OED with human exposure values for known chemicals, and the application of the method to screen food-related chemicals. In this study we correctly flagged 7/11 chemicals known to induce lipid accumulation, and also identified two chemi- cals, metformin and WY-14693, as having the potential to induce lipid accumulation, that were unexpected (Fig. 8). A potential rationale for this observation, in the case of metformin, is that it induced lipid accumulation only at extremely high, physiologically irrelevant concentrations (25,000 and 50,000 µM), whereas data suggesting that met- formin reduces lipid accumulation was observed for con- centrations ranging from 100 to 2 mM (Zare et al. 2019; Kim et al. 2020). Furthermore, four chemicals, lomitapide, etomoxir, fialuridine and β-naphthoflavone did not give rise to lipid accumulation as we would have expected based on their established mechanism of action and previous clini- cal observations (Supplementary table 2) (McKenzie et al. 1995; Cuchel et al. 2013). A possible explanation for these Ta bl e 1 IV IV E of se le ct ed fo od -r el at ed c he m ic al s w hi ch sh ow ed a re sp on se a bo ve th e ca lc ul at ed th re sh ol d in th e im ag in g as sa y LA li pi d ac cu m ul at io n, M M P m ito ch on dr ia l m em br an e po te nt ia l, O S ox id at iv e str es s, N M n uc le ar m or ph ol og y, 1 B en ch m ar k co nc en tra tio n (µ M ) p er im ag in g en dp oi nt c al cu la te d w ith B M D S so ftw ar e, 2 D er iv ed fr om in  v itr o to in  v iv o ex tra po la tio n, 3 In v itr o to in  v iv o ex tra po la tio n de riv ed m ed ia n C he m ic al LA 1 M M P1 O S1 N M 1 M ed ia n C ss 2 (µ M ) In v iv o bl oo d pl as m a co nc . (µ M ) O ED 3 (m g/ kg B W /d ay ) N O A EL /L O A EL fr om li te ra - tu re (m g/ kg B W /d ay ) Ex po su re es tim at e fro m lit er at ur e M ax im um re si du e le ve ls (m g/ kg ) Re fe re nc e O ro tic a ci d – 42 3 – – 7. 04 0. 15 –0 .3 61 .2 (9 5% C I 2 1. 1– 13 4) 50 /1 00 2– 10 0  m g/ kg B W /d ay – (A gu ila r e t a l. 20 09 ; D ’A po lit o et  a l. 20 12 ) U ric a ci d 59 1 – – – 6. 84 21 7– 40 2 87 .3 (9 5% C I 6 7. 3– 15 0. 6) – – – (L i e t a l. 20 09 ) A tra zi ne 19 – – – 40 – 0. 5 (9 5% C I 0 .2 –1 .7 ) 5/ 50 (i nc re as e re la tiv e liv er w ei gh t) – 0. 05 –0 .1 (A ge nc y fo r t ox ic su bs ta nc es an d di se as e re gi str y A TS D R 20 03 ) C ar bo fu ra n – 49 75 – – 9 – 56 5 (9 5% C I 9 1. 7– in fin ity 0. 03 –0 .1 /0 .1 –0 .3 – 0. 01 –0 .5 (F oo d an d A ut ho rit y 20 09 a) C ar bo su lfa n 77 – – – 8. 5 – 10 .4 (9 5% C I 7 –1 9. 3) 0. 5– 1. 2/ 5– 62 – 0. 00 5 (F oo d an d A ut ho rit y 20 09 b) Table 2 Comparison of predicted oral equivalent doses (OED) with published exposure values a In vitro to in vivo extrapolation derived median b Exposure values include human therapeutic ranges or doses given to rat and rabbits c (United States Environmental Protection Agency (USEPA) 2007) d (Kleiner et al. 1997) e (Siddoway 2003) f (Bienenfeld et al. 2017) g (Pollinger and Merk 2017) h (Koch-Weser and Browne 1980) i (Hu et al. 1996) j (Panno et al. 2014) k (Garber et al. 1997) Chemical OEDa (mg/kg/d) In vivo exposure from literatureb (mg/ kg/d) Rotenone 0.3 0.01c Fialuridine 5.6 0.07–0.25d Amiodarone 15.8 2.9–22.9e Tetracycline 30.2 5.4–14.3f WY14643 1.7 0.4–0.7 g Valproic acid 92.7 3.6–30 h Menadione 55.4 1.1–36.6i Lomitapide 1333 0.07–0.86j Metformin 6778 7.1–35.7 k 1715Archives of Toxicology (2023) 97:1701–1721 1 3 apparent false negative results is that the 24-h duration of exposure may be too short for these chemicals to induce lipid accumulation. Indeed, lipid accumulation in patients was observed only after 13 weeks for fialuridine and 26 weeks for lomitapide (McKenzie et al. 1995; Cuchel et al. 2013). This limitation suggests that future research using the assay reported here, but beyond the scope of the present work, should aim to evaluate how longer duration in vitro expo- sures, i.e. on the order of one to two weeks may resolve false negative results. Encouragingly, in previous studies, hepato- toxicity was effectively characterized in HepaRG cells after up to 14 days of chemical exposure (Anthérieu et al. 2011; Dietrich et al. 2020; Donato et al. 2022). In addition to lipid accumulation data, we could assess changes in mitochondrial membrane potential, oxidative stress and nuclear morphology with varying degrees of accuracy (Fig. 8). Thus, in the case of the 10 chemicals expected to disrupt MMP, 9 were indeed positive in our assay, whereas one chemical, valproic acid, did not affect MMP as we would have expected based on its mechanism of action (Supplementary table 2) (Fig. 8). This is poten- tially expected since valproic acid has been characterized to induce mitochondrial dysfunction at very high concentra- tions, i.e. above 15 mM, and after prolonged exposure up to 72 h (Caiment et al. 2020) whereas we tested concentrations ranging from 94 to 6,000 µM and exposed the cells for 24 h. Interestingly, etomoxir and fialuridine also increased MMP, which to our knowledge has not been reported previously. However, the biological relevance of this observation is still unclear and needs further investigation. In the case of oxidative stress, sensitivity was poor; of the 9 chemicals expected to induce oxidative stress, only amiodarone, arsenite and tetracycline actually did, while six were unexpected negative outcomes, including the well- known ROS-inducer menadione (Fig. 8). This apparent low sensitivity could also be due to potential duration differ- ences and a capacity for adaptive responses since menadione induced oxidative stress in embryo chick cardiomyocytes Fig. 8 Quantification of the predictivity of the in vitro assay for A lipid accumulation B mitochondrial membrane potential C oxidative stress and D nuclear morphology for 16 reference chemicals using confusion matrices. Sensitiv- ity was calculated as “number of true positives/(number of true positives + number of false negatives)” and specificity was calculated as “number of true negatives/(number of true negatives + number of false positives)”. Overall sensitivity and specificity were calculated across all 4 endpoints. Of the 16 reference compounds tested, the outcomes were catego- rized as negative or positive, excluding amongst positive outcomes instances where no dose–response relationship was observed 7 2 4 3 Yes No Yes No Ac tu al Predicted Sensitivity: 64% Specificity: 60% 2 3 4 5 6 7 9 1 1 5 Yes No Yes No Ac tu al Predicted Sensitivity: 90% Specificity: 83% 2 4 6 8 3 0 6 7 Yes No Yes No Ac tu al Predicted Sensitivity: 33% Specificity: 100% 0 2 4 6 2 7 0 7 Yes No Yes No Ac tu al Predicted Sensitivity: 100% Specificity: 50% 0 2 4 6 A C B D Overall sensitivity: 66% Overall specificity: 69% 1716 Archives of Toxicology (2023) 97:1701–1721 1 3 after exposure durations as short as 25 min (Loor et al. 2010). Thus, in addition to longer than 24 h durations men- tioned previously, also shorter exposure durations than 24 h could be important to provide a more dynamic picture of oxidative stress responses in hepatocytes. While the simul- taneous detection of various markers is highly appealing, there is also a limitation that the dynamics of various pro- cesses arise following different durations. Finally, nuclear morphology appeared to be an over-sensitive read-out in this assay. Whereas only two chemicals were expected to induce such changes (arsenite and menadione), altered nuclear mor- phology was actually observed for a total of 8 chemicals. Nuclear morphology response mostly was above the calcu- lated threshold at the highest tested concentrations, consist- ent with cell viability data (Supplementary Fig. 4). Following the characterization of the various reference chemicals, we applied the assay to test 14 food-related chem- icals, 6 of which induced steatosis. Since the test chemi- cals were selected from a large number of chemicals with a focus on predictions of receptor-binding potential, it is not surprising that so many induced steatosis, however, at high concentrations relative to real exposures. Interestingly, orotic acid did not induce steatosis but disrupted mitochon- drial membrane potential at extrapolated doses overlapping with actual exposure values (Table 2). Indeed, orotic acid is known to induce fatty liver in rats, whereas other animal species including mice and monkeys did not develop fatty liver upon orotic acid exposure (Durschlag and Robinson 1980; Löffler et al. 2015). These findings are consistent with a concern that orotic acid might be a potential risk factor for hepatotoxicity and steatosis in humans, however further investigations are needed. For the other two end- points, namely oxidative stress and apoptosis, none of the tested compounds appeared to induce responses. This may be related to the selection criteria, which was not directly related to associated molecular initiating events, or may be related to the requirement for longer exposure duration to stimulate these further cellular effects. In vitro response data typically describe cell population average responses, missing potentially relevant information about how sub-populations might react to chemical exposure (Singh et al. 2014). In our study, we performed single cell analysis and investigated how hepatocytes and cholangio- cyte-like cells responded on a single cell level to tetracycline exposure. For tetracycline in particular, we could readily observe from preliminary measurements that responses of individual cells were heterogeneous, so we quantified the distribution. To our knowledge, this is the first time that chemically exposed HepaRG cells were analyzed on a sin- gle cell level, discriminating cellular responses in hepato- cyte and cholangiocyte-like cells. Interestingly, hepatocytes always responded more than cholangiocyte-like cells, pos- sibly because hepatocytes may take up chemicals more efficiently than cholangiocyte-like cells do. Furthermore, hepatocytes with larger lipid droplets were rather protected from a decrease in mitochondrial membrane potential com- pared to hepatocytes with small- and medium-sized lipid droplets. Likewise, hepatocytes with larger lipid droplets tended to have lower oxidative stress levels than hepato- cytes with smaller ones suggesting a protective effects of larger sized lipid droplets (Jarc and Petan 2019). Further mechanistic evidence showed that lipid droplet formation may protect cells against fatty acid-induced lipid toxicity. For example, Listenberger et al. showed that excess palmi- tate, which could not be stored in lipid droplets, induced apoptosis (CHO cells) whereas cells were protected from apoptosis when palmitate was integrated into triglycerides and stored in lipid droplets (Listenberger et al. 2003). Fur- thermore palmitic acid decreased mitochondrial membrane potential in HepG2 and THP-1 cells (Li et al. 2018; Alnahdi et al. 2019). Tetracycline is known to decrease β-oxidation and free fatty acids efflux via MTP inhibition. This could lead to palmitic acid accumulating in the cytoplasm and thus further increase ROS and decrease mitochondrial membrane potential. Hepatocytes, which appear to have the capacity to shuttle excess palmitic acids into lipid droplets, may be protected from such stress events. This kind of analysis can be extended and used as a basis to formulate new hypothesis about mechanisms of toxicity. PBPK models can be used to perform in vitro to in vivo extrapolation and predict possible human exposure sce- narios that connect observed in vitro effects with human exposure levels via predicted OED values (Blaauboer 2010). In our study, predicted Css values overlapped with Css values observed in vivo for all chemicals except BNF and WY14643 (Fig. 7A). Importantly, for compounds with available human exposure data, estimated OEDs were within the range of published exposure values. Neverthe- less, nominal in vitro concentrations for selected chemicals were higher compared to in vivo Css (e.g., metformin or valproic acid), leading to a high estimate of the OED. The following assumptions were made during IVIVE performed in this study: (1) Constant dose rate of 1 mg/kg/d with com- plete absorption for every chemical and excretion was lim- ited to the renal route; (2) Plasma protein binding was not considered when in vitro assay concentrations were used for deriving BMC; (3) Metabolism of the chemicals was not accounted for. Aside from the quantitative aspects, there are also uncertainties associated with IVIVE regarding transla- tion of whether cellular effects observed in vitro also occur in vivo, or whether adaptive or further detrimental responses at the organ level are relevant. Quantitative relationships between higher level key events in the AOP are needed to further link the cellular phenotypes with human disease. This would need to be further validated in more complex and physiologically relevant models and compared to relevant 1717Archives of Toxicology (2023) 97:1701–1721 1 3 human exposures, including in susceptible populations, to inform risk assessment. While the use of HepaRG cells offers benefits of meta- bolic competency and reproducibility, there are limitations. For example, they were derived from a single donor, mak- ing it impossible to study population variability. In addition, despite being metabolically active, they do not express high levels of CYP2D6, an enzyme responsible for about 25% of drug metabolism, due to being derived from a poor CYP2D6 metabolizer patient (Guillouzo et al. 2007). Finally, they are differentiated using DMSO, a histone deacetylase inhibitor that may interfere with drug metabolism (Wang et al. 2019). In further studies, therefore, the approach may be adapted for more complex in vitro models such as primary human hepat- ocytes from donor pools or from multiple single donors, sandwich cultures, liver spheroids from primary human cells (including hepatocytes, Kupffer cells and hepatic stel- late cells), liver-on-a-chip or body-on-a-chip models which combine liver spheroids with additional tissues such as adi- pose tissue or the gut with/without bacterial co-cultivation to potentially improve the prediction accuracy and biologi- cal relevance of the hepatic in vitro models to progressive NAFLD. The results presented here comprise a new method com- bining high-content imaging of chemically exposed human liver cells and IVIVE to prioritize chemicals for further evaluation for their potential to induce steatosis and hepa- totoxicity. Further refinement of the specificity and sensi- tivity of the in vitro model and methodology may include measuring additional key events related to the liver steatosis AOP such as ER stress, fatty acid influx/efflux or triglyc- eride synthesis, adapting the approach for high-throughput screening, testing more types of chemicals (including mix- tures of chemicals), and evaluating the temporal dynamics of shorter/longer exposure scenarios. Of 14 food-relevant exposures and pesticides tested with this method, a num- ber could induce lipid accumulation, but at levels outside the range of realistic exposures, and none induced oxida- tive stress or cell death, suggesting the need to screen larger numbers of compounds, at varying durations. Nonetheless, of these compounds, orotic acid was flagged for its capac- ity to disrupt mitochondrial membrane potential. Thus, aside from expanding the scope of screening, these apical endpoint investigations could be combined with molecular mechanistic data derived from more in-depth transcriptomic and proteomic analyses to identify additional predictive key events and refine mechanistic understanding. Future testing of existing chemicals of concern and emerging chemicals would support a first step in a tiered approach for next gener- ation chemical risk assessment. Furthermore, this approach also could be used in a drug discovery pipeline to screen out molecules that could induce steatosis and/or unwanted changes in the additional hepatotoxicity markers, as well as to gain mechanistic insights for candidate molecules with preclinical in vivo liver steatosis findings. Supplementary Information The online version contains supplemen- tary material available at https:// doi. org/ 10. 1007/ s00204- 023- 03490-8. Acknowledgements The authors acknowledge the financial support from the Swiss Federal Food and Veterinary Office and the support of the Scientific Center for Optical and Electron Microscopy ScopeM of the ETH Zürich. The authors would also like to thank Dr. Katherine Hurley and Shuhuan Zhai for their critical feedback during manuscript writing Funding Open access funding provided by Swiss Federal Institute of Technology Zurich. This work was supported by the Swiss Federal Food and Veterinary Office (Grant Number 4.17.01). Data availability Data are available from the corresponding author upon request. Declarations Conflict of interest Marianna Stamou is an employee of AstraZeneca and has stock ownership and/or stock options or interests in the com- pany. Open Access This article is licensed under a Creative Commons Attri- bution 4.0 International License, which permits use, sharing, adapta- tion, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. References Agency for Toxic Substances and Disease Registry (ATSDR) (2003) Toxicological Profile for Atrazine. US Department of Health and Human Services, Public Health Service, Atlanta, GA Aguilar F, Charrondiere U, Dusemund B et al (2009) Orotic acid salts as sources of orotic acid and various minerals added for nutri- tional purposes to food supplements. EFSA J 7:1–25. https:// doi. org/ 10. 2903/j. efsa. 2009. 1187 Akaike H (1973) Information theory and an extension of the maximum likelihood principle. In: Proceeding of the second international symposium on information theory, BN Petrov and F Caski, eds, Akademiai Kiado, Budapest 267–281. Doi: https:// doi. org/ 10. 1007/ 978-1- 4612- 1694-0_ 15 Allard J, Bucher S, Massart J et al (2020) Drug-induced hepatic stea- tosis in absence of severe mitochondrial dysfunction in HepaRG cells: proof of multiple mechanism-based toxicity. Cell Biol Toxicol. https:// doi. org/ 10. 1007/ s10565- 020- 09537-1 Allen TEH, Goodman JM, Gutsell S, Russell PJ (2014) Defining molec- ular initiating events in the adverse outcome pathway framework for risk assessment. Chem Res Toxicol 27:2100–2112. https:// doi. org/ 10. 1021/ TX500 345J Alnahdi A, John A, Raza H (2019) Augmentation of glucotoxicity, oxi- dative stress, apoptosis and mitochondrial dysfunction in HepG2 https://doi.org/10.1007/s00204-023-03490-8 http://creativecommons.org/licenses/by/4.0/ https://doi.org/10.2903/j.efsa.2009.1187 https://doi.org/10.2903/j.efsa.2009.1187 https://doi.org/10.1007/978-1-4612-1694-0_15 https://doi.org/10.1007/978-1-4612-1694-0_15 https://doi.org/10.1007/s10565-020-09537-1 https://doi.org/10.1021/TX500345J https://doi.org/10.1021/TX500345J 1718 Archives of Toxicology (2023) 97:1701–1721 1 3 cells by palmitic acid. Nutrients. https:// doi. org/ 10. 3390/ NU110 91979 Andersson TB, Kanebratt KP, Kenna JG (2012) The HepaRG cell line: a unique in vitro tool for understanding drug metabolism and tox- icology in human. Expert Opin Drug Metab Toxicol 8:909–920. https:// doi. org/ 10. 1517/ 17425 255. 2012. 685159 Angrish MM, Mcqueen CA, Cohen-Hubal E et al (2017) Mechanistic toxicity tests based on an adverse outcome pathway network for hepatic steatosis. Toxicol Sci 159:159–169. https:// doi. org/ 10. 1093/ toxsci/ kfx121 Aninat C, Piton A, Glaise D et al (2006) Expression of cytochromes P450, conjugating enzymes and nuclear receptors in human hepatoma HepaRG cells. Drug Metab Dispos 34:75–83. https:// doi. org/ 10. 1124/ dmd. 105. 006759 Anthérieu S, Rogue A, Fromenty B et al (2011) Induction of vesicu- lar steatosis by amiodarone and tetracycline is associated with up-regulation of lipogenic genes in heparg cells. Hepatology 53:1895–1905. https:// doi. org/ 10. 1002/ hep. 24290 Anthérieu S, Chesné C, Li R et al (2012) Optimization of the Hep- aRG cell model for drug metabolism and toxicity studies. Toxi- col in Vitro 26:1278–1285. https:// doi. org/ 10. 1016/J. TIV. 2012. 05. 008 Bai X, Hong W, Cai P et al (2017) Valproate induced hepatic steatosis by enhanced fatty acid uptake and triglyceride synthesis. Toxi- col Appl Pharmacol 324:12–25. https:// doi. org/ 10. 1016/j. taap. 2017. 03. 022 Begriche K, Massart J, Robin M-A et al (2011) Drug-induced toxicity on mitochondria and lipid metabolism: mechanistic diversity and deleterious consequences for the liver. J Hepatol 54:773–794. https:// doi. org/ 10. 1016/J. JHEP. 2010. 11. 006 Benz R, McLaughlin S (1983) The molecular mechanism of action of the proton ionophore FCCP (carbonylcyanide p-trifluorometh- oxyphenylhydrazone). Biophys J 41:381–398. https:// doi. org/ 10. 1016/ S0006- 3495(83) 84449-X Bienenfeld A, Nagler AR, Orlow SJ (2017) Oral antibacterial therapy for acne vulgaris: an evidence-based review. Am J Clin Dermatol 18:469–490 Blaauboer BJ (2010) Biokinetic modeling and in Vitro–in vivo extrapo- lations. J Toxicol Environm Health Part B 13:242–252. https:// doi. org/ 10. 1080/ 10937 404. 2010. 483940 Bucher S, Tête A, Podechard N et al (2018) Co-exposure to benzo[a] pyrene and ethanol induces a pathological progression of liver steatosis in vitro and in vivo. Sci Rep 8:5963. https:// doi. org/ 10. 1038/ s41598- 018- 24403-1 Caiment F, Wolters J, Smit E et al (2020) Valproic acid promotes mito- chondrial dysfunction in primary human hepatocytes in vitro; impact of C/EBPα-controlled gene expression. Archives of Toxi- col. https:// doi. org/ 10. 1007/ S00204- 020- 02835-X Cha JY, Repa JJ (2007) The Liver X Receptor (LXR) and hepatic lipo- genesis: The carbohydrate-response element-binding protein is a target gene of LXR. J Biol Chem 282:743–751. https:// doi. org/ 10. 1074/ jbc. M6050 23200 Chalasani N, Younossi Z, Lavine JE et al (2012) The diagnosis and management of non-alcoholic fatty liver disease: practice guide- line by the american association for the study of liver diseases, american college of gastroenterology, and the american gastro- enterological association. Hepatology 55:2005–2023. https:// doi. org/ 10. 1002/ hep. 25762 Cuchel M, Meagher EA, Theron HDT et al (2013) Efficacy and safety of a microsomal triglyceride transfer protein inhibitor in patients with homozygous familial hypercholesterolaemia: a single-arm, open-label, phase 3 study. The Lancet 381:40–46. https:// doi. org/ 10. 1016/ S0140- 6736(12) 61731-0 Cuykx M, Claes L, Rodrigues RM et al (2018) Metabolomics profil- ing of steatosis progression in HepaRG ® cells using sodium valproate. Toxicol Lett 286:22–30. https:// doi. org/ 10. 1016/j. toxlet. 2017. 12. 015 D’Apolito O, Garofalo D, la Marca G et al (2012) Reference inter- vals for orotic acid in urine, plasma and dried blood spot using hydrophilic interaction liquid chromatography-tandem mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci 883–884:155–160. https:// doi. org/ 10. 1016/j. jchro mb. 2011. 09. 054 D’Souza AA, Devarajan PV (2015) Asialoglycoprotein receptor medi- ated hepatocyte targeting - Strategies and applications. J Control Release 203:126–139 Dietrich J, Schindler M, Lampen A et al (2020) Comparison of long- term versus short-term effects of okadaic acid on the apoptotic status of human HepaRG cells. Chem Biol Interact 317:108937. https:// doi. org/ 10. 1016/J. CBI. 2020. 108937 Donato MT, Tolosa L, Jiménez N et al (2012) High-content imaging technology for the evaluation of drug-induced steatosis using a multiparametric cell-based assay. J Biomol Screen 17:394–400. https:// doi. org/ 10. 1177/ 10870 57111 427586 Donato MT, Jiménez N, Pelechá M, Tolosa L (2022) Oxidative- stress and long-term hepatotoxicity: comparative study in Upcyte human hepatocytes and hepaRG cells. Archives Toxicol 96:1021–1037. https:// doi. org/ 10. 1007/ S00204- 022- 03236-Y Durschlag RP, Robinson JL (1980) Species specificity in the metabolic consequences of orotic acid consumption. J Nutr 110:822–828. https:// doi. org/ 10. 1093/ JN/ 110.4. 822 Food E, Authority S (2009a) Conclusion on pesticide peer review regarding the risk assessment of the active substance carbofuran. EFSA J. https:// doi. org/ 10. 2903/j. efsa. 2009. 310r Food E, Authority S (2009b) Conclusion on the peer review of the pes- ticide risk assessment of the active substance carbosulfan. EFSA J 7:1354. https:// doi. org/ 10. 2903/j. efsa. 2009. 1354 Garber AJ, Duncan TG, Goodman AM et al (1997) Efficacy of met- formin in type II diabetes: Results of a double-blind, placebo- controlled, dose-response trial. Am J Med 103:491–497. https:// doi. org/ 10. 1016/ S0002- 9343(97) 00254-4 Gómez-Lechón MJ, Donato MT, Martínez-Romero A et al (2007) A human hepatocellular in vitro model to investigate steatosis. Chem Biol Interact 165:106–116. https:// doi. org/ 10. 1016/j. cbi. 2006. 11. 004 Graffmann N, Ring S, Kawala M-A et al (2016) Modeling nonalco- holic fatty liver disease with human pluripotent stem cell-derived immature hepatocyte-like cells reveals activation of PLIN2 and confirms regulatory functions of peroxisome proliferator-acti- vated receptor alpha. Stem Cells Dev 25:1119–1133. https:// doi. org/ 10. 1089/ scd. 2015. 0383 Guillouzo A, Corlu A, Aninat C et al (2007) The human hepatoma HepaRG cells: a highly differentiated model for studies of liver metabolism and toxicity of xenobiotics. Chem Biol Interact 168:66–73. https:// doi. org/ 10. 1016/j. cbi. 2006. 12. 003 Hartung T (2018) Perspectives on In Vitro to In Vivo Extrapolations. Appl in Vitro Toxicol 4:305. https:// doi. org/ 10. 1089/ AIVT. 2016. 0026 Holubarsch CJF, Rohrbach M, Karrasch M et al (2007) A double-blind randomized multicentre clinical trial to evaluate the efficacy and safety of two doses of etomoxir in comparison with placebo in patients with moderate congestive heart failure: the ERGO (etomoxir for the recovery of glucose oxidation) study. Clin Sci 113:205–212. https:// doi. org/ 10. 1042/ CS200 60307 Hu OY-P, Wu C-Y, Chan W-K et al (1996) A pharmacokinetic study with the high-dose anticancer agent menadione in rabbits. Biop- harm Drug Dispos 17:493–499. https:// doi. org/ 10. 1002/ (SICI) 1099- 081X(199608) 17:6% 3c493:: AID- BDD972% 3e3.0. CO;2-6 Jarc E, Petan T (2019) Lipid droplets and the management of cellular stress. Yale J Biol Med 92:435–452 https://doi.org/10.3390/NU11091979 https://doi.org/10.3390/NU11091979 https://doi.org/10.1517/17425255.2012.685159 https://doi.org/10.1093/toxsci/kfx121 https://doi.org/10.1093/toxsci/kfx121 https://doi.org/10.1124/dmd.105.006759 https://doi.org/10.1124/dmd.105.006759 https://doi.org/10.1002/hep.24290 https://doi.org/10.1016/J.TIV.2012.05.008 https://doi.org/10.1016/J.TIV.2012.05.008 https://doi.org/10.1016/j.taap.2017.03.022 https://doi.org/10.1016/j.taap.2017.03.022 https://doi.org/10.1016/J.JHEP.2010.11.006 https://doi.org/10.1016/S0006-3495(83)84449-X https://doi.org/10.1016/S0006-3495(83)84449-X https://doi.org/10.1080/10937404.2010.483940 https://doi.org/10.1080/10937404.2010.483940 https://doi.org/10.1038/s41598-018-24403-1 https://doi.org/10.1038/s41598-018-24403-1 https://doi.org/10.1007/S00204-020-02835-X https://doi.org/10.1074/jbc.M605023200 https://doi.org/10.1074/jbc.M605023200 https://doi.org/10.1002/hep.25762 https://doi.org/10.1002/hep.25762 https://doi.org/10.1016/S0140-6736(12)61731-0 https://doi.org/10.1016/S0140-6736(12)61731-0 https://doi.org/10.1016/j.toxlet.2017.12.015 https://doi.org/10.1016/j.toxlet.2017.12.015 https://doi.org/10.1016/j.jchromb.2011.09.054 https://doi.org/10.1016/j.jchromb.2011.09.054 https://doi.org/10.1016/J.CBI.2020.108937 https://doi.org/10.1177/1087057111427586 https://doi.org/10.1007/S00204-022-03236-Y https://doi.org/10.1093/JN/110.4.822 https://doi.org/10.2903/j.efsa.2009.310r https://doi.org/10.2903/j.efsa.2009.1354 https://doi.org/10.1016/S0002-9343(97)00254-4 https://doi.org/10.1016/S0002-9343(97)00254-4 https://doi.org/10.1016/j.cbi.2006.11.004 https://doi.org/10.1016/j.cbi.2006.11.004 https://doi.org/10.1089/scd.2015.0383 https://doi.org/10.1089/scd.2015.0383 https://doi.org/10.1016/j.cbi.2006.12.003 https://doi.org/10.1089/AIVT.2016.0026 https://doi.org/10.1089/AIVT.2016.0026 https://doi.org/10.1042/CS20060307 https://doi.org/10.1002/(SICI)1099-081X(199608)17:6%3c493::AID-BDD972%3e3.0.CO;2-6 https://doi.org/10.1002/(SICI)1099-081X(199608)17:6%3c493::AID-BDD972%3e3.0.CO;2-6 1719Archives of Toxicology (2023) 97:1701–1721 1 3 Jennings P, Schwarz M, Landesmann B et  al (2014) SEURAT-1 liver gold reference compounds: a mechanism-based review. Arch Toxicol 88:2099–2133. https:// doi. org/ 10. 1007/ s00204- 014- 1410-8 Jones TR, Carpenter A, Golland P (2005) Voronoi-based segmentation of cells on image manifolds. In: Liu Y, Jiang T, Zhang C (eds) Computer vision for biomedical image applications. Springer, Berlin, Heidelberg Jones TR, Kang IH, Wheeler DB et al (2008) Cell Profiler Analyst: data exploration and analysis software for complex image- based screens. BMC Bioinform 9:482. https:// doi. org/ 10. 1186/ 1471- 2105-9- 482 Jossé R, Aninat C, Glaise D et al (2008) Long-term functional stability of human HepaRG hepatocytes and use for chronic toxicity and genotoxicity studies. Drug Metab Dispos 36:1111–1118. https:// doi. org/ 10. 1124/ dmd. 107. 019901 Kaiser JP, Lipscomb JC, Wesselkamper SC (2012) Putative mecha- nisms of environmental chemical-induced steatosis. Int J Toxicol 31:551–563. https:// doi. org/ 10. 1177/ 10915 81812 466418 Kim HS, Ren G, Kim T et al (2020) Metformin reduces saturated fatty acid-induced lipid accumulation and inflammatory response by restoration of autophagic flux in endothelial cells. Sci Rep. https:// doi. org/ 10. 1038/ s41598- 020- 70347-w Klein S, Maggioni S, Bucher J et al (2016) In Silico modeling for the prediction of dose and pathway-related adverse effects in humans from in vitro repeated-dose studies. Toxicol Sci 149:55–66. https:// doi. org/ 10. 1093/ toxsci/ kfv218 Kleiner DE, Gaffey MJ, Sallie R et  al (1997) Histopathologic changes associated with fialuridine hepatotoxicity. Mod Pathol 10:192–199 Koch-Weser J, Browne TR (1980) Valproic acid. N Engl J Med 302:661–666. https:// doi. org/ 10. 1056/ NEJM1 98003 20302 1204 Lambert CB, Spire C, Claude N, Guillouzo A (2009) Dose- and time- dependent effects of phenobarbital on gene expression profil- ing in human hepatoma HepaRG cells. Toxicol Appl Pharmacol 234:345–360. https:// doi. org/ 10. 1016/j. taap. 2008. 11. 008 Lasch A, Marx-Stoelting P, Al B, Lichtenstein D (2021) More than additive effects on liver triglyceride accumulation by com- binations of steatotic and non-steatotic pesticides in HepaRG cells. Arch Toxicol 95:1397–1411. https:// doi. org/ 10. 1007/ S00204- 021- 02997-2 Lee JH, Wada T, Febbraio M et al (2010) A novel role for the dioxin receptor in fatty acid metabolism and hepatic steatosis. Gastro- enterology 139:653–663. https:// doi. org/ 10. 1053/j. gastro. 2010. 03. 033 Lewis W, Day BJ, Copeland WC (2003) Mitochondrial toxicity of NRTI antiviral drugs: an integrated cellular perspective. Nat Rev Drug Discov 2:812–822 Li Y, Xu C, Yu C et al (2009) Association of serum uric acid level with non-alcoholic fatty liver disease: a cross-sectional study. J Hepa- tol 50:1029–1034. https:// doi. org/ 10. 1016/j. jhep. 2008. 11. 021 Li H, Xiao Y, Tang L et al (2018) Adipocyte fatty acid-binding pro- tein promotes palmitate-induced mitochondrial dysfunction and apoptosis in macrophages. Front Immunol 9:81. https:// doi. org/ 10. 3389/ FIMMU. 2018. 00081 Lichtenstein D, Mentz A, Schmidt FF et al (2020) Transcript and pro- tein marker patterns for the identification of steatotic compounds in human HepaRG cells. Food and Chem Toxicol. https:// doi. org/ 10. 1016/J. FCT. 2020. 111690 Lin M, Zhao S, Shen L, Xu D (2014) Potential approaches to amelio- rate hepatic fat accumulation seen with MTP inhibition. Drug Saf 37:213–224 Listenberger LL, Han X, Lewis SE et al (2003) Triglyceride accu- mulation protects against fatty acid-induced lipotoxicity. Proc Natl Acad Sci 100:3077–3082. https:// doi. org/ 10. 1073/ PNAS. 06305 88100 Löffler M, Carrey EA, Zameitat E (2015) Orotic acid, more than just an intermediate of pyrimidine de novo synthesis. J Genet Genomics 42:207–219. https:// doi. org/ 10. 1016/J. JGG. 2015. 04. 001 Loor G, Kondapalli J, Schriewer JM et al (2010) Menadione triggers cell death through ROS-dependent mechanisms involving PARP activation without requiring apoptosis. Free Radic Biol Med 49:1925–1936. https:// doi. org/ 10. 1016/j. freer adbio med. 2010. 09. 021 Lübberstedt M, Müller-Vieira U, Mayer M et al (2011) HepaRG human hepatic cell line utility as a surrogate for primary human hepato- cytes in drug metabolism assessment in vitro. J Pharmacol Toxi- col Methods 63:59–68. https:// doi. org/ 10. 1016/j. vascn. 2010. 04. 013 Luckert C, Braeuning A, de Sousa G et al (2018) Adverse outcome pathway-driven analysis of liver steatosis in vitro : a case study with cyproconazole. Chem Res Toxicol 31:784–798. https:// doi. org/ 10. 1021/ acs. chemr estox. 8b001 12 Lv Y, Hu Q, Shi M et al (2020) The role of PSMB5 in sodium arsenite– induced oxidative stress in L-02 cells. Cell Stress Chaperones. https:// doi. org/ 10. 1007/ S12192- 020- 01104-1 Manning FJ, Swartz M (1995) Review of the fialuridine (FIAU) clinical trials. The National Academies Press, Washington, DC Masarone M, Rosato V, Dallio M et al (2018) Role of oxidative stress in pathophysiology of nonalcoholic fatty liver disease. Oxid Med Cell Longev. https:// doi. org/ 10. 1155/ 2018/ 95476 13 McKenzie R, Fried MW, Sallie R et al (1995) Hepatic failure and lactic acidosis due to fialuridine (FIAU), an investigational nucleoside analogue for chronic hepatitis B. N Engl J Med 333:1099–1105. https:// doi. org/ 10. 1056/ NEJM1 99510 26333 1702 McQuin C, Goodman A, Chernyshev V et al (2018) Cell profiler 30: next-generation image processing for biology. PLoS Biol 16:e2005970. https:// doi. org/ 10. 1371/ journ al. pbio. 20059 70 Mellor CL, Steinmetz FP, Cronin MTD (2015) The identification of nuclear receptors associated with hepatic steatosis to develop and extend adverse outcome pathways. Crit Rev Toxicol 8444:1–15. https:// doi. org/ 10. 3109/ 10408 444. 2015. 10894 71 Mellor CL, Steinmetz FP, Cronin MTD (2016) Using molecular initiat- ing events to develop a structural alert based screening workflow for nuclear receptor ligands associated with hepatic steatosis. Chem Res Toxicol 29:203–212. https:// doi. org/ 10. 1021/ acs. chemr estox. 5b004 80 Merrill CL, Ni H, Yoon LW et al (2002) Etomoxir-induced oxidative stress in HepG2 cells detected by differential gene expression is confirmed biochemically. Toxicol Sci 68:93–101. https:// doi. org/ 10. 1093/ toxsci/ 68.1. 93 Michaut A, le Guillou D, Moreau C et al (2016) A cellular model to study drug-induced liver injury in nonalcoholic fatty liver dis- ease: application to acetaminophen. Toxicol Appl Pharmacol 292:40–55. https:// doi. org/ 10. 1016/j. taap. 2015. 12. 020 Mitro N, Vargas L, Romeo R et al (2007) T0901317 is a potent PXR ligand: implications for the biology ascribed to LXR. FEBS Lett 581:1721–1726. https:// doi. org/ 10. 1016/j. febsl et. 2007. 03. 047 Moya M, José Gómez-Lechón M, Castell JV, Jover R (2010) Enhanced steatosis by nuclear receptor ligands: A study in cultured human hepatocytes and hepatoma cells with a characterized nuclear receptor expression profile. Chem Biol Interact 184:376–387. https:// doi. org/ 10. 1016/j. cbi. 2010. 01. 008 Müller K, Wickham H, James DA, Falcon S (2020) RSQLite: SQLite Interface for R. https:// rsqli te.r- dbi. org, https:// github. com/r- dbi/ RSQLi te Musso G, Cassader M, Gambino R (2016) Non-alcoholic steatohepa- titis: emerging molecular targets and therapeutic strategies. Nat https://doi.org/10.1007/s00204-014-1410-8 https://doi.org/10.1007/s00204-014-1410-8 https://doi.org/10.1186/1471-2105-9-482 https://doi.org/10.1186/1471-2105-9-482 https://doi.org/10.1124/dmd.107.019901 https://doi.org/10.1124/dmd.107.019901 https://doi.org/10.1177/1091581812466418 https://doi.org/10.1038/s41598-020-70347-w https://doi.org/10.1093/toxsci/kfv218 https://doi.org/10.1056/NEJM198003203021204 https://doi.org/10.1016/j.taap.2008.11.008 https://doi.org/10.1007/S00204-021-02997-2 https://doi.org/10.1007/S00204-021-02997-2 https://doi.org/10.1053/j.gastro.2010.03.033 https://doi.org/10.1053/j.gastro.2010.03.033 https://doi.org/10.1016/j.jhep.2008.11.021 https://doi.org/10.3389/FIMMU.2018.00081 https://doi.org/10.3389/FIMMU.2018.00081 https://doi.org/10.1016/J.FCT.2020.111690 https://doi.org/10.1016/J.FCT.2020.111690 https://doi.org/10.1073/PNAS.0630588100 https://doi.org/10.1073/PNAS.0630588100 https://doi.org/10.1016/J.JGG.2015.04.001 https://doi.org/10.1016/j.freeradbiomed.2010.09.021 https://doi.org/10.1016/j.freeradbiomed.2010.09.021 https://doi.org/10.1016/j.vascn.2010.04.013 https://doi.org/10.1016/j.vascn.2010.04.013 https://doi.org/10.1021/acs.chemrestox.8b00112 https://doi.org/10.1021/acs.chemrestox.8b00112 https://doi.org/10.1007/S12192-020-01104-1 https://doi.org/10.1155/2018/9547613 https://doi.org/10.1056/NEJM199510263331702 https://doi.org/10.1371/journal.pbio.2005970 https://doi.org/10.3109/10408444.2015.1089471 https://doi.org/10.1021/acs.chemrestox.5b00480 https://doi.org/10.1021/acs.chemrestox.5b00480 https://doi.org/10.1093/toxsci/68.1.93 https://doi.org/10.1093/toxsci/68.1.93 https://doi.org/10.1016/j.taap.2015.12.020 https://doi.org/10.1016/j.febslet.2007.03.047 https://doi.org/10.1016/j.cbi.2010.01.008 https://rsqlite.r-dbi.org https://github.com/r-dbi/RSQLite https://github.com/r-dbi/RSQLite 1720 Archives of Toxicology (2023) 97:1701–1721 1 3 Rev Drug Discov 15:249–274. https:// doi. org/ 10. 1038/ nrd. 2015.3 Panno MD, Cefalù AB, Averna MR (2014) Lomitapide: a novel drug for homozygous familial hypercholesterolemia. Clin Lipidol 9:19–32. https:// doi. org/ 10. 2217/ clp. 13. 74 Parent R, Marion MJ, Furio L et al (2004) Origin and characterization of a human bipotent liver progenitor cell line. Gastroenterology 126:1147–1156. https:// doi. org/ 10. 1053/j. gastro. 2004. 01. 002 Pavlik L, Regev A, Ardayfio PA, Chalasani NP (2019) Drug-induced steatosis and steatohepatitis: the search for novel serum biomark- ers among potential biomarkers for non-alcoholic fatty liver dis- ease and non-alcoholic steatohepatitis. Drug Saf 42(6):701–711. https:// doi. org/ 10. 1007/ S40264- 018- 00790-2 Pearce RG, Setzer RW, Strope CL et al (2017) httk : r package for high-throughput toxicokinetics. J Stat Softw. https:// doi. org/ 10. 18637/ jss. v079. i04 Pernicova I, Korbonits M (2014) Metformin-mode of action and clini- cal implications for diabetes and cancer. Nat Rev Endocrinol 10:143–156 Persson M, Løye AF, Mow T, Hornberg JJ (2013) A high content screening assay to predict human drug-induced liver injury dur- ing drug discovery. J Pharmacol Toxicol Methods 68:302–313. https:// doi. org/ 10. 1016/j. vascn. 2013. 08. 001 Pollinger J, Merk D (2017) Therapeutic applications of the versatile fatty acid mimetic WY14643. Expert Opin Ther Pat 27:517–525 Richard AM, Judson RS, Houck KA et al (2016) ToxCast chemical landscape: paving the road to 21st century toxicology. Chem Res Toxicol 29:1225–1251 Ritz C, Baty F, Streibig JC, Gerhard D (2015) Dose-response analysis using R. PLoS ONE 10:e0146021. https:// doi. org/ 10. 1371/ journ al. pone. 01460 21 Rogue A, Anthérieu S, Vluggens A et al (2014) PPAR agonists reduce steatosis in oleic acid-overloaded HepaRG cells. Toxicol Appl Pharmacol 276:73–81. https:// doi. org/ 10. 1016/j. taap. 2014. 02. 001 Rotroff DM, Wetmore BA, Dix DJ et al (2010) Incorporating human dosimetry and exposure into high-throughput in vitro toxicity screening. Toxicol Sci 117:348–358. https:// doi. org/ 10. 1093/ toxsci/ kfq220 Saito J, Okamura A, Takeuchi K et al (2016) High content analysis assay for prediction of human hepatotoxicity in HepaRG and HepG2 cells. Toxicol in Vitro 33:63–70. https:// doi. org/ 10. 1016/j. tiv. 2016. 02. 019 Sakamuru S, Attene-Ramos MS, Xia M (2016) Mitochondrial mem- brane potential assay. High-throughput screening assays in toxi- cology. Methods Molec Biol 1473:17–22. https:// doi. org/ 10. 1007/ 978-1- 4939- 6346-1_2 Sankur B (2004) Survey over image thresholding techniques and quantitative performance evaluation. J Electron Imaging 13:146. https:// doi. org/ 10. 1117/1. 16313 15 Schindelin J, Arganda-Carreras I, Frise E et al (2012) Fiji: an open- source platform for biological-image analysis. Nat Methods 9:676–682 Schumacher JD, Guo GL (2015) Mechanistic review of drug-induced steatohepatitis. Toxicol Appl Pharmacol 289:40–47. https:// doi. org/ 10. 1016/j. taap. 2015. 08. 022 ThermoFisher Scientific (2011) Probes for reactive oxygen species, including nitric oxide. in: molecular probes handbook-a guide to fluorescent probes and labeling technologies. //www. therm ofish er. com/ uk/ en/ home/ refer ences/ molec ular- probes- the- handb ook/ mp- handb ook- downl oad. html. Accessed 8 Dec 2020 Shah I, Antonijevic T, Chambers B et al (2021) Estimating hepatotoxic doses using high-content imaging in primary hepatocytes. Toxi- col Sci. https:// doi. org/ 10. 1093/ TOXSCI/ KFAB0 91 Sharma S, Mells JE, Fu PP et al (2011) GLP-1 analogs reduce hepato- cyte steatosis and improve survival by enhancing the unfolded protein response and promoting macroautophagy. PLoS ONE 6:e25269. https:// doi. org/ 10. 1371/ journ al. pone. 00252 69 Siddiqui MA, Ahmad J, Farshori NN et  al (2013) Rotenone- induced oxidative stress and apoptosis in human liver HepG2 cells. Mol Cell Biochem 384:59–69. https:// doi. org/ 10. 1007/ s11010- 013- 1781-9 Siddoway LA (2003) Amiodarone: guidelines for use and monitoring. Am Fam Physician 68(11):2189–2196 Singh S, Carpenter AE, Genovesio A (2014) Increasing the content of high-content screening: an overview. J Biomol Screen 19:640– 650. https:// doi. org/ 10. 1177/ 10870 57114 528537 Sipes NS, Wambaugh JF, Pearce R et al (2017) An intuitive approach for predicting potential human health risk with the Tox21 10k library. Environ Sci Technol 51:10786–10796. https:// doi. org/ 10. 1021/ acs. est. 7b006 50 Soret P-A, Magusto J, Housset C, Gautheron J (2020) In vitro and in vivo models of non-alcoholic fatty liver disease: a critical appraisal. J Clin Med 10:1–18. https:// doi. org/ 10. 3390/ JCM10 010036 Tolosa L, Gómez-Lechón MJ, Jiménez N et al (2016) Advantageous use of HepaRG cells for the screening and mechanistic study of drug-induced steatosis. Toxicol Appl Pharmacol 302:1–9. https:// doi. org/ 10. 1016/j. taap. 2016. 04. 007 United States Environmental Protection Agency (USEPA) (2007) Reregistration Eligibility Decision for Rotenone. US EPA, Pes- ticides, and Toxic substances, Special review and Reregistration Division 259. Doi: EPA738-R-07-005 United States Environmental Protection Agency (USEPA) (2012) Benchmark Dose Technical Guidance. US EPA, Washinton, DC 1–99 United States Environmental Protection Agency (2023) Benchmark Dose Software (BMDS) (Build 3.3.2; Model Library Version 2023.03.1) [Computer Software]. Available from https:// www. epa. gov/ bmds/ downl oad- bmds Vinken M (2015) Adverse outcome pathways and drug-induced liver injury testing. Chem Res Toxicol 28:1391–1397. https:// doi. org/ 10. 1021/ acs. chemr estox. 5b002 08 Wahlang B, Beier JI, Clair HB et al (2013) Toxicant-associated Stea- tohepatitis. Toxicol Pathol. https:// doi. org/ 10. 1177/ 01926 23312 468517 Wang ZY, Li WJ, Li QG et al (2019) A DMSO-free hepatocyte matura- tion medium accelerates hepatic differentiation of HepaRG cells in vitro. Biomed Pharmacotherapy 116:109010. https:// doi. org/ 10. 1016/J. BIOPHA. 2019. 109010 Weininger D (1988) SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules. J Chem Inf Comput Sci 28:31–36. https:// doi. org/ 10. 1021/ CI000 57A005 Wetmore BA (2015) Quantitative in vitro-to-in vivo extrapolation in a high-throughput environment. Toxicology 332:94–101. https:// doi. org/ 10. 1016/j. tox. 2014. 05. 012 Wetmore BA, Wambaugh JF, Allen B et  al (2015) Incorporating high-throughput exposure predictions with dosimetry-adjusted in vitro bioactivity to inform chemical toxicity testing. Toxicol Sci 148:121–136. https:// doi. org/ 10. 1093/ toxsci/ kfv171 Wickham H (2016) ggplot2: elegant graphics for data analysis. Springer-Verlag, New York Willebrords J, Pereira IVA, Maes M et al (2015) Strategies, models and biomarkers in experimental non-alcoholic fatty liver disease research. Prog Lipid Res 59:106–125. https:// doi. org/ 10. 1016/j. plipr es. 2015. 05. 002 Woods CG, Burns AM, Bradford BU et al (2007) WY-14,643-induced cell proliferation and oxidative stress in mouse liver are inde- pendent of NADPH oxidase. Toxicol Sci 98:366–374. https:// doi. org/ 10. 1093/ toxsci/ kfm104 https://doi.org/10.1038/nrd.2015.3 https://doi.org/10.1038/nrd.2015.3 https://doi.org/10.2217/clp.13.74 https://doi.org/10.1053/j.gastro.2004.01.002 https://doi.org/10.1007/S40264-018-00790-2 https://doi.org/10.18637/jss.v079.i04 https://doi.org/10.18637/jss.v079.i04 https://doi.org/10.1016/j.vascn.2013.08.001 https://doi.org/10.1371/journal.pone.0146021 https://doi.org/10.1371/journal.pone.0146021 https://doi.org/10.1016/j.taap.2014.02.001 https://doi.org/10.1016/j.taap.2014.02.001 https://doi.org/10.1093/toxsci/kfq220 https://doi.org/10.1093/toxsci/kfq220 https://doi.org/10.1016/j.tiv.2016.02.019 https://doi.org/10.1016/j.tiv.2016.02.019 https://doi.org/10.1007/978-1-4939-6346-1_2 https://doi.org/10.1007/978-1-4939-6346-1_2 https://doi.org/10.1117/1.1631315 https://doi.org/10.1016/j.taap.2015.08.022 https://doi.org/10.1016/j.taap.2015.08.022 http://www.thermofisher.com/uk/en/home/references/molecular-probes-the-handbook/mp-handbook-download.html http://www.thermofisher.com/uk/en/home/references/molecular-probes-the-handbook/mp-handbook-download.html http://www.thermofisher.com/uk/en/home/references/molecular-probes-the-handbook/mp-handbook-download.html https://doi.org/10.1093/TOXSCI/KFAB091 https://doi.org/10.1371/journal.pone.0025269 https://doi.org/10.1007/s11010-013-1781-9 https://doi.org/10.1007/s11010-013-1781-9 https://doi.org/10.1177/1087057114528537 https://doi.org/10.1021/acs.est.7b00650 https://doi.org/10.1021/acs.est.7b00650 https://doi.org/10.3390/JCM10010036 https://doi.org/10.3390/JCM10010036 https://doi.org/10.1016/j.taap.2016.04.007 https://doi.org/10.1016/j.taap.2016.04.007 https://www.epa.gov/bmds/download-bmds https://www.epa.gov/bmds/download-bmds https://doi.org/10.1021/acs.chemrestox.5b00208 https://doi.org/10.1021/acs.chemrestox.5b00208 https://doi.org/10.1177/0192623312468517 https://doi.org/10.1177/0192623312468517 https://doi.org/10.1016/J.BIOPHA.2019.109010 https://doi.org/10.1016/J.BIOPHA.2019.109010 https://doi.org/10.1021/CI00057A005 https://doi.org/10.1021/CI00057A005 https://doi.org/10.1016/j.tox.2014.05.012 https://doi.org/10.1016/j.tox.2014.05.012 https://doi.org/10.1093/toxsci/kfv171 https://doi.org/10.1016/j.plipres.2015.05.002 https://doi.org/10.1016/j.plipres.2015.05.002 https://doi.org/10.1093/toxsci/kfm104 https://doi.org/10.1093/toxsci/kfm104 1721Archives of Toxicology (2023) 97:1701–1721 1 3 Younossi Z, Anstee QM, Marietti M et al (2018) Global burden of NAFLD and NASH: trends, predictions, risk factors and preven- tion. Nat Rev Gastroenterol Hepatol 15:11–20 Zare M, Panahi G, Koushki M et al (2019) Metformin reduces lipid accumulation in HepG2 cells via downregulation of miR-33b. Arch Physiol Biochem. https:// doi. org/ 10. 1080/ 13813 455. 2019. 16807 00 Zhao T, Li X, Sun D, Zhang Z (2019) Oxidative stress: one potential factor for arsenite-induced increase of N6-methyladenosine in human keratinocytes. Environ Toxicol Pharmacol 69:95–103. https:// doi. org/ 10. 1016/J. ETAP. 2019. 04. 005 Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. https://doi.org/10.1080/13813455.2019.1680700 https://doi.org/10.1080/13813455.2019.1680700 https://doi.org/10.1016/J.ETAP.2019.04.005 In vitro to in vivo extrapolation and high-content imaging for simultaneous characterization of chemically induced liver steatosis and markers of hepatotoxicity Abstract Introduction Materials and methods Chemicals and reagents Cell culture and chemical exposure Cell viability assessment High-content imaging assay Automated quantitative image analysis Data analysis of single cell data Machine learning approach to quantify apoptotic cells Dose–response data and model fitting Benchmark concentration modeling In vitro to in vivo extrapolation (IVIVE) Results Multiparametric analysis of reference chemicals Analysis of food-related chemicals and pesticides Chemical exposures induce distinct cellular endophenotypes within cell populations IVIVE Discussion Anchor 21 Acknowledgements References