Non-alcoholic fatty liver disease (NAFLD) has emerged as a public health concern as reflected in its widespread distribution in the general population. Yet, treatment options are scarce which is at least in part due to lack of reliable human in vitro disease models. Here, we report a human hepatic 3D spheroid system cultured under defined chemical conditions that has the potential to mimic steatotic conditions in a reversible manner, useful for identification of novel drug treatment conditions. Primary human hepatocytes (PHH) from different donors were cultured as spheroid microtissues in physiological in vivo -like culture conditions. Hepatic steatosis was induced over the course of three weeks in culture by supplementing the culture medium with pathophysiological concentrations of free fatty acids, carbohydrates and insulin. Effects of steatosis in the 3D system were evaluated on transcriptional, metabolomic and lipidomic levels. Free fatty acids on one hand as well as a combination of insulin and monosaccharides, promoted lipid accumulation in hepatocytes and increased expression of lipogenic genes, such as fatty acid synthase. This milieu also promoted development of insulin resistance within 2 weeks as manifested by an increase in gluconeogenic and insulin resistance markers, which are observed in type 2 diabetes mellitus and metabolic syndrome. Induced steatosis was reversible after withdrawal of lipogenic substrates and a further reduction in cellular fat content was observed following treatment with different antisteatotic compounds, such as metformin, glucagon, olaparib and antioxidants. Taken together, these results demonstrate that the 3D hepatic spheroids can serve as a valuable, HTS compatible model for the study of liver steatosis and facilitate translational discovery of novel drug targets.

Drug-induced cholestasis (DIC) is poorly understood and its preclinical prediction is mainly limited to assessing the compound's potential to inhibit the bile salt export pump (BSEP). Here, we evaluated two 3D spheroid models, one from primary human hepatocytes (PHH) and one from HepaRG cells, for the detection of compounds with cholestatic liability. By repeatedly co-exposing both models to a set of compounds with different mechanisms of hepatotoxicity and a non-toxic concentrated bile acid (BA) mixture for 8 days we observed a selective synergistic toxicity of compounds known to cause cholestatic or mixed cholestatic/hepatocellular toxicity and the BA mixture compared to exposure to the compounds alone, a phenomenon that was more pronounced after extending the exposure time to 14 days. In contrast, no such synergism was observed after both 8 and 14 days of exposure to the BA mixture for compounds that cause non-cholestatic hepatotoxicity. Mechanisms behind the toxicity of the cholestatic compound chlorpromazine were accurately detected in both spheroid models, including intracellular BA accumulation, inhibition of ABCB11 expression and disruption of the F-actin cytoskeleton. Furthermore, the observed synergistic toxicity of chlorpromazine and BA was associated with increased oxidative stress and modulation of death receptor signalling. Combined, our results demonstrate that the hepatic spheroid models presented here can be used to detect and study compounds with cholestatic liability.

Non-alcoholic steatohepatitis (NASH) is a major health problem leading to liver fibrosis and hepatocellular carcinoma, among other diseases, and for which there is still no approved drug treatment. Previous studies in animal models and in LX-2 cells have indicated a role for serotonin (5-HT) and 5-HT receptors in stellate cell activation and the development of NASH. In the current study, we investigated the extent to which these findings are applicable to a human NASH in vitro model consisting of human liver spheroids containing hepatocytes and non-parenchymal cells. Treatment of the spheroids with 5-HT or free fatty acids (FFA) induced fibrosis, whereas treatment of the spheroids with the 5-HT receptor antagonists ketanserin, pimavanserin, sarpogrelate, and SB269970 inhibited FFA-induced fibrosis via a reduction in stellate cell activation as determined by the expression of vimentin, TGF-β1 and COL1A1 production. siRNA-based silencing of 5-HT2A receptor expression reduced the anti-fibrotic properties of ketanserin, suggesting a role for 5-HT receptors in general and 5-HT2A receptors in particular in the FFA-mediated increase in fibrosis in the human liver spheroid model. The results suggest a contribution of the 5-HT receptors in the development of FFA-induced human liver fibrosis with implications for further efforts in drug development.

Liver biology and function, drug-induced liver injury (DILI) and liver diseases are difficult to study using current in vitro models such as primary human hepatocyte (PHH) monolayer cultures, as their rapid de-differentiation restricts their usefulness substantially. Thus, we have developed and extensively characterized an easily scalable 3D PHH spheroid system in chemically-defined, serum-free conditions. Using whole proteome analyses, we found that PHH spheroids cultured this way were similar to the liver in vivo and even retained their inter-individual variability. Furthermore, PHH spheroids remained phenotypically stable and retained morphology, viability, and hepatocyte-specific functions for culture periods of at least 5 weeks. We show that under chronic exposure, the sensitivity of the hepatocytes drastically increased and toxicity of a set of hepatotoxins was detected at clinically relevant concentrations. An interesting example was the chronic toxicity of fialuridine for which hepatotoxicity was mimicked after repeated-dosing in the PHH spheroid model, not possible to detect using previous in vitro systems. Additionally, we provide proof-of-principle that PHH spheroids can reflect liver pathologies such as cholestasis, steatosis and viral hepatitis. Combined, our results demonstrate that the PHH spheroid system presented here constitutes a versatile and promising in vitro system to study liver function, liver diseases, drug targets and long-term DILI.

The transcription factor Nrf2 coordinates an adaptive response to chemical and oxidative stress characterised by the upregulated expression of cytoprotective target genes. In order to understand the mechanistic relevance of Nrf2 as a marker of drug-induced stress it is important to know if this adaptive response is truly localised in the context of organ-specific drug toxicity. Here, we address this knowledge gap through real-time bioluminescence imaging of transgenic Nrf2-luciferase (Nrf2-luc) reporter mice following administration of the metabolism-dependent hepatotoxin acetaminophen (APAP) or the direct nephrotoxin cisplatin. We detected localised bioluminescence in the liver (APAP) and kidneys (cisplatin) in vivo and ex vivo, whilst qPCR, Taqman low-density array and immunoblot analysis of these tissues further revealed increases in the expression level of several endogenous Nrf2-regulated genes/proteins, including heme oxygenase 1 (Hmox1). Consistent with the toxic effects of APAP in the liver and cisplatin in the kidney, immunohistochemical analysis revealed the elevated expression of luciferase and Hmox1 in centrilobular hepatocytes and in tubular epithelial cells, respectively. In keeping with the role of reactive metabolite formation in APAP-induced chemical stress, both the hepatotoxicity and localised Nrf2-luc response were ameliorated by the cytochrome P450 inhibitor aminobenzotriazole. Together, these findings show that Nrf2 can reflect highly-localised cellular perturbations associated with relevant toxicological mechanisms.

There is a high unmet need for developing treatments for nonalcoholic fatty liver disease (NAFLD), for which there are no approved drugs today. Here, we used a human in vitro disease model to understand mechanisms linked to genetic risk variants associated with NAFLD. The model is based on 3D spheroids from primary human hepatocytes from five different donors. Across these donors, we observed highly reproducible differences in the extent of steatosis induction, demonstrating that inter-donor variability is reflected in the in vitro model. Importantly, our data indicates that the genetic variant TM6SF2 E167K, previously associated with increased risk for NAFLD, induces increased hepatocyte fat content by reducing APOB particle secretion. Finally, differences in gene expression pathways involved in cholesterol, fatty acid and glucose metabolism between wild type and TM6SF2 E167K mutation carriers (N = 125) were confirmed in the in vitro model. Our data suggest that the 3D in vitro spheroids can be used to investigate the mechanisms underlying the association of human genetic variants associated with NAFLD. This model may also be suitable to discover new treatments against NAFLD.

Maximising the use of preclinical murine models of progressive kidney disease as test beds for therapies ideally requires kidney function to be measured repeatedly in a safe, minimally invasive manner. To date, most studies of murine nephropathy depend on unreliable markers of renal physiological function, exemplified by measuring blood levels of creatinine and urea, and on various end points necessitating sacrifice of experimental animals to assess histological damage, thus counteracting the principles of Replacement, Refinement and Reduction. Here, we applied two novel minimally invasive techniques to measure kidney function in SCID mice with adriamycin-induced nephropathy. We employed i) a transcutaneous device that measures the half-life of intravenously administered FITC-sinistrin, a molecule cleared by glomerular filtration; and ii) multispectral optoacoustic tomography, a photoacoustic imaging device that directly visualises the clearance of the near infrared dye, IRDye 800CW carboxylate. Measurements with either technique showed a significant impairment of renal function in experimental animals versus controls, with significant correlations with the proportion of scarred glomeruli five weeks after induction of injury. These technologies provide clinically relevant functional data and should be widely adopted for testing the efficacies of novel therapies. Moreover, their use will also lead to a reduction in experimental animal numbers.

Chronic lymphocytic leukaemia (CLL) exhibits variable clinical course and response to therapy, but the molecular basis of this variability remains incompletely understood. Data independent acquisition (DIA)-MS technologies, such as SWATH (Sequential Windowed Acquisition of all THeoretical fragments), provide an opportunity to study the pathophysiology of CLL at the proteome level. Here, a CLL-specific spectral library (7736 proteins) is described alongside an analysis of sample replication and data handling requirements for quantitative SWATH-MS analysis of clinical samples. The analysis was performed on 6 CLL samples, incorporating biological (IGHV mutational status), sample preparation and MS technical replicates. Quantitative information was obtained for 5169 proteins across 54 SWATH-MS acquisitions: the sources of variation and different computational approaches for batch correction were assessed. Functional enrichment analysis of proteins associated with IGHV mutational status showed significant overlap with previous studies based on gene expression profiling. Finally, an approach to perform statistical power analysis in proteomics studies was implemented. This study provides a valuable resource for researchers working on the proteomics of CLL. It also establishes a sound framework for the design of sufficiently powered clinical proteomics studies. Indeed, this study shows that it is possible to derive biologically plausible hypotheses from a relatively small dataset.