Toxins in Botanical Drugs and Plant-derived Food and Feed - from Science to Regulation: A Workshop Review.

In September 2022, the 3rd International Workshop on pyrrolizidine alkaloids (PAs) and related phytotoxins was held on-line, entitled 'Toxins in botanical drugs and plant-derived food and feed - from science to regulation'. The workshop focused on new findings about the occurrence, exposure, toxicity, and risk assessment of PAs. In addition, new scientific results related to the risk assessment of alkenylbenzenes, a distinct class of herbal constituents, were presented. The presence of PAs and alkenylbenzenes in plant-derived food, feed, and herbal medicines has raised health concerns with respect to their acute and chronic toxicity but mainly related to the genotoxic and carcinogenic properties of several congeners. The compounds are natural constituents of a variety of plant families and species widely used in medicinal, food, and feed products. Their individual occurrence, levels, and toxic properties, together with the broad range of congeners present in nature, represent a striking challenge to modern toxicology. This review tries to provide an overview of the current knowledge on these compounds and indicates needs and perspectives for future research.

Dose addition in mixtures of compounds with dissimilar endocrine modes of action in in vitro receptor activation assays and the zebrafish sexual development test.

BACKGROUND: Human exposure to pesticides is being associated with feminisation for which a decrease of the anogenital distance (AGD) is a sensitive endpoint. Dose addition for the cumulative risk assessment of pesticides in food is considered sufficiently conservative for combinations of compounds with both similar and dissimilar modes of action (MoA). OBJECTIVE: The present study was designed to test the dose addition hypothesis in a binary mixture of endocrine active compounds with a dissimilar mode of action for the endpoint feminisation. METHODS: Compounds were selected from a list of chemicals of which exposure is related to a decrease of the AGD in rats and completed with reference compounds. These chemicals were characterised using specific in vitro transcriptional activation (TA) assays for estrogenic and androgenic properties, leading to a final selection of dienestrol as an ER-agonist and flutamide, linuron, and deltamethrin as AR-antagonists. These compounds were then tested in an in vivo model, i.e. in zebrafish (Danio rerio), using sex ratio in the population as an endpoint in order to confirm their feminising effect and MoA. Ultimately, the fish model was used to test a binary mixture of flutamide and dienestrol. RESULTS: Statistical analysis of the binary mixture of flutamide and dienestrol in the fish sexual development tests (FSDT) with zebrafish supported dose addition.

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

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

New approach methodologies: A quantitative in vitro to in vivo extrapolation case study with PFASs.

PER: and polyfluoroalkyl substances (PFASs) have been associated with increased blood lipids in humans. Perfluorooctanoic acid (PFOA) has been also linked with elevated alanine transferase (ALT) serum levels in humans, and in rodents the liver is a main target organ for many PFASs. With the focus on New Approach Methodologies, the chronic oral equivalent effect doses were calculated for PFOA, PFNA (perfluorononanoic acid), PFHxS (perfluorohexanesulfonic acid) and PFOS (perfluorooctane sulfonic acid) based on in vitro effects measured in the HepaRG cell line. Selected in vitro readouts were considered biomarkers for lipid disturbances and hepatotoxicity. Concentration-response data obtained from HepaRG cells on triglyceride (TG) accumulation and expression changes of 12 selected genes (some involved in cholesterol homeostasis) were converted into corresponding human dose-response data, using physiologically based kinetic (PBK) model-facilitated reverse dosimetry. Next to this, the biokinetics of the chemicals were studied in the cell system. The current European dietary PFASs exposure overlaps with the calculated oral equivalent effect doses, indicating that the latter may lead to interference with hepatic gene expression and lipid metabolism. These findings illustrate an in vitro-in silico methodology, which can be applied for more PFASs, to select those that should be prioritized for further hazard characterization.

Perfluoroalkyl substances (PFASs) decrease the expression of recombination-activating genes (RAG1 and RAG2) in human B lymphoma Namalwa cells.

Per- and polyfluoroalkyl substances (PFASs) are omnipresent and have been shown to induce a wide range of adverse effects, including hepatotoxicity, developmental toxicity and immunotoxicity. So far, little information is available about the mechanisms underlying the toxicity of PFASs, including those related to their immunotoxicity. Reported immunotoxic effects of PFASs include decreased antibody responses in experimental animals and humans, indicating that PFASs may, among others, affect B cell function. In the present study, we first assessed the effects of PFOA on the transcriptome of the human Namalwa B cell line using RNA seq analysis. Gene expression changes, analyzed using Ingenuity Pathway Analysis, pointed to various cellular processes affected by PFOA, including 'B cell development' and 'Primary immunodeficiency signaling'. Interestingly, PFOA decreased the expression of RAG1 and RAG2, genes involved in immunoglobulin and T cell receptor V(D)J recombination. As a next step, time- and concentration-dependent changes in the expression of RAG1 and RAG2 upon exposure to PFOA, PFNA, PFHxS and PFOS were studied through RT-qPCR analysis. Analysis with the concentration-response modeling software PROAST resulted in the following potency ranking: PFNA > PFOA > PFOS > PFHxS. Altogether, the present in vitro study provides insights into the effects of selected PFASs on B cells, identifying RAG1 and RAG2 expression as possible relevant targets that may play a role in the immunotoxicity of PFASs.

Bioassay-directed analysis-based identification of relevant pyrrolizidine alkaloids.

Pyrrolizidine alkaloids (PAs) are produced by various plant species and have been detected as contaminants in food and feed. Monitoring programmes should include PAs that are present in relevant matrices and that exhibit a high toxic potential. The aim of the present study was to use a bioassay-directed analysis approach to identify relevant PAs not yet included in monitoring programmes. To that end, extracts of Heliotropium europaeum and H. popovii were prepared and analysed with LC-MS/MS for the presence of 35 PAs included in monitoring programmes, as well as for genotoxic activity in the HepaRG/γH2AX assay. Europine, heliotrine and lasiocarpine were found to be the most abundant PAs. The extracts showed a higher γH2AX activity than related artificial mixtures of quantified known PAs, which might point to the presence of unknown toxic PAs. The H. europaeum extract was fractionated and γH2AX activities of individual fractions were determined. Fractions were further analysed applying LC-Orbitrap-MS analysis and Compound Discoverer software, identifying various candidate PAs responsible for the non-explained genotoxic activity. Altogether, the results obtained show that bioassay-directed analysis allows identification of candidate PAs that can be included in monitoring programmes.

Providing Biological Plausibility for Exposure-Health Relationships for the Mycotoxins Deoxynivalenol (DON) and Fumonisin B1 (FB1) in Humans Using the AOP Framework.

Humans are chronically exposed to the mycotoxins deoxynivalenol (DON) and fumonisin B1 (FB1), as indicated by their widespread presence in foods and occasional exposure in the workplace. This exposure is confirmed by human biomonitoring (HBM) studies on (metabolites of) these mycotoxins in human matrices. We evaluated the exposure-health relationship of the mycotoxins in humans by reviewing the available literature. Since human studies did not allow the identification of unequivocal chronic health effects upon exposure to DON and FB1, the adverse outcome pathway (AOP) framework was used to structure additional mechanistic evidence from in vitro and animal studies on the identified adverse effects. In addition to a preliminary AOP for DON resulting in the adverse outcome (AO) 'reduced body weight gain', we developed a more elaborated AOP for FB1, from the molecular initiating event (MIE) 'inhibition of ceramide synthases' leading to the AO 'neural tube defects'. The mechanistic evidence from AOPs can be used to support the limited evidence from human studies, to focus FB1- and DON-related research in humans to identify related early biomarkers of effect. In order to establish additional human exposure-health relationships in the future, recommendations are given to maximize the information that can be obtained from HBM.

Novel Insights into Pyrrolizidine Alkaloid Toxicity and Implications for Risk Assessment: Occurrence, Genotoxicity, Toxicokinetics, Risk Assessment-A Workshop Report.

This paper reports on the major contributions and results of the 2nd International Workshop of Pyrrolizidine Alkaloids held in September 2020 in Kaiserslautern, Germany. Pyrrolizidine alkaloids are among the most relevant plant toxins contaminating food, feed, and medicinal products of plant origin. Hundreds of PA congeners with widespread occurrence are known, and thousands of plants are assumed to contain PAs. Due to certain PAs' pronounced liver toxicity and carcinogenicity, their occurrence in food, feed, and phytomedicines has raised serious human health concerns. This is particularly true for herbal teas, certain food supplements, honey, and certain phytomedicinal drugs. Due to the limited availability of animal data, broader use of in vitro data appears warranted to improve the risk assessment of a large number of relevant, 1,2-unsaturated PAs. This is true, for example, for the derivation of both toxicokinetic and toxicodynamic data. These efforts aim to understand better the modes of action, uptake, metabolism, elimination, toxicity, and genotoxicity of PAs to enable a detailed dose-response analysis and ultimately quantify differing toxic potencies between relevant PAs. Accordingly, risk-limiting measures comprising production, marketing, and regulation of food, feed, and medicinal products are discussed.

Systemic PFOS and PFOA exposure and disturbed lipid homeostasis in humans: what do we know and what not?

Associations between per- and polyfluoroalkyl substances (PFASs) and increased blood lipids have been repeatedly observed in humans, but a causal relation has been debated. Rodent studies show reverse effects, i.e. decreased blood cholesterol and triglycerides, occurring however at PFAS serum levels at least 100-fold higher than those in humans. This paper aims to present the main issues regarding the modulation of lipid homeostasis by the two most common PFASs, PFOS and PFOA, with emphasis on the underlying mechanisms relevant for humans. Overall, the apparent contrast between human and animal data may be an artifact of dose, with different molecular pathways coming into play upon exposure to PFASs at very low versus high levels. Altogether, the interpretation of existing rodent data on PFOS/PFOA-induced lipid perturbations with respect to the human situation is complex. From a mechanistic perspective, research on human liver cells shows that PFOS/PFOA activate the PPARα pathway, whereas studies on the involvement of other nuclear receptors, like PXR, are less conclusive. Other data indicate that suppression of the nuclear receptor HNF4α signaling pathway, as well as perturbations of bile acid metabolism and transport might be important cellular events that require further investigation. Future studies with human-relevant test systems would help to obtain more insight into the mechanistic pathways pertinent for humans. These studies shall be designed with a careful consideration of appropriate dosing and toxicokinetics, so as to enable biologically plausible quantitative extrapolations. Such research will increase the understanding of possible perturbed lipid homeostasis related to PFOS/ PFOA exposure and the potential implications for human health.

Cytochrome P450 expression, induction and activity in human induced pluripotent stem cell-derived intestinal organoids and comparison with primary human intestinal epithelial cells and Caco-2 cells.

Human intestinal organoids (HIOs) are a promising in vitro model consisting of different intestinal cell types with a 3D microarchitecture resembling native tissue. In the current study, we aimed to assess the expression of the most common intestinal CYP enzymes in a human induced pluripotent stem cell (hiPSC)-derived HIO model, and the suitability of that model to study chemical-induced changes in CYP expression and activity. We compared this model with the commonly used human colonic adenocarcinoma cell line Caco-2 and with a human primary intestinal epithelial cell (IEC)-based model, closely resembling in vivo tissue. We optimized an existing protocol to differentiate hiPSCs into HIOs and demonstrated that obtained HIOs contain a polarized epithelium with tight junctions consisting of enterocytes, goblet cells, enteroendocrine cells and Paneth cells. We extensively characterized the gene expression of CYPs and activity of CYP3A4/5, indicating relatively high gene expression levels of the most important intestinal CYP enzymes in HIOs compared to the other models. Furthermore, we showed that CYP1A1 and CYP1B1 were induced by ß-naphtoflavone in all three models, whereas CYP3A4 was induced by phenobarbital and rifampicin in HIOs, in the IEC-based model (although not statistically significant), but not in Caco-2 cells. Interestingly, CYP2B6 expression was not induced in any of the models by the well-known liver CYP2B6 inducer phenobarbital. In conclusion, our study indicates that hiPSC-based HIOs are a useful in vitro intestinal model to study biotransformation of chemicals in the intestine.

Development of a Web-Based Toolbox to Support Quantitative In-Vitro-to-In-Vivo Extrapolations (QIVIVE) within Nonanimal Testing Strategies.

The goal of the present study was to develop an online web-based toolbox that contains generic physiologically based kinetic (PBK) models for rats and humans, including underlying calculation tools to predict plasma protein binding and tissue:plasma distribution, to be used for quantitative in-vitro-to-in-vivo extrapolations (QIVIVE). The PBK models within the toolbox allow first estimations of internal plasma and tissue concentrations of chemicals to be made, based on the logP and pKa of the chemicals and values for intestinal uptake and intrinsic hepatic clearance. As a case study, the toolbox was used to predict oral equivalent doses of in vitro ToxCast bioactivity data for the food additives methylparaben, propyl gallate, octyl gallate, and dodecyl gallate. These oral equivalent doses were subsequently compared with human exposure estimates, as a low tier assessment allowing prioritization for further assessment. The results revealed that daily intake levels of especially propyl gallate can lead to internal plasma concentrations that are close to in vitro biological effect concentrations, particularly with respect to the inhibition of human thyroid peroxidase (TPO). Estrogenic effects were not considered likely to be induced by the food additives, as daily exposure levels of the different compounds remained 2 orders of magnitude below the oral equivalent doses for in vitro estrogen receptor activation. Overall, the results of the study show how the toolbox, which is freely accessible through www.qivivetools.wur.nl, can be used to obtain initial internal dose estimates of chemicals and to prioritize chemicals for further assessment, based on the comparison of oral equivalent doses of in vitro biological activity data with human exposure levels.

Perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid (PFOS), and perfluorononanoic acid (PFNA) increase triglyceride levels and decrease cholesterogenic gene expression in human HepaRG liver cells.

Per- and polyfluoroalkyl substances (PFASs) are omnipresent in the environment, food chain, and humans. Epidemiological studies have shown a positive association between serum levels of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS), and increased serum cholesterol and, in some cases, also triglyceride levels. However, causality has been questioned, as animal studies, as well as a human trial, showed a decrease in serum cholesterol and no effects or a decrease in plasma triglycerides. To obtain more insight into the effects of PFASs on these processes, the present study investigated the effects of PFOA, PFOS, and perfluorononanoic acid (PFNA) on intracellular triglyceride and cholesterol levels in human HepaRG liver cells. DNA microarray analyses were performed to provide insight into underlying mechanisms. All PFASs induced an increase in cellular triglyceride levels, but had no effect on cholesterol levels. Gene set enrichment analysis (GSEA) of the microarray data indicated that gene sets related to cholesterol biosynthesis were repressed by PFOA, PFOS, and PFNA. Other gene sets commonly affected by all PFAS were related to PERK/ATF4 signaling (induced), tRNA amino-acylation (induced), amino acid transport (induced), and glycolysis/gluconeogenesis (repressed). Moreover, numerous target genes of peroxisome proliferator-activated receptor α (PPARα) were found to be upregulated. Altogether, the present study shows that PFOA, PFOS, and PFNA increase triglyceride levels and inhibit cholesterogenic gene expression in HepaRG cells. In addition, the present study indicates that PFASs induce endoplasmic reticulum stress, which may be an important mechanism underlying some of the toxic effects of these chemicals.

New approach methodologies (NAMs) for human-relevant biokinetics predictions. Meeting the paradigm shift in toxicology towards an animal-free chemical risk assessment

For almost fifteen years, the availability and regulatory acceptance of new approach methodologies (NAMs) to assess the absorption, distribution, metabolism and excretion (ADME/biokinetics) in chemical risk evaluations are a bottleneck. To enhance the field, a team of 24 experts from science, industry, and regulatory bodies, including new generation toxicologists, met at the Lorentz Centre in Leiden, The Netherlands. A range of possibilities for the use of NAMs for biokinetics in risk evaluations were formulated (for example to define species differences and human variation or to perform quantitative in vitro-in vivo extrapolations). To increase the regulatory use and acceptance of NAMs for biokinetics for these ADME considerations within risk evaluations, the development of test guidelines (protocols) and of overarching guidance documents is considered a critical step. To this end, a need for an expert group on biokinetics within the Organisation of Economic Cooperation and Development (OECD) to supervise this process was formulated. The workshop discussions revealed that method development is still required, particularly to adequately capture transporter mediated processes as well as to obtain cell models that reflect the physiology and kinetic characteristics of relevant organs. Developments in the fields of stem cells, organoids and organ-on-a-chip models provide promising tools to meet these research needs in the future.

Hepatotoxicity of the pesticides imazalil, thiacloprid and clothianidin - Individual and mixture effects in a 28-day study in female Wistar rats.

Humans are exposed to pesticide residues through various food products. As these residues can occur in mixtures, there is a need to investigate possible mixture effects on human health. Recent exposure studies revealed the preponderance of imazalil, thiacloprid, and clothianidin in food diets. In this study, we assessed their toxicity alone and in binary mixtures in a 28-day gavage study in female Wistar rats. Five dose levels (up to 350 mg/kg bw/day) ranging from a typical toxicological reference value to a clear effect dose were applied. Data show that the liver was a target organ of all pesticides and their mixtures. Increases in liver weight were observed and histopathological examination revealed centrilobular hepatocellular hypertrophy and cytoplasm degeneration for all treatment conditions. No accumulation of hepatic triglycerides was reported. Tissue residue analysis showed altered pesticide residues in the liver and the kidney when being in mixture as compared to the levels of pesticide residues for the single compound treatment, indicating possible toxicokinetic interactions. Overall, all mixtures appeared to follow the additivity concept, even though quantitative analysis was limited for some endpoints due to the semi-quantitative nature of the data, raising no specific concern for the risk assessment of the examined pesticides.

In vitro toxicological characterisation of the antifungal compound soybean toxin (SBTX).

Soybean toxin (SBTX) is a protein isolated from soybean seeds and composed of two polypeptide subunits (17 and 27 kDa). SBTX has in vitro activity against phytopathogenic fungi such as Cercospora sojina, Aspergillus niger, and Penicillium herguei, and yeasts like Candida albicans, C. parapsilosis, Kluyveromyces marxiannus, and Pichia membranifaciens. The present study aimed to analyze in vitro whether SBTX causes any side effects on non-target bacterial and mammalian cells that could impede its potential use as a novel antifungal agent. SBTX at 100 µg/mL and 200 µg/mL did not hinder the growth of the bacteria Salmonella enterica (subspecies enterica serovar choleraesuis), Bacillus subtilis (subspecies spizizenii) and Staphylococcus aureus. Moreover, SBTX at concentrations up to 500 µg/mL did not significantly affect the viability of erythrocytes, neutrophils, and human intestinal Caco-2 cells. To study whether SBTX could induce relevant alterations in gene expression, in vitro DNA microarray experiments were conducted in which differentiated Caco-2 cells were exposed for 24 h to 100 µg/mL or 200 µg/mL SBTX. SBTX up-regulated genes involved in cell cycle and immune response pathways, but down-regulated genes that play a role in cholesterol biosynthesis and platelet degranulation pathways. Thus, although SBTX did not affect bacteria, nor induced cytotoxity in mammalian cells, it affected some biological pathways in the human Caco-2 cell line that warrants further investigation.

An adverse outcome pathway-based approach to assess steatotic mixture effects of hepatotoxic pesticides in vitro.

Exposure to complex chemical mixtures requires a tiered strategy for efficient mixture risk assessment. As a part of the EuroMix project we developed an adverse outcome pathway (AOP)-based assay toolbox to investigate the combined effects of the liver steatosis-inducing compounds imazalil, thiacloprid, and clothianidin in human HepaRG hepatocarcinoma cells. Compound-specific relative potency factors were determined using a benchmark dose approach. Equipotent mixtures were tested for nuclear receptor activation, gene and protein expression, and triglyceride accumulation, according to the molecular initiating events and key events proposed in the steatosis AOP. All three compounds affected the activity of nuclear receptors, but not key genes/proteins as proposed. Triglyceride accumulation was observed with three different methods. Mixture effects were in agreement with the assumption of dose additivity for all the combinations and endpoints tested. Compound-specific RPFs remained similar over the different endpoints studied downstream the AOP. Therefore, it might be possible to reduce testing to a smaller battery of key tests. The results demonstrate the suitability of our in vitro assay toolbox, integrated within an AOP framework and combined with the RPF approach, for the analysis of steatotic effects of chemical mixtures. However, mRNA results suggest that the steatosis AOP still needs improvement.

Current Insights into Monitoring, Bioaccumulation, and Potential Health Effects of Microplastics Present in the Food Chain.

Microplastics (MPs) are considered an emerging issue as environmental pollutants and a potential health threat. This review will focus on recently published data on concentrations in food, possible effects, and monitoring methods. Some data are available on concentrations in seafood (fish, bivalves, and shrimps), water, sugar, salt, and honey, but are lacking for other foods. Bottled water is a considerable source with numbers varying between 2600 and 6300 MPs per liter. Particle size distributions have revealed an abundance of particles smaller than 25 µm, which are considered to have the highest probability to pass the intestinal border and to enter the systemic circulation of mammals. Some studies with mice and zebrafish with short- or medium-term exposure (up to 42 days) have revealed diverse results with respect to both the type and extent of effects. Most notable modifications have been observed in gut microbiota, lipid metabolism, and oxidative stress. The principal elements of MP monitoring in food are sample preparation, detection, and identification. Identified data gaps include a lack of occurrence data in plant- and animal-derived food, a need for more data on possible effects of different types of microplastics, a lack of in silico models, a lack of harmonized monitoring methods, and a further development of quality assurance.

Pyrrolizidine alkaloids in food and phytomedicine: Occurrence, exposure, toxicity, mechanisms, and risk assessment - A review.

Among naturally occurring plant constituents, the 1,2-unsaturated pyrrolizidine alkaloids (in the following termed 'PAs') play a distinct role because of the large number of congeners occurring in nature and the pronounced toxicity of some congeners. Several PAs are hepatotoxic in humans, experimental and farm animals and were shown to be potent hepatocarcinogens in laboratory rodents. Although the general mode of action leading to toxicity has been elucidated, i.e., being mediated by metabolic conversion of the parent molecule into a highly reactive electrophile capable of attacking cellular target molecules, major questions related to the risk assessment of PAs remain unresolved. It was the aim of a workshop held in September 2018 to shed more light on the occurrence, exposure, mode of action, toxicokinetics and -dynamics of PAs to improve the scientific basis for an advanced toxicological risk assessment. The contributions in nine chapters describe the scientific progress using advanced analytical methods, studies in subcellular fractions, cell culture, experimental animals and humans and the use of PBPK modeling and structure-activity relationship considerations aiming at a better understanding of PA toxicity and genotoxicity. Since PAs differ considerably in their toxic potencies and substantial species differences in sensitivity towards PA exposure exist, a special emphasis was placed on these issues.

Towards harmonization of test methods for in vitro hepatic clearance studies.

Non-animal methods for toxicokinetics, such as in vitro hepatic metabolic clearance studies, play an important role in chemical risk evaluations. To gain regulatory acceptance of such clearance data, the development of a test guideline for performing in vitro hepatic clearance studies is crucial. The aim of the present study was to obtain insight in the experimental conditions of clearance studies that influence obtained intrinsic clearance (CLint) values. To that end, in vitro hepatic CLint data obtained with rat or human hepatocytes and methodological aspects of the experiments, were collected from 42 different suitable studies published between 1995 and 2018. The CLint values for the majority of chemicals differed by more than one order of magnitude. We estimated the systematic effect of different experimental setups on the CLint values using a random forest regression analysis, revealing that 'hepatocyte concentration', 'species' (rat or human hepatocytes) and 'culture medium' have the largest impact. Calculating unbound CLint (CLint,u) values slightly reduced the variation for most chemicals. Given that in vivo clearance is in general underpredicted based on in vitro clearance data, a harmonized protocol is preferably based on a protocol that provides relatively high in vitro CLint values.