1. Introduction
Drug-induced liver injury (DILI) is the most common cause of liver injury affecting approximately 14–19 cases per 100,000 patient-years in Western patients. Among the patients diagnosed with DILI in a French study, 6% died and 12% required hospitalization [
1]. DILI is also the most common reason for failures during the drug development process [
2] and for safety-related withdrawal of drugs from the pharmaceutical market [
3]. Therefore, having tools and techniques that can detect hepatotoxic properties in drug candidates at an early stage of discovery is highly desirable. Traditionally, this has been accomplished by performing in vivo studies in mice and rats by measuring the serum alanine aminotransferase levels and by histological analysis of liver sections. However, these models are time demanding, costly, and in addition interfere with the welfare of the animals. Hence there is an urgent need to identify fast and cost-efficient methods using alternative models with high human predictability that in the preclinical phase identify and exclude drugs with hepatotoxic properties.
Zebrafish (
Danio rerio) are a promising small vertebrate animal model for use in biomedical research and represent a powerful alternative model to implement the 3R principle of humane animal research (replacement, reduction, and refinement) [
4,
5,
6]. In the context of compound hepatotoxicity testing the model offers several advantages. First and foremost, zebrafish show a high similarity to mammals in hepatocyte metabolism and liver functionality [
7] and display similar responses to hepatotoxicants [
4]. In addition, the zebrafish liver is fully functional at 72 hpf [
5,
8], so hepatic responses after exposure of zebrafish larvae to hepatotoxicants can be monitored from an early stage onwards. Secondly, zebrafish are highly amenable to medium/high throughput screening. One adult pair of fish can produce approximately 200–300 eggs/week, resulting in access to large numbers of eggs/embryos. Moreover, the developing zebrafish larvae have a size of 1–4 mm, thus experiments can easily be performed in microtiter plates using only microgram amounts of compound.
Current hepatotoxicity testing using zebrafish includes endpoints such as liver degeneration and yolk retention [
6], expression of liver-specific fatty acid binding protein (
lfabp10a) in hepatocytes, and liver morphology using a Tg(
fabp10a:DsRed) zebrafish transgenic line that specifically expresses fluorescent DsRed in hepatocytes [
9]. Unfortunately, all these methods require blinding of the investigators to minimize bias and maximize the validity of the results. Besides, the 2D image-based results do not necessarily reflect the actual 3D changes in the liver.
As hepatocytes account for 70%–75% of hepatic cells and carry out most of the liver functions [
10], it is anticipated that the reduction in number of hepatocytes after drug exposure is a key indicator of hepatotoxicity. Presently, fluorescence-activated cell sorting (FACS) is the only method used to estimate ex vivo the number of hepatocytes present in zebrafish [
11]. However, as the procedure requires single sample filtering and processing, the method is inherently slow and hence not suitable for medium/high throughput purposes.
In this study cell imaging counting (CIC) was used to quantify the amount of DsRed-labeled hepatocytes recovered from Tg(
fabp10a:DsRed) transgenic zebrafish in a fast and straightforward way. Two compounds known to cause hepatotoxicity in humans, i.e., paracetamol and tetracycline [
6], as well as two non-hepatotoxicants (sucrose and rutin) were selected for the study. The outcome was compared with results found with commonly used hepatotoxicity-testing methods that are based on assessing liver morphology.
The studies published to date only used immersion as the route of exposing zebrafish to drugs investigating their hepatotoxic effects. This situation can result in false-negatives due to poor absorption as observed when investigating cardiotoxicity [
12]. We therefore also explored the applicability of using yolk and pericardial injection as alternative administration routes for testing hepatoxicants. Finally, we also studied the pharmacokinetics of paracetamol and tetracycline in zebrafish larvae using the concerned administration routes.
3. Discussion
The first objective of this study was to evaluate cell imaging counting (CIC) as a new and efficient approach for ex vivo hepatotoxicity assessment in zebrafish larvae. CIC is a high-content imaging and quantification tool that is used for detecting cells and even changes in cells at a subcellular level [
13,
14]. Using a transgenic zebrafish line that specifically expresses fluorescent DsRed in the hepatocytes, we were able to quantify the amount of recovered hepatocytes in a reproducible way, even after a simple protease-digestion of the larvae in combination with a short centrifugation step. Indeed, in contrast to FACS analysis of cells recovered from whole organisms that needs sample filtering in order to prevent lines from clogging [
11], our CIC-based method allows the straightforward detection of a limited number of fluorescent cells in the midst of non-fluorescent debris, without the need for time-consuming purification steps. This implies that the method can be used to test larger collections of compounds in a medium-throughput way.
To evaluate the sensitivity of CIC to detect drug-induced hepatotoxicity, visual analysis of liver morphology by fluorescence microscopy and size analysis of fluorescent 2D liver images were employed as reference methods. The data show that all three methods were able to identify the hepatotoxic effects of paracetamol and tetracycline, at least when larvae were immersed in the non-lethal maximal concentration (MC) (
Table 4). As compared to the other methods, CIC was more sensitive in detecting the effects induced by tetracycline at concentrations lower than MC. CIC also identified hepatotoxicity induced by paracetamol at MD via both yolk injection and pericardial injection, whereas the other methods were only successful for yolk injections. However, liver morphology assessment turned out to be more sensitive than CIC, e.g., in case of pericardial or yolk injections of tetracycline (
Table 4). In general, this study also suggests that by combining visual assessment of liver morphology and CIC, a proper strategy is obtained for a very successful screening of the hepatotoxicity of compounds. As a matter of fact, by combining these methods, both paracetamol and tetracycline are recognized as hepatotoxic compounds via any of the three administration routes.
The reason underlying the differences in sensitivity of the three methodological approaches is presently not clear. Obviously, further studies including detailed microscopic analyses of liver sections are needed that would allow us to correlate the outcome, as observed in this study, with histological changes and cellular viability.
Herein we also explored the applicability of using injections for the delivery of compounds investigated for their hepatotoxic effects. In case of paracetamol a clear effect was detected after yolk injection at MD, as identified by all three methods. Thus hepatotoxicity can also be induced by a single overdose of paracetamol in zebrafish larvae, as observed for humans [
15]. Tetracycline induced pronounced effects, both after yolk and pericardial injections, although the hepatotoxicity was detected only by visual analysis of liver morphology, and not by the other methods (
Table 4).
Overall, the outcome therefore shows that also injections, especially in case of administering compounds into the yolk, can produce significant results in the framework of hepatotoxicity testing. However, as the pharmacokinetic data indicate that both paracetamol and tetracycline are absorbed by immersion, leading to obvious liver effects, it is not clear from the present data that yolk injections could be advantageous over simple bathing of larvae in compound solutions. Obviously, the yolk administration route that can also be performed automatically [
16] is particularly of interest when testing hepatotoxicity of compounds that are poorly absorbed, i.e., in case of highly hydrophilic or large compounds [
12].
The pharmacokinetic profile of paracetamol shows that almost no compound could be recovered from the larvae 24 h after yolk injection, although the hepatotoxic effects detected at 6 dpf are similar to the immersion results. This seems to indicate that the injurious effects take place particularly during the first 24 h after administration (3–4 dpf). Alternatively, it is possible that metabolites like
N-acetyl-p-benzoquinone imine (NAPQI) are formed (not tracked down in this study) that exert their effect over time. Of interest, NAPQI-GSH adduct (GSH is the reduced form of glutathion) has been found when paracetamol was incubated with liver microsomes of adult zebrafish, possibly as the result of the expression of cytochromes P450 (CYPs), CYP3A65 with contributions of CYP3C1 [
17]. Of interest, the genes for these CYPs are expressed in embryonal zebrafish as well [
18], although it is not known whether the protein concentrations reach functional levels.
Interestingly, the pharmacokinetic profiles after yolk or pericardial injection of paracetamol are somewhat similar, although the hepatotoxic effects induced were different. Likely, both administration routes result in a dissimilar overall distribution of the compound in organs and tissue, affecting the final hepatotoxic effect.
In case of injected tetracycline, different pharmacokinetic profiles were observed as compared to the ones seen after injection of paracetamol. The concentrations remained somewhat constant over the first 24 h after administration, and it looks like the compound is not readily excreted from the larval body. It is possible that tetracycline binds to a large extent and with high affinity to local tissue without being distributed over the rest of the larvae. Anyhow, the final hepatotoxicity induced by the different administration routes as detected by liver morphology analysis is surprisingly quite similar. In contrast, when using CIC to assess hepatotoxicity, immersion led to pronounced effects (all concentrations active), whereas no impact was observed after injection. Again, the overall outcome regarding tetracycline is likely the consequence of a complex interplay between distribution, metabolism and detection of hepatotoxicity. These first intriguing results regarding pharmacokinetics and organ toxicity in zebrafish larvae definitely demand further investigation, and it looks like detailed metabolism data are required to fully understand the results. However, this experimental work was beyond the scope of the present study.
4. Materials and Methods
4.1. Chemicals
Paracetamol, tetracycline, sucrose and rutin were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Protease inhibitor cocktail (1.5× solution) was made by dissolving 1 tablet (cOmplete mini Protease Inhibitor Cocktail Tablets, Roche, Basel, Switzerland) in 6.5 mL Dulbecco’s phosphate-buffered saline (DPBS) (Life Technologies, Ledeberg, Belgium) without calcium and magnesium.
Dispase II (Sigma-Aldrich) was dissolved in DPBS (Life Technologies, Ledeberg, Belgium) without calcium and magnesium to make a 20 mg/mL stock solution.
UHPLC experiments were performed using ammonium formate purchased from Fluka (Buchs, Switzerland). Formic acid was obtained from Acros Organics (Geel, Belgium), acetonitrile and methanol was purchased from Fischer Scientific (Loughborough, UK). All reagents were of analytical grade or better. Ultrapure water was produced using a Millipore Milli-Q Gradient system (Milford, MA, USA). All pH modifications were done with a Metrohm 691 pH meter (Antwerp, Belgium).
4.2. Zebrafish
All zebrafish experiments carried out were approved by the Ethics Committee of the University of Leuven (approval number P154/2015, the start date: 1 October 2015, the end date: 30 September 2019) and by the Belgian Federal Department of Public Health, Food Safety & Environment (approval number LA1210261).
A transgenic zebrafish line expressing DsRed in hepatocytes under the
fabp10a promoter (herein referred to as Tg(
fabp10a:DsRed)) was created and raised in our laboratory [
9]. Tg(
fabp10a:DsRed) zebrafish (larvae and adults) were kept in standard conditions for zebrafish (28 °C, 14/10-h light/dark cycle). Tg(
fabp10a:DsRed) adults from the third and fourth generation were mated (two males and three females per cage). Only embryos of good quality (fertilized, clear cytoplasm and symmetric cleavage) were selected for experiments and kept in petri dishes containing Danieau’s solution until compound exposure. All larvae were derived from the same spawns of eggs for statistical comparison between the control and treated groups.
4.3. Compound Administration
Compounds tested were dissolved in DMSO (stock solution) and then further diluted in Danieau’s solution with a final concentration of 1% DMSO or diluted with 0.9% NaCl solution (1:1 ratio) to obtain a working solution for immersion and yolk/pericardial injections, respectively. At 3 dpf, hatched Tg(
fabp10a:DsRed) larvae were transferred to 24-well plates (Falcon
®, Corning, Lasne, Belgium) with each well containing 10 larvae in 1 mL Danieau’s solution and the dissolved drug or the vehicle control for immersion. For yolk injections, and pericardial injections into the bloodstream via the common cardinal veins [
19,
20], 12 larvae were used for each condition, and each larva was injected with 1 nL of solution at 3 dpf.
4.4. Hepatotoxicity Assessed by Visual Analysis of Liver Morphology Using Fluorescence Microscopy and Size Analysis of Fluorescent 2D Liver Images
At 6 dpf treated larvae were anesthetized by 0.5 mM tricaine in a 6 well-plate (Falcon®) and arranged so that they were lying on their right sides with their heads facing left. This enabled us to view the left lateral side of the larvae and the larger lobe of the liver that was scored in a blinded manner under a fluorescence microscope (Leica MZ10F-Dsred channel, 40× magnification, Leica Microsystems, Wetzlar, Germany). Livers were considered abnormal if their typical globular structure underwent gross morphological changes. Images of the livers (including a scale bar) were imported to ImageJ software (National Institutes of Health, Bethesda, MD, USA) for determining the liver size. To that end, the picture was set to scale (functions: analyze, set scale), the liver area extracted (functions: split channel, threshold) and the area calculated (function: analyze particles, settings: size = 0−infinity, circularity = 0.00−1.00).
4.5. Quantification of Hepatocytes Using Cell Imaging Counting (CIC)
4.5.1. Larvae Digestion and Recovery of Hepatocytes
In order to digest the larvae and separate hepatocytes from the liver, larvae were transferred to 96 well-plates (Falcon®, Corning, Lasne, Belgium) with one larva per well and incubated with 200 µL 2.5 mg dispase/mL DPBS for 4 h. The proteolytic activity was stopped by adding 100 µL of 1.5× protease inhibitors cocktail solution. The hepatocytes were dispersed by pipetting up and down 20 times with a 1 mL automatic pipette. After digestion, hepatocytes were positioned at the bottom of the well by centrifugation at 200× g (Jouan CR412 centrifuge, Jouan, Saint-Herblain, France) for 15 min.
4.5.2. Cell Imaging Counting (CIC) of Hepatocytes
Fluorescent hepatocytes present in the wells were imaged and counted on an IN Cell Analyzer 2000 using software IN Cell Investigator Developer Toolbox (GE Healthcare, Machelen, Belgium), with a 2×/0.1 NA lens, DsRed fluorescence channel, a 0.1-s exposure time, and the intensity segmentation (172-4095) set to ensure hepatocytes were recognized selectively.
4.6. Pharmacokinetics Study of Hepatotoxicants
4.6.1. Sample Preparation
At given time points (0, 2, 6, 12 and 24 h), treated larvae were collected and transferred to 1.5 mL Eppendorf
® tubes pre-filled with acid washed glass beads (diameter: 710–1180 µm, Sigma-Aldrich) and 270 µL of extraction medium (1:2 water:methanol, each containing 0.1% formic acid) was added (10 larvae per Eppendorf
® tube). Samples were centrifuged for 15 s at 14,100×
g (Eppendorf Mini Spin Plus, Hamburg, Germany) to insure high homogenization efficiency. Subsequently, the samples were homogenized with ultrasonication (Diagenode Bioruptor Plus, Seraing, Belgium) at +4 °C and high energy input setting to ensure the highest possible disruption of the cells. The overall treatment time was 15 min spread over 30 cycles of 30 s treatment followed by pauses of 30 s in-between. Six samples were homogenized during each treatment cycle. Subsequently, treated samples were centrifuged at 14,100×
g for 15 min and 200 µL of supernatant was transferred to an empty Eppendorf tube (1.5 mL) and stored on ice until all samples were processed. All samples were then placed into a vacuum oven (Model 1410, Sheldon Manufacturing Inc., Cornelius, OR, USA) and evaporated until dry at +45 °C during 45 min. The samples were then reconstituted with a mobile phase-like solution (6% acetonitrile in water), vortexed for 5 s, resuspended at least three times with a micropipette and sonicated (Branson 3510, Danbury, CT, USA) for 15 min. The reconstitution volume was chosen such that the linearity range of the method was not exceeded and is given in
Table 5. Finally, samples were transferred into amber glass vials (Supelco, Bellefonte, PA, USA) with 250 µL glass inserts (Supelco) for further analysis on UHPLC-UV.
4.6.2. Instrumental Analysis
An Agilent 1290 UHPLC system (Agilent, Waldbronn, Germany) consisting of an autosampler, quaternary pump, degasser and DAD-detector (diode array detector) was used for all uptake measurements. Data acquisition and peak processing was performed using OpenLAB CDS Chemstation Edition 01.04 software (Agilent). All calculations concerning the evaluation of recorded data were made in MS Excel (Microsoft Corporation, Seattle, WA, USA).
Compounds of interest were separated chromatographically on an Acquity BEH C18 column (100 mm × 2.1 mm, dp (particle diameter) = 1.7 µm) from Waters Inc. (Milford, MA, USA) at a flow rate of 0.4 mL/min. Gradient elution was performed starting at 97:3 (v/v) 10 mM ammonium formate (pH 2.8 modified with formic acid):acetonitrile, and changed to 18:82 (v/v) 10 mM ammonium formate (pH 2.8):acetonitrile in 10.5 min. The injection volume (Vinj) of each sample was 2 µL. After elution of the compounds, a cleaning step at 100% acetonitrile was applied for 1.3 min after which the column was re-equilibrated for 7.0 min at initial conditions. Paracetamol was detected at a wavelength of 254 nm, while a wavelength of 274 nm was used for the measurement of tetracycline.
The larval amount of hepatotoxicant (mg/kg) was calculated assuming the weight of a 3 or 4 dpf Tg(fabp10a:DsRed) zebrafish larva to be 0.25 mg.
4.7. Statistical Analysis
Differences between drug-treated groups and control vehicle control (VHC) regarding numbers of hepatocytes and liver sizes were analyzed for statistical significance using an ANOVA test. Regarding the liver morphology, Fisher’s exact test was used for the analysis. p < 0.05 was considered statistically significant.
5. Conclusions
We successfully applied CIC ex vivo to detect hepatotoxicity induced by paracetamol and tetracycline in our transgenic zebrafish model. The outcome was in general comparable with the results obtained using two reference methods, i.e., visual analysis of liver morphology by fluorescence microscopy and size analysis of fluorescent 2D liver images, although clear differences exist depending on the administration routes investigated, i.e., immersion, yolk and pericardial injections. In general, immersion of larvae resulted in more consistent results among the different hepatotoxicity tests.
Unlike the other methods, our CIC-dependent approach measures directly and in an unbiased way the effect of compounds on the number of hepatocytes in the liver, is quick and reliable, and therefore provides a novel and straightforward tool for screening hepatotoxicants in early stages of drug development. As CIC allows for high-content imaging, the system is also amenable to more detailed analysis of recovered cells, e.g., multiparametric cytotoxicity testing. In addition, the method is also suitable for testing of other organ-related toxicities subject to the organs and tissue expressing fluorescent proteins in transgenic zebrafish lines.