Inhibition of the renal apical sodium dependent bile acid transporter prevents cholemic nephropathy in mice with obstructive cholestasis.

BACKGROUND & AIMS: Cholemic nephropathy (CN) is a severe complication of cholestatic liver diseases for which there is no specific treatment. We revisited its pathophysiology with the aim of identifying novel therapeutic strategies. METHODS: Cholestasis was induced by bile duct ligation (BDL) in mice. Bile flux in kidneys and livers was visualized by intravital imaging, supported by MALDI mass spectrometry imaging and liquid chromatography-tandem mass spectrometry. The effect of AS0369, a systemically bioavailable apical sodium-dependent bile acid transporter (ASBT) inhibitor, was evaluated by intravital imaging, RNA-sequencing, histological, blood, and urine analyses. Translational relevance was assessed in kidney biopsies from patients with CN, mice with a humanized bile acid (BA) spectrum, and via analysis of serum BAs and KIM-1 (kidney injury molecule 1) in patients with liver disease and hyperbilirubinemia. RESULTS: Proximal tubular epithelial cells (TECs) reabsorbed and enriched BAs, leading to oxidative stress and death of proximal TECs, casts in distal tubules and collecting ducts, peritubular capillary leakiness, and glomerular cysts. Renal ASBT inhibition by AS0369 blocked BA uptake into TECs and prevented kidney injury up to 6 weeks after BDL. Similar results were obtained in mice with humanized BA composition. In patients with advanced liver disease, serum BAs were the main determinant of KIM-1 levels. ASBT expression in TECs was preserved in biopsies from patients with CN, further highlighting the translational potential of targeting ASBT to treat CN. CONCLUSIONS: BA enrichment in proximal TECs followed by oxidative stress and cell death is a key early event in CN. Inhibiting renal ASBT and consequently BA enrichment in TECs prevents CN and systemically decreases BA concentrations. IMPACT AND IMPLICATIONS: Cholemic nephropathy (CN) is a severe complication of cholestasis and an unmet clinical need. We demonstrate that CN is triggered by the renal accumulation of bile acids (BAs) that are considerably increased in the systemic blood. Specifically, the proximal tubular epithelial cells of the kidney take up BAs via the apical sodium-dependent bile acid transporter (ASBT). We developed a therapeutic compound that blocks ASBT in the kidneys, prevents BA overload in tubular epithelial cells, and almost completely abolished all disease hallmarks in a CN mouse model. Renal ASBT inhibition represents a potential therapeutic strategy for patients with CN.

Acetaminophen overdose causes a breach of the blood-bile barrier in mice but not in rats.

Acetaminophen (APAP) is known to cause a breach of the blood-bile barrier in mice that, via a mechanism called futile bile acid (BA) cycling, increases BA concentrations in hepatocytes above cytotoxic thresholds. Here, we compared this mechanism in mice and rats, because both species differ massively in their susceptibility to APAP and compared the results to available human data. Dose and time-dependent APAP experiments were performed in male C57BL6/N mice and Wistar rats. The time course of BA concentrations in liver tissue and in blood was analyzed by MALDI-MSI and LC-MS/MS. APAP and its derivatives were measured in the blood by LC-MS. APAP-induced liver damage was analyzed by histopathology, immunohistochemistry, and by clinical chemistry. In mice, a transient increase of BA in blood and in peri-central hepatocytes preceded hepatocyte death. The BA increase coincided with oxidative stress in liver tissue and a compromised morphology of bile canaliculi and immunohistochemically visualized tight junction proteins. Rats showed a reduced metabolic activation of APAP compared to mice. However, even at very high doses that caused cell death of hepatocytes, no increase of BA concentrations was observed neither in liver tissue nor in the blood. Correspondingly, no oxidative stress was detectable, and the morphology of bile canaliculi and tight junction proteins remained unaltered. In conclusion, different mechanisms cause cell death in rats and mice, whereby oxidative stress and a breach of the blood-bile barrier are seen only in mice. Since transient cholestasis also occurs in human patients with APAP overdose, mice are a clinically relevant species to study APAP hepatotoxicity but not rats.

Integrated data from intravital imaging and HPLC-MS/MS analysis reveal large interspecies differences in AFB1 metabolism in mice and rats.

Large interspecies differences between rats and mice concerning the hepatotoxicity and carcinogenicity of aflatoxin B1 (AFB1) are known, with mice being more resistant. However, a comprehensive interspecies comparison including subcellular liver tissue compartments has not yet been performed. In this study, we performed spatio-temporal intravital analysis of AFB1 kinetics in the livers of anesthetized mice and rats. This was supported by time-dependent analysis of the parent compound as well as metabolites and adducts in blood, urine, and bile of both species by HPLC-MS/MS. The integrated data from intravital imaging and HPLC-MS/MS analysis revealed major interspecies differences between rats and mice: (1) AFB1-associated fluorescence persisted much longer in the nuclei of rat than mouse hepatocytes; (2) in the sinusoidal blood, AFB1-associated fluorescence was rapidly cleared in mice, while a time-dependent increase was observed in rats in the first three hours after injection followed by a plateau that lasted until the end of the observation period of six hours; (3) this coincided with a far stronger increase of AFB1-lysine adducts in the blood of rats compared to mice; (4) the AFB1-guanine adduct was detected at much higher concentrations in bile and urine of rats than mice. In both species, the AFB1-glutathione conjugate was efficiently excreted via bile, where it reached concentrations at least three orders of magnitude higher compared to blood. In conclusion, major differences between mice and rats were observed, concerning the nuclear persistence, formation of AFB1-lysine adducts, and the AFB1-guanine adducts.

Guided interactive image segmentation using machine learning and color-based image set clustering.

MOTIVATION: Over the last decades, image processing and analysis have become one of the key technologies in systems biology and medicine. The quantification of anatomical structures and dynamic processes in living systems is essential for understanding the complex underlying mechanisms and allows, i.e. the construction of spatio-temporal models that illuminate the interplay between architecture and function. Recently, deep learning significantly improved the performance of traditional image analysis in cases where imaging techniques provide large amounts of data. However, if only a few images are available or qualified annotations are expensive to produce, the applicability of deep learning is still limited. RESULTS: We present a novel approach that combines machine learning-based interactive image segmentation using supervoxels with a clustering method for the automated identification of similarly colored images in large image sets which enables a guided reuse of interactively trained classifiers. Our approach solves the problem of deteriorated segmentation and quantification accuracy when reusing trained classifiers which is due to significant color variability prevalent and often unavoidable in biological and medical images. This increase in efficiency improves the suitability of interactive segmentation for larger image sets, enabling efficient quantification or the rapid generation of training data for deep learning with minimal effort. The presented methods are applicable for almost any image type and represent a useful tool for image analysis tasks in general. AVAILABILITY AND IMPLEMENTATION: The presented methods are implemented in our image processing software TiQuant which is freely available at tiquant.hoehme.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Interruption of bile acid uptake by hepatocytes after acetaminophen overdose ameliorates hepatotoxicity.

BACKGROUND & AIMS: Acetaminophen (APAP) overdose remains a frequent cause of acute liver failure, which is generally accompanied by increased levels of serum bile acids (BAs). However, the pathophysiological role of BAs remains elusive. Herein, we investigated the role of BAs in APAP-induced hepatotoxicity. METHODS: We performed intravital imaging to investigate BA transport in mice, quantified endogenous BA concentrations in the serum of mice and patients with APAP overdose, analyzed liver tissue and bile by mass spectrometry and MALDI-mass spectrometry imaging, assessed the integrity of the blood-bile barrier and the role of oxidative stress by immunostaining of tight junction proteins and intravital imaging of fluorescent markers, identified the intracellular cytotoxic concentrations of BAs, and performed interventions to block BA uptake from blood into hepatocytes. RESULTS: Prior to the onset of cell death, APAP overdose causes massive oxidative stress in the pericentral lobular zone, which coincided with a breach of the blood-bile barrier. Consequently, BAs leak from the bile canaliculi into the sinusoidal blood, which is then followed by their uptake into hepatocytes via the basolateral membrane, their secretion into canaliculi and repeated cycling. This, what we termed 'futile cycling' of BAs, led to increased intracellular BA concentrations that were high enough to cause hepatocyte death. Importantly, however, the interruption of BA re-uptake by pharmacological NTCP blockage using Myrcludex B and Oatp knockout strongly reduced APAP-induced hepatotoxicity. CONCLUSIONS: APAP overdose induces a breach of the blood-bile barrier which leads to futile BA cycling that causes hepatocyte death. Prevention of BA cycling may represent a therapeutic option after APAP intoxication. LAY SUMMARY: Only one drug, N-acetylcysteine, is approved for the treatment of acetaminophen overdose and it is only effective when given within ∼8 hours after ingestion. We identified a mechanism by which acetaminophen overdose causes an increase in bile acid concentrations (to above toxic thresholds) in hepatocytes. Blocking this mechanism prevented acetaminophen-induced hepatotoxicity in mice and evidence from patients suggests that this therapy may be effective for longer periods after ingestion compared to N-acetylcysteine.

Intravital Dynamic and Correlative Imaging of Mouse Livers Reveals Diffusion-Dominated Canalicular and Flow-Augmented Ductular Bile Flux.

BACKGROUND AND AIMS: Small-molecule flux in tissue microdomains is essential for organ function, but knowledge of this process is scant due to the lack of suitable methods. We developed two independent techniques that allow the quantification of advection (flow) and diffusion in individual bile canaliculi and in interlobular bile ducts of intact livers in living mice, namely fluorescence loss after photoactivation and intravital arbitrary region image correlation spectroscopy. APPROACH AND RESULTS: The results challenge the prevailing "mechano-osmotic" theory of canalicular bile flow. After active transport across hepatocyte membranes, bile acids are transported in the canaliculi primarily by diffusion. Only in the interlobular ducts is diffusion augmented by regulatable advection. Photoactivation of fluorescein bis-(5-carboxymethoxy-2-nitrobenzyl)-ether in entire lobules demonstrated the establishment of diffusive gradients in the bile canalicular network and the sink function of interlobular ducts. In contrast to the bile canalicular network, vectorial transport was detected and quantified in the mesh of interlobular bile ducts. CONCLUSIONS: The liver consists of a diffusion-dominated canalicular domain, where hepatocytes secrete small molecules and generate a concentration gradient and a flow-augmented ductular domain, where regulated water influx creates unidirectional advection that augments the diffusive flux.

Inhibition of cytochrome P450 enhances the nephro- and hepatotoxicity of ochratoxin A.

The mycotoxin ochratoxin A (OTA) is a contaminant in food that causes nephrotoxicity and to a minor degree hepatotoxicity. Recently, we observed that OTA induces liver damage preferentially to the cytochrome P450 (CYP)-expressing pericentral lobular zone, similar to hepatotoxic substances known to be metabolically toxified by CYP, such as acetaminophen or carbon tetrachloride. To investigate whether CYP influences OTA toxicity, we used a single dose of OTA (7.5 mg/kg; intravenous) with and without pre-treatment with the pan CYP-inhibitor 1-aminobenzotriazole (ABT) 2 h before OTA administration. Blood, urine, as well as liver and kidney tissue samples were collected 24 h after OTA administration for biochemical and histopathological analyses. Inhibition of CYPs by ABT strongly increased the nephro- and hepatotoxicity of OTA. The urinary kidney damage biomarkers kidney injury molecule-1 (KIM-1) and neutrophil gelatinase-associated lipocalin (NGAL) were increased > 126-fold and > 20-fold, respectively, in mice treated with ABT and OTA compared to those receiving OTA alone. The blood biomarkers of liver damage, alanine transaminase (ALT) and aspartate transaminase (AST) both increased > 21- and 30-fold, respectively, when OTA was administered to ABT pre-treated mice compared to the effect of OTA alone. Histological analysis of the liver revealed a pericentral lobular damage induced by OTA despite CYP-inhibition by ABT. Administration of ABT alone caused no hepato- or nephrotoxicity. Overall, the results presented are compatible with a scenario where CYPs mediate the detoxification of OTA, yet the mechanisms responsible for the pericental liver damage pattern still remain to be elucidated.

Colchicine overdose impairs the capacity of Kupffer cells to clear foreign particles and endotoxins.

Colchicine is an anti-inflammatory drug with a narrow therapeutic index. Its binding to tubulin prevents microtubule polymerization; however, little is known about how depolymerization of microtubules interferes with the phagocytosis function of Kupffer cells (KC). Here, we applied functional intravital imaging techniques to investigate the influence of microtubule disruption by colchicine on KC morphology, as well as its capacity to clear foreign particles and bacterial lipopolysaccharide (LPS) in anesthetized mice. Intravital imaging of KC in healthy mice showed the typical elongated morphology, localization at the luminal side of the sinusoidal endothelial cells, and moving cell protrusions. In contrast, at colchicine doses of 1 mg/kg and higher (intraperitoneal), KC appeared roundish with strongly reduced protrusions and motility. To study the functional consequences of these alterations, we analyzed the capacity of KC to phagocytose fluorescent nanospheres (100 nm-size) and LPS. After tail vein injection, the nanospheres formed aggregates of up to ~ 5 µm moving along the sinusoidal bloodstream. In controls, the nanosphere aggregates were rapidly captured by the Kupffer cell protrusions, followed by an internalization process that lasted up to 10 min. Similar capture events and internalization processes were observed after the administration of fluorescently labeled LPS. In contrast, capture and internalization of both nanospheres and LPS by KC were strongly reduced in colchicine-treated mice. Reduced phagocytosis of LPS was accompanied by aggravated production of inflammatory cytokines. Since 0.4 mg/kg colchicine in mice has been reported to be bio-equivalent to human therapeutic doses, the here-observed adverse effects on KC occurred at doses only slightly above those used clinically, and may be critical for patients with endotoxemia due to a leaky gut-blood barrier.

Hypoalbuminemia affects the spatio-temporal tissue distribution of ochratoxin A in liver and kidneys: consequences for organ toxicity.

Hypoalbuminemia (HA) is frequently observed in systemic inflammatory diseases and in liver disease. However, the influence of HA on the pharmacokinetics and toxicity of compounds with high plasma albumin binding remained insufficiently studied. The 'lack-of-delivery-concept' postulates that HA leads to less carrier mediated uptake of albumin bound substances into hepatocytes and to less glomerular filtration; in contrast, the 'concept-of-higher-free-fraction' argues that increased concentrations of non-albumin bound compounds facilitate hepatocellular uptake and enhance glomerular filtration. To address this question, we performed intravital imaging on livers and kidneys of anesthetized mice to quantify the spatio-temporal tissue distribution of the mycotoxin ochratoxin A (OTA) based on its auto-fluorescence in albumin knockout and wild-type mice. HA strongly enhanced the uptake of OTA from the sinusoidal blood into hepatocytes, followed by faster secretion into bile canaliculi. These toxicokinetic changes were associated with increased hepatotoxicity in heterozygous albumin knockout mice for which serum albumin was reduced to a similar extent as in patients with severe hypoalbuminemia. HA also led to a shorter half-life of OTA in renal capillaries, increased glomerular filtration, and to enhanced uptake of OTA into tubular epithelial cells. In conclusion, the results favor the 'concept-of-higher-free-fraction' in HA; accordingly, HA causes an increased tissue uptake of compounds with high albumin binding and increased organ toxicity. It should be studied if this concept can be generalized to all compounds with high plasma albumin binding that are substrates of hepatocyte and renal tubular epithelial cell carriers.

Subcellular spatio-temporal intravital kinetics of aflatoxin B1 and ochratoxin A in liver and kidney.

Local accumulation of xenobiotics in human and animal tissues may cause adverse effects. Large differences in their concentrations may exist between individual cell types, often due to the expression of specific uptake and export carriers. Here we established a two-photon microscopy-based technique for spatio-temporal detection of the distribution of mycotoxins in intact kidneys and livers of anesthetized mice with subcellular resolution. The mycotoxins ochratoxin A (OTA, 10 mg/kg b.w.) and aflatoxin B1 (AFB1, 1.5 mg/kg b.w.), which both show blue auto-fluorescence, were analyzed after intravenous bolus injections. Within seconds after administration, OTA was filtered by glomeruli, and enriched in distal tubular epithelial cells (dTEC). A striking feature of AFB1 toxicokinetics was its very rapid uptake from sinusoidal blood into hepatocytes (t1/2 ~ 4 min) and excretion into bile canaliculi. Interestingly, AFB1 was enriched in the nuclei of hepatocytes with zonal differences in clearance. In the cytoplasm of pericentral hepatocytes, the half-life (t1/2~ 63 min) was much longer compared to periportal hepatocytes of the same lobules (t1/2 ~ 9 min). In addition, nuclear AFB1 from periportal hepatocytes cleared faster compared to the pericentral region. These local differences in AFB1 clearance may be due to the pericentral expression of cytochrome P450 enzymes that activate AFB1 to protein- and DNA-binding metabolites. In conclusion, the present study shows that large spatio-temporal concentration differences exist within the same tissues and its analysis may provide valuable additional information to conventional toxicokinetic studies.

Quantitative analysis of hepatic macro- and microvascular alterations during cirrhogenesis in the rat.

Cirrhosis represents the end-stage of any persistent chronically active liver disease. It is characterized by the complete replacement of normal liver tissue by fibrosis, regenerative nodules, and complete fibrotic vascularized septa. The resulting angioarchitectural distortion contributes to an increasing intrahepatic vascular resistance, impeding liver perfusion and leading to portal hypertension. To date, knowledge on the dynamically evolving pathological changes of the hepatic vasculature during cirrhogenesis remains limited. More specifically, detailed anatomical data on the vascular adaptations during disease development is lacking. To address this need, we studied the 3D architecture of the hepatic vasculature during induction of cirrhogenesis in a rat model. Cirrhosis was chemically induced with thioacetamide (TAA). At predefined time points, the hepatic vasculature was fixed and visualized using a combination of vascular corrosion casting and deep tissue microscopy. Three-dimensional reconstruction and data-fitting enabled cirrhogenic features to extracted at multiple scales, portraying the impact of cirrhosis on the hepatic vasculature. At the macrolevel, we noticed that regenerative nodules severely compressed pliant venous vessels from 12 weeks of TAA intoxication onwards. Especially hepatic veins were highly affected by this compression, with collapsed vessel segments severely reducing perfusion capabilities. At the microlevel, we discovered zone-specific sinusoidal degeneration, with sinusoids located near the surface being more affected than those in the middle of a liver lobe. Our data shed light on and quantify the evolving angioarchitecture during cirrhogenesis. These findings may prove helpful for future targeted invasive interventions.

Model Prediction and Validation of an Order Mechanism Controlling the Spatiotemporal Phenotype of Early Hepatocellular Carcinoma.

Recently, hepatocyte-sinusoid alignment (HSA) has been identified as a mechanism that supports the coordination of hepatocytes during liver regeneration to reestablish a functional micro-architecture (Hoehme et al. in Proc Natl Acad Sci 107(23):10371-10376, 2010). HSA means that hepatocytes preferentially align along the closest micro-vessels. Here, we studied whether this mechanism is still active in early hepatocellular tumors. The same agent-based spatiotemporal model that previously correctly predicted HSA in liver regeneration was further developed to simulate scenarios in early tumor development, when individual initiated hepatocytes gain increased proliferation capacity. The model simulations were performed under conditions of realistic liver micro-architectures obtained from 3D reconstructions of confocal laser scanning micrographs. Interestingly, the established model predicted that initiated hepatocytes at first arrange in elongated patterns. Only when the tumor progresses to cell numbers of approximately 4000, does it adopt spherical structures. This prediction may have relevant consequences, since elongated tumors may reach critical structures faster, such as larger vessels, compared to a spherical tumor of similar cell number. Interestingly, this model prediction was confirmed by analysis of the spatial organization of initiated hepatocytes in a rat liver tumor initiation study using single doses of 250 mg/kg of the genotoxic carcinogen N-nitrosomorpholine (NNM). Indeed, small clusters of GST-P positive cells induced by NNM were elongated, almost columnar, while larger GDT-P positive foci of approximately the size of liver lobuli adopted spherical shapes. From simulations testing numerous possible mechanisms, only HSA could explain the experimentally observed initial deviation from spherical shape. The present study demonstrates that the architecture of small cell clusters of hepatocytes early after initiation is still controlled by physiological mechanisms. However, this coordinating influence is lost when the tumor grows to approximately 4000 cells, leading to further growth in spherical shape. Our findings stress the potential importance of organ micro-architecture in understanding tumor phenotypes.

Model-guided identification of a therapeutic strategy to reduce hyperammonemia in liver diseases.

BACKGROUND & AIMS: Recently, spatial-temporal/metabolic mathematical models have been established that allow the simulation of metabolic processes in tissues. We applied these models to decipher ammonia detoxification mechanisms in the liver. METHODS: An integrated metabolic-spatial-temporal model was used to generate hypotheses of ammonia metabolism. Predicted mechanisms were validated using time-resolved analyses of nitrogen metabolism, activity analyses, immunostaining and gene expression after induction of liver damage in mice. Moreover, blood from the portal vein, liver vein and mixed venous blood was analyzed in a time dependent manner. RESULTS: Modeling revealed an underestimation of ammonia consumption after liver damage when only the currently established mechanisms of ammonia detoxification were simulated. By iterative cycles of modeling and experiments, the reductive amidation of alpha-ketoglutarate (α-KG) via glutamate dehydrogenase (GDH) was identified as the lacking component. GDH is released from damaged hepatocytes into the blood where it consumes ammonia to generate glutamate, thereby providing systemic protection against hyperammonemia. This mechanism was exploited therapeutically in a mouse model of hyperammonemia by injecting GDH together with optimized doses of cofactors. Intravenous injection of GDH (720 U/kg), α-KG (280 mg/kg) and NADPH (180 mg/kg) reduced the elevated blood ammonia concentrations (>200 µM) to levels close to normal within only 15 min. CONCLUSION: If successfully translated to patients the GDH-based therapy might provide a less aggressive therapeutic alternative for patients with severe hyperammonemia.

TiQuant: software for tissue analysis, quantification and surface reconstruction.

MOTIVATION: TiQuant is a modular software tool for efficient quantification of biological tissues based on volume data obtained by biomedical image modalities. It includes a number of versatile image and volume processing chains tailored to the analysis of different tissue types which have been experimentally verified. TiQuant implements a novel method for the reconstruction of three-dimensional surfaces of biological systems, data that often cannot be obtained experimentally but which is of utmost importance for tissue modelling in systems biology. AVAILABILITY AND IMPLEMENTATION: TiQuant is freely available for non-commercial use at msysbio.com/tiquant. Windows, OSX and Linux are supported. CONTACT: hoehme@uni-leipzig.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integrated metabolic spatial-temporal model for the prediction of ammonia detoxification during liver damage and regeneration.

UNLABELLED: The impairment of hepatic metabolism due to liver injury has high systemic relevance. However, it is difficult to calculate the impairment of metabolic capacity from a specific pattern of liver damage with conventional techniques. We established an integrated metabolic spatial-temporal model (IM) using hepatic ammonia detoxification as a paradigm. First, a metabolic model (MM) based on mass balancing and mouse liver perfusion data was established to describe ammonia detoxification and its zonation. Next, the MM was combined with a spatial-temporal model simulating liver tissue damage and regeneration after CCl4 intoxication. The resulting IM simulated and visualized whether, where, and to what extent liver damage compromised ammonia detoxification. It allowed us to enter the extent and spatial patterns of liver damage and then calculate the outflow concentrations of ammonia, glutamine, and urea in the hepatic vein. The model was validated through comparisons with (1) published data for isolated, perfused livers with and without CCl4 intoxication and (2) a set of in vivo experiments. Using the experimentally determined portal concentrations of ammonia, the model adequately predicted metabolite concentrations over time in the hepatic vein during toxin-induced liver damage and regeneration in rodents. Further simulations, especially in combination with a simplified model of blood circulation with three ammonia-detoxifying compartments, indicated a yet unidentified process of ammonia consumption during liver regeneration and revealed unexpected concomitant changes in amino acid metabolism in the liver and at extrahepatic sites. CONCLUSION: The IM of hepatic ammonia detoxification considerably improves our understanding of the metabolic impact of liver disease and highlights the importance of integrated modeling approaches on the way toward virtual organisms.

How predictive quantitative modelling of tissue organisation can inform liver disease pathogenesis.

From the more than 100 liver diseases described, many of those with high incidence rates manifest themselves by histopathological changes, such as hepatitis, alcoholic liver disease, fatty liver disease, fibrosis, and, in its later stages, cirrhosis, hepatocellular carcinoma, primary biliary cirrhosis and other disorders. Studies of disease pathogeneses are largely based on integrating -omics data pooled from cells at different locations with spatial information from stained liver structures in animal models. Even though this has led to significant insights, the complexity of interactions as well as the involvement of processes at many different time and length scales constrains the possibility to condense disease processes in illustrations, schemes and tables. The combination of modern imaging modalities with image processing and analysis, and mathematical models opens up a promising new approach towards a quantitative understanding of pathologies and of disease processes. This strategy is discussed for two examples, ammonia metabolism after drug-induced acute liver damage, and the recovery of liver mass as well as architecture during the subsequent regeneration process. This interdisciplinary approach permits integration of biological mechanisms and models of processes contributing to disease progression at various scales into mathematical models. These can be used to perform in silico simulations to promote unravelling the relation between architecture and function as below illustrated for liver regeneration, and bridging from the in vitro situation and animal models to humans. In the near future novel mechanisms will usually not be directly elucidated by modelling. However, models will falsify hypotheses and guide towards the most informative experimental design.

Protocols for staining of bile canalicular and sinusoidal networks of human, mouse and pig livers, three-dimensional reconstruction and quantification of tissue microarchitecture by image processing and analysis.

Histological alterations often constitute a fingerprint of toxicity and diseases. The extent to which these alterations are cause or consequence of compromised organ function, and the underlying mechanisms involved is a matter of intensive research. In particular, liver disease is often associated with altered tissue microarchitecture, which in turn may compromise perfusion and functionality. Research in this field requires the development and orchestration of new techniques into standardized processing pipelines that can be used to reproducibly quantify tissue architecture. Major bottlenecks include the lack of robust staining, and adequate reconstruction and quantification techniques. To bridge this gap, we established protocols employing specific antibody combinations for immunostaining, confocal imaging, three-dimensional reconstruction of approximately 100-µm-thick tissue blocks and quantification of key architectural features. We describe a standard procedure termed 'liver architectural staining' for the simultaneous visualization of bile canaliculi, sinusoidal endothelial cells, glutamine synthetase (GS) for the identification of central veins, and DAPI as a nuclear marker. Additionally, we present a second standard procedure entitled 'S-phase staining', where S-phase-positive and S-phase-negative nuclei (stained with BrdU and DAPI, respectively), sinusoidal endothelial cells and GS are stained. The techniques include three-dimensional reconstruction of the sinusoidal and bile canalicular networks from the same tissue block, and robust capture of position, size and shape of individual hepatocytes, as well as entire lobules from the same tissue specimen. In addition to the protocols, we have also established image analysis software that allows relational and hierarchical quantifications of different liver substructures (e.g. cells and vascular branches) and events (e.g. cell proliferation and death). Typical results acquired for routinely quantified parameters in adult mice (C57Bl6/N) include the hepatocyte volume (5,128.3 ± 837.8 µm(3)) and the fraction of the hepatocyte surface in contact with the neighbouring hepatocytes (67.4 ± 6.7 %), sinusoids (22.1 ± 4.8 %) and bile canaliculi (9.9 ± 3.8 %). Parameters of the sinusoidal network that we also routinely quantify include the radius of the sinusoids (4.8 ± 2.25 µm), the branching angle (32.5 ± 11.2°), the length of intersection branches (23.93 ± 5.9 µm), the number of intersection nodes per mm(3) (120.3 × 103 ± 42.1 × 10(3)), the average length of sinusoidal vessel per mm(3) (5.4 × 10(3) ± 1.4 × 10(3)mm) and the percentage of vessel volume in relation to the whole liver volume (15.3 ± 3.9) (mean ± standard deviation). Moreover, the provided parameters of the bile canalicular network are: length of the first-order branches (7.5 ± 0.6 µm), length of the second-order branches (10.9 ± 1.8 µm), length of the dead-end branches (5.9 ± 0.7 µm), the number of intersection nodes per mm(3) (819.1 × 10(3) ± 180.7 × 10(3)), the number of dead-end branches per mm(3) (409.9 × 10(3) ± 95.6 × 10(3)), the length of the bile canalicular network per mm(3) (9.4 × 10(3) ± 0.7 × 10(3) mm) and the percentage of the bile canalicular volume with respect to the total liver volume (3.4 ± 0.005). A particular strength of our technique is that quantitative parameters of hepatocytes and bile canalicular as well as sinusoidal networks can be extracted from the same tissue block. Reconstructions and quantifications performed as described in the current protocols can be used for quantitative mathematical modelling of the underlying mechanisms. Furthermore, protocols are presented for both human and pig livers. The technique is also applicable for both vibratome blocks and conventional paraffin slices.

Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME.

This review encompasses the most important advances in liver functions and hepatotoxicity and analyzes which mechanisms can be studied in vitro. In a complex architecture of nested, zonated lobules, the liver consists of approximately 80 % hepatocytes and 20 % non-parenchymal cells, the latter being involved in a secondary phase that may dramatically aggravate the initial damage. Hepatotoxicity, as well as hepatic metabolism, is controlled by a set of nuclear receptors (including PXR, CAR, HNF-4α, FXR, LXR, SHP, VDR and PPAR) and signaling pathways. When isolating liver cells, some pathways are activated, e.g., the RAS/MEK/ERK pathway, whereas others are silenced (e.g. HNF-4α), resulting in up- and downregulation of hundreds of genes. An understanding of these changes is crucial for a correct interpretation of in vitro data. The possibilities and limitations of the most useful liver in vitro systems are summarized, including three-dimensional culture techniques, co-cultures with non-parenchymal cells, hepatospheres, precision cut liver slices and the isolated perfused liver. Also discussed is how closely hepatoma, stem cell and iPS cell-derived hepatocyte-like-cells resemble real hepatocytes. Finally, a summary is given of the state of the art of liver in vitro and mathematical modeling systems that are currently used in the pharmaceutical industry with an emphasis on drug metabolism, prediction of clearance, drug interaction, transporter studies and hepatotoxicity. One key message is that despite our enthusiasm for in vitro systems, we must never lose sight of the in vivo situation. Although hepatocytes have been isolated for decades, the hunt for relevant alternative systems has only just begun.

Prediction and validation of cell alignment along microvessels as order principle to restore tissue architecture in liver regeneration.

Only little is known about how cells coordinately behave to establish functional tissue structure and restore microarchitecture during regeneration. Research in this field is hampered by a lack of techniques that allow quantification of tissue architecture and its development. To bridge this gap, we have established a procedure based on confocal laser scans, image processing, and three-dimensional tissue reconstruction, as well as quantitative mathematical modeling. As a proof of principle, we reconstructed and modeled liver regeneration in mice after damage by CCl(4), a prototypical inducer of pericentral liver damage. We have chosen the regenerating liver as an example because of the tight link between liver architecture and function: the complex microarchitecture formed by hepatocytes and microvessels, i.e. sinusoids, ensures optimal exchange of metabolites between blood and hepatocytes. Our model captures all hepatocytes and sinusoids of a liver lobule during a 16 days regeneration process. The model unambiguously predicted a so-far unrecognized mechanism as essential for liver regeneration, whereby daughter hepatocytes align along the orientation of the closest sinusoid, a process which we named "hepatocyte-sinusoid alignment" (HSA). The simulated tissue architecture was only in agreement with the experimentally obtained data when HSA was included into the model and, moreover, no other likely mechanism could replace it. In order to experimentally validate the model of prediction of HSA, we analyzed the three-dimensional orientation of daughter hepatocytes in relation to the sinusoids. The results of this analysis clearly confirmed the model prediction. We believe our procedure is widely applicable in the systems biology of tissues.

Digital twin demonstrates significance of biomechanical growth control in liver regeneration after partial hepatectomy.

Partial liver removal is an important therapy option for liver cancer. In most patients within a few weeks, the liver is able to fully regenerate. In some patients, however, regeneration fails with often severe consequences. To better understand the control mechanisms of liver regeneration, experiments in mice were performed, guiding the creation of a spatiotemporal 3D model of the regenerating liver. The model represents cells and blood vessels within an entire liver lobe, a macroscopic liver subunit. The model could reproduce the experimental data only if a biomechanical growth control (BGC)-mechanism, inhibiting cell cycle entrance at high compression, was taken into account and predicted that BGC may act as a short-range growth inhibitor minimizing the number of proliferating neighbor cells of a proliferating cell, generating a checkerboard-like proliferation pattern. Model-predicted cell proliferation patterns in pigs and mice were found experimentally. The results underpin the importance of biomechanical aspects in liver growth control.